/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_diags.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_diags.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by gforget, Fri Jul 1 18:03:43 2011 UTC revision 1.5 by gforget, Tue Jul 29 16:47:50 2014 UTC
# Line 1  Line 1 
1  function []=insitu_diags(dirModel,doComp,varargin);  function []=insitu_diags(dirMat,doComp,dirTex,nameTex);
2  %object:        driver for insitu_misfit and insitu_cost  %object:     driver for insitu_misfit and insitu_cost
3  %inputs:        dirModel is the model run directory      %inputs:     dirMat is the directory where diagnozed .mat files will be saved
4  %               doComp states whether to compute (1) or display (0)  %                     -> 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  if nargin==3; addToTex={'addToTex',varargin{1}}; else; addToTex={'addToTex',0}; end;  gcmfaces_global;
13    
14  if doComp;  dirMat=[dirMat '/'];
15  insitu_misfit(dirModel,1);  
16  insitu_cost(dirModel,1);  if isfield(myenv,'profiles');
17      dirData=myenv.profiles;
18  else;  else;
19  insitu_misfit(dirModel,0,addToTex);    dirData=[dirMat '/profiles/output/'];
20  insitu_cost(dirModel,0,addToTex);  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;  end;
25    
26  return  dirMat={'dirMat',dirMat};
27    dirData={'dirData',dirData};
28    
29  dirModel='./apr1alpha_it0/';  if isempty(who('dirTex'));
30  listData={'listData','WOD09_XBT*'};    addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
31  listVar={'listVar','T'};  else;
32  suffMat={'suffMat','XBT'};    if ~ischar(dirTex); error('mis-specified dirTex'); end;
33  addToTex={'addToTex',1};    if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path
34      addToTex={'addToTex',1}; dirTex={'dirTex',[dirTex '/']};
35  insitu_misfit(dirModel,1,listData,listVar,suffMat);    if isempty(who('nameTex')); nameTex='myPlots'; end;
36  insitu_cost(dirModel,1,listData,listVar,suffMat);    nameTex={'nameTex',nameTex};
37    end;
38  mkdir([dirModel 'tex']);  
39  write2tex([dirModel 'tex/myPlots.tex'],0);  if doComp;
40  insitu_misfit(dirModel,0,suffMat,addToTex);  insitu_misfit(1,dirData,dirMat);
41  insitu_cost(dirModel,0,suffMat,addToTex);  insitu_cost(1,dirData,dirMat);
42  write2tex([dirModel 'tex/myPlots.tex'],4);  else;
43    insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
44    insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
45    end;
46    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22