% script newazmetatsts.m
% metaSonarExample - an example script for processing a single ADV data file.
% This program sets up metadata and runs the sonar processing appropriate
% to the kind of sonar employwd.
%
% you should run this script from the directory in which the data
% files reside
%

%%% START USGS BOILERPLATE -------------
% These programs are intened for us in adjusting metadata terms in
% netCDF files from the USGS CMGP Oceanographic time-series data
% archive
%
% Program written Matlab 7.6.0.342 (R2008a)
% Program ran on Linux PC with RHEL4 and on Windows XP PC.
%
% "Although this program has been used by the USGS, no warranty,
% expressed or implied, is made by the USGS or the United States
% Government as to the accuracy and functioning of the program
% and related program material nor shall the fact of distribution
% constitute any such warranty, and no responsibility is assumed
% by the USGS in connection therewith."
%%% END USGS BOILERPLATE --------------

% fname must contain a file of Global attributes
fname='glob_att_unh.txt';
userMeta=read_globalatts(fname);  % provides values to insert below

% Written by Marinna Martini for the U.S. Geological Survey
% Coastal & Marine Program Woods Hole Field Center, Woods Hole, MA
% http://woodshole.er.usgs.gov/ Please report bugs to mmartini@usgs.gov

% add more EPIC fields
userMeta.INST_TYPE='Imagenex Model 881a Profiling sonar with azimuth drive';
userMeta.history='Sonar data downloaded from instrument';
userMeta.DATA_TYPE='TIME';
userMeta.DELTA_T='180 min';
userMeta.metadata_author='Ellyn Montgomery, USGS';
% these are specific to the sonar
% this one is key- says which kind of sonar to process, and may affect which 
% additional fields are available choose from: 'fan' | 'pen' | 'azm'
userMeta.SonartoAnimate='azm';
% these are configuration items not found in the instrument file - most of the 
%  instrument set-up is obtained from header or switches 
userMeta.Height=1.055;     %  mab at deployment
userMeta.FirstSonarDay='09-Dec-2009';
userMeta.LastSonarDay='10-Dec-2009';    % make sure there's more than 1 day
userMeta.RootDataDir='C:\home\data\sonar_tst\unh_tank\dives';
userMeta.sweeps=1;  % there's nothing in the header that has this info
userMeta.note='blue and green legs up- run with new targets'
%
% the ones below come into play if processing is to take place
% the adcp related ones help determine how much to rotate the fan sonar image
procMeta.SonartoAnimate='azm';
procMeta.tidx=1;
procMeta.plottype='scat_frm_img';
% pencil tilt is used to flatten the profile of the pencil sweeps
procMeta.Pencil_tilt=0.0;
procMeta.rot2compass=0.0;
% dxy is the grid size images are interpolated onto
% .01 maintains resolution, but takes a long time and makes large files
% .05 is faster and makes smaller files, so is useful for quick-looks to verify rotaion
procMeta.dxy=0.005;
%
% save the output information
diary(sprintf('run%s',datestr(now,30))) 
outFile='dummy';   % not used for asimuth data- otherwise use a real name

% create the raw cdf- repeat once for each type
if 0
    mkrawcdf(userMeta, outFile);
end
%
% make processed images
if 1
    if strcmp(userMeta.SonartoAnimate,'azm')
        plotrange_cdf(['az2009-12-10_raw.cdf'],procMeta.tidx,procMeta)
    else
        disp ('could not interpret command- try again')
    end
end
diary off