/[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.11 by gforget, Sun Jan 17 16:00:24 2016 UTC
# Line 21  vecCaxis=[[3 3 1];[0.5 0.5 0.2]]; Line 21  vecCaxis=[[3 3 1];[0.5 0.5 0.2]];
21    
22  %load paths and grid:  %load paths and grid:
23  %--------------------  %--------------------
24  global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end;  gcmfaces_global; if isempty(mygrid); grid_load; end;
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);
# 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.11
  ViewVC Help
Powered by ViewVC 1.1.22