LAST CHANGED : 2014/02/05 02:01:15 CRUISE NAME(S) : ka0903 CRUISE DATES : 2009/08/25 to 2009/09/24 SHIP NAME : Ka`imimoana PORTS : Honolulu, HI to Honolulu, HI CHIEF SCIENTIST : DATABASE NAME : aship DATA FILES : ka2009_236_04332.raw to ka2009_266_79200.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 2009/08/25 01:17:12 off 300 60 16 5 8 44.40 0.00 0001 POSITIONS : gps (NorthStar952, POSMV positions had trouble) CALIBRATION : original heading alignment: 44.4 final transducer angle (original transducer angle) - (rotate_angle) 44.4 - (-0.11) = 44.51 original scale factor 1 additional rotation or scale factor:1.002 COMMENTS : (1) gaps in posmv positions; (reprocess with other gps) (2) "gyro" feed is from the comparitor, not a gyro (3) secondary gps wandered a bit, causing problems with ships speed in some cases (deleted various 30min chunks) - this is going to persist for the rest of 2009 and possibly into 2010. (4) correcting for gps-transducer offset dx=-6, dy=29 (5) single-ping processing *included* gps-transducer offsets, and in postprocessing "uvship" was used as refuv source (improves ship speeds when PG is bad) PROCESSOR : JHummon ============================================ # # written 2014/02/04 15:50:33 # 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 badbeam None beamangle 30 configtype python cruisename ka0903 datatype uhdas dbname aship ens_len 300 fixfile aship.gps hcorr_inst posmv pgmin 50 proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source uvship sonar os75nb txy_file aship.agt xducer_dx -6 xducer_dy 29 yearbase 2009 ============================================ background notes on this cruise: What a mess. POSMV : gaps of up to 30 sec in position. headings are often bad gpsnav : positions aren't great gyro : has large gaps gyro1 : seems OK, but then failed gyro2 : seems OK solution: process using - gpsnav for positions - gyro12 (links of gyro1 and gyro2, with gyro filling in the gap) - posmv for "heading correction" What a mess. (0) This cruise came in many parts ka0903_legs/ ka_09_03/ ka_09_03a/ ka_09_03b/ ka_09_03c/ ka_09_03d/ cd datadir mkdir ka0903 for dir in ka0903_legs do link_uhdaslegs.py ka0903_legs/$dir ka0903 done ==== now try processing ========================== cd ../config uhdas_proc_gen.py -s ka mv ka_proc.py ka0903_proc.py ## edit, add yearbase = 2009 # usually year of first data logged uhdas_dir = '/home/data/ka0903' # path to uhdas data directory shipname = 'Kaimimoana' # for documentation cruiseid = 'ka0903' # for titles #adjust to use # gpsnav # gyro # posmv cd .. (2) now we can use that file. Set up the processing directory by typing: mkdir ka0903 cd ka0903 adcptree.py os75nb --datatype uhdas --cruisename ka0903 --configpath ../config NOTE: Now there is no error message about missing configuration files (3) (a) change directories to ADCP processing directory just created cd os75nb NOTE: there is a "config" directory and it has a copy of the file ka0903_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 2009 --cruisename ka0903 # 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!! # --py_gbindirbase gbin # (will put them adjacent to q_py.cnt) --configtype python ## <=== USE THE NEW FILE WE CREATED --sonar os75nb --dbname aship --datatype uhdas --ens_len 300 --xducer_dx -6 --xducer_dy 29 --refuv_source uvship --ping_headcorr ## applies heading correction. ## settings found in config files --max_search_depth 3000 #if the topography says the ocean is deeper # than this, do not autodetect the bottom. # (reduces likelihood of shallow scattering # layers being identified as the bottom) 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 gaps in heading correction: cd cal/rotate patch_hcorr.py NOW: (i) inspect the newhcorr*.png files: figview.py *.png (ii) remove earlier time-dependent heading correction rotate unrotate.tmp (iii) apply new heading correction rotate rotate_fixed.tmp (iv) run navigation steps and inspect calibrations cd ../.. quick_adcp.py --steps2rerun navsteps:calib --auto (b) check calibration: **watertrack** ------------ Number of edited points: 47 out of 53 amp = 1.0033 + 0.0009 (t - 248.5) phase = -0.11 + 0.0150 (t - 248.5) median mean std amplitude 1.0010 1.0033 0.0125 phase -0.1290 -0.1093 0.4848 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS fixfile is aship.gps calculation done at 2014/02/05 02:09:15 xducer_dx = 0.192962 xducer_dy = 0.155303 signal = 5582.593979 ------------ #=========== ## apply phase calib quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle -0.11 --rotate_amplitude 1.002 --auto (c) look at the data: dataviewer.py (5) edit "gautoedit" cd edit gautoedit.py # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto # check editing -- looks OK **watertrack** ------------ Number of edited points: 44 out of 53 amp = 1.0012 + 0.0009 (t - 248.7) phase = 0.07 + 0.0066 (t - 248.7) median mean std amplitude 0.9990 1.0012 0.0128 phase -0.0135 0.0661 0.4251 ------------ 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.2013/webpy/sectinfo.txt webpy quick_web.py --redo --vecscale 1 - view with a browser, look at webpy/index.html (8) extract data adcp_nc.py adcpdb contour/os75nb ka0903 os75nb quick_adcp.py --steps2rerun matfiles --auto ==== done