/[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.6 by gforget, Tue May 10 21:38:51 2016 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  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  dirModel='./apr1alpha_it0/';  if doComp;
43  listData={'listData','WOD09_XBT*'};  fprintf('starting insitu_misfit\n'); clock
44  listVar={'listVar','T'};  insitu_misfit(1,dirData,dirMat,listData);
45  suffMat={'suffMat','XBT'};  fprintf('starting insitu_cost\n'); clock
46  addToTex={'addToTex',1};  insitu_cost(1,dirData,dirMat,listData);
47    fprintf('done with insitu_diags\n'); clock
48  insitu_misfit(dirModel,1,listData,listVar,suffMat);  else;
49  insitu_cost(dirModel,1,listData,listVar,suffMat);  insitu_misfit(0,dirMat,addToTex,dirTex,nameTex);
50    insitu_cost(0,dirMat,addToTex,dirTex,nameTex);
51  mkdir([dirModel 'tex']);  end;
 write2tex([dirModel 'tex/myPlots.tex'],0);  
 insitu_misfit(dirModel,0,suffMat,addToTex);  
 insitu_cost(dirModel,0,suffMat,addToTex);  
 write2tex([dirModel 'tex/myPlots.tex'],4);  
52    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22