LAST CHANGED    : Wed Sep  7 16:54:15 HST 2005
CRUISE NAME(S)  : ka0401, gp104_ens
CRUISE DATES    : 2004/03/26  to 2004/04/22  
SHIP NAME       : Ka`imimoana
PORTS           : San Diego, CA to Manzanillo, Mexico
CHIEF SCIENTIST : Dave Zimmerman
DB_NAME         : a104
DATA FILES      : ADCP087*.ENS
STATUS          : averaged, calibrated , edited, done.  ready for JASADCP

INSTRUMENT      : OS75 instrument #1

ACQUISITION
     PROGRAM    : VmDAS
     PING TYPE  : narrowband
      PARAMS    : blank=            8m 
                  pulse =           16 
                  bin length=       16 
                  number of bins=   60 
HEADING
     PRIMARY    : gyro
     CORRECTION : posmv 

POSITIONS       : ?? primary gps, whatever that is

CALIBRATION     : phase calibration = 2.4
                : final transducer angle = 56.5 - (-2.4) = 58.9 relative to posmv

COMMENTS        : had to redo processing from ens files; apply 
                          posmv correction manually.
                  lots of bad serial feeds ("bad navigation" in ENS files)
                  only $PRDID messages; time series had junms (see heading_correction.txt)
                  often computer time was reset to local time
PROCESSED BY    : Jules Hummon, University of Hawaii

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

  ### q_py.cnt would be:
   --yearbase 2004                       # for decimal day conversion
   --dbname a104                         # database name; in adcpdb
   --ducer_depth 5                       # ENS needs transducer depth
   --ens_halign  56.5                    ## h_align for gyro
   --use_refsm                           # could use --use_smoothr
   --instclass os                        # instrument type
   --datatype ens                        # datafile type
   --pingtype  nb                        # specify nb or bb
   --datafile_glob *.ENS                 # filename extension
   --tktail
   --datadir /mnt/ulili2/vmdas_data/osdata/kaimi075/gp104
   ### end of q_py.cnt


quick_adcp.py --auto --cntfile q_py.cnt


(6b) cal/watertrk indicates additional rotation needed:

Number of edited points:  67 out of 111
   amp   = 1.0029  + 0.0000 (t - 101.0)
   phase =  -3.05  + -0.0945 (t - 101.0)
            median     mean      std
amplitude   1.0030   1.0029   0.0099
phase      -3.1930  -3.0493   1.2725


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

now get posmv correction
see heading_correction.txt.  ugh.  


(1a) identify the data directory:
/mnt/removable/vmdas_data_archive/osdata/kaimi075/gp104

(1b) set the stage:
cd cal/rotate
mkdir rbin
cd rbin

(2) make the rbin files from N2R 
serasc2bin.py  -m pmv -c last -y 2004 -v -t vmdas \
    /mnt/removable/vmdas_data_archive/osdata/kaimi075/gp104/*.N2R

(3) fire up matlab, get the headings

%%%
matlab
datadir = '/mnt/removable/vmdas_data_archive/osdata/kaimi075/gp104';
filelist = dirs(fullfile(datadir, '*.ENS'),'fullfile',1);
for ifile = 1:length(filelist)
  ensfile = filelist(ifile).name
  ensgyro2rbin(os, ensfile);
end
%%% MANY BAD navigations



apply the correction: 
    cal/rotate: " rotate rotate1.tmp"
    quick_adcp.py --use_refsm --yearbase 2004 --instname os75 \
       --steps2rerun navsteps:calib --auto 

cat cal/watertrk/*out


Number of edited points:  93 out of 112
   amp   = 1.0027  + 0.0000 (t - 100.7)
   phase =  -2.29  + 0.0474 (t - 100.7)
            median     mean      std
amplitude   1.0030   1.0027   0.0099
phase      -2.3340  -2.2894   0.7263


    quick_adcp.py --use_refsm --yearbase 2004 --instname os75 \
       --steps2rerun rotate:navsteps:calib:matfiles --auto \
       --rotate_angle -2.3

cat cal/watertrk/*out


Number of edited points:  92 out of 112
   amp   = 1.0024  + -0.0000 (t - 100.6)
   phase =  -0.01  + 0.0431 (t - 100.6)
            median     mean      std
amplitude   1.0030   1.0024   0.0094
------------------



(7) check transducer offset from gps:

best_xducerxy('yearbase',2004,'xoffsets',[-5],'yoffsets',[0:5:70],...
           'enl',6000000,'fixfile','a0401.gps')
best_xducerxy('yearbase',2004,'xoffsets',[-10:2:10],'yoffsets',[30],...
            'enl',6000000,'fixfile','a0401.gps')


NOTE:run this to get the correct agt file:
 xducerxy_core('dx',6,'dy',30','yearbase',2004,'fixfile','a0401.gps','enl',6000000)
THEN edit asetup.m to point to this file.


(8) gps is not directly over transducer?? (from gp0401b_ens)


    quick_adcp.py --use_refsm --yearbase 2004 --instname os75 \
       --steps2rerun apply_edit:navsteps:calib:matfiles \
       --xducer_dx 6 --xducer_dy 30 --auto


Number of edited points:  80 out of  95
   amp   = 1.0006  + 0.0000 (t - 100.5)
   phase =  -0.11  + 0.0257 (t - 100.5)
            median     mean      std
amplitude   1.0010   1.0006   0.0063
phase      -0.1630  -0.1061   0.5542


    quick_adcp.py --use_refsm --yearbase 2004 --instname os75 \
       --steps2rerun rotate:navsteps:calib:matfiles --auto \
       --xducer_dx 6 --xducer_dy 30 --auto --rotate_angle -.1


----------