/[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.2 by gforget, Wed Apr 18 20:24:45 2012 UTC revision 1.3 by gforget, Sun Dec 30 00:35:53 2012 UTC
# Line 1  Line 1 
1  function []=insitu_diags(dirModel,dirMat,doComp,dirTex);  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
 %            dirMat is the directory where diagnozed .mat files will be saved  
4  %                     -> set it to '' to use the default [dirModel 'mat/']  %                     -> set it to '' to use the default [dirModel 'mat/']
5  %            doComp states whether to compute (1) or display (0)  %            doComp states whether to compute (1) or display (0)
6  %            dirTex is the directory where tex and fig files will be created  %            dirTex is the directory where tex and fig files will be created
7    %            nameTex is the tex file name (default : 'myPlots')
8    
9  if dirMat(1)~='/'; dirMat=[pwd '/' dirMat]; end; %make full path  dirMat=[dirMat '/'];
 subdirMat={'subdirMat',[dirMat '/']};  
10    
11  subdirData={'subdirData',[dirMat '/profiles/output/']};  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    dirMat={'dirMat',dirMat};
18    dirData={'dirData',dirData};
19    
20  if isempty(who('dirTex'));  if isempty(who('dirTex'));
21    addToTex={'addToTex',0}; subdirTex={'subdirTex',''};    addToTex={'addToTex',0}; dirTex={'dirTex',''}; nameTex={'nameTex',''};
22  else;  else;
23    if ~ischar(dirTex); error('mis-specified dirTex'); end;    if ~ischar(dirTex); error('mis-specified dirTex'); end;
24    if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path    if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path
25    addToTex={'addToTex',1}; subdirTex={'subdirTex',[dirTex '/']};    addToTex={'addToTex',1}; dirTex={'dirTex',[dirTex '/']};
26      if isempty(who('nameTex')); nameTex='myPlots'; end;
27      nameTex={'nameTex',nameTex};
28  end;  end;
29    
30  if doComp;  if doComp;
31  insitu_misfit('',1,subdirData,subdirMat);  insitu_misfit(1,dirData,dirMat);
32  insitu_cost('',1,subdirData,subdirMat);  insitu_cost(1,dirData,dirMat);
33  else;  else;
34  insitu_misfit('',0,subdirData,subdirMat,addToTex,subdirTex);  insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
35  insitu_cost('',0,subdirData,subdirMat,addToTex,subdirTex);  insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
36  end;  end;
37    
 return  
 %old example that may or may not work ...  
   
 dirModel='./apr1alpha_it0/';  
 listData={'listData','WOD09_XBT*'};  
 listVar={'listVar','T'};  
 suffMat={'suffMat','XBT'};  
 addToTex={'addToTex',1};  
   
 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);  
   
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22