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')
%
%notes : MITprof files will be used from the myenv.profiles directory
%            or [dirMat '/profiles/output/'] is myenv.profiles is missing

gcmfaces_global;

dirMat=[dirMat '/'];

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

listData=dir([dirData '*.nc']);

dirMat={'dirMat',dirMat};
dirData={'dirData',dirData};
listData={'listData',listData(:).name};

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;
fprintf('starting insitu_misfit\n'); clock
insitu_misfit(1,dirData,dirMat,listData);
fprintf('starting insitu_cost\n'); clock
insitu_cost(1,dirData,dirMat,listData);
fprintf('done with insitu_diags\n'); clock
else;
insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
end;

1	gforget	1.3	function []=insitu_diags(dirMat,doComp,dirTex,nameTex);
2	gforget	1.2	%object: driver for insitu_misfit and insitu_cost
3	gforget	1.3	%inputs: dirMat is the directory where diagnozed .mat files will be saved
4	gforget	1.2	% -> 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	gforget	1.3	% nameTex is the tex file name (default : 'myPlots')
8	gforget	1.5	%
9			%notes : MITprof files will be used from the myenv.profiles directory
10			% or [dirMat '/profiles/output/'] is myenv.profiles is missing
11
12			gcmfaces_global;
13	gforget	1.1
14	gforget	1.3	dirMat=[dirMat '/'];
15	gforget	1.2
16	gforget	1.5	if isfield(myenv,'profiles');
17			dirData=myenv.profiles;
18			else;
19			dirData=[dirMat '/profiles/output/'];
20			end;
21	gforget	1.4	while ~isdir(dirData)&doComp;
22	gforget	1.3	fprintf(['directory : ' dirData '\n'])
23			dirData=input(['does not exist. Specify directory of nc file : \n']);
24			end;
25
26	gforget	1.6	listData=dir([dirData '*.nc']);
27
28	gforget	1.3	dirMat={'dirMat',dirMat};
29			dirData={'dirData',dirData};
30	gforget	1.6	listData={'listData',listData(:).name};
31	gforget	1.2
32			if isempty(who('dirTex'));
33	gforget	1.3	addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
34	gforget	1.2	else;
35			if ~ischar(dirTex); error('mis-specified dirTex'); end;
36			if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path
37	gforget	1.3	addToTex={'addToTex',1}; dirTex={'dirTex',[dirTex '/']};
38			if isempty(who('nameTex')); nameTex='myPlots'; end;
39			nameTex={'nameTex',nameTex};
40	gforget	1.2	end;
41	gforget	1.1
42			if doComp;
43	gforget	1.6	fprintf('starting insitu_misfit\n'); clock
44			insitu_misfit(1,dirData,dirMat,listData);
45			fprintf('starting insitu_cost\n'); clock
46			insitu_cost(1,dirData,dirMat,listData);
47			fprintf('done with insitu_diags\n'); clock
48	gforget	1.1	else;
49	gforget	1.3	insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
50			insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
51	gforget	1.1	end;
52