/[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.3 by gforget, Sun Dec 30 00:35:53 2012 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    dirMat=[dirMat '/'];
10    
11    dirData=[dirMat '/profiles/output/'];
12    while ~isdir(dirData);
13      fprintf(['directory : ' dirData '\n'])
14      dirData=input(['does not exist. Specify directory of nc file : \n']);
15    end;
16    
17  if nargin==3; addToTex={'addToTex',varargin{1}}; else; addToTex={'addToTex',0}; end;  dirMat={'dirMat',dirMat};
18    dirData={'dirData',dirData};
19    
20  if doComp;  if isempty(who('dirTex'));
21  insitu_misfit(dirModel,1);    addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
 insitu_cost(dirModel,1);  
22  else;  else;
23  insitu_misfit(dirModel,0,addToTex);    if ~ischar(dirTex); error('mis-specified dirTex'); end;
24  insitu_cost(dirModel,0,addToTex);    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;  end;
29    
30  return  if doComp;
31    insitu_misfit(1,dirData,dirMat);
32  dirModel='./apr1alpha_it0/';  insitu_cost(1,dirData,dirMat);
33  listData={'listData','WOD09_XBT*'};  else;
34  listVar={'listVar','T'};  insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
35  suffMat={'suffMat','XBT'};  insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
36  addToTex={'addToTex',1};  end;
   
 insitu_misfit(dirModel,1,listData,listVar,suffMat);  
 insitu_cost(dirModel,1,listData,listVar,suffMat);  
   
 mkdir([dirModel 'tex']);  
 write2tex([dirModel 'tex/myPlots.tex'],0);  
 insitu_misfit(dirModel,0,suffMat,addToTex);  
 insitu_cost(dirModel,0,suffMat,addToTex);  
 write2tex([dirModel 'tex/myPlots.tex'],4);  
37    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22