/[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.1 by gforget, Fri Jul 1 18:03:43 2011 UTC revision 1.4 by gforget, Fri Jul 8 12:48:59 2011 UTC
# Line 26  global mygrid; if isempty(mygrid); grid_ Line 26  global mygrid; if isempty(mygrid); grid_
26    
27  %set more optional paramaters to default values  %set more optional paramaters to default values
28  year0=1992; year1=2008;  year0=1992; year1=2008;
29  %listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'};  listData={'WOD09_XBT*','WOD09_CTD*','argo_in*','argo_at*','argo_pa*','seals*'};
30  %listVar={'T','S'};  listVar={'T','S'};
31  listData={'argo_in*'}; listVar={'T'};  %shorter test case: listData={'argo_in*'}; listVar={'T'};
32  subdirData='/profiles/output/';  subdirData='/profiles/output/';
33  subdirMat='mat/'; suffMat='_all';  subdirMat='mat/'; suffMat='_all';
34  subdirTex='tex/'; addToTex=0;  subdirTex='tex/'; addToTex=0;
# Line 81  else; Line 81  else;
81  eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']);  eval(['load ' dirData '/' subdirMat 'insitu_misfit' suffMat '.mat;']);
82    
83  nk=size(vecDepth,1);  nk=size(vecDepth,1);
84  figureL; orient tall;  figureL;
85  for vv=1:length(listVar);  for vv=1:length(listVar);
86      varCur=listVar{vv};      varCur=listVar{vv};
87      for kkk=1:size(vecDepth,1);      for kkk=1:size(vecDepth,1);
# Line 95  for vv=1:length(listVar); Line 95  for vv=1:length(listVar);
95      end;      end;
96  end;  end;
97    
98  myCaption={'model minus in situ profiles'};  myCaption={'Time mean misfit (model-data) for in situ profiles, ',...
99             'at various depths (rows), for T (left; in K) and S (right; in psu).'};
100  if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end;  if addToTex; write2tex([dirData '/' subdirTex '/myPlots.tex'],2,myCaption,gcf); end;
101    
102  end;  end;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22