--- MITgcm_contrib/jmc_script/grph_StD.m 2015/09/06 18:15:09 1.8 +++ MITgcm_contrib/jmc_script/grph_StD.m 2017/08/23 22:17:40 1.10 @@ -5,7 +5,7 @@ Nexp=size(namA,2); %- -% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/grph_StD.m,v 1.8 2015/09/06 18:15:09 jmc Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/grph_StD.m,v 1.10 2017/08/23 22:17:40 jmc Exp $ % $Name: $ nItMx=1e10*ones(1,Nexp); %nItMx(3)=11; @@ -108,8 +108,8 @@ %titall='Dyncore test-case 5 (cs-32)' ; %========================================================= -ng=0; fxb=100; fyb=60; -%fyb=160; fxb=-2600; +ng=0; fxb=100; fyb=60; fdx=100; fdy=40; fsc=1.; +%fyb=-360; fxb=-2600; fdy=60; fsc=1.5; for jv=1:nbV, %------------------- @@ -127,7 +127,7 @@ ng=ng+1; %- reset "ng" to jv for fix fig number (independent of which one is in list_on): %ng=jv; - 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; var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp); for n=1:Nexp, tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0; @@ -150,7 +150,8 @@ axis(AA); grid ; if nv == 1, title(['Avr ',titv,' ',ttmn]); end if nv == 2, title(['Std-Dev ',titv,' ',ttav]); end - if nv == 3, title(['min ',titv,' ',ttav]); legend(namLg,0); end + if nv == 3, title(['min ',titv,' ',ttav]); + legend(namLg,'Location','best'); end if nv == 4, title(['Max ',titv,' ',ttav]); end end ; xlabel(titT); %--