/[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.9 by gforget, Sat Jan 18 23:42:50 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 114  eval(['save ' dirMat 'insitu_cost' suffM Line 115  eval(['save ' dirMat 'insitu_cost' suffM
115    
116  else;%display result  else;%display result
117    
118    if isdir([dirMat 'cost/']); dirMat=[dirMat 'cost/']; end;
119    
120  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);
121    
122  figureL;  figureL;
# Line 143  myCaption={'Cost function (top) for in s Line 146  myCaption={'Cost function (top) for in s
146          '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).'};
147  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;
148    
149    if ~isempty(listBas);
150    
151  ii=0;  ii=0;
152  ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=3;  ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=3;
153  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 186  if addToTex; write2tex([dirTex '/' nameT
186    
187  end;  end;
188    
189    end;
190    
191  end;  end;
192    
193    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22