--- MITgcm_contrib/jmc_script/grph_StD.m	2015/03/05 20:57:26	1.6
+++ MITgcm_contrib/jmc_script/grph_StD.m	2016/09/10 22:05:35	1.9
@@ -1,14 +1,15 @@
  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.6 2015/03/05 20:57:26 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;
 %-----------
@@ -32,7 +33,6 @@
   [ntA(n),rList,tim,vv1,listV] = ...
     read_StD(prefix,namA(n),listV);
   nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
-  vv1(find(vv1==undef))=NaN;
 %- set global dims: & load vvA --> vvB
   nbV=size(listV,2);
   nrec=nIt; n3d=nk; nReg=nRg;
@@ -43,7 +43,6 @@
   [ntA(n),rList,tim,vv1,listV] = ...
     read_StD(prefix,namA(n),listV);
   nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
-  vv1(find(vv1==undef))=NaN;
   if (nrec < nIt),
     fprintf('\n');
     error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]);
@@ -84,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];
@@ -109,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;
@@ -127,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,