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

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

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

revision 1.8 by gforget, Sat Jan 18 23:42:50 2014 UTC revision 1.12 by gforget, Mon Feb 15 18:07:32 2016 UTC
# Line 19  vecDepth=[[75 125];[250 350];[900 1100]] Line 19  vecDepth=[[75 125];[250 350];[900 1100]]
19  vecCaxis=[[3 3 1];[0.5 0.5 0.2]];  vecCaxis=[[3 3 1];[0.5 0.5 0.2]];
20  %if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;  %if nargin>6; vecCaxis=varargin{1}; else; vecCaxis=[]; end;
21    
22  %load paths and grid:  gcmfaces_global;
23  %--------------------  global myparms;
24  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_*'};  
25  listVar={'T','S'};  listVar={'T','S'};
26  %shorter test case: listData={'argo_in*'}; listVar={'T'};  %shorter test case: listData={'argo_in*'}; listVar={'T'};
27  dirData='mat/profiles/output/';  dirData='mat/profiles/output/';
# Line 50  for ii=1:nargin-1; Line 44  for ii=1:nargin-1;
44     end;     end;
45  end;  end;
46    
47    if isempty(myparms)|isempty(mygrid);
48      load([dirMat 'diags_grid_parms.mat']);
49    end;
50    
51    if isempty(whos('year0'));
52      year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);
53    end;
54    
55    if isempty(whos('listData'));
56      listData=dir([dirMat 'profiles/output/*.nc']);
57      listData={listData(:).name};
58      for ff=1:length(listData); listData{ff}=[listData{ff}(1:end-3) '*']; end;
59    end;
60    
61  if doComp;  if doComp;
62    
63  %initialize delaunay triangulation:  %initialize delaunay triangulation:
# Line 76  for vv=1:length(listVar); Line 84  for vv=1:length(listVar);
84      end;      end;
85  end;  end;
86    
87  eval(['save ' dirMat 'insitu_misfit' suffMat '.mat misfit_map_* vecDepth  year* listData listVar;']);  if ~isdir([dirMat 'cost/']); mkdir([dirMat 'cost/']); end;
88    eval(['save ' dirMat 'cost/insitu_misfit' suffMat '.mat misfit_map_* vecDepth  year* listData listVar;']);
89    
90  else;  else;
91    
# Line 94  for vv=1:length(listVar); Line 103  for vv=1:length(listVar);
103    
104          if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end;          if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end;
105          subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));          subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));
106          m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc});          m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc},{'myShading','flat'});
107          title(sprintf('%s misfits at %dm   (%d to %d)',varCur,dep,year0,year1));          title(sprintf('%s misfits at %dm   (%d to %d)',varCur,dep,year0,year1));
108      end;      end;
109  end;  end;
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22