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

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

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

revision 1.7 by gforget, Wed Apr 3 14:30:12 2013 UTC revision 1.10 by gforget, Sun Jan 19 15:00:17 2014 UTC
# Line 26  listData={'argo_*_1992_to_2007*','argo_* Line 26  listData={'argo_*_1992_to_2007*','argo_*
26            'ctd_*','itp_*','seals_*','xbt_*'};            'ctd_*','itp_*','seals_*','xbt_*'};
27  listVar={'T','S'};  listVar={'T','S'};
28  listBas={'atlExt','pacExt','indExt','arct'};  listBas={'atlExt','pacExt','indExt','arct'};
29    if sum([90 1170]~=mygrid.ioSize)>0; listBas={}; end;
30  listBasTxt=''; for bb=1:length(listBas); listBasTxt=[listBasTxt ' ' listBas{bb} '*']; end;  listBasTxt=''; for bb=1:length(listBas); listBasTxt=[listBasTxt ' ' listBas{bb} '*']; end;
31  %shorter test case: listData={'argo_in*'}; listVar={'T'};  %shorter test case: listData={'argo_in*'}; listVar={'T'};
32  dirData='/profiles/output/';  dirData='/profiles/output/';
# Line 109  for vv=1:length(listVar); Line 110  for vv=1:length(listVar);
110    
111  end;  end;
112    
113  eval(['save ' dirMat 'insitu_cost' suffMat '.mat costAve* misfitDistrib* depthTimeCost* year* ' ...  if ~isdir([dirMat 'cost/']); mkdir([dirMat 'cost/']); end;
114    eval(['save ' dirMat 'cost/insitu_cost' suffMat '.mat costAve* misfitDistrib* depthTimeCost* year* ' ...
115        'listData listVar listBas ' listBasTxt ';']);        'listData listVar listBas ' listBasTxt ';']);
116    
117  else;%display result  else;%display result
118    
119    if isdir([dirMat 'cost/']); dirMat=[dirMat 'cost/']; end;
120    
121  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);
122    
123  figureL;  figureL;
# Line 143  myCaption={'Cost function (top) for in s Line 147  myCaption={'Cost function (top) for in s
147          'Distribution of normalized misfits (bottom) as a function of latitude. For T (left) and S (right).'};          'Distribution of normalized misfits (bottom) as a function of latitude. For T (left) and S (right).'};
148  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;
149    
150    if ~isempty(listBas);
151    
152  ii=0;  ii=0;
153  ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=3;  ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=3;
154  ii=ii+1; listPanels(ii).bas='pacExt'; listPanels(ii).ll=3;  ii=ii+1; listPanels(ii).bas='pacExt'; listPanels(ii).ll=3;
# Line 181  if addToTex; write2tex([dirTex '/' nameT Line 187  if addToTex; write2tex([dirTex '/' nameT
187    
188  end;  end;
189    
190    end;
191    
192  end;  end;
193    
194    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22