LAST CHANGED : 2017/04/07 12:10:58 CRUISE NAME(S) : P18_leg2 (rb1606_leg2) CRUISE DATES : 2017/01/01 to 2017/02/02 SHIP NAME : Ronald H. Brown PORTS : Hanga Roa, Chile to Punta Arenas, Chile CHIEF SCIENTIST : R. Sonnerup DATABASE NAME : aship DATA FILES : rb2017_000_00000.raw to rb2017_032_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 : os75 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 2016/01/01 00:04:59 off 300 80 8 5 8 51.40 0.00 0001 2016/01/02 01:04:57 off 300 80 8 5 8 51.40 0.00 0001 2016/01/03 02:04:57 off 300 80 8 5 8 51.40 0.00 0001 2016/01/04 06:51:52 off 300 80 8 5 8 51.40 0.00 0001 2016/01/05 23:30:19 off 300 80 8 5 8 51.40 0.00 0001 2016/01/06 19:58:30 off 300 80 8 5 8 51.40 0.00 0001 2016/01/07 20:58:30 off 300 80 8 5 8 51.40 0.00 0001 2016/01/08 01:42:20 off 300 80 8 5 8 51.40 0.00 0001 2016/01/09 02:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/10 03:42:22 off 300 80 8 5 8 51.40 0.00 0001 2016/01/11 04:42:22 off 300 80 8 5 8 51.40 0.00 0001 2016/01/13 06:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/14 07:42:23 off 300 80 8 5 8 51.40 0.00 0001 2016/01/15 08:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/16 09:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/18 11:42:23 off 300 80 8 5 8 51.40 0.00 0001 2016/01/20 13:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/21 14:42:22 off 300 80 8 5 8 51.40 0.00 0001 2016/01/22 15:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/23 16:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/24 17:42:21 off 300 80 8 5 8 51.40 0.00 0001 2016/01/25 18:42:23 off 300 80 8 5 8 51.40 0.00 0001 2016/01/26 00:04:59 off 300 80 8 5 8 51.40 0.00 0001 2016/01/26 14:05:00 off 300 80 8 5 8 51.40 0.00 0001 2016/01/26 16:04:58 off 300 80 8 5 8 51.40 0.00 0001 2016/01/28 04:04:58 off 300 80 8 5 8 51.40 0.00 0001 2016/01/29 05:04:58 off 300 80 8 5 8 51.40 0.00 0001 2016/01/30 06:04:57 off 300 80 8 5 8 51.40 0.00 0001 2016/01/31 07:05:01 off 300 80 8 5 8 51.40 0.00 0001 2016/02/01 08:04:58 off 300 80 8 5 8 51.40 0.00 0001 2016/02/02 09:04:58 off 300 80 8 5 8 51.40 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 : Heading alignment: 51.4 additional rotation: -0.1085 final transducer angle is: 51.5085 (original transducer angle) - (rotate_angle) applied scale factor: 1.0 additional scale factor: 1.001 COMMENTS : reprocessed data in order to clean it up gaps in heading correction numerous scattering layers # There was a bug (oversight) in our treatment of zmq_publisher.py: # (1) zmq_publisher.py is started at boot and its yearbase was determined at boot # (2) therefore if zmq rolls over to a new year, the $UNIXD lines would have # dates higher than 365 # # eg. ./RB1606_Leg2 # ./RB1606_Leg2b # ./RB1606_Leg2c # ./RB1606_Leg2d # # And because the cruise was started in the same year, everything should have been OK # # (3) But then they rebooted and the UNIXD in gpsnav (zmq_publisher) was reset to the # subsequent year, but the cruise was not, so the base year changed # # The original solution is below, but that resulted in the whole cruise # having time range 2016/01/01 00:04:59 to 2016/02/02 16:04:59 # but it should have been 2017/01/01 # # # solution : add 366*86400 seconds PROCESSOR : Dr. Thomas Roc --- processing parameters ---------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 75 ## (determined from "sonar"): instname = os75 ## (determined from "sonar"): pingtype = bb badbeam None beamangle 30 configtype python cruisename RB1606_Leg2 datatype uhdas dbname aship ens_len 300 fixfile aship.gps hcorr_inst posmv pgmin 50 pingpref None proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source nav sonar os75bb txy_file aship.gps xducer_dx 0 xducer_dy 0 yearbase 2016 ======================================================================== CRUISE=RB1606_Leg2 sonar=os75bb (1) config dir mkdir config cd config cp $datadir/RB1606_Leg2_proc.py/config - copied config directory from at-sea processing - edit to point to data source - NOTE: gpsnav has UNIXD off by a year so use POSMV instead. cd .. (2) now we can use that file. Set up the processing directory by typing: adcptree.py os75bb --datatype uhdas --cruisename RB1606_Leg2 cd os75bb cat << EOF > q_py.cnt ####----- begin q_py.cnt------------ ## all lines after the first "#" sign are ignored ## python processing --yearbase 2016 --cruisename RB1606_Leg2 # used to identify configuration files # *must* match prefix of files in config dir --configtype python ## <=== USE THE NEW FILE WE CREATED --sonar os75bb --dbname aship --datatype uhdas --ens_len 300 --ping_headcorr ## applies heading correction. ## settings found in config files --max_search_depth -1 ## -1 = "never search for the bottom" EOF (c) run quick_adcp.py: quick_adcp.py --cntfile q_py.cnt --auto (4) review the data (a) check calibration:heading correction device: figview.py cal/rotate/*png conclude: gaps, run patch_hcorr.py (b) check calibration: catwt **watertrack** ------------ Number of edited points: 189 out of 189 amp = 1.0010 + -0.0000 (t - 12.8) phase = -0.11 + 0.0023 (t - 12.8) median mean std amplitude 1.0010 1.0010 0.0048 phase -0.1060 -0.1114 0.2957 ------------ catbt...no bottom track catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2017/04/08 03:45:23 xducer_dx = 0.786464 xducer_dy = 7.526452 signal = 2281.250212 ------------ (c) look at the data: dataviewer.py (5) edit "gautoedit" cd edit gautoedit.py -n6 # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto # check editing -- looks OK dataviewer.py (6) check calibrations again and apply correction(s) catwt **watertrack** ------------ Number of edited points: 184 out of 188 amp = 1.0011 + -0.0000 (t - 12.6) phase = -0.11 + 0.0024 (t - 12.6) median mean std amplitude 1.0010 1.0011 0.0048 phase -0.1065 -0.1097 0.2962 ------------ catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2017/04/21 21:35:05 xducer_dx = 0.896718 xducer_dy = 6.859967 signal = 2281.250212 ------------ # Apply scale factor, angle and transducer offset correction: quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_amplitude 1.001 --rotate_angle -0.1085 --xducer_dx 1 --xducer_dy 7 --auto catwt **watertrack** ------------ Number of edited points: 184 out of 188 amp = 1.0001 + -0.0000 (t - 12.6) phase = -0.00 + 0.0025 (t - 12.6) median mean std amplitude 1.0000 1.0001 0.0048 phase 0.0000 -0.0010 0.2962 ------------ # OK, done (7) make plots: quick_web.py --interactive (8) extract data adcp_nc.py adcpdb contour/os75bb RB1606_Leg2 os75bb quick_adcp.py --steps2rerun matfiles --auto