/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_cost.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/profilesMatlabProcessing/ecco_v4/insitu_cost.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.6 by gforget, Tue Apr 2 13:55:53 2013 UTC revision 1.9 by gforget, Sat Jan 18 23:42:50 2014 UTC
# Line 25  listData={'argo_*_1992_to_2007*','argo_* Line 25  listData={'argo_*_1992_to_2007*','argo_*
25            'argo_*_2011_to_2012*','climode_*',...            'argo_*_2011_to_2012*','climode_*',...
26            'ctd_*','itp_*','seals_*','xbt_*'};            'ctd_*','itp_*','seals_*','xbt_*'};
27  listVar={'T','S'};  listVar={'T','S'};
28    listBas={'atlExt','pacExt','indExt','arct'};
29    if sum([90 1170]~=mygrid.ioSize)>0; listBas={}; end;
30    listBasTxt=''; for bb=1:length(listBas); listBasTxt=[listBasTxt ' ' listBas{bb} '*']; end;
31  %shorter test case: listData={'argo_in*'}; listVar={'T'};  %shorter test case: listData={'argo_in*'}; listVar={'T'};
32  dirData='/profiles/output/';  dirData='/profiles/output/';
33  dirMat='mat/'; suffMat='_all';  dirMat='mat/'; suffMat='_all';
# Line 82  for vv=1:length(listVar); Line 85  for vv=1:length(listVar);
85      eval(['misfitDistrib' varCur '.x=x;']); eval(['misfitDistrib' varCur '.y=y;']);      eval(['misfitDistrib' varCur '.x=x;']); eval(['misfitDistrib' varCur '.y=y;']);
86      eval(['misfitDistrib' varCur '.z=z;']); eval(['misfitDistrib' varCur '.n=n;']);      eval(['misfitDistrib' varCur '.z=z;']); eval(['misfitDistrib' varCur '.n=n;']);
87    
88        %upper ocean temporal distribution
89        for bb=1:length(listBas);
90        for ll=1:3;
91          bbb=listBas{bb};
92          if ll==1;     Lmin=-90; Lmax=-25; txt=[bbb '_90S25S_misDis' varCur];
93          elseif ll==2; Lmin=-25; Lmax=25; txt=[bbb '_25S25N_misDis' varCur];
94          elseif ll==3; Lmin=25; Lmax=90; txt=[bbb '_25N90N_misDis' varCur];
95          end;
96          msk=v4_basin(bbb);
97          msk(mygrid.YC<Lmin|mygrid.YC>Lmax)=0;
98          %
99          msk=convert2array(msk);
100          MITprof.prof_msk=msk(MITprof.prof_point);    
101          MITprofSub=MITprof_subset(MITprof,'msk',1,'depth',[0 700]);      
102          %
103          tmp_date=(MITprofSub.prof_date-date0)/365;
104          [x,y,z,n]=MITprof_stats(tmp_date,[0:1/4:1+(year1-year0)],MITprofSub.prof,[-5:0.25:5]);
105          z(n==0)=NaN; n(n==0)=NaN;
106          eval([txt '.x=x;']); eval([txt '.y=y;']);
107          eval([txt '.z=z;']); eval([txt '.n=n;']);
108        end;
109        end;
110    
111  end;  end;
112    
113  eval(['save ' dirMat 'insitu_cost' suffMat '.mat costAve* misfitDistrib* depthTimeCost* year* listData listVar;']);  eval(['save ' dirMat 'insitu_cost' suffMat '.mat costAve* misfitDistrib* depthTimeCost* year* ' ...
114          'listData listVar listBas ' listBasTxt ';']);
115    
116  else;%display result  else;%display result
117    
118    if isdir([dirMat 'cost/']); dirMat=[dirMat 'cost/']; end;
119    
120  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);  eval(['load ' dirMat 'insitu_cost' suffMat '.mat;']);
121    
122  figureL;  figureL;
# Line 99  for vv=1:length(listVar); Line 128  for vv=1:length(listVar);
128    
129      eval(['x=depthTimeCost' varCur '.x;']); eval(['y=depthTimeCost' varCur '.y;']);      eval(['x=depthTimeCost' varCur '.x;']); eval(['y=depthTimeCost' varCur '.y;']);
130      eval(['z=depthTimeCost' varCur '.z;']); eval(['n=depthTimeCost' varCur '.n;']);      eval(['z=depthTimeCost' varCur '.z;']); eval(['n=depthTimeCost' varCur '.n;']);
131      z(n<1e2)=NaN;      z(n<1e2)=NaN; y=y+year0-1;
132      subplot(2,length(listVar),vv); depthStretchPlot('pcolor',{y,x,z},[0:200:1000 1500:500:6000],[0 1000 4000]);      subplot(2,length(listVar),vv); depthStretchPlot('pcolor',{y,x,z},[0:200:1000 1500:500:6000],[0 1000 4000]);
133      caxis([0 8]); shading flat; colorbar; ylabel('depth (in m)');      caxis([0 8]); shading flat; colorbar; ylabel('depth (in m)');
     xlabel(['date (in years since ' num2str(year0) ')']);  
134      title(sprintf('mean cost for %s: %0.3g (median: %0.3g)',varCur,costCur));      title(sprintf('mean cost for %s: %0.3g (median: %0.3g)',varCur,costCur));
135    
136      %latitudinal distribution:      %latitudinal distribution:
# Line 118  myCaption={'Cost function (top) for in s Line 146  myCaption={'Cost function (top) for in s
146          'Distribution of normalized misfits (bottom) as a function of latitude. For T (left) and S (right).'};          'Distribution of normalized misfits (bottom) as a function of latitude. For T (left) and S (right).'};
147  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;  if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;
148    
149    if ~isempty(listBas);
150    
151    ii=0;
152    ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=3;
153    ii=ii+1; listPanels(ii).bas='pacExt'; listPanels(ii).ll=3;
154    ii=ii+1; listPanels(ii).bas='arct'; listPanels(ii).ll=3;
155    ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=2;
156    ii=ii+1; listPanels(ii).bas='pacExt'; listPanels(ii).ll=2;
157    ii=ii+1; listPanels(ii).bas='indExt'; listPanels(ii).ll=2;
158    ii=ii+1; listPanels(ii).bas='atlExt'; listPanels(ii).ll=1;
159    ii=ii+1; listPanels(ii).bas='pacExt'; listPanels(ii).ll=1;
160    ii=ii+1; listPanels(ii).bas='indExt'; listPanels(ii).ll=1;
161    
162    for vv=1:length(listVar);
163      varCur=listVar{vv};
164      figureL;
165      for ii=1:9;
166          bbb=listPanels(ii).bas; ll=listPanels(ii).ll;
167          if ll==1;     Lmin=-90; Lmax=-25; txt=[bbb '_90S25S_misDis' varCur];
168          elseif ll==2; Lmin=-25; Lmax=25; txt=[bbb '_25S25N_misDis' varCur];
169          elseif ll==3; Lmin=25; Lmax=90; txt=[bbb '_25N90N_misDis' varCur];
170          end;
171    
172        eval(['x=' txt '.x;']); eval(['y=' txt '.y;']);
173        eval(['z=' txt '.z;']); eval(['n=' txt '.n;']);
174        x=x+year0-1;
175        tmp1=(nansum(n,2)>1e4)*ones(1,size(n,2)); z(tmp1==0)=NaN;
176        subplot(3,3,ii); contourf(x,y,z,[0:0.05:0.5]); grid on;
177        caxis([0 0.5]); colorbar; title(txt,'Interpreter','none');
178        %title([txt ', normalized misfit pdf'],'Interpreter','none');
179        if mod(ii,3)==1; ylabel('normalized misfit'); end;
180        %if ii<=6; set(gca,'XTick',[]); end;
181      end;
182    
183    myCaption={'Distribution of normalized misfits per basin (panel) as ',...
184            'a function of latitude, for ',varCur};
185    if addToTex; write2tex([dirTex '/' nameTex '.tex'],2,myCaption,gcf); end;
186    
187    end;
188    
189    end;
190    
191  end;  end;
192    
193    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22