--- MITgcm_contrib/jmc_script/grph_StD.m	2014/08/25 21:35:12	1.5
+++ MITgcm_contrib/jmc_script/grph_StD.m	2015/03/07 14:45:49	1.7
@@ -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.7 2015/03/07 14:45:49 jmc Exp $
 % $Name:  $
 
-
 nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;
 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,
@@ -128,9 +129,12 @@
   figure(ng); set(ng,'position',[100+100*ng 60+40*ng 500 700]);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,