LAST CHANGED : Tue Dec 28 06:28:21 HST 2004 CRUISE NAME(S) : ka0407, gp804_ens CRUISE DATES : 2004/11/23 to 2004/11/23 PORTS : Kwajalein, RMI to Honolulu, HI SHIP NAME : Ka`imimoana CHIEF SCIENTIST : Ben Moore DB_NAME : a804 DATA FILES : ADCP087*.ENS STATUS : 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 = -.145 : final transducer angle = 58.73 - (-.145) = 58.8750 COMMENTS : had to redo processing from ens files; apply posmv correction manually. bad weather rendered several days useless unless on station, with problems including - very low PG - underway bias gps feed is aft of the transducer; (compensated) averages have teh wrong depth: use "depthcng" to decrease depths by 5 m and change transducer depth to 5m PROCESSED BY : Jules Hummon University of Hawaii ------------------------------------------------------- Times: 4 hours to average+edit the data few minutes for the rest of the codas processing several minutes for pmv.rbin (from N2R) 1 hour for hdg.rbin (from ens) ### q_py.cnt is --yearbase 2004 # for decimal day conversion --dbname a804 # database name; in adcpdb --ducer_depth 5 # ENS needs transducer depth --ens_halign 58.73 ## 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/gp804 ### end of q_py.cnt quick_adcp.py --auto --cntfile q_py.cnt ------- get first calibration estimate (from gyro) ------------ ----------------------------- now get posmv correction (1a) identify the data directory: /mnt/ulili2/vmdas_data/osdata/kaimi075/gp804 (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/ulili2/vmdas_data/osdata/kaimi075/gp804/*.N2R (3) fire up matlab, get the headings %%% matlab datadir = '/mnt/ulili2/vmdas_data/osdata/kaimi075/gp804'; filelist = dirs(fullfile(datadir, '*.ENS'),'fullfile',1); for ifile = 1:length(filelist) ensfile = filelist(ifile).name ensgyro2rbin(os, ensfile); end %%% (4) go back to cal/rotate, in matlab, edit and run mk_ensgyrodh: cd .. % rotate by ang file rotate rotate1.tmp cd ../.. quick_adcp.py --use_refsm --yearbase 2004 \ --steps2rerun apply_edit:navsteps:calib:matfiles --auto cat cal/watertrk/*out Number of edited points: 80 out of 97 amp = 1.0031 + -0.0003 (t - 335.5) phase = -0.14 + -0.0312 (t - 335.5) median mean std amplitude 1.0030 1.0031 0.0090 phase -0.1110 -0.1449 0.5699 (6b) cal/watertrk indicates additional rotation needed: quick_adcp.py --use_refsm --yearbase 2004 --instname os75 \ --steps2rerun rotate:navsteps:calib:matfiles --auto \ --rotate_angle -.145 -------------- (7) check transducer offset from gps: best_xducerxy('yearbase',2004,'xoffsets',[-5],'yoffsets',[0:5:70],'fixfile','a804.gps') best_xducerxy('yearbase',2004,'xoffsets',[-10:2:10],'yoffsets',[50],'fixfile','a804.gps') NOTE:run this to get the correct agt file: xducerxy_core('dx',0,'dy',50','yearbase',2004,'fixfile','a804.gps') THEN edit asetup.m to point to this file. (7) checked editing, fixed depth ("depthcng", subtract 7, fix xducer depth at 5) (8) quick_adcp.py --use_refsm --yearbase 2004 --instname os75 \ --steps2rerun apply_edit:navsteps:calib:matfiles --auto \ --xducer_dx 0 --xducer_dy 50 ** watertrk ** Number of edited points: 82 out of 97 amp = 1.0018 + 0.0001 (t - 335.7) phase = 0.04 + -0.0349 (t - 335.7) median mean std amplitude 1.0010 1.0018 0.0071 phase 0.0175 0.0355 0.5053