% script do_cvt_seabird 
% do_cvt_seabird - script to convert ascii seabird data to netcdf
% 
% modify the in_info structure and th gatt_name as appropriate for
% the instrument to be processed
%   The tinyurl that describes this process more fully is: 
%    http://tinyurl.com/y2cr4q7.


%%% START USGS BOILERPLATE -------------% There is no published documentation for these programs- use the 
%  internal comments and help function.
% Programs written in Matlab- as early as 7.2.0.232 (R2006a) and 
% most recent modifications in  v7.8.0 (R2009a)
% Programs ran on PC with Windows XP Professional OS, and RHEL4 linux.
%
% "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 --------------

in_info.instno='MC-0727';    % Needs MC or SC prefix
in_info.type='mc';
in_info.hab=3.46;           % meters
in_info.pos_order=1;        % from filename
in_info.ascfile	='/home/data/proc/ctd/8371mc/8371mc.asc';
gatt_nm='/home/data/proc/ctd/globatt_837.txt';
[nnam,vi]=ipf_seabird(gatt_nm,in_info);
% now convert to epic
asc2epic(vi,nnam)

% add salinity using em_dolly
%  You need to provide a reference pressure to do this, so have it ready
%of=em_dolly('8371mc/8371mc.cdf','8371mc/8371mc1.cdf','salinity','y')

% add density
%of=em_dolly('8371mc/8371mc1.cdf','8371mc/8371mc2.cdf','sigma_theta','y')

% and finally truncate
%trunc_cdf('8371mc2.cdf','8371mc-trm.cdf',[15 18120])