/[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.10 by jmc, Wed Aug 23 22:17:40 2017 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;
# Line 150  for jv=1:nbV, Line 150  for jv=1:nbV,
150      axis(AA); grid ;      axis(AA); grid ;
151      if nv == 1, title(['Avr ',titv,'  ',ttmn]); end      if nv == 1, title(['Avr ',titv,'  ',ttmn]); end
152      if nv == 2, title(['Std-Dev ',titv,'  ',ttav]); end      if nv == 2, title(['Std-Dev ',titv,'  ',ttav]); end
153      if nv == 3, title(['min ',titv,'  ',ttav]); legend(namLg,0); end      if nv == 3, title(['min ',titv,'  ',ttav]);
154                    legend(namLg,'Location','best'); end
155      if nv == 4, title(['Max ',titv,'  ',ttav]); end      if nv == 4, title(['Max ',titv,'  ',ttav]); end
156    end ; xlabel(titT);    end ; xlabel(titT);
157  %--  %--
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22