--- MITgcm_contrib/jmc_script/grph_StD.m 2015/03/07 14:45:49 1.7 +++ MITgcm_contrib/jmc_script/grph_StD.m 2016/09/10 22:05:35 1.9 @@ -5,11 +5,11 @@ Nexp=size(namA,2); %- -% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/grph_StD.m,v 1.7 2015/03/07 14:45:49 jmc Exp $ +% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/grph_StD.m,v 1.9 2016/09/10 22:05:35 jmc Exp $ % $Name: $ nItMx=1e10*ones(1,Nexp); %nItMx(3)=11; -nItMx=2400*ones(1,Nexp); +%nItMx=2400*ones(1,Nexp); namLg=namA ; namLg=strrep(namLg,'_','\_'); undef=123456.7; %----------- @@ -83,7 +83,7 @@ %-- list_on=zeros(1,nbV); -nbG=8; +nbG=18; nbG=min(nbG,nbV); list_on(1:nbG)=1 ; %if nbG < nbV, list_on=0; list_on(nbG:nbV)=1; end %- to get 2nd set of plots %list_on(1:6)=[1 1 1 1 1 1]; @@ -108,8 +108,9 @@ %titall='Dyncore test-case 5 (cs-32)' ; %========================================================= +ng=0; fxb=100; fyb=60; fdx=100; fdy=40; fsc=1.; +%fyb=-360; fxb=-2600; fdy=60; fsc=1.5; -ng=0; for jv=1:nbV, %------------------- flag=list_on(jv); kl=0; @@ -126,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',[100+100*ng 60+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;