/[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.2 by gforget, Fri Jul 1 18:15:16 2011 UTC revision 1.5 by gforget, Wed Apr 18 20:24:45 2012 UTC
# Line 25  vecCaxis=[[3 3 1];[0.5 0.5 0.2]]; Line 25  vecCaxis=[[3 3 1];[0.5 0.5 0.2]];
25  global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end;  global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end;
26    
27  %set more optional paramaters to default values  %set more optional paramaters to default values
28  year0=1992; year1=2008;  global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2);
29    %year0=1992; year1=2008;
30  listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'};  listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'};
31  listVar={'T','S'};  listVar={'T','S'};
32  %shorter test case: listData={'argo_in*'}; listVar={'T'};  %shorter test case: listData={'argo_in*'}; listVar={'T'};
# Line 81  else; Line 82  else;
82  eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']);  eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']);
83    
84  nk=size(vecDepth,1);  nk=size(vecDepth,1);
85  figureL; orient tall;  figureL;
86  for vv=1:length(listVar);  for vv=1:length(listVar);
87      varCur=listVar{vv};      varCur=listVar{vv};
88      for kkk=1:size(vecDepth,1);      for kkk=1:size(vecDepth,1);
# Line 95  for vv=1:length(listVar); Line 96  for vv=1:length(listVar);
96      end;      end;
97  end;  end;
98    
99  myCaption={'model minus in situ profiles'};  myCaption={'Time mean misfit (model-data) for in situ profiles, ',...
100             'at various depths (rows), for T (left; in K) and S (right; in psu).'};
101  if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end;  if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end;
102    
103  end;  end;
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22