LAST CHANGED : Mon Feb 10 11:21:11 HST 2014
CRUISE NAME(S) : ps0918
CRUISE DATES : 2009/09/22 14:35:59 to 2009/09/22 23:13:33
SHIP NAME : Point Sur
PORTS : Monterey, CA to Monterey, CA
CHIEF SCIENTIST : Curt Collins
DATABASE NAME : aship
DATA FILES : *LTA
STATUS : to do done
------ -----------
averaged [ x ] LTA files
loaded [ x ]
NOTE: No automated time-dependent
heading correction exists
check heading correction [ x ]
calibration [ x ]
edited [ x ]
re-check heading correction [ x ]
check editing [ x ]
figures [ x ]
INSTRUMENT : os75
ACQUISITION :
PROGRAM : lta
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/09/22 13:40:59 on 297 80 8 3 8 1.18 0.00 0001
2009/09/22 16:12:50 on 298 80 8 3 8 1.18 0.00 0001
2009/09/22 22:18:30 on 297 80 8 3 8 1.18 0.00 0001
HEADING :
PRIMARY : heading from gyro?? (specify primary heading device here)
CORRECTION : specify heading correction device (if it exists)
POSITIONS : gps
CALIBRATION : original heading angle 1.18
additional rotation 1.5
final transducer angle is:
(original transducer angle) - (rotate_angle)
applied scale factor 1
additional scale factor (none)
COMMENTS : looks OK
PROCESSOR : I.M. Persistent
--- processing parameters ----------
## (determined from "sonar"): model = os
## (determined from "sonar"): frequency = 75
## (determined from "sonar"): pingtype = bb
## (determined from "sonar"): instname = os75
badbeam None
beamangle 30
configtype python
cruisename ps0918
datatype lta
dbname aship
ens_len 300
fixfile aship.gps
pgmin 50
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 2009
This 'help' refers to the qpy_demos ps0918 cruise.
=======================================================
(1) find out information about the data
Make a **project** directory (ps0918_vmdas)
From here:
vmdas_info.py --logfile lta_manual_info.txt os ../../vmdas_data/ps0918/*LTA
vmdas_info.py --logfile enr_manual_info.txt os ../../vmdas_data/ps0918/*ENR
The output says:
- 3 files, all have data
- averaging length: 300sec
- instrument frequency 75
- pingtypes: bb pings
- bottom track was on for part of the cruise
- beam angle 30 (normal for OS75)
- transducer angle (EA) 1.18
- messages available and what was used:
N1R HEHDT <-- VmDAS primary heading
N2R GPGGA
N2R GPHDT
N2R PASHR,AT2 <-- VmDAS backup heading
N2R PASHR,ATT
N2R PASHR,POS
N3R GPGGA <-- VmDAS primary position
N3R GPGLL
Look at the info files and note:
- what kind of pings
- averaging length
- are filenames sorted in chronological order?
- what was the heading source?
(2) Start a text file "os75bb_enr.txt" adjacent to the processing directory,
"os75bb_lta" and put this information in it:
--> these are 'bb' pings
--> averaged over 300sec
--> frequency is 75kHz, so "sonar" os "os75bb"
--> primary heading was HEHDT
--> there is an Ashtech, used as secondary heading
(3) convert vmdas data to uhdas-style data
mkdir config
reform_vmdas.py --project_dir_path ./ --vmdas_dir_path ../../vmdas_data/ps0918 --uhdas_style_dir ../../uhdas_style_data --cruisename ps0918_manual --start ..
cd config
python3 vmdas2uhdas.py
(4) create the processing configuration file
proc_starter.py reform_defs.py
cd ..
(5a) set up processing directory:
adcptree.py os75bb_manual --datatype uhdas --cruisename ps0918_manual
cd os75bb_manual
# make q_py.cnt
quick_adcp.py --cntfile q_py.cnt --auto
(5b) create a quick_adcp.py control file
create a quick_adcp.py control file
windows:
- you can create q_py.cnt with an editor to have the
parts between (or including) the comments
- it may be easier to call it "q_py.txt" to avoid
fighting with your editor, which may try to set the suffix
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
##### q_py.cnt is
--yearbase 2009 ## required, for decimal day conversion
## (year of first data)
--cruisename ps0918_manual ## *must* match prefix in config dir
--dbname aship ## database name; in adcpdb. eg. a0918
##
--datatype uhdas ## datafile type
--sonar os75bb ## specify instrument letters, frequency,
## (and ping type for ocean surveyors)
--ens_len 300 ## averages of 300sec duration
##
--update_gbin ## required for this kind of processing
--configtype python ## file used in config/ dirctory is python
##
--ping_headcorr ## ps0918_manual_proc.py says use HDT first,
## correct to ashtech
--max_search_depth 1500 ## try to identify the bottom and eliminate
## data below the bottom IF topo says
## the bottom is shallower than 1000m
#### end of q_py.cnt
EOF
(6) run quick_adcp.py:
quick_adcp.py --cntfile q_py.cnt --auto
(7) Check calibrations
**watertrack**
------------
Number of edited points: 2 out of 2
amp = 1.0045 + 0.1842 (t - 264.8)
phase = -0.10 + 1.8246 (t - 264.8)
median mean std
amplitude 1.0045 1.0045 0.0148
phase -0.1000 -0.1000 0.1471
------------
codaspy:(os75bb_manual.calibrated)-PY3$ catbt
**bottomtrack**
------------
unedited: 20 points
edited: 16 points, 2.0 min speed, 2.5 max dev
median mean std
amplitude 0.9997 1.0000 0.0025
phase 0.2755 0.2877 0.1113
------------
(8) apply rotation correction (no scale factor necessary)
(all in one line)
quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle 0.28 --auto
**bottomtrack**
------------
unedited: 20 points
edited: 18 points, 2.0 min speed, 2.5 max dev
median mean std
amplitude 1.0010 1.0005 0.0044
phase 0.0042 0.0563 0.1798
------------
(9) edit bad data
cd edit
gautoedit.py
# apply editing:
cd ..
quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto
(10) make web plots
quick_web.py --interactive
- view with a browser, look at webpy/index.html
(11) extract data
quick_adcp.py --steps2rerun matfiles --auto
adcp_nc.py adcpdb contour/os75bb ps0918_demo os75bb
