LAST CHANGED : 2017/06/14 21:03:39 CRUISE NAME(S) : kn197_8 CRUISE DATES : 2010/05/22 to 2010/06/24 SHIP NAME : Knorr PORTS : Barbados to Barbados CHIEF SCIENTIST : P. Yager DATABASE NAME : a_kn DATA FILES : kn2010_141_46063.raw to kn2010_174_28800.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 2010/05/22 12:52:49 on 300 110 8 5 8 -44.05 0.00 0001 2010/05/26 08:48:57 on 300 110 8 5 8 -44.05 0.00 0001 2010/05/27 09:48:59 on 300 110 8 5 8 -44.05 0.00 0001 2010/05/30 12:48:54 on 300 110 8 5 8 -44.05 0.00 0001 2010/06/01 14:48:57 on 300 110 8 5 8 -44.05 0.00 0001 2010/06/13 00:34:08 on 300 110 8 5 8 -44.05 0.00 0001 2010/06/14 01:34:08 on 300 110 8 5 8 -44.05 0.00 0001 2010/06/22 00:07:57 on 300 110 8 5 8 -44.05 0.00 0001 2010/06/23 01:07:57 on 300 110 8 5 8 -44.05 0.00 0001 HEADING : PRIMARY : heading from gyro CORRECTION : heading correction from phins 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) additional rotation 0.02 final transducer angle is: (original transducer angle) - (rotate_angle) applied scale factor 1 additional scale factor 1.004 dx/dy correction was appleied: --xducer_dx 3 --xducer_dy 15 COMMENTS : No gaps in the heading correction. PROCESSOR : U. Pinho --- processing parameters ---------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 75 ## (determined from "sonar"): instname = os75 ## (determined from "sonar"): pingtype = bb badbeam None beamangle 30 configtype python cruisename kn197_8 datatype uhdas dbname a_kn ens_len 300 fixfile a_kn.gps hcorr_inst phins pgmin 50 pingpref None proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source nav sonar os75bb txy_file a_kn.gps xducer_dx 0 xducer_dy 0 yearbase 2010 ============================================= Cruise kn197_8 single ping and post processing: (1a) get some information about this cruise: uhdas_info.py --overview --logfile uhdas_info.txt /home/data/kn197_8/ We see: - positions from gpsnav, gyro and posmv - headings from gpsnav,time,phins,posmv,gyro - 383 files; all present gbin: INFO This cruise was probably originally processed using Matlab gbin: INFO Must regenerate gbins to use CODAS Python processing proc: os75bb 141.537 - 174.377 (2010/05/22 to 2010/06/24) proc: os75nb 141.537 - 174.377 (2010/05/22 to 2010/06/24) proc: wh300 141.535 - 174.398 (2010/05/22 to 2010/06/24) (1b)# in kn197_8 make a directory called 'config' # and run: mkdir config cd config uhdas_proc_gen.py -s kn # This creates a file called "kn_proc.py" # copy to a new version for our cruise (use the cruisename "kn197_8"): cp kn_proc.py kn197_8_proc.py # Edit the file and add these five (5) # lines to the top of kn197_8_proc.py # full path to uhdas data directory (do not use relative path) shipname = 'Knorr' cruiseid = 'kn197_8' yearbase = 2010 uhdas_dir = '/home/data/kn197_8' # go back one level cd .. (2) now we can use that file. Set up the processing directory by typing: adcptree.py os75bb --datatype uhdas --cruisename kn197_8 (3) (a) change directories to ADCP processing directory just created cd os75bb NOTE: there is a "config" directory and it has a copy of the file kn197_8_proc.py that we created in #2 (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 2010 --cruisename kn197_8 # 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!! --configtype python ## <=== USE THE NEW FILE WE CREATED --sonar os75bb --dbname a_kn --datatype uhdas --ens_len 300 --ping_headcorr ## applies heading correction. ## settings found in config files --max_search_depth 2000 ## 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 and this one in q_py.cnt file: --cruisename kn197_8 # used to identify configuration files # *must* match prefix of files in config dir mv /home/data/kn197_8/gbin /home/data/kn197_8/gbin.original (c) run quick_adcp.py: quick_adcp.py --cntfile q_py.cnt --auto # POST-PROCESSING: ------------------------------------------------------------ cd os75bb plot_nav.py nav/a_kn.gps - This is a cruise from south Caribean Sea (Barbados) to the Guianas and north of Brazil, and cruise back to Barbados - Caribbean. --------- (4) review the data (a) check calibration:heading correction device: figview.py cal/rotate/*png conclude: no action needed figview.py . (b) check calibration: catwt **watertrack** ------------ Number of edited points: 119 out of 129 amp = 1.0041 + -0.0001 (t - 157.4) phase = -0.02 + -0.0076 (t - 157.4) median mean std amplitude 1.0040 1.0041 0.0081 phase 0.0130 -0.0184 0.5936 ------------ catbt **bottomtrack** ------------ unedited: 1151 points edited: 588 points, 2.0 min speed, 2.5 max dev median mean std amplitude 1.0014 1.0016 0.0033 phase 0.0739 0.0893 0.2626 ------------ catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_kn.agt calculation done at 2017/06/14 21:03:37 xducer_dx = 2.628176 xducer_dy = 16.744705 signal = 3986.773865 ------------ quick_adcp.py --steps2rerun navsteps:calib --xducer_dx 3 --xducer_dy 17 --auto **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2017/06/13 20:57:31 xducer_dx = 0.119873 xducer_dy = -1.804224 signal = 3064.367030 ------------ quick_adcp.py --steps2rerun navsteps:calib --xducer_dy -2 --auto **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_kn.agt calculation done at 2017/06/14 21:14:47 xducer_dx = -0.335252 xducer_dy = -0.117725 signal = 3986.773865 ------------ (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 catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_kn.agt calculation done at 2017/06/15 01:25:15 xducer_dx = -0.699188 xducer_dy = -0.245500 signal = 3986.773865 ------------ cd edit gautoedit.py -n5 **watertrack** ------------ Number of edited points: 109 out of 125 amp = 1.0040 + -0.0001 (t - 157.2) phase = -0.02 + -0.0065 (t - 157.2) median mean std amplitude 1.0040 1.0040 0.0083 phase 0.0090 -0.0175 0.6168 ------------ applied a small amplitude and phase corrections: quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.004 --rotate_angle -0.02 --auto **watertrack** ------------ Number of edited points: 109 out of 125 amp = 1.0001 + -0.0001 (t - 157.2) phase = 0.00 + -0.0066 (t - 157.2) median mean std amplitude 1.0000 1.0001 0.0083 phase 0.0270 0.0027 0.6177 ------------ 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 **watertrack** ------------ Number of edited points: 102 out of 117 amp = 1.0001 + -0.0001 (t - 157.2) phase = 0.00 + -0.0040 (t - 157.2) median mean std amplitude 1.0000 1.0001 0.0090 phase 0.0245 0.0043 0.6366 ------------ OK, done (7) make plots: quick_web.py --interactive ff webpy/index.html (8) extract data adcp_nc.py adcpdb contour/os75bb_short kn197_8 os75bb adcp_nc.py --long adcpdb contour/os75bb_long kn197_8 os75bb quick_adcp.py --steps2rerun matfiles --auto ncdump -h contour/os75bb_short.nc ---- done with processing. Add os75bb/cruise_info.txt to the top of this file. Edit with correct info