LAST CHANGED : Wed May 12 18:09:13 HST 2010 CRUISE NAME(S) : at15_63 CRUISE DATES : 2010/03/15 to 2010/04/14 PORTS : Galapagos, Equador to Galapagos, Equador SHIP NAME : Atlantis CHIEF SCIENTIST : John Sinton DB_NAME : a_at DATA FILES : at2010_073_57600. - at2010_103_28800.raw STATUS : averaged, loaded, edited, heading corrected, done. POSITIONS : gps STATUS : to do done ------ ----------- averaged [ x ] loaded [ x ] NOTE: time-dependent heading corrections applied IN the ensembles check heading correction [ x ] calibration [ x ] edited [ x ] re-check heading correction [ x ] check editing [ x ] figures [ x ] ready for 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 2010/03/07 17:38:49 off 300 60 16 5 8 66.90 0.00 0001 HEADING : PRIMARY : phins CALIBRATION : original heading alignment: 66.9 : need it to be 67.30 , so -0.4 calibration : (consistent with watertrack) : : tweak by -0.4 : amplitude: 1.00 amplitude, tweaked by 1.006 COMMENTS : no heading correction: phins was coming in the gyro port. PROCESSOR : J.Hummon ======================================================= # # written 2012/07/14 16:04:05 # this file contains CODAS quick_adcp.py info #name, value pairs #-------------------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 75 ## (determined from "sonar"): pingtype = nb ## (determined from "sonar"): instname = os75 ## ping_headcorr is True ## ping_headcorr is False beamangle 30 configtype python cruisename at15_63 datatype uhdas ens_len 300 fixfile None pgmin 50 proc_engine python ref_method refsm sonar os75nb yearbase 2010 ======================================================= (1) mkdir config cd config uhdas_proc_gen.py -s at mv at_proc.py at15_63_proc.py # edit # add these yearbase = 2010 # usually year of first data logged uhdas_dir = '/home/data/at15_63' # path to uhdas data directory shipname = 'Atlantis' # for documentation cruiseid = 'at15_63' # for titles ## change this hdg_inst = 'phins' # same as gyro, as it turns out... # disable this hcorr_inst = '' # disable with '' (not None) hcorr_msg = '' # disable with '' (not None) ## cd .. (2) adcptree.py os75nb --datatype uhdas --cruisename at15_63 (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 2010 --cruisename at15_63 # 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 matlab ## python processing can read either type --sonar os75nb --dbname aship --datatype uhdas --ens_len 300 --ping_headcorr ## applies heading correction. ## settings found in config files --proc_engine python ## use python to do calculations --configtype python --max_search_depth 500 ## do not look for the bottom in open ocean 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 conclude: no action needed (b) check calibration: **watertrack** ------------ Number of edited points: 109 out of 113 amp = 1.0065 + 0.0002 (t - 88.8) phase = -0.41 + 0.0012 (t - 88.8) median mean std amplitude 1.0060 1.0065 0.0088 phase -0.3700 -0.4051 0.4187 ------------ **watertrack** (cal/watertrk/adcpcal.out) ------------------- (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 **watertrack** ------------ Number of edited points: 109 out of 111 amp = 1.0065 + 0.0002 (t - 88.8) phase = -0.41 + 0.0011 (t - 88.8) median mean std amplitude 1.0060 1.0065 0.0088 phase -0.3700 -0.4058 0.4197 ------------ dataviewer.py - check calibrations again EG. still need scale factor: quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.006 --rotate_angle -0.4 --auto **watertrack** ------------ Number of edited points: 109 out of 111 amp = 1.0006 + 0.0002 (t - 88.8) phase = -0.01 + 0.0009 (t - 88.8) median mean std amplitude 1.0000 1.0006 0.0087 phase 0.0300 -0.0056 0.4192 ------------ make plots: quick_web.py --interactive - view with a browser, look at webpy/index.html