LAST CHANGED    : 2021/12/08 00:27:23
CRUISE NAME(S)  : km2118, hot334
CRUISE DATES    : 2021/12/02 to 2021/12/07
SHIP NAME       : Kilo Moana
PORTS           : Honolulu, HI to Honolulu, HI
CHIEF SCIENTIST : Karin Björkman
DATABASE NAME   : a_km
DATA FILES      : km2021_335_68882.raw to km2021_341_00000.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      : wh300

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/12/02 19:10:03  off  120  70   2   7  2  45.45   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.113
                  final transducer angle is: 45.563
                       (original transducer angle) - (rotate_angle) 45.45 - (-0.124) - (-0.012) - 0.023 =  45.563
                  applied scale factor 1.006
                  additional scale factor (none)

COMMENTS        : gaps in heading correction? One small one at beginning of cruise, ignored it.
                  scattering layers? No
                  bubbles? Yes, rough weather in this cruise. Upper water column data ugly.
                  underway bias? Minimal
                 
PROCESSOR       : T. Rohrer


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


## (determined from "sonar"): model = wh
## (determined from "sonar"): frequency = 300
## (determined from "sonar"): instname = wh300
## (determined from "sonar"): pingtype = bb
           beamangle   20
          cruisename   km2118
            datatype   uhdas
             ens_len   120
             fixfile   a_km.gps
           frequency   300
          hcorr_inst   posmv
            instname   wh300
               model   wh
            pingtype   bb
         proc_engine   python
          ref_method   refsm
     refuv_smoothwin   3
        refuv_source   nav
               sonar   wh300
            txy_file   a_km.agt
           xducer_dx   -2
           xducer_dy   1
            yearbase   2021
(unless otherwise stated, comments are relative
to root processing directory)

codaspy:(wh300)$ 

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

(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!
One gap right at beginning of cruise. I'm just going to delete that data.

Did NOT run patch_hcorr.py.


(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)

Number of edited points:  16 out of  16
   amp   = 0.9997  + 0.0004 (t - 338.3)
   phase =  -0.06  + -0.0912 (t - 338.3)
            median     mean      std
amplitude   1.0045   0.9997   0.0149
phase      -0.1110  -0.0611   0.5077
nav - pc   -1.0000  -0.8125   1.5586
var         0.0010   0.0009   0.0009
min var     0.0010   0.0008   0.0009
delta-u    -0.7200  -0.3737   3.4778
delta-v     0.1650  -0.3206   3.0367

Good enough to do edits first.

(3) edit: gautoedit.py
    from gautoedit, for every screen do
	- mark seabed if necessary
	- additional manual editing if necessary
	- apply edits (creates/appends new a*.asc files with new flags)

Minor edits.

(4) apply editing

   quick_adcp.py --steps2rerun apply_edit:navsteps:calib --proc_engine python --datatype uhdas --auto
	(if you don't want to sit and hit y for all questions, add "yes | " before the command (no quotes))

Number of edited points:  15 out of  16
   amp   = 1.0025  + 0.0001 (t - 338.3)
   phase =  -0.01  + -0.0906 (t - 338.3)
            median     mean      std
amplitude   1.0060   1.0025   0.0117
phase      -0.1240  -0.0074   0.4823
nav - pc    0.0000  -0.3333   1.3973
var         0.0010   0.0007   0.0007
min var     0.0010   0.0006   0.0006
delta-u    -0.7100  -0.1660   3.4944
delta-v     0.2100  -0.3507   3.1415

Tweak it!

quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.006 --rotate_angle -0.124  --auto

Number of edited points:  16 out of  17
   amp   = 0.9987  + -0.0009 (t - 338.2)
   phase =  -0.08  + 0.0101 (t - 338.2)
            median     mean      std
amplitude   1.0005   0.9987   0.0139
phase      -0.0120  -0.0791   0.9147
nav - pc    0.0000   0.1250   1.3601
var         0.0010   0.0009   0.0008
min var     0.0010   0.0008   0.0007
delta-u    -0.7300  -0.2544   3.3944
delta-v     0.2400  -0.2519   3.0606


(5) apply calibration

correct to third sig fig for amplitude; second for phase 
if both are positive values, enter a positive correction

Try fixing the phase.
quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.00 --rotate_angle -0.012  --auto

Number of edited points:  15 out of  17
   amp   = 0.9966  + 0.0001 (t - 338.3)
   phase =   0.13  + -0.0901 (t - 338.3)
            median     mean      std
amplitude   1.0000   0.9966   0.0118
phase       0.0230   0.1315   0.4830
nav - pc    0.0000  -0.1333   0.9155
var         0.0010   0.0008   0.0006
min var     0.0010   0.0008   0.0006
delta-u    -0.7100  -0.1660   3.4944
delta-v     0.2100  -0.3507   3.1415

Made it WORSE. WTF?

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

Number of edited points:  16 out of  17
   amp   = 0.9987  + -0.0009 (t - 338.2)
   phase =  -0.09  + 0.0092 (t - 338.2)
            median     mean      std
amplitude   1.0005   0.9987   0.0141
phase      -0.0120  -0.0882   0.9112
nav - pc    0.0000   0.1250   1.3601
var         0.0010   0.0009   0.0008
min var     0.0010   0.0008   0.0007
delta-u    -0.7300  -0.2544   3.3944
delta-v     0.2400  -0.2519   3.0606

Back to where it was. Guess I'm going to leave well enough alone.


(6) web plots:

quick_web.py --interactive
cd ../
cat wh300/cruise_info.txt dm_notes_wh300.txt > wh300.txt

ALL done