LAST CHANGED    : 2021/11/02 17:10:01
CRUISE NAME(S)  : km2116, hot333
CRUISE DATES    : 2021/10/27 to 2021/11/2
SHIP NAME       : Kilo Moana
PORTS           : Honolulu, HI to Honolulu, HI
CHIEF SCIENTIST : Karin Björkman
DATABASE NAME   : a_km
DATA FILES      : km2021_299_74961.raw to km2021_305_57600.raw

STATUS          :  to do                           done
                  ------                       -----------
                  averaged                       [  x   ]
                  loaded                         [  x   ]

                  NOTE: heading correction instrument exists
                        NOTE: time-dependent heading corrections
                        applied IN the ensembles
                        (see   cal/rotate/ens_hcorr.ang)

                  check heading correction       [  x   ]
                  calibration                    [  x   ]
                  edited                         [  x   ]
                  re-check heading correction    [  x   ]
                  check editing                  [  x   ]
                  figures                        [  x   ]

INSTRUMENT      : os38

ACQUISITION     : 
     PROGRAM    : uhdas
     PROCESSING:    : python

LOGGING         : 

     PARAMETERS : BT : bottom track mode (on or off)
                  SI : sampling interval or averaging period for ensemble (sec)
                  NB : number of bins
                  BL : bin length (m)
                  TD : transducer depth (m)
                  BK : blanking length (m)
                  HO : heading offset applied by DAS (deg)
                  HB : heading bias (deg)
                  CRPH : compensation for roll-pitch-heading, 1:on, 0:off)
                  
                  yy/mm/dd  hh:mm:ss   BT   SI  NB  BL  TD BK     HO     HB CRPH
                  2021/10/27 20:54:22  off  300  75  24   7 16  44.55   0.00 0001

HEADING         : 
     PRIMARY    : heading from gyro
     CORRECTION : heading correction from posmv
                        NOTE: time-dependent heading corrections
                        applied IN the ensembles
                        (see   cal/rotate/ens_hcorr.ang)

POSITIONS       : gps positions from posmv

CALIBRATION     : (check original processing parameters)
                  additional rotation -0.119
                  final transducer angle is:
                       (original transducer angle) - (rotate_angle) 44.55 - (-0.14) - 0.15 - (-0.129) = 44.669
                  applied scale factor 1.003
                  additional scale factor (none)

COMMENTS        : gaps in heading correction? 1 small set
                  scattering layers? Not bad
                  bubbles? Not bad
                  underway bias? Minimal
                  
PROCESSOR       : T. Rohrer


--- processing parameters ----------


## (determined from "sonar"): model = os
## (determined from "sonar"): frequency = 38
## (determined from "sonar"): instname = os38
## (determined from "sonar"): pingtype = nb
           beamangle   30
          cruisename   km2116
            datatype   uhdas
             ens_len   300
             fixfile   a_km.gps
           frequency   38
          hcorr_inst   posmv
            instname   os38
               model   os
            pingtype   nb
         proc_engine   python
          ref_method   refsm
     refuv_smoothwin   3
        refuv_source   nav
               sonar   os38nb
            txy_file   a_km.agt
           xducer_dx   16
           xducer_dy   6
            yearbase   2021


(unless otherwise stated, comments are relative
to root processing directory)

codaspy:(os38nb)$ 

===================================


(0) remake all the plots

quick_mplplots.py --plots2run all --yearbase 2021

(1) look at heading correction: are there gaps?

    figview.py cal/rotate/*png

[] There were no gaps!
[x] There WERE some gaps!

Just a few small ones.

patch_hcorr.py


- Water Track Calibration:
            median     mean      std
amplitude   1.0030   1.0024   0.0074
phase      -0.1400  -0.1018   0.3275



---New Calibration---
- Water Track Calibration:
            median     mean      std
amplitude   1.0030   1.0026   0.0074
phase      -0.1180  -0.1000   0.3269




 quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.003 --rotate_angle -0.14 --auto


(2) see what calibration looks like:

                amp  in [0.995, 1.005] : good enough to proceed 
                                         to the next step.  tweak later
                phase in (-.5, .5) : good enough to proceed 
                                         to the next step.  tweak later

    rule of thumb: errors when underway:

         scale factor of 1.02 (i.e. 2%) gives 10cm/s alongtrack error
         phase error of 1deg gives 10cm crosstrack error 

         (goal is 1cm/s)

tail -20 cal/watertrk/adcpcal.out

Number of edited points:  32 out of  32
   amp   = 0.9999  + -0.0004 (t - 197.7)
   phase =  -0.05  + 0.0459 (t - 197.7)
            median     mean      std
amplitude   1.0000   0.9999   0.0067
phase       0.0045  -0.0501   0.3977
nav - pc    1.0000   0.6875   1.4906
var         0.0010   0.0010   0.0006
min var     0.0010   0.0009   0.0006
delta-u    -1.2400  -0.3797   3.2409
delta-v    -0.0200  -0.2150   3.9792


Looking good, despite median/mean offsets.


(3) gautoedit.py

Minor edits.

(4) apply editing

    quick_adcp.py --steps2rerun apply_edit:navsteps:calib --proc_engine python --datatype uhdas --auto


(5) rechecking calibration

tail -20 cal/watertrk/adcpcal.out

** watertrack ** 

Number of edited points:  49 out of  50
   amp   = 0.9996  + 0.0001 (t - 302.9)
   phase =   0.04  + -0.0270 (t - 302.9)
            median     mean      std
amplitude   1.0000   0.9996   0.0072
phase       0.0150   0.0367   0.3314
nav - pc    0.0000   0.1020   1.9605
var         0.0000   0.0005   0.0005
min var     0.0000   0.0004   0.0005
delta-u    -0.5900  -0.2382   3.3926
delta-v    -0.0300  -0.1343   3.4020

A little bit left to tweak.

quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.00 --rotate_angle 0.15 --auto

Goddamn it. Tweaked it 0.15 instead of 0.015.

Number of edited points:  49 out of  50
   amp   = 0.9996  + 0.0001 (t - 302.9)
   phase =  -0.11  + -0.0264 (t - 302.9)
            median     mean      std
amplitude   1.0000   0.9996   0.0072
phase      -0.1290  -0.1130   0.3313
nav - pc    0.0000  -0.0204   1.9737
var         0.0000   0.0005   0.0005
min var     0.0000   0.0004   0.0005
delta-u    -0.5900  -0.2382   3.3926
delta-v    -0.0300  -0.1343   3.4020

quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.00 --rotate_angle -0.129 --auto

Number of edited points:  49 out of  50
   amp   = 0.9995  + 0.0001 (t - 302.9)
   phase =   0.02  + -0.0265 (t - 302.9)
            median     mean      std
amplitude   1.0000   0.9995   0.0073
phase       0.0030   0.0164   0.3301
nav - pc    0.0000   0.0816   1.9773
var         0.0010   0.0005   0.0005
min var     0.0000   0.0004   0.0005
delta-u    -0.5900  -0.2382   3.3926
delta-v    -0.0300  -0.1343   3.4020

There we go.

(6) web plots (make sure to save this file now):

Making web plots.

mkdir webpy
cp ../wh300/webpy/sectinfo.txt webpy/
quick_web.py --redo
cd ../
cat os38nb/cruise_info.txt dm_notes_os38nb.txt > os38nb.txt

ALL done
