--- MITgcm_contrib/jmc_script/read_StD.m	2008/10/08 18:18:04	1.3
+++ MITgcm_contrib/jmc_script/read_StD.m	2014/05/14 16:53:04	1.5
@@ -15,11 +15,12 @@
 %  vvA      = 5 dims output array:
 %           ( kLev, time_rec, region_rec, [ave,std,min,max,vol], var_rec )
 
-% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/read_StD.m,v 1.3 2008/10/08 18:18:04 jmc Exp $
+% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/read_StD.m,v 1.5 2014/05/14 16:53:04 jmc Exp $
 % $Name:  $
 
 %- Remove insignificant whitespace:
-sufx=strtrim(char(sufx));
+%sufx=strtrim(char(sufx)); % <-- only with matlab-7 or more recent
+sufx=strrep(char(sufx),' ','');
 namfhd=[namF,'_head','.',sufx];
 namfil=[namF,'_Iter','.',sufx];
  fprintf(['read ',sufx,' :']);
@@ -128,7 +129,7 @@
 
 %- build time: 
 delT=0;
-if nIt > 1, delT=tim(2,1)-tim(1,1); delT=max(0,delT); end
+if nIt > 1, delT=(tim(nIt,1)-tim(1,1))/(nIt-1); delT=max(0,delT); end
 if delT > 0,
   delT=abs(frq)/delT;
   delta=delT-round(delT);