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	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	%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
14	dirMat=[dirMat '/'];
15
16	if isfield(myenv,'profiles');
17	dirData=myenv.profiles;
18	else;
19	dirData=[dirMat '/profiles/output/'];
20	end;
21	while ~isdir(dirData)&doComp;
22	fprintf(['directory : ' dirData '\n'])
23	dirData=input(['does not exist. Specify directory of nc file : \n']);
24	end;
25
26	listData=dir([dirData '*.nc']);
27
28	dirMat={'dirMat',dirMat};
29	dirData={'dirData',dirData};
30	listData={'listData',listData(:).name};
31
32	if isempty(who('dirTex'));
33	addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
34	else;
35	if ~ischar(dirTex); error('mis-specified dirTex'); end;
36	if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path
37	addToTex={'addToTex',1}; dirTex={'dirTex',[dirTex '/']};
38	if isempty(who('nameTex')); nameTex='myPlots'; end;
39	nameTex={'nameTex',nameTex};
40	end;
41
42	if doComp;
43	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	else;
49	insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
50	insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
51	end;
52