LAST CHANGED : 2016/01/27 21:35:11 CRUISE NAME(S) : ARC01 (hly1502) CRUISE DATES : 2015/09/11 to 2015/08/23 SHIP NAME : Healy PORTS : Dutch Harbor, AK to Dutch Harbor, AK CHIEF SCIENTIST : D.Kadko DATABASE NAME : aship DATA FILES : hly2015_220_77865.raw to hly2015_284_00000.raw STATUS : to do done ------ ----------- averaged [ x ] loaded [ x ] NOTE: heading correction instrument exists NOTE: time-dependent heading corrections staged but NOT applied (see cal/rotate/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 2015/08/09 21:42:46 off 300 80 8 8 8 43.40 0.00 0001 HEADING : PRIMARY : heading from gyro39 CORRECTION : heading correction from ashpaq5 NOTE: time-dependent heading corrections staged but NOT applied (see cal/rotate/hcorr.ang) --> check if applied using "rotate" POSITIONS : gps positions from gpsnav CALIBRATION : (check original processing parameters) additional rotation 0 final transducer angle is: (original transducer angle) - (rotate_angle) 43.4 - 0.53 = 42.87 applied scale factor 1 additional scale factor (1.006) COMMENTS : see notes below: POSMV and seapath were broken, ashtech had holes, gyro was MK39 (very good) but fell apart at 90N Lots of ice, poor scattering, some electrical bias PROCESSOR : D.Cardoso, J.Hummon --- processing parameters ---------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 75 ## (determined from "sonar"): instname = os75 ## (determined from "sonar"): pingtype = nb badbeam None beamangle 30 configtype python cruisename hly1502 datatype uhdas dbname aship ens_len 300 fixfile aship.gps hcorr_inst ashpaq5 pgmin 50 pingpref None proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source nav sonar os75nb txy_file aship.agt xducer_dx 0 xducer_dy 17 yearbase 2015 Process UHDAS data using Python ============================================= In order to use all-python processing, we will generate the configuration file necessary from modern ship's settings and edit as needed. This file will be (a) consistent with all-python processing (b) guarantees newest file format (c) is more flexible ***Special alteration of 'best' gbins in /home/data/hly1502/gbin/os75/time/*** - wrote a script to create a list of certain 'best' gbin files that have bad gyro data for 2 instruments os75/150 - 247.61, 253.1 are the dday times when Gyro data is bad,note Unix day is very similar to dd_day on this ship - replaced this range of bad gyro data with Ashteck data - checked the ashteck data for bad flags which show up in the gbins as Nans and Interpolated to remove the Nans - wrote new gbins 'best' files - replace the bad data near the pole region [247.675, 249.172] with a constant heading 90 **************************************************** (1) copy config directory from /home/data/hly1502/proc/os75nb/config edit hly1502_proc.py cd .. (2) now we can use that file. Set up the processing directory by typing: adcptree.py os75nb --datatype uhdas --cruisename hly1502 # using 'hly1502" will specify config/hly1502_proc.py) (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 hly1502_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 2015 --cruisename hly1502 # used to identify configuration files # *must* match prefix of files in config dir # --update_gbin ## see notes above: we are not updating these --configtype python ## <=== USE THE NEW FILE WE CREATED --sonar os75nb --dbname aship --datatype uhdas --ens_len 300 --xducer_dx 0 # not --xducer_dx = 0 --xducer_dy 27 # --ping_headcorr # do not apply heading correction --max_search_depth 3000 EOF (c) run quick_adcp.py: quick_adcp.py --cntfile q_py.cnt --auto switch over to navigation method "uvship" quick_adcp.py --steps2rerun navsteps:calib --refuv_source uvship --auto (4) review the data (a) check calibration:heading correction device: figview.py cal/rotate/*png conclude: no action needed did not apply heading correction (b) check calibration: ** ADCP watertrack calibration ** Number of edited points: 80 out of 100 amp = 1.0054 + -0.0001 (t - 249.9) phase = 0.56 + 0.0089 (t - 249.9) median mean std amplitude 1.0065 1.0054 0.0120 phase 0.5155 0.5558 0.8695 (5) edit "gautoedit" on the os75nb_uvship data 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; apply correction quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_amplitude 1.006 --rotate_angle 0.53 --auto **watertrack** ------------ Number of edited points: 80 out of 102 amp = 0.9992 + -0.0001 (t - 249.9) phase = 0.00 + 0.0083 (t - 249.9) median mean std amplitude 1.0000 0.9992 0.0121 phase -0.0280 0.0041 0.8760 ------------ ** guess_xducerxy** ------------ guessing ADCP (dx=startboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2016/02/03 20:13:31 xducer_dx = -0.097757 xducer_dy = 8.960362 signal = 8431.788284 ------------ (7) compare the os75 and os150 dataviewer.py -m5 -c os150nb_uvship os75nb_uvship (8)edit "gautoedit" on the os75nb_uvship data cd edit gautoedit.py -n5 # compare with os150 (9) - check calibrations again **watertrack** ------------ Number of edited points: 80 out of 101 amp = 0.9992 + -0.0001 (t - 250.0) phase = -0.01 + 0.0074 (t - 250.0) median mean std amplitude 1.0000 0.9992 0.0121 phase -0.0630 -0.0089 0.8500 ------------ OK, done (10) make plots: mkdir webpy cp ../ os75nb_uvship/webpy/sectinfo.txt webpy quick_web.py --redo - view with a browser, look at webpy/index.html (11) extract data adcp_nc.py adcpdb contour/os75nb hly1502_newgbins os75nb ## dbpath outfilebase cruisetitle sonar test your file by running adcp_ncdump.py contour/os75nb.nc quick_adcp.py --steps2rerun matfiles --auto ---- done with processing. Add os75nb/cruise_info.txt to the top of this file. Edit with correct info