profilesMatlabProcessing/profiles_IO_v2/MITprof_read.m

%function:      MITprof_read
%object:        read netcdf data files in the "MIT format"
%author:        Gael Forget (gforget@mit.edu)
%date:          june 21th, 2006
%
%usage:    [MITprof]=MITprof_read(fileIn); 
%               ---> loads full data set
%          [MITprof]=MITprof_read(fileIn,list_vars); 
%               ---> loads only the files listed in list_vars cell
%               array (e.g. list_vars={'prof_T','prof_Tweight'}) 
%               plus the one dimensional information (prof_lon etc.)
%
%inputs:        fileIn          data file name
%               list_vars       variables list (optional)
%
%outputs:       MITprof structure containing the various fields/vectors

function [MITprof]=MITprof_read(fileIn,varargin);

if nargin>1; list_vars=varargin{1}; 
else; list_vars={'prof_T','prof_Tweight','prof_Tmask','prof_Testim','prof_Tflag',...
              'prof_S','prof_Sweight','prof_Smask','prof_Sestim','prof_Sflag',...
              'prof_U','prof_Uweight','prof_Umask',...
              'prof_V','prof_Vweight','prof_Vmask',...
              'prof_ptr','prof_ptrweight','prof_ptrmask',...
              'prof_ssh','prof_sshweight','prof_sshmask'};
end;

list_vars_plus=[{'prof_depth','prof_date','prof_YYYYMMDD','prof_HHMMSS',...
        'prof_lon','prof_lat','prof_basin','prof_point','prof_descr'}...
         list_vars];

%get directory name:
tmp1=strfind(fileIn,'/'); 
if ~isempty(tmp1); dirIn=fileIn(1:tmp1(end)); else; dirIn='./'; end;
%check that file exists:
tmp1=dir(fileIn); 
if isempty(tmp1); tmp1=dir([fileIn '.nc']); end;
if isempty(tmp1); error([fileIn ' file not found']); end;
%open file:
nc=netcdf([dirIn tmp1.name],'read');

for ii=1:length(list_vars_plus);
        tmp1=nc{list_vars_plus{ii}}(:);
        if ii==1; MITprof=struct('prof_depth',tmp1); 
        elseif ~isempty(tmp1); eval(['MITprof.' list_vars_plus{ii} '=tmp1;']);
        end;
end;

close(nc);

1	gforget	1.1	%function: MITprof_read
2			%object: read netcdf data files in the "MIT format"
3			%author: Gael Forget (gforget@mit.edu)
4			%date: june 21th, 2006
5			%
6	gforget	1.3	%usage: [MITprof]=MITprof_read(fileIn);
7	gforget	1.1	% ---> loads full data set
8	gforget	1.3	% [MITprof]=MITprof_read(fileIn,list_vars);
9	gforget	1.1	% ---> loads only the files listed in list_vars cell
10			% array (e.g. list_vars={'prof_T','prof_Tweight'})
11			% plus the one dimensional information (prof_lon etc.)
12			%
13	gforget	1.3	%inputs: fileIn data file name
14			% list_vars variables list (optional)
15	gforget	1.1	%
16			%outputs: MITprof structure containing the various fields/vectors
17
18	gforget	1.3	function [MITprof]=MITprof_read(fileIn,varargin);
19	gforget	1.1
20	gforget	1.3	if nargin>1; list_vars=varargin{1};
21	gforget	1.2	else; list_vars={'prof_T','prof_Tweight','prof_Tmask','prof_Testim','prof_Tflag',...
22			'prof_S','prof_Sweight','prof_Smask','prof_Sestim','prof_Sflag',...
23	gforget	1.1	'prof_U','prof_Uweight','prof_Umask',...
24			'prof_V','prof_Vweight','prof_Vmask',...
25			'prof_ptr','prof_ptrweight','prof_ptrmask',...
26			'prof_ssh','prof_sshweight','prof_sshmask'};
27			end;
28
29			list_vars_plus=[{'prof_depth','prof_date','prof_YYYYMMDD','prof_HHMMSS',...
30			'prof_lon','prof_lat','prof_basin','prof_point','prof_descr'}...
31			list_vars];
32
33	gforget	1.3	%get directory name:
34			tmp1=strfind(fileIn,'/');
35			if ~isempty(tmp1); dirIn=fileIn(1:tmp1(end)); else; dirIn='./'; end;
36			%check that file exists:
37			tmp1=dir(fileIn);
38			if isempty(tmp1); tmp1=dir([fileIn '.nc']); end;
39	gforget	1.1	if isempty(tmp1); error([fileIn ' file not found']); end;
40	gforget	1.3	%open file:
41	gforget	1.1	nc=netcdf([dirIn tmp1.name],'read');
42
43			for ii=1:length(list_vars_plus);
44			tmp1=nc{list_vars_plus{ii}}(:);
45			if ii==1; MITprof=struct('prof_depth',tmp1);
46			elseif ~isempty(tmp1); eval(['MITprof.' list_vars_plus{ii} '=tmp1;']);
47			end;
48			end;
49
50			close(nc);
51