LAST CHANGED : 2013/10/12 00:23:18 CRUISE NAME(S) : A16N_2013 (rb1304) CRUISE DATES : 2013/09/18 to 2013/08/12 SHIP NAME : Ronald H Brown PORTS : Reykjavik, Iceland to Madeira to Natal, Brazil CHIEF SCIENTIST : M.Baringer, J.Bullister DATABASE NAME : aship DATA FILES : rb2013_214_32084.raw to rb2013_274_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 ] OK to release 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 2013/08/03 08:59:48 off 300 60 16 5 8 22.00 0.00 0001 HEADING : PRIMARY : heading from mahrs CORRECTION : heading correction from posmv NOTE: time-dependent heading corrections applied IN the ensembles (see cal/rotate/ens_hcorr.ang) POSITIONS : gps positions from gpsnav CALIBRATION : original transducer orientation: 22.00 transducer depth: 5.0 (check original processing parameters) additional rotation 0.47 final transducer angle is: (original transducer angle) - (rotate_angle) applied scale factor 1.003 additional scale factor (none) COMMENTS : strong scattering layers ==> MAHRS is poorly calibrated/maintained and needs a latitude-dependent correction ==> POSMV is doing very badly: very little of the data are trustworthy, but in a gross sense it can provide the latitudinal correction for the MAHRS some kind of oscillation in deeper bins when on station 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 badbeam None beamangle 30 configtype python cruisename rb1304 datatype uhdas dbname aship ens_len 300 fixfile aship.gps hcorr_inst posmv pgmin 50 proc_engine python ref_method refsm sonar os75nb yearbase 2013 ================================================================== cd /home/data ln -s /home/ipu1/uhdas_data_archive/ron_brown/RB1304_legs/shipboard_adcp/RB1304_all rb1304 rsync -axl /home/data/rb1304/proc/os75nb/config . cd config # edit cd .. adcptree.py os75nb --datatype uhdas --cruisename rb1304 cd os75nb cat << EOF > q_py.cnt ####----- begin q_py.cnt------------ ## all lines after the first "#" sign are ignored --yearbase 2013 --cruisename rb1304 # 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 --ping_headcorr ## applies heading correction. ## settings found in config files --proc_engine python ## use python to do calculations --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 #small gaps in mahrs; go ahead and correct cd cal/rotate patch_hcorr.py - inspect the newhcorr*.png files: figview.py *.png - remove earlier time-dependent heading correction rotate unrotate.tmp - apply new heading correction rotate rotate_fixed.tmp - run navigation steps and inspect calibrations cd ../.. quick_adcp.py --steps2rerun navsteps:calib --auto (b) check calibration: adjust phase before editing: **watertrack** ------------ Number of edited points: 301 out of 308 amp = 1.0027 + 0.0000 (t - 243.7) phase = 0.47 + 0.0060 (t - 243.7) median mean std amplitude 1.0030 1.0027 0.0052 phase 0.4890 0.4714 0.3720 ------------ quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle 0.47 --rotate_amplitude 1.003 --auto **watertrack** ------------ Number of edited points: 301 out of 308 amp = 0.9996 + 0.0000 (t - 243.7) phase = 0.00 + 0.0060 (t - 243.7) median mean std amplitude 1.0000 0.9996 0.0051 phase 0.0220 0.0017 0.3736 ------------ (5) edit "gautoedit" (go to "edit" directory with command line, type "gautoedit.py" # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto make plots: quick_web.py --interactive - view with a browser, look at webpy/index.html (8) extract data adcp_nc.py adcpdb contour/os75nb rb1304 os75nb quick_adcp.py --steps2rerun matfiles --auto