LAST CHANGED    : 2010/04/05 03:36:15
CRUISE NAME(S)  : P6_leg1 (mv0911)
CRUISE DATES    : 2009/11/21 to 2010/01/02
SHIP NAME       : Melville
PORTS           : Brisbane, Australia to Papeete, Tahiti
CHIEF SCIENTIST : Ruth Curry
DATABASE NAME   : aship
DATA FILES      : mv2009_324_21341.raw to mv2010_002_79200.raw

STATUS          :  to do                           done
                  ------                       -----------
                  averaged                       [  x   ]
                  loaded                         [  x   ]
                  check heading correction       [  x   ]
                  calibration                    [  x   ]
                  edited                         [  x   ]
                  re-check heading correction    [  x   ]
                  check editing                  [  x   ]
                  figures                        [  x   ]

ready for JASADCP


INSTRUMENT      : os150

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
                  2009/11/21 06:00:51   on  300  70   8   5  4   0.00   0.00 0001
                  2009/11/22 07:49:48  off  300  70   8   5  4   0.00   0.00 0001

HEADING         :
     PRIMARY    : gyro
     CORRECTION : ashtech


POSITIONS       : gps (specify more detail if available)

CALIBRATION     : CALIBRATION     : original heading alignment: 0.0
                  additional rotation 0.41
                  final transducer angle estimate is 0.0 - (0.41) = -0.41
                         (original transducer angle) - (rotate_angle)
                  original scale factor 1
                  additional tweak of 1.007 was applied (a little high)

COMMENTS        : (1) heading:

                  long ashtech gap 330.05-340.2;
                  no ashtech data at all, 345.05-345.3
                  very long ashtech gap 351.95-354.05
                  very long ashtech gap 364.45 - 368

                  Fortunately, the heading correction appears to be locked in
                  at about 0.8 +-0.1, so we're not in bad shape

                  (2) : on-station

                  The transducer well evidently has the potential for
                  trapping air, so sometimes on station it simply cannot
                  find any water to ping into. About half of the
                  "On station" it often does not have any good velocities.

PROCESSOR       : J.Hummon


#===========================================

reprocessing from scratch to regain data in regions of low scattering.

# rsync config from the data directory to here;
# modify to have useful paths

(1) # make processing directory

adcptree.py MV0911_os150nb  --datatype uhdas --configbase MV0911 --configpath config

(2) # set up control file

cd $procdir

cat << EOF > q_py.cnt

--yearbase    2009
--pingtype    nb
--cfgbase     MV1001
--cruisename  MV0911
--instclass   os
--instname    os150
--dbname      aship
--datatype    uhdas
--datadir     path_to_data/MV0911/raw/os150
--cfgpath     $basedir/config
--ens_len     300
--ping_headcorr
--hcorr_inst ashtech
--use_topo4edit
--delblks     1
--auto

EOF


(3)

quick_adcp.py --cntfile q_py.cnt --mplplots all --npycalc  all


(4) check rotation

Ashtech gaps -- fix

        cd cal/rotate
        # copy file in, edit for this case
        python fix_heading.py # makes newhcorr.ang, newhcorr.asc
        cp Hcorrplot_script.py Hcorrplot_newhcorr.py
        # edit, run, look at pictures -- all OK?
        rotate unrotate.tmp
        rotate rotate1.tmp


        cd ../..

        quick_adcp.py --yearbase 2009 --use_refsm --steps2rerun \
              calib:matfiles --auto  --mplplots all --npycalc  all


** watertrack **

Number of edited points: 188 out of 201
   amp   = 0.0001  + 1.0072 (t - 342.7)
   phase =   0.00  + 0.4229 (t - 342.7)
            median     mean      std
amplitude   1.0070   1.0072   0.0058
phase       0.4120   0.4229   0.3301



===>  tweak calibration before editing

 quick_adcp.py --yearbase 2009 --use_refsm --steps2rerun \
   rotate:navsteps:calib:matfiles --auto  --mplplots all --npycalc  all \
   --rotate_angle -.41 --rotate_amplitude 1.007


** watertrack **

Number of edited points: 188 out of 201
   amp   = 0.0001  + 1.0002 (t - 342.7)
   phase =   0.00  + 0.8348 (t - 342.7)
            median     mean      std
amplitude   1.0000   1.0002   0.0058
phase       0.8185   0.8348   0.3345

OOPS -- got the sign wrong on the phase -- do it again:


 quick_adcp.py --yearbase 2009 --use_refsm --steps2rerun \
   rotate:navsteps:calib:matfiles --auto  --mplplots all --npycalc  all \
   --rotate_angle .83


** watertrack **

Number of edited points: 188 out of 201
   amp   = 0.0001  + 1.0002 (t - 342.7)
   phase =   0.00  + 0.0050 (t - 342.7)
            median     mean      std
amplitude   1.0000   1.0002   0.0058
phase      -0.0110   0.0050   0.3341


That's better; now edit


(5) edit

 quick_adcp.py --yearbase 2009 --use_refsm --steps2rerun \
   apply_edit:navsteps:calib:matfiles --auto  --mplplots all --npycalc  all


 ** watertrack **



review editing

=================================

(7) web plots


mkdir webpy
cp ../MV0911_os75nb/webpy/sectinfo.txt webpy/
quick_web.py --instping os150nb --cruiseid MV0911_os150nb --redo