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 60 16 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.133 final transducer angle is: 51.533 (original transducer angle) - (rotate_angle) applied scale factor:1.0 additional scale factor: none COMMENTS : reprocessed data in order to clean it up gaps in heading correction few scattering layers transducer offset correction # 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 = nb 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 os75nb txy_file aship.gps xducer_dx 0 xducer_dy 0 yearbase 2016 ======================================================================== CRUISE=RB1606_Leg2 sonar=os75nb - copied config directory from at-sea processing - edit to point to data source cd .. (2) now we can use that file. Set up the processing directory by typing: adcptree.py os75nb --datatype uhdas --cruisename RB1606_Leg2 cd os75nb 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 --update_gbin --configtype python ## <=== USE THE NEW FILE WE CREATED --sonar os75nb --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 ==> gaps in heading correction; run patch_hcorr.py NOW: (b) inspect the newhcorr*.png files: figview.py *.png (c) remove earlier time-dependent heading correction rotate unrotate.tmp (d) apply new heading correction rotate rotate_fixed.tmp (e) run navigation steps and inspect calibrations cd ../.. quick_adcp.py --steps2rerun navsteps:calib --auto (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 catwt **watertrack** ------------ Number of edited points: 185 out of 186 amp = 1.0004 + -0.0000 (t - 12.6) phase = -0.13 + 0.0090 (t - 12.6) median mean std amplitude 1.0000 1.0004 0.0040 phase -0.1420 -0.1258 0.2796 ------------ catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2017/04/21 19:48:12 xducer_dx = 0.412187 xducer_dy = 7.065594 signal = 2281.250212 ------------ # applying correction quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle -0.1339 --xducer_dx 0 --xducer_dy 7 --auto catwt **watertrack** ------------ Number of edited points: 185 out of 186 amp = 1.0004 + -0.0000 (t - 12.6) phase = 0.01 + 0.0090 (t - 12.6) median mean std amplitude 1.0000 1.0004 0.0040 phase -0.0080 0.0079 0.2798 ------------ # OK, done (7) make plots: quick_web.py --interactive (8) extract data adcp_nc.py adcpdb contour/os75nb RB1606_Leg2 os75nb quick_adcp.py --steps2rerun matfiles --auto