profilesMatlabProcessing/ecco_v4/insitu_diags.m

function []=insitu_diags(dirMat,doComp,dirTex,nameTex);
%object:     driver for insitu_misfit and insitu_cost
%inputs:     dirMat is the directory where diagnozed .mat files will be saved
%                     -> set it to '' to use the default [dirModel 'mat/']
%            doComp states whether to compute (1) or display (0)
%            dirTex is the directory where tex and fig files will be created
%            nameTex is the tex file name (default : 'myPlots')

dirMat=[dirMat '/'];

dirData=[dirMat '/profiles/output/'];
while ~isdir(dirData)&doComp;
  fprintf(['directory : ' dirData '\n'])
  dirData=input(['does not exist. Specify directory of nc file : \n']);
end;

dirMat={'dirMat',dirMat};
dirData={'dirData',dirData};

if isempty(who('dirTex'));
  addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
else;
  if ~ischar(dirTex); error('mis-specified dirTex'); end;
  if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path
  addToTex={'addToTex',1}; dirTex={'dirTex',[dirTex '/']};
  if isempty(who('nameTex')); nameTex='myPlots'; end;
  nameTex={'nameTex',nameTex};
end;

if doComp;
insitu_misfit(1,dirData,dirMat);
insitu_cost(1,dirData,dirMat);
else;
insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
end;

1	function []=insitu_diags(dirMat,doComp,dirTex,nameTex);
2	%object: driver for insitu_misfit and insitu_cost
3	%inputs: dirMat is the directory where diagnozed .mat files will be saved
4	% -> set it to '' to use the default [dirModel 'mat/']
5	% doComp states whether to compute (1) or display (0)
6	% dirTex is the directory where tex and fig files will be created
7	% nameTex is the tex file name (default : 'myPlots')
8
9	dirMat=[dirMat '/'];
10
11	dirData=[dirMat '/profiles/output/'];
12	while ~isdir(dirData)&doComp;
13	fprintf(['directory : ' dirData '\n'])
14	dirData=input(['does not exist. Specify directory of nc file : \n']);
15	end;
16
17	dirMat={'dirMat',dirMat};
18	dirData={'dirData',dirData};
19
20	if isempty(who('dirTex'));
21	addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
22	else;
23	if ~ischar(dirTex); error('mis-specified dirTex'); end;
24	if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path
25	addToTex={'addToTex',1}; dirTex={'dirTex',[dirTex '/']};
26	if isempty(who('nameTex')); nameTex='myPlots'; end;
27	nameTex={'nameTex',nameTex};
28	end;
29
30	if doComp;
31	insitu_misfit(1,dirData,dirMat);
32	insitu_cost(1,dirData,dirMat);
33	else;
34	insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
35	insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
36	end;
37