LAST CHANGED : 2016/04/21 06:37:40 CRUISE NAME(S) : I08S (rr1603) CRUISE DATES : 2016/02/08 to 2016/03/15 SHIP NAME : R.Revelle PORTS : Fremantle, Australia to Fremantle, Australia CHIEF SCIENTIST : Alison MacDonald DATABASE NAME : aship DATA FILES : rr2016_038_48356.raw to rr2016_074_72000.raw STATUS : to do done ------ ----------- averaged [ x ] redone from scratch 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/02/08 13:30:58 off 300 40 8 5 4 -5.00 0.00 0001 2016/02/13 23:05:32 off 300 55 8 5 4 -5.00 0.00 0001 2016/02/26 16:58:52 off 300 60 8 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: -7.7 additional rotation 0.12 final transducer angle is: (original transducer angle) - (rotate_angle) (-7.7) - (0.12) = -7.82 (2) scale factor: applied scale factor 1.009 (large, but seems to work) (3) Estimate of offset between ADCP transducer and gps values are: ADCP (dx=startboard, dy=fwd) meters from GPS original offsets: --xducer_dx 2 --xducer_dy 26 residual (close enough) xducer_dx = -3 xducer_dy = 1 COMMENTS : - no gaps in heading correction - few scattering layers; island wake at 32N? - internal waves dday 54.6-54.9? - bad weather : used navigation from uvship NOTE : some strange stripe at 200m early in the cruise (lost velocities but not visible in PG : early part of the cruise, looks like the top 4 bins were biased in some way. : had to NOT look for the bottom at all, because some other acoustic device pinged at *almost* the same rate and was interpreted as the bottom. required another pass with processing-from-scratch. PROCESSOR : J.Hummon --- processing parameters ---------- ============= os150nb dbinfo ================ #this file was automatically generated. DO NOT EDIT # # written 2016/07/03 16:42:56 # 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 = nb badbeam None beamangle 30 configtype python cruisename RR1603 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 os150nb txy_file aship.agt xducer_dx 2 xducer_dy 30 yearbase 2016 ===================== processing notes ================================= (1a) uhdas_info.py --overview --logfile uhdas_info.txt ../../uhdas_data/RR1603 We see: - positions from gpsnav, ashtech and seapath - headings from phins, seapath and ashtech - 439 files; all present (1b) copy proc/os150nb/config dir to here edit RR1603_proc.py #Edit the file and change the uhdas_dir uhdas_dir = '/Volumes/BIGDATA/uhdas_data_archive/revelle/RR1603' # go back one level cd .. (2) now we can use that file. Set up the processing directory by typing: adcptree.py os150nb_seapath_dxdy --datatype uhdas --cruisename RR1603 (3) (a) change directories to ADCP processing directory just created cd os150nb_seapath_dxdy 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 RR1603 # 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 os150nb --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 -1 ## 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: 162 out of 165 amp = 1.0089 + 0.0001 (t - 60.2) phase = 0.11 + 0.0008 (t - 60.2) median mean std amplitude 1.0090 1.0089 0.0062 phase 0.1270 0.1057 0.4686 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2016/07/04 02:52:39 xducer_dx = -3.931838 xducer_dy = 1.411197 signal = 2604.607394 ------------ (c) look at the data: dataviewer.py ## apply alternate navigation quick_adcp.py --steps2rerun navsteps:calib --refuv_source uvship --auto (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 still need scale factor: pretty high quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_amplitude 1.009 --rotate_angle -.11 --auto **watertrack** ------------ Number of edited points: 161 out of 165 amp = 0.9998 + 0.0001 (t - 60.2) phase = 0.22 + 0.0005 (t - 60.2) median mean std amplitude 1.0000 0.9998 0.0057 phase 0.2390 0.2217 0.4719 ------------ Whoops, wrong sign on phase quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle 0.23 --auto phase correction is very small: **watertrack** ------------ Number of edited points: 161 out of 165 amp = 1.0002 + 0.0002 (t - 60.2) phase = 0.01 + 0.0006 (t - 60.2) median mean std amplitude 1.0000 1.0002 0.0058 phase 0.0270 0.0120 0.4807 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2016/07/04 01:59:12 xducer_dx = -3.208209 xducer_dy = 1.146567 signal = 2604.660137 ------------ **watertrack** ------------ Number of edited points: 161 out of 165 amp = 0.9998 + 0.0001 (t - 60.2) phase = -0.01 + 0.0005 (t - 60.2) median mean std amplitude 1.0000 0.9998 0.0057 phase 0.0040 -0.0083 0.4721 ------------ 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/os150nb RR1603 os150nb quick_adcp.py --steps2rerun matfiles --auto ---- done with processing. Add os150nb/cruise_info.txt to the top of this file. Edit with correct info