LAST CHANGED : 2017/04/07 11:49:32 CRUISE NAME(S) : P18_leg1 (rb1606_leg1) CRUISE DATES : 2016/11/12 to 2016/12/22 SHIP NAME : Ronald H. Brown PORTS : San Diego, CA, USA to Hanga Roa, Chile CHIEF SCIENTIST : B. Carter DATABASE NAME : aship DATA FILES : rb2016_316_80678.raw to rb2016_356_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 : 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/11/12 22:29:42 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.13 final transducer angle is: 51.53 (original transducer angle) - (rotate_angle) applied scale factor: 1.0 additional scale factor: 1.0025 COMMENTS : gaps in heading correction yearbase was wrong. Had to reprocess data transducer offset had to be corrected PROCESSOR : Dr. J. Hummon --- 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_Leg1 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_Leg1 sonar=os75nb (1) config dir mkdir config cd config cp $datadir/RB1606_Leg1_proc.py/config - 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_Leg1 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_Leg1 # 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: (1) inspect the newhcorr*.png files: figview.py *.png (2) remove earlier time-dependent heading correction rotate unrotate.tmp (3) apply new heading correction rotate rotate_fixed.tmp (4) run navigation steps and inspect calibrations cd ../.. quick_adcp.py --steps2rerun navsteps:calib --auto (b) check calibration: **watertrack** ------------ Number of edited points: 239 out of 241 amp = 1.0026 + -0.0001 (t - 341.6) phase = -0.13 + 0.0024 (t - 341.6) median mean std amplitude 1.0030 1.0026 0.0036 phase -0.1340 -0.1265 0.2253 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2017/04/17 22:58:22 xducer_dx = 0.422084 xducer_dy = 6.266173 signal = 2892.628782 ------------ (c) look at the data: dataviewer.py (5) edit "gautoedit" cd edit gautoedit.py -n5 # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto # check editing -- looks OK dataviewer.py (6) - check calibrations again; apply quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_amplitude 1.0025 --rotate_angle -0.13 --auto **watertrack** ------------ Number of edited points: 239 out of 241 amp = 1.0001 + -0.0001 (t - 341.6) phase = 0.00 + 0.0023 (t - 341.6) median mean std amplitude 1.0000 1.0001 0.0036 phase -0.0060 0.0038 0.2277 ------------ quick_adcp.py --steps2rerun navsteps:calib --xducer_dx 0 --xducer_dy 6 --auto **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2017/04/18 01:10:58 xducer_dx = 0.480640 xducer_dy = 0.460840 signal = 2892.628782 ------------ OK, done (7) make plots: quick_web.py --interactive - or, to use the same sections as were used in the postprocessing demo mkdir webpy cp ../os75nb/webpy/sectinfo.txt webpy quick_web.py --redo - view with a browser, look at webpy/index.html (8) extract data adcp_nc.py adcpdb contour/os75nb RB1606_Leg1 os75nb quick_adcp.py --steps2rerun matfiles --auto