/[MITgcm]/MITgcm/verification/offline_exf_seaice/input/gendata.m
ViewVC logotype

Diff of /MITgcm/verification/offline_exf_seaice/input/gendata.m

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

revision 1.3 by jmc, Wed Jan 9 22:10:20 2013 UTC revision 1.4 by jmc, Fri Mar 7 22:30:26 2014 UTC
# Line 1  Line 1 
1    
2  kwr=1;  kwr=1; kprt=0;
3    kwr=0; kprt=0;
4  nx=80; ny=42; nr=3; nt=1;  nx=80; ny=42; nr=3; nt=1;
5    
6  xc=[1:nx]; xc=xc-mean(xc);  xc=[1:nx]; xc=xc-mean(xc);
# Line 195  dewPt=(qa_x*atmrho)/cvapor_fac; Line 196  dewPt=(qa_x*atmrho)/cvapor_fac;
196  dewPt=-cvapor_exp./log(dewPt);  dewPt=-cvapor_exp./log(dewPt);
197    
198  figure(2);clf;  figure(2);clf;
199  subplot(211)  subplot(311)
200  plot(xc,ta_x-cel2K,'r-'); hold on;  P(1)=plot(xc,ta_x-cel2K,'r-'); hold on;
201  plot(xc,dewPt-cel2K,'b-');  P(2)=plot(xc,dewPt-cel2K,'b-');
202  plot(xc,tfreeze*ones(nx,1),'k-');  P(3)=plot(xc,tfreeze*ones(nx,1),'k-');
203  hold off;  set(P,'LineWidth',1);
204  AA=axis; axis([-nx/2 nx/2 AA(3:4)]);  hold off; AA=axis;
205    axis([-nx/2 nx/2 AA(3:4)]);
206  legend('ta','dew');  legend('ta','dew');
207  grid  grid
208  title(['del-Temp-X= ',int2str(dtx),' ; RH= ',int2str(rh),' ; Air Temp (^oC)']);  xlabel('X')
209  subplot(212)  title(['Air Temp (^oC): del-Temp-X = ',int2str(dtx),' , RH= ',int2str(rh)]);
210  plot(yc,to_y,'b-'); hold on;  subplot(312)
211  plot(yc,tfreeze*ones(ny,1),'k-');  P(1)=plot(yc,to_y,'b-'); hold on;
212  hold off;  P(2)=plot(yc,tfreeze*ones(ny,1),'k-');
213  AA=axis; axis([0 ny AA(3:4)]);  set(P,'LineWidth',1);
214    hold off; AA=axis;
215    L=line([1 1],AA(3:4)); set(L,'LineWidth',2.,'Color',[0 0 0]);
216    axis([0 ny AA(3:4)]);
217  grid  grid
218    xlabel('Y')
219  title('Ocean Temp ^oC');  title('Ocean Temp ^oC');
220    
221  %--  subplot(313)
222    var=iceConc(1,:);
223    P(1)=semilogy(yc,var,'b-x'); hold on;
224    %plot(yc,var,'b-x'); hold on;
225    var=iceVol(1,:);
226    P(2)=semilogy(yc,var,'r-x');
227    %plot(yc,var,'r-x');
228    set(P,'LineWidth',1);
229    hold off; AA=axis;
230    L=line([1 1],AA(3:4)); set(L,'LineWidth',2.,'Color',[0 0 0]);
231    axis([0 ny [0 2]*iceC0]);
232    grid
233    xlabel('Y')
234    legend('iceC','hEff','Location','South');
235    title('Initial ice in Channel : y-section');
236    %-----
237    if kprt == 1, f=2;
238     namfig=sprintf('forcing_%2.2i',f);
239     fprintf([' print fig= %2i to file: ',namfig,' '],f);
240     set(f,'PaperOrientation','portrait')
241    %set(f,'PaperPosition',[0.25 1.5 6. 8.]);
242     set(f,'PaperPosition',[0.25 1.5 5.25 7.]);
243     print(f,'-depsc2',namfig); fprintf('\n');
244    end
245    %-----
246    
247  figure(3);clf;  figure(3);clf;
248  subplot(311)  subplot(311)
# Line 242  AA=axis; axis([0 ny [0 2]*iceC0]); Line 272  AA=axis; axis([0 ny [0 2]*iceC0]);
272  grid  grid
273  legend('iceC','hEff','Location','South');  legend('iceC','hEff','Location','South');
274  title('Initial ice in Channel : y-section');  title('Initial ice in Channel : y-section');
275  %--  %-----
276    
277    return
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22