/[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.6 by gforget, Sun Dec 30 00:35:53 2012 UTC revision 1.10 by gforget, Sat Feb 1 18:06:35 2014 UTC
# Line 25  global mygrid; if isempty(mygrid); grid_ Line 25  global mygrid; if isempty(mygrid); grid_
25    
26  %set more optional paramaters to default values  %set more optional paramaters to default values
27  global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);  global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);
28  listData={'argo_june2012_1992_to_2007*','argo_june2012_2008_to_2010*',...  listData={'argo_*_1992_to_2007*','argo_*_2008_to_2010*',...
29            'argo_june2012_2011_to_2012*','climode_june2012*',...            'argo_*_2011_to_2012*','climode_*',...
30            'ctd_june2012*','itp_june2012*','seals_june2012*','xbt_june2012*'};            'ctd_*','itp_*','seals_*','xbt_*'};
31  listVar={'T','S'};  listVar={'T','S'};
32  %shorter test case: listData={'argo_in*'}; listVar={'T'};  %shorter test case: listData={'argo_in*'}; listVar={'T'};
33  dirData='mat/profiles/output/';  dirData='mat/profiles/output/';
# Line 76  for vv=1:length(listVar); Line 76  for vv=1:length(listVar);
76      end;      end;
77  end;  end;
78    
79  eval(['save ' dirMat 'insitu_misfit' suffMat '.mat misfit_map_* vecDepth  year* listData listVar;']);  if ~isdir([dirMat 'cost/']); mkdir([dirMat 'cost/']); end;
80    eval(['save ' dirMat 'cost/insitu_misfit' suffMat '.mat misfit_map_* vecDepth  year* listData listVar;']);
81    
82  else;  else;
83    
84    if isdir([dirMat 'cost/']); dirMat=[dirMat 'cost/']; end;
85    
86  eval(['load ' dirMat 'insitu_misfit' suffMat '.mat;']);  eval(['load ' dirMat 'insitu_misfit' suffMat '.mat;']);
87    
88  nk=size(vecDepth,1);  nk=size(vecDepth,1);
# Line 92  for vv=1:length(listVar); Line 95  for vv=1:length(listVar);
95    
96          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;
97          subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));          subplot(nk,length(listVar),vv+(kkk-1)*length(listVar));
98          m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc});          m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc},{'myShading','flat'});
99          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));
100      end;      end;
101  end;  end;
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22