LAST CHANGED : 2016/05/18 02:19:44 CRUISE NAME(S) : I09N (rr1604) CRUISE DATES : 2016/03/21 to 2016/04/27 SHIP NAME : R.Revelle PORTS : Fremantle, Australia to Phuket, Thailand CHIEF SCIENTIST : Leticia Barbero DATABASE NAME : aship DATA FILES : rr2016_080_20895.raw to rr2016_117_50400.raw STATUS : to do done ------ ----------- averaged [ x ] loaded [ x ] from scratch 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 ] ready to submit to JASADCP 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/03/21 05:53:17 off 300 60 16 5 8 28.00 0.00 0001 HEADING : PRIMARY : heading from phins3 CORRECTION : heading correction from seapath NOTE: time-dependent heading corrections applied IN the ensembles (see cal/rotate/ens_hcorr.ang) POSITIONS : gps positions from gpsnav CALIBRATION : Check original processing parameters; Note changes during postprocessing (1) transducer angle: original processing: 28 additional rotation: 3.27 final transducer angle is: (original transducer angle) - (rotate_angle) 28 - (3.27) = 24.73 (2) scale factor: applied scale factor 1.0035 (pretty typical) (3) Estimate of offset between ADCP transducer and gps values are: ADCP (dx=startboard, dy=fwd) meters from GPS original offsets: correciton final estimate --xducer_dx 2 xducer_dx = -2 2 + (-2) = 0 --xducer_dy 26 xducer_dy = 3 26 + 3 = 29 COMMENTS All of the following reprocessing components helped: - seapath instead of ashtech (reduces wobble; ashtech isn't perfect) - transducer-gps offset (removes artifacts at ons-station/underway transitions) - navigation from 'uvship' (clears up a couple of egregious outliers) In addition: - strong small-scale internal wave field from 5N-7N (dday 103-106) ?? - little eddy (200m-300m) at 11.5S PROCESSOR : J.Hummon ============= os75nb dbinfo ================ #this file was automatically generated. DO NOT EDIT # # written 2016/05/17 16:05:19 # this file contains CODAS quick_adcp.py info #name, value pairs #-------------------- ## (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 RR1604 datatype uhdas dbname aship ens_len 300 fixfile aship.gps hcorr_inst seapath pgmin 50 pingpref None proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source uvship sonar os75nb txy_file aship.agt xducer_dx 0 xducer_dy 29 yearbase 2016 ================================= processing notes ================================= (1a) cruise was split into 2 segments; glue together : mkdir RR1604_legs mv RR1604* RR1604_legs mkdir RR1604_merged for dir in RR1604_legs/* do link_uhdaslegs.py $dir RR1604_merged done (1b) stage for single-ping processing: - mkdir config; cd config - # copy RR1604_legs/RR1604/proc/os75nb/config/RR1604_proc.py to processing directory - # edit RR1604_proc.py: make heading correction device be seapath - cd .. (2) now we can use that file. Set up the processing directory by typing: adcptree.py os75nb --datatype uhdas --cruisename RR1604 (3) (a) change directories to ADCP processing directory just created cd os75nb NOTE: there is a "config" directory and it has a copy of the file RR1603_proc.py that we created in #1 (b) create a quick_adcp.py control file cat << EOF > q_py.cnt ####----- begin q_py.cnt------------ ## all lines after the first "#" sign are ignored ## python processing --yearbase 2016 --cruisename RR1604 # used to identify configuration files # *must* match prefix of files in config dir --update_gbin ## NOTE: You should generally remake gbins ## - you are not sure ## - if parameters for averaging changed ## - various other reasons. ## ==> MAKE SURE you move the original gbin directory ## to another name first!! --sonar os75nb --dbname aship --datatype uhdas --ens_len 300 --ping_headcorr ## applies heading correction. ## settings found in config files --xducer_dx 2 --xducer_dy 26 --max_search_depth 3000 ## use topography for editing? ## 0 = "always use amplitude to guess the bottom; ## flag data below the bottom as bad ## -1 = "never search for the bottom" ## positive integer: use ADCP amp to autodetect ## the bottom. Only do this in "deep water", ## i.e. topo says bottom is deeper than this 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: no action needed (b) check calibration: **watertrack** ------------ Number of edited points: 227 out of 228 amp = 1.0037 + 0.0000 (t - 99.9) phase = 3.28 + -0.0019 (t - 99.9) median mean std amplitude 1.0030 1.0037 0.0047 phase 3.2710 3.2793 0.2454 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2016/05/18 02:19:42 xducer_dx = -2.557637 xducer_dy = 3.496350 signal = 3529.059243 ------------ # apply calibrations quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle 3.27 --rotate_amplitude 1.0035 --xducer_dx -2 --xducer_dy 3 --auto (c) look at the data: dataviewer.py # use navigation from uvship quick_adcp.py --steps2rerun navsteps:calib --refuv_source uvship --auto 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 **watertrack** ------------ Number of edited points: 228 out of 228 amp = 1.0002 + 0.0000 (t - 99.9) phase = 0.01 + -0.0019 (t - 99.9) median mean std amplitude 1.0000 1.0002 0.0046 phase 0.0005 0.0098 0.2441 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2016/07/04 22:55:53 xducer_dx = 0.212436 xducer_dy = 0.959606 signal = 3529.059243 OK, done (7) make plots: quick_web.py --interactive - view with a browser, look at webpy/index.html (8) extract data adcp_nc.py adcpdb contour/os75nb RR1604 os75nb quick_adcp.py --steps2rerun matfiles --auto Done. Ready to submit to JASADCP