/[MITgcm]/MITgcm_contrib/jmc_script/grph_StD.m
ViewVC logotype

Diff of /MITgcm_contrib/jmc_script/grph_StD.m

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

revision 1.8 by jmc, Sun Sep 6 18:15:09 2015 UTC revision 1.9 by jmc, Sat Sep 10 22:05:35 2016 UTC
# Line 108  titall='Global Ocean, Cubic-G (32x32) , Line 108  titall='Global Ocean, Cubic-G (32x32) ,
108  %titall='Dyncore test-case 5 (cs-32)' ;  %titall='Dyncore test-case 5 (cs-32)' ;
109    
110  %=========================================================  %=========================================================
111  ng=0; fxb=100; fyb=60;  ng=0; fxb=100; fyb=60; fdx=100; fdy=40; fsc=1.;
112  %fyb=160; fxb=-2600;  %fyb=-360; fxb=-2600; fdy=60; fsc=1.5;
113    
114  for jv=1:nbV,  for jv=1:nbV,
115  %-------------------  %-------------------
# Line 127  for jv=1:nbV, Line 127  for jv=1:nbV,
127    ng=ng+1;    ng=ng+1;
128  %- reset "ng" to jv for fix fig number (independent of which one is in list_on):  %- reset "ng" to jv for fix fig number (independent of which one is in list_on):
129    %ng=jv;    %ng=jv;
130    figure(ng); set(ng,'position',[fxb+100*ng fyb+40*ng 500 700]);clf;    figure(ng); set(ng,'position',[fxb+fdx*ng fyb+fdy*ng [500 700]*fsc]);clf;
131    var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp);    var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp);
132    for n=1:Nexp,    for n=1:Nexp,
133     tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0;     tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0;
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22