7.7.37. quick_postproc_help.txtΒΆ

===========================================
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/