SHIP NAME       : Thomas G. Thompson 
LAST CHANGED    : Fri Apr  7 16:16:36 HST 2006
CRUISE          : P16N (TN191_leg1)
CRUISE DATES    : 2006/02/13 to 2006/03/03
PORTS           : Papeete, Tahiti to Honolulu, HI
CHIEF SCIENTIST : Chris Sabine
DB_NAME         : a_tt
DATA FILES      : tt2006_042_82400.raw to tt2006_066_79205.raw
STATUS          : processed, edited, calibrated (klooged); 

INSTRUMENTs     : OS75

ACQUISITION
     PROGRAM    : UHDAS

     PING TYPE  : narrowband

      PARAMS    : blank=            8  (BK) (see adcpdb/\*cnh after first pass)
                  pulse =          16  (PL)
                  bin length=      16  (BL)
                  number of bins=  70  (NB)
HEADING
     PRIMARY    : gyro
     CORRECTION : posmv

POSITIONS       : gps

CALIBRATION     : original heading alignment: 58.67, amp = 1.02
                : final phase calibration 58.67 - (-.31)
                : further amplitude correction necessary; for now, use
                :       soundspeed to get finer-grained control.
                : Final amplitude is corrected as:
                :       amp = 1.0318 - .0006*temp
                : see calibration section, below


COMMENTS        : data collected at sea; processing directory cleaned
                :     up here for final product

PROCESSOR       : Jules Hummon

-------------------------

(1) Clean up in general:

 (i)

    Files to remove               why
    ---------------            -------------
    adcpdb/*blk                remove the database; we will load it again
    cal/rotate/rotate.log      A new database requires a new "rotate.log"
                               file, or quick_adcp.py will be confused
    cal/watertrk/*.out         clean calibrations are coming; not incremental
    cal/botmtrk/*.out          clean calibrations are coming; not incremental
    edit/a*.asc                clean editing is coming, not just using defaults
    edit/a*.log                clean editing is coming, make a clean log file

    ** do not remove any other files!!

 (ii)

    - in load, edit ldcodas.tmp so it lists all cmd files
    - run "ldcodas ldcodas.tmp"

 (iii)

    - in nav/ run "putnav putnav.tmp"

  (iv)

   - in processing directory, run quick_adcp with nav steps:

   quick_adcp.py --yearbase 2005 --use_refsm --steps2rerun  navsteps --auto

  (v) mark as bad, everything below 50 percent good
      in edit/, run this:

      setflags setflags.tmp

(2) Make nice overview plots of heading correction, and look at quality

    (a) create figures and statistics for the whole cruise
        (the original directory came from underway processing and so
        has some fossils)

    - change directories to cal/rotate
    - copy the diagnostic programs to more descriptive names:

         copy from           to this name
         ---------          ---------------
       plot_hcorrstats_tmp.m      plot_hcorrstats_all.m
       print_hcorrstats_tmp.m     print_hcorrstats_all.m

    - edit the two new files:
        (i)  (plot_hcorrstats_all.m only:)
             remove the first three lines of code starting "set(0,..."

        (ii) change 'plot_ddrange' and 'print_ddrange' from 1 to 360

    - in matlab, run these two programs

    (b) load ens_hcorr.asc into matlab and look at the fields.  Are
        there gaps that ned to be interpolated over? spikes?
        (answer = 'yes' for this dataset).  See
        cal/manual_rotate/improve_dh.m for an example of a cleaner
        correction)

        Steps to follow:
           (i)  in cal/manual_rotate, create a better heading  correction file
          (ii)  in cal/rotate, remove the heaading corrections that are embedded
                   in the ensembles 'rotate unrotate.tmp'
                   (for the control file, use time range from adcpdb/*.tr)
         (iii)  in cal/rotate, apply the new heading correction 'rotate rotate1.tmp'
                   (for the control file, use time range from adcpdb/*.tr)
         (iii)  rerun nav steps:
           quick_adcp.py --yearbase 2005 --use_refsm --steps2rerun  navsteps --auto


(3) calibration:

watertrack:



Number of edited points:  90 out of  90
   amp   = 0.9956  + 0.0004 (t -  53.0)
   phase =  -0.34  + -0.0082 (t -  53.0)
            median     mean      std
amplitude   0.9960   0.9956   0.0056
phase      -0.3430  -0.3397   0.2294




%% newvel = oldvel*amp = oldvel*(newsoundspeed/oldsoundspeed)
%% ==> amp = newsoundspeed/oldsoundspeed
%% ==> newsoundspeed = amp*oldsoundspeed;


Apply time-dependent amplitude correction
(this actually comes in through fix_temp
fix_temp fix_temp_ss.tmp
cd ../..
  quick_adcp.py --yearbase 2006 --use_refsm --steps2rerun  \
         navsteps:calib:matfiles --auto

%%  watertrack calib %%%


Number of edited points:  90 out of  90
   amp   = 0.9998  + 0.0002 (t -  53.0)
   phase =  -0.34  + -0.0083 (t -  53.0)
            median     mean      std
amplitude   1.0000   0.9998   0.0055
phase      -0.3385  -0.3377   0.2281



%%% rotate by -.31 deg:

quick_adcp.py --yearbase 2006 --use_refsm --steps2rerun \
  rotate:navsteps:calib:matfiles --auto \
  --rotate_angle -.31

%%  watertrack calib %%%


Number of edited points:  90 out of  90
   amp   = 0.9998  + 0.0002 (t -  53.0)
   phase =  -0.03  + -0.0082 (t -  53.0)
            median     mean      std
amplitude   0.9995   0.9998   0.0056
phase      -0.0250  -0.0272   0.2278


END
Final phase calibration  58.67 - (-.31) from watertrack
Final amplitude calibration roughly 1.02, but temperature -dependent