/[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.7 by gforget, Tue Apr 2 13:55:53 2013 UTC revision 1.10 by gforget, Sat Feb 1 18:06:35 2014 UTC
# 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.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22