profilesMatlabProcessing/profiles_misc/MITprof_struct.m

function [MITprof]=MITprof_struct(nProf,dataset,varargin);
%   [MITprof]=MITprof_struct;
%   create an empty struct variable using the format MITprof
%       nProf: number of profiles
%       prof_depth: list of depth levels
%       list_var: specify name of created fields
%
%   list of fields always created:
%       prof_depth      [nDepth x 1]
%       prof_descr      {nProf x 1}
%       prof_date       [nProf x 1]
%       prof_YYYYMMDD   [nProf x 1]
%       prof_HHMMSS     [nProf x 1]
%       prof_lon        [nProf x 1]
%       prof_lat        [nProf x 1]
%       prof_basin      [nProf x 1]
%       prof_point      [nProf x 1]
%       np, nr, list_descr, nd
%
%   list of fields created by default, if list_var not specified
%       prof_T, prof_Tflag, prof_Terr, prof_Tweight, prof_Testim
%       prof_S, prof_Sflag, prof_Serr, prof_Sweight, prof_Sestim
%           array size: [nProf x nLev]

prof_depth=dataset.z_std;

if isfield(dataset,'var_tmp');
    var_tmp=dataset.var_tmp;
    var_out=dataset.var_out;
    jj=find(~strcmp(var_tmp,'z'));
    var_tmp={var_tmp{jj}};
    var_out={var_out{jj}};
    list_vars={};
    for ii=1:length(var_tmp);
        list_vars={list_vars{:},['prof_' var_out{ii} ]};
        list_vars={list_vars{:},['prof_' var_out{ii} 'weight']};
        list_vars={list_vars{:},['prof_' var_out{ii} 'estim']};
        list_vars={list_vars{:},['prof_' var_out{ii} 'err']};
        list_vars={list_vars{:},['prof_' var_out{ii} 'flag']};
    end;
else;
    list_vars={'prof_T','prof_Tweight','prof_Testim','prof_Terr','prof_Tflag',...
        'prof_S','prof_Sweight','prof_Sestim','prof_Serr','prof_Sflag'};
    if ~strcmp(dataset.coord,'depth');
        list_vars={list_vars{:},'prof_D','prof_Destim'};
    end;
end;
if nargin>2,
    list_vars=varargin{1};
end


MITprof=[];

nLev=length(prof_depth);
prof_depth=reshape(prof_depth,length(prof_depth),1);
MITprof.prof_depth=prof_depth;

MITprof.prof_date=zeros(nProf,1);
MITprof.prof_YYYYMMDD=zeros(nProf,1);
MITprof.prof_HHMMSS=zeros(nProf,1);
MITprof.prof_lon=zeros(nProf,1);
MITprof.prof_lat=zeros(nProf,1);
MITprof.prof_basin=zeros(nProf,1);
MITprof.prof_point=zeros(nProf,1);

MITprof.prof_descr=cell(nProf,1);

for ii=1:length(list_vars),
    MITprof=setfield(MITprof,list_vars{ii},zeros(nProf,nLev));
end

MITprof.np=nProf;
MITprof.nr=nLev;
MITprof.nd=0;
MITprof.list_descr={};


1	gforget	1.4	function [MITprof]=MITprof_struct(nProf,dataset,varargin);
2	gforget	1.1	% [MITprof]=MITprof_struct;
3			% create an empty struct variable using the format MITprof
4			% nProf: number of profiles
5	roquet	1.3	% prof_depth: list of depth levels
6			% list_var: specify name of created fields
7	gforget	1.5	%
8	roquet	1.3	% list of fields always created:
9			% prof_depth [nDepth x 1]
10			% prof_descr {nProf x 1}
11			% prof_date [nProf x 1]
12			% prof_YYYYMMDD [nProf x 1]
13			% prof_HHMMSS [nProf x 1]
14			% prof_lon [nProf x 1]
15			% prof_lat [nProf x 1]
16			% prof_basin [nProf x 1]
17			% prof_point [nProf x 1]
18			% np, nr, list_descr, nd
19			%
20			% list of fields created by default, if list_var not specified
21			% prof_T, prof_Tflag, prof_Terr, prof_Tweight, prof_Testim
22			% prof_S, prof_Sflag, prof_Serr, prof_Sweight, prof_Sestim
23			% array size: [nProf x nLev]
24
25	gforget	1.4	prof_depth=dataset.z_std;
26	roquet	1.3
27	gforget	1.5	if isfield(dataset,'var_tmp');
28			var_tmp=dataset.var_tmp;
29			var_out=dataset.var_out;
30			jj=find(~strcmp(var_tmp,'z'));
31			var_tmp={var_tmp{jj}};
32			var_out={var_out{jj}};
33			list_vars={};
34			for ii=1:length(var_tmp);
35			list_vars={list_vars{:},['prof_' var_out{ii} ]};
36			list_vars={list_vars{:},['prof_' var_out{ii} 'weight']};
37			list_vars={list_vars{:},['prof_' var_out{ii} 'estim']};
38			list_vars={list_vars{:},['prof_' var_out{ii} 'err']};
39			list_vars={list_vars{:},['prof_' var_out{ii} 'flag']};
40			end;
41			else;
42			list_vars={'prof_T','prof_Tweight','prof_Testim','prof_Terr','prof_Tflag',...
43	roquet	1.3	'prof_S','prof_Sweight','prof_Sestim','prof_Serr','prof_Sflag'};
44	gforget	1.5	if ~strcmp(dataset.coord,'depth');
45			list_vars={list_vars{:},'prof_D','prof_Destim'};
46			end;
47	gforget	1.4	end;
48	roquet	1.3	if nargin>2,
49			list_vars=varargin{1};
50			end
51	roquet	1.2
52	gforget	1.1
53			MITprof=[];
54	roquet	1.3
55			nLev=length(prof_depth);
56			prof_depth=reshape(prof_depth,length(prof_depth),1);
57			MITprof.prof_depth=prof_depth;
58
59	gforget	1.1	MITprof.prof_date=zeros(nProf,1);
60			MITprof.prof_YYYYMMDD=zeros(nProf,1);
61			MITprof.prof_HHMMSS=zeros(nProf,1);
62			MITprof.prof_lon=zeros(nProf,1);
63			MITprof.prof_lat=zeros(nProf,1);
64			MITprof.prof_basin=zeros(nProf,1);
65			MITprof.prof_point=zeros(nProf,1);
66	roquet	1.3
67	gforget	1.1	MITprof.prof_descr=cell(nProf,1);
68	roquet	1.3
69			for ii=1:length(list_vars),
70			MITprof=setfield(MITprof,list_vars{ii},zeros(nProf,nLev));
71			end
72
73	gforget	1.1	MITprof.np=nProf;
74			MITprof.nr=nLev;
75			MITprof.nd=0;
76			MITprof.list_descr={};
77
78
79