LAST CHANGED : Sat Mar 9 23:45:53 UTC 2019
CRUISE NAME(S) : km1001c
CRUISE DATES : 2010/01/30 17:41:17 to 2010/01/30 18:11:47
SHIP NAME : Kilo Moana
PORTS : Coos Bay, OR to Coos Bay, OR
CHIEF SCIENTIST : Brian Taylor
DATABASE NAME : a_km
DATA FILES : km2010_029_64800.raw to km2010_029_63671.raw
STATUS : to do done
------ -----------
averaged [ at sea ]
loaded [ at sea ]
------------------------
check heading correction [ x ]
calibration [ x ]
edited [ x ]
re-check heading correction [ x ]
check editing [ x ]
figures [ x ]
--> OK to release to JASADCP (PI gave permission)
INSTRUMENT : os38
ACQUISITION :
PROGRAM : uhdas
PING TYPE : nb
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
2010/01/30 17:46:23 off 300 70 24 7 16 44.55 0.00 0001
HEADING :
PRIMARY : gyro
CORRECTION : posmv
POSITIONS : gps (specify more detail if available)
CALIBRATION : original heading alignment: 44.55
additional rotation (no change)
additional scale factor (1.0065)
no gaps in heading correction
PROCESSOR : I.M. Persistent
===========================================
dbinfo.txt
#--------------------
## (determined from "sonar"): model = os
## (determined from "sonar"): frequency = 38
## (determined from "sonar"): pingtype = nb
## (determined from "sonar"): instname = os38
badbeam None
beamangle 30
configtype python
cruisename km1001c
datatype uhdas
dbname a_km
ens_len 300
fixfile a_km.gps
pgmin 50
proc_engine python
ref_method refsm
refuv_smoothwin 3
refuv_source nav
sonar os38nb
txy_file a_km.gps
xducer_dx 0
xducer_dy 0
yearbase 2010
===========================================
Cruise km1001c processing:
The os38nb at-sea processing directory was copied from
the "proc" directory of the km1001c cruise disk.
All postprocessing is done with full codas+python (i.e. using
numpy+matplotlib)
All commands are run from "os38nb" unless otherwise noted.
(1) Make directory km1001c for post-procssing
- copy proc/os38nb directory from cruise disk
(i) copy the os38nb directory to "os38nb_postproc"
(ii) make another copy called "os38nb_postproc_orig" for later comparison
- start editing a text file called os38nb.txt, with the
contents of os38nb/cruise_info.txt at the top.
- keep notes down below, and fill in above later
Look at the cruise track to see what we're expecting.
plot_nav.py nav/a_km.gps
- This is a cruise with a "patch test" to calibrate the
multibeam. There should be ample watertrack calibration
points from the patch test.
(2) check accurate heading device (POSMV):
- look at figures in cal/rotate *hcorr.png for gaps
(gaps would be red "+" signs -- there are none)
figview.py cal/rotate/*png
conclude: no action needed
(3) check calibration:
(a) look at heading correction: are there gaps? does it need to be patched?
(b) look at watertrack and bottom track calibrations: what can we expect?
(a)run this command, then add part (or all) of the
record displayed, to the documetation:
tail -20 cal/watertrk/adcpcal.out
**watertrack**
Number of edited points: 32 out of 36
amp = 1.0080 + -0.0027 (t - 30.3)
phase = -0.01 + 0.0589 (t - 30.3)
median mean std
amplitude 1.0060 1.0080 0.0085 <--- slight scale factor
phase 0.0775 -0.0058 0.4343 <--- no phase adjustment
-----------------
===> look at the watertrack figures:
figview.py cal/watertrk/*png
(b) cal/botmtrk/btcaluv.out (no bottom track)
- action: watertrack suggests a slight scale factor
might be applied after editing
- this is "close enough" to allow editing at this point
(4) view the data
dataviewer.py
(5)
- update the navigations steps and redo calibration -- get
the directory in line with modern tools and metadata
This will fail:
quick_adcp.py --steps2rerun calib --auto
#
# - NOTE: quick_adcp.py will complain about missing some information.
# At each error, include the information in the quick_adcp.py
# command line until it runs.
#
# ERROR solution: add to commands
#
#
# ERROR -- must select "datatype" "--datatype uhdas"
# ERROR -- must set "sonar" "--sonar os38nb"
# ERROR -- must set "yearbase" "--yearbase 2010"
# ERROR -- must set "ens_len" "--ens_len 300"
#
# Look in the file adcpdb/a_km.cnh to determine the number of
# seconds used for the averaging. We need the column "SI" for
# "sampling interval"
#
# ERROR -- must set "cruisename" "--cruisename km1001c"
#
# For post-processing, the cruise name is just used for titles
#
#
# ERROR -- must set "beamangle" "--beamangle 30"
#
# Usually, beamangles are:
# os38, os75, os150 30
# wh300, wh600, wh1200 20
# bb75, bb150, bb300 30
# nb150, nb300 30
#
#
#
# - now try again
#
quick_adcp.py --steps2rerun navsteps:calib --datatype uhdas --sonar os38nb --beamangle 30 --yearbase 2010 --ens_len 300 --cruisename km1001c --auto
Now make sure all the edits have been applied and look at the calibration
cd ..
quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto
tail -20 cal/watertrk/adcpcal.out
**watertrack**
------------
Number of edited points: 32 out of 35
amp = 1.0080 + -0.0027 (t - 30.3)
phase = -0.01 + 0.0589 (t - 30.3)
median mean std
amplitude 1.0060 1.0080 0.0085
phase 0.0775 -0.0058 0.4343
------------
(6) Edit out any bad points or profiles, or the data below the bottom using
dataviewer.py -e
(7) Apply final calibration
There are 32 points, enough for reasonable statistics
# A phase correction is not warranted because the
# mean and median are under 0.1 degree. If the phase
# values above had said X.YY, then we would include
# this in the quick_adcp.py command:
#
# --rotate_angle X.YY
#
# But a scale factor ("amplitude") should be applied:
quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.007 --auto
**watertrack**
------------
Number of edited points: 32 out of 36
amp = 1.0011 + -0.0028 (t - 30.3)
phase = -0.01 + 0.0589 (t - 30.3)
median mean std
amplitude 0.9990 1.0011 0.0084
phase 0.0810 -0.0081 0.4342
------------
Looks good
Look at the figures too
figview.py cal/watertrk/wtcal1.png
(6) make web plots
quick_web.py --interactive
- view with a browser, look at webpy/index.html
(7) extract data
quick_adcp.py --steps2rerun matfiles --auto
adcp_nc.py adcpdb contour/os38nb km1001c_demo os38nb --ship_name Kilo Moana
Submit to JASADCP (Joint Archive for Shipboard ADCP)
http://ilikai.soest.hawaii.edu/sadcp/
