LAST CHANGED : 2016/05/18 03:53:27 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_43200.raw STATUS : to do done ------ ----------- averaged [ x ] reprocess 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 ] ready to submit to JASADCP INSTRUMENT : os150 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 85 4 5 4 -5.00 0.00 0001 2016/03/23 07:53:17 off 300 85 4 5 4 -5.00 0.00 0001 2016/03/24 08:53:17 off 300 85 4 5 4 -5.00 0.00 0001 2016/03/26 10:53:15 off 300 85 4 5 4 -5.00 0.00 0001 2016/03/29 13:53:17 off 300 85 4 5 4 -5.00 0.00 0001 2016/04/02 17:53:16 off 300 85 4 5 4 -5.00 0.00 0001 2016/04/20 10:53:16 off 300 85 4 5 4 -5.00 0.00 0001 2016/04/21 11:53:17 off 300 85 4 5 4 -5.00 0.00 0001 2016/04/27 10:57:31 off 300 85 4 5 4 -5.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: -5 additional rotation: 2.77 final transducer angle is: (original transducer angle) - (rotate_angle) -5 - (2.77) - -7.77 (2) scale factor: applied scale factor 1.008 (pretty large) (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 = -4 2 + (-4) = -2 --xducer_dy 30 xducer_dy = 1 30 + (1) = 31 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 --- processing parameters ---------- #this file was automatically generated. DO NOT EDIT # # written 2016/05/17 17:39:31 # this file contains CODAS quick_adcp.py info #name, value pairs #-------------------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 150 ## (determined from "sonar"): instname = os150 ## (determined from "sonar"): pingtype = bb 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 os150bb txy_file aship.agt xducer_dx -2 xducer_dy 31 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 os150bb --datatype uhdas --cruisename RR1604 (3) (a) change directories to ADCP processing directory just created cd os150bb (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 os150bb --dbname aship --datatype uhdas --ens_len 300 --ping_headcorr ## applies heading correction. ## settings found in config files --xducer_dx 2 --xducer_dy 30 --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: 226 out of 228 amp = 1.0089 + 0.0000 (t - 99.9) phase = 2.77 + -0.0023 (t - 99.9) median mean std amplitude 1.0090 1.0089 0.0061 phase 2.7780 2.7664 0.3854 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2016/05/18 03:53:26 xducer_dx = -4.262983 xducer_dy = 1.010890 signal = 3527.821511 ------------ # apply calibrations quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle 2.77 --rotate_amplitude 1.008 --xducer_dx -4 --xducer_dy 1 --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: 226 out of 228 amp = 1.0007 + -0.0000 (t - 99.9) phase = -0.00 + -0.0023 (t - 99.9) median mean std amplitude 1.0010 1.0007 0.0062 phase 0.0105 -0.0034 0.3849 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2016/07/04 23:02:07 xducer_dx = 0.415240 xducer_dy = 0.308106 signal = 3527.821511 ------------ (7) make plots: mkdir webpy cp ../os75bb/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/os150bb RR1604 os150bb quick_adcp.py --steps2rerun matfiles --auto Done. Ready to submit to JASADCP