LAST CHANGED : 2021/05/19 16:49:17 CRUISE NAME(S) : km2105, hot330 CRUISE DATES : 2021/05/15 to 2021/05/19 SHIP NAME : Kilo Moana PORTS : Honolulu, HI to Honolulu, HI CHIEF SCIENTIST : Tully Rohrer DATABASE NAME : a_km DATA FILES : km2021_134_70981.raw to km2021_138_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 : 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/05/15 19:45:01 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.118, -0.092 final transducer angle is: (original transducer angle) - (rotate_angle) 45.45 - -0.118 - -0.092 = 45.66 applied scale factor 1.007 additional scale factor (none) COMMENTS : gaps in heading correction? A few small ones. Good corrections scattering layers? No bubbles? Minimal underway bias? No 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 km2105 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! Small gaps at 135.6, 135.6, 137 Ran patch_hcorr.py ---Original Calibration--- - Water Track Calibration: median mean std amplitude 1.0060 1.0063 0.0077 phase -0.2220 -0.2689 0.4881 ---New Calibration--- - Water Track Calibration: median mean std amplitude 1.0070 1.0065 0.0076 phase -0.1880 -0.2644 0.4735 (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) Applying correction after heading correction: quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.0070 --rotate_angle -0.118 --auto Number of edited points: 50 out of 54 amp = 0.9995 + -0.0002 (t - 137.2) phase = -0.15 + -0.0566 (t - 137.2) median mean std amplitude 1.0000 0.9995 0.0076 phase -0.0820 -0.1544 0.4743 nav - pc 0.0000 0.0400 0.9467 var 0.0010 0.0007 0.0009 min var 0.0010 0.0006 0.0007 delta-u -0.6050 -0.2210 3.2743 delta-v -0.9350 -0.1972 3.5222 could tweak phase again, but going to edit 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: 51 out of 54 amp = 0.9995 + -0.0003 (t - 137.2) phase = -0.15 + -0.0566 (t - 137.2) median mean std amplitude 1.0000 0.9995 0.0075 phase -0.0920 -0.1540 0.4695 nav - pc 0.0000 0.0588 0.9468 var 0.0010 0.0230 0.1591 min var 0.0010 0.0006 0.0007 delta-u -0.5600 -0.1302 3.3056 delta-v -0.8400 -0.1649 3.4945 One more tweak needed. (5) apply calibration correct to third sig fig for amplitude; second for phase if both are positive values, enter a positive correction quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.00 --rotate_angle -0.092 --auto Number of edited points: 51 out of 54 amp = 0.9995 + -0.0003 (t - 137.2) phase = -0.06 + -0.0578 (t - 137.2) median mean std amplitude 1.0000 0.9995 0.0075 phase 0.0030 -0.0629 0.4687 nav - pc 0.0000 0.0392 0.9583 var 0.0010 0.0230 0.1590 min var 0.0010 0.0006 0.0007 delta-u -0.5600 -0.1302 3.3056 delta-v -0.8400 -0.1649 3.4945 Very nice. Split mean/median, as has been typical lately. (6) web plots: quick_web.py --interactive cd ../ cat wh300/cruise_info.txt dm_notes_wh300.txt > wh300.txt ALL done