CRUISE NAME(S) : at15_42 CRUISE DATES : 2009/01/09 to 2009/02/05 PORTS : Puerto Ayora, Ecuador to Puntarenas, Costa Rica SHIP NAME : Atlantis CHIEF SCIENTIST : J. McGuire DATABASE NAME : aship DATA FILES : at2009_008_58461.raw to at2009_035_64800.raw STATUS : to do done ------ ----------- averaged [ x ] loaded [ x ] NOTE: heading correction instrument exists NOTE: time-dependent heading corrections staged (cal/rotate/ens_hcorr.ang) --> check if applied using "rotate" 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 2009/01/09 16:19:22 off 300 60 16 5 8 66.90 0.00 0001 HEADING : PRIMARY : phins POSITIONS : gps (specify more detail if available) CALIBRATION : check parameters for original heading alignment additional rotation 0 final transducer angle is: (original transducer angle) - (-0.39) original scale factor 1 additional scale factor (1.003) COMMENTS : headings are a little strange: 'gyro' port: up to dday 13.5 : phins at 1Hz, subsampled by 10 (so 1 pt every 10sec) but phins failed 13.5-15.5 after 13.5 : gyro 'phins' port: phins heading used: phins port; deleted those two days. PROCESSOR : J.Hummon --- processing parameters ---------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 75 ## (determined from "sonar"): pingtype = nb ## (determined from "sonar"): instname = os75 ## ping_headcorr is False ## hcorr_inst is phins beamangle 30 configtype python cruisename at15_42 datatype uhdas ens_len 300 fixfile None pgmin 50 proc_engine python ref_method refsm sonar os75nb yearbase 2009 ======================================================= redo from scratch (1) mkdir config cd config uhdas_proc_gen.py -s at mv at_proc.py at15_42_proc.py # Edit the file and add these lines to the top yearbase = 2009 # usually year of first data logged uhdas_dir = '/home/data/at15_42' # path to uhdas data directory shipname = 'Atlantis' cruiseid = 'at15_42' # for titles (2) now we can use that file. Set up the processing directory by typing: cd .. adcptree.py os75nb --datatype uhdas --cruisename at15_42 (3) change directories to ADCP processing directory just created cd os75nb (4) create a quick_adcp.py control file unix: this example is written as a bash "heredoc" for cut-and-past ease, but you can create q_py.cnt with an editor to have the parts between (or including) the comments cat << EOF > q_py.cnt ####----- begin q_py.cnt------------ ## all lines after the first "#" sign are ignored --yearbase 2009 --cruisename at15_42 # used to identify configuration files # *must* match prefix of files in config dir --update_gbin ## NOTE: You must remake gbins for python processing if ## - you are not sure ## - if parameters for averaging changed ## - various other reasons. ## ==> just do it, and put them somewhere else ## (see next argument) # --py_gbindirbase gbin # (will put them locally) --configtype python ## <=== USE THE NEW FILE WE CREATED --sonar os75nb --dbname aship --datatype uhdas --ens_len 300 # --ping_headcorr ## DO NOT APPLY --proc_engine python ## use matlab to do calculations --max_search_depth 1000 EOF ---------- (5) run quick_adcp.py: quick_adcp.py --cntfile q_py.cnt --auto ---------- (6) (a) check calibration:heading correction device: figview.py cal/rotate/*png Must delete a little over 2 days around 13.5-15.5 because of bad headings (b) check calibration: ------------------- (7) edit "gautoedit" (go to "edit" directory with command line, type "gautoedit.py" # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto # check editing -- looks OK dataviewer.py - check calibrations again **watertrack** ------------ Number of edited points: 115 out of 115 amp = 1.0033 + 0.0000 (t - 21.3) phase = -0.39 + 0.0054 (t - 21.3) median mean std amplitude 1.0030 1.0033 0.0045 phase -0.3870 -0.3945 0.2489 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS fixfile is aship.gps calculation done at 2013/01/29 23:12:23 xducer_dx = 3.447210 xducer_dy = 5.901256 signal = 2556.137637 ------------ still need scale factor: quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.003 --rotate_angle -0.39 --auto **watertrack** ------------ Number of edited points: 115 out of 115 amp = 1.0003 + 0.0000 (t - 21.3) phase = -0.00 + 0.0053 (t - 21.3) median mean std amplitude 1.0000 1.0003 0.0045 phase 0.0020 -0.0045 0.2493 ------------ OK, done make plots: quick_web.py --interactive - view with a browser, look at webpy/index.html ==================================================================