===========================================
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
# pn45 20
# 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/