LAST CHANGED : 2017/04/13 21:50:38 CRUISE NAME(S) : tn188_transit CRUISE DATES : 2005/11/26 to 2005/11/29 SHIP NAME : Thomas G. Thompson PORTS : unspecified CHIEF SCIENTIST : unspecified DATABASE NAME : aship DATA FILES : tt2005_325_69562.raw to tt2005_336_45341.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 2005/11/22 19:24:23 off 300 80 8 5 16 58.90 0.00 0001 2005/11/23 20:24:24 off 300 80 8 5 16 58.90 0.00 0001 2005/11/24 21:24:27 off 300 80 8 5 16 58.90 0.00 0001 2005/11/25 22:24:25 off 300 80 8 5 16 58.90 0.00 0001 2005/11/26 23:24:24 off 300 80 8 5 16 58.90 0.00 0001 2005/11/28 00:24:23 off 300 80 8 5 16 58.90 0.00 0001 2005/11/29 01:24:24 off 300 80 8 5 16 58.90 0.00 0001 2005/11/30 02:24:25 off 300 80 8 5 16 58.90 0.00 0001 2005/12/01 03:24:26 off 300 80 8 5 16 58.90 0.00 0001 2005/12/02 04:24:26 off 300 80 8 5 16 58.90 0.00 0001 2005/12/03 05:24:26 off 300 80 8 5 16 58.90 0.00 0001 HEADING : PRIMARY : heading from gyro 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 : heading alignment: 58.9 additional rotation: 0.203 final transducer angle is: 58.697 (original transducer angle) - (rotate_angle) applied scale factor: 1.0 additional scale factor: 1.014 COMMENTS : gaps in heading correction scattering layers PROCESSOR : Dr. Thomas Roc --- processing parameters ---------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 75 ## (determined from "sonar"): instname = os75 ## (determined from "sonar"): pingtype = bb badbeam None beamangle 30 configtype python cruisename tn188_transit datatype uhdas dbname aship ens_len 300 fixfile aship.gps hcorr_inst posmv pgmin 50 pingpref None proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source nav sonar os75bb txy_file aship.gps xducer_dx 0 xducer_dy 0 yearbase 2005 ========================================================================== 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 (1a) get some information about this cruise: uhdas_info.py --overview --logfile uhdas_info.txt ./TN188_transit We see: - positions from gpsnav and posmv - headings from gyro and posmv - 132 files; all present (1b) no need to create config/tn188_transit_proc.py because it had been created earlier while processing os75nb (2) now we can use that file. Set up the processing directory by typing: adcptree.py os75bb --datatype uhdas --cruisename tn188_transit (3) (a) change directories to ADCP processing directory just created cd os75bb (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 2005 --cruisename tn188_transit # 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) # The other option is to move the # original gbins out of the way, # to a new directory (eg. # gbin.orgimat) --configtype python ## <=== USE THE NEW FILE WE CREATED --sonar os75bb --dbname aship --datatype uhdas --ens_len 300 --ping_headcorr ## applies heading correction. ## settings found in config files --max_search_depth 3000 ## use topography for editing? ## 0 = "always use amplitude to guess the bottom; ## flag data below the bottom as bad ## -1 = "never search for the bottom" ## positive integer: use ADCP amp to autodetect ## the bottom. Only do this in "deep water", ## i.e. topo says bottom is deeper than this EOF (c) run quick_adcp.py: quick_adcp.py --cntfile q_py.cnt (4) review the data (a) check calibration:heading correction device: figview.py cal/rotate/*png conclude: there are gaps in heading thus patch_hcoor.py needs to be ran: cd cal/rotate patch_hcorr.py (1) inspect the newhcorr*.png files: figview.py *.png (2) remove earlier time-dependent heading correction rotate unrotate.tmp (3) apply new heading correction rotate rotate_fixed.tmp (4) run navigation steps and inspect calibrations cd ../.. quick_adcp.py --steps2rerun navsteps:calib --auto (b) check calibration: ------------------- ADCP watertrack calibration ## Number of edited points: 3 out of 4 amp = 1.0117 + -0.0017 (t - 333.8) phase = -0.03 + -0.1684 (t - 333.8) median mean std amplitude 1.0130 1.0117 0.0042 phase 0.1860 -0.0293 0.4028 ------------------- (c) look at the data: dataviewer.py (5) edit "gautoedit" cd edit gautoedit.py -n5 Notes: TR: Many scattering layers TR: Odd v surface currents around 329.0 days TR: same between 329.5 and 330.0 days TR: and again around 331.7 days # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto # check editing cd .. dataviewer.py -c os75nb os75bb # Notes: TR: Comparison between os75nb and os75bb TR: Solid difference between them both TR: hopefully will disapear with corrections (6) check calibrations again cd os75bb tail -20 cal/watertrk/adcpcal.out ------------ ADCP watertrack calibration ## Number of edited points: 2 out of 2 amp = 1.0140 + 0.0107 (t - 332.4) phase = 0.20 + 0.1818 (t - 332.4) median mean std amplitude 1.0140 1.0140 0.0014 phase 0.2030 0.2030 0.0240 ------------ TOO FEW POINTS TO USE; instead use values from the previous cruise (c) look at the data: dataviewer.py (5) edit "gautoedit" cd edit gautoedit.py -n5 #Notes: TR: Many scattering layers # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto # check editing -- looks OK dataviewer.py cd .. dataviewer.py -c os75nb os75bb # Notes: TR: Still solid difference at the beginning of the cruise TR: yet over all better now with amplitude and angle TR: corrections applied. TR: Editing is considered done (6) Reviews from Jules cd os75nb/cal/rotate \rm newhcorr_mpl_*.png patch_hcorr.png mv rotate.log rotate.log.old \rm newhcorr* patch_hcorr.py - nav - zoom - delete (select) one point -- target = -0.9 NOW: (1) inspect the newhcorr*.png files: figview.py *.png (2) remove earlier time-dependent heading correction rotate unrotate.tmp (3) apply new heading correction rotate rotate_fixed.tmp (4) run navigation steps and inspect calibrations cd ../.. quick_adcp.py --steps2rerun navsteps:calib --auto Apply calibrations from previous cruise: quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_amplitude 1.018 --auto catwt (7) make plots: quick_web.py --interactive (8) extract data adcp_nc.py adcpdb contour/os75bb tn188_transit os75bb quick_adcp.py --steps2rerun matfiles --auto # done with processing. (9) Add cruise_info.txt to the top of this file. # Edit with correct info **watertrack** ------------ Number of edited points: 2 out of 2 amp = 0.9975 + 0.0160 (t - 332.4) phase = 0.11 + 1.0856 (t - 332.4) median mean std amplitude 0.9975 0.9975 0.0021 phase 0.1085 0.1085 0.1435 ------------ cp ../os75nb/webpy/sectinfo.txt webpy/ quick_web.py --redo adcp_nc.py adcpdb contour/os75bb_short tn188_transit os75bb adcp_nc.py adcpdb contour/os75bb_long tn188_transit os75bb --long quick_adcp.py --steps2rerun matfiles --auto TR: Notes for Jules - done and reviewed