/[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.12 by gforget, Mon Feb 15 18:07:31 2016 UTC
# Line 15  function []=insitu_cost(doComp,varargin) Line 15  function []=insitu_cost(doComp,varargin)
15  %  %
16  %example:         insitu_cost('./',1,{'listData','argo_in*'});  %example:         insitu_cost('./',1,{'listData','argo_in*'});
17    
18  %load paths and grid:  gcmfaces_global;
19  %--------------------  global myparms;
20  global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end;  
   
 %set more optional paramaters to default values  
 global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);  
 listData={'argo_*_1992_to_2007*','argo_*_2008_to_2010*',...  
           'argo_*_2011_to_2012*','climode_*',...  
           'ctd_*','itp_*','seals_*','xbt_*'};  
21  listVar={'T','S'};  listVar={'T','S'};
22  listBas={'atlExt','pacExt','indExt','arct'};  listBas={'atlExt','pacExt','indExt','arct'};
23    if sum([90 1170]~=mygrid.ioSize)>0; listBas={}; end;
24  listBasTxt=''; for bb=1:length(listBas); listBasTxt=[listBasTxt ' ' listBas{bb} '*']; end;  listBasTxt=''; for bb=1:length(listBas); listBasTxt=[listBasTxt ' ' listBas{bb} '*']; end;
25  %shorter test case: listData={'argo_in*'}; listVar={'T'};  %shorter test case: listData={'argo_in*'}; listVar={'T'};
26  dirData='/profiles/output/';  dirData='/profiles/output/';
# Line 48  for ii=1:nargin-1; Line 43  for ii=1:nargin-1;
43     end;     end;
44  end;  end;
45    
46    if isempty(myparms)|isempty(mygrid);
47      load([dirMat 'diags_grid_parms.mat']);
48    end;
49    
50    if isempty(whos('year0'));
51      year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);
52    end;
53    
54    if isempty(whos('listData'));
55      listData=dir([dirMat 'profiles/output/*.nc'])
56      listData={listData(:).name};
57      for ff=1:length(listData); listData{ff}=[listData{ff}(1:end-3) '*']; end;
58    end;
59    
60  if doComp;  if doComp;
61    
62  %time limits:  %time limits:
# Line 109  for vv=1:length(listVar); Line 118  for vv=1:length(listVar);
118    
119  end;  end;
120    
121  eval(['save ' dirMat 'insitu_cost' suffMat '.mat costAve* misfitDistrib* depthTimeCost* year* ' ...  if ~isdir([dirMat 'cost/']); mkdir([dirMat 'cost/']); end;
122    eval(['save ' dirMat 'cost/insitu_cost' suffMat '.mat costAve* misfitDistrib* depthTimeCost* year* ' ...
123        'listData listVar listBas ' listBasTxt ';']);        'listData listVar listBas ' listBasTxt ';']);
124    
125  else;%display result  else;%display result
126    
127    if isdir([dirMat 'cost/']); dirMat=[dirMat 'cost/']; end;
128    
129  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);
130    
131  figureL;  figureL;
# Line 143  myCaption={'Cost function (top) for in s Line 155  myCaption={'Cost function (top) for in s
155          '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).'};
156  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;
157    
158    if ~isempty(listBas);
159    
160  ii=0;  ii=0;
161  ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=3;  ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=3;
162  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 195  if addToTex; write2tex([dirTex '/' nameT
195    
196  end;  end;
197    
198    end;
199    
200  end;  end;
201    
202    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22