--- MITgcm_contrib/jmc_script/grph_StD.m	2014/08/25 21:35:12	1.5
+++ MITgcm_contrib/jmc_script/grph_StD.m	2017/08/23 22:17:40	1.10
@@ -1,16 +1,17 @@
  prefix='dynStD';
+%prefix='iceStD';
  namA={'g02','b02'};
 %namA={'r17a'};
  Nexp=size(namA,2);
 %-
 
-% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/grph_StD.m,v 1.5 2014/08/25 21:35:12 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;
-nItMx=2400*ones(1,Nexp);
+%nItMx=2400*ones(1,Nexp);
 namLg=namA ; namLg=strrep(namLg,'_','\_');
+undef=123456.7;
 %-----------
 %- test if the variable krd is define :
 if size(who('krd'),1) > 0,
@@ -82,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];
@@ -107,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;
@@ -125,12 +127,15 @@
   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,
-   dd(:,n)=max(var(isA(n):ieA(n),:,n))-min(var(isA(n):ieA(n),:,n));
-   av(:,n)=mean(var(isA(n):ieA(n),:,n));
+   tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0;
+   av(:,n)=sum(tmp);
+   if length(I) < size(tmp,1), av(:,n)=av(:,n)/(size(tmp,1)-length(I)); end
+   tmp(I,:)=NaN; dd(:,n)=max(tmp)-min(tmp);
   end
+  var(find(var==undef))=NaN;
   for nv=1:4,
     subplot(410+nv); ttmn=' Mx-mn:'; ttav=' Av:';
     for n=1:Nexp,
@@ -145,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);
 %--