--- MITgcm_contrib/jmc_script/read_StD.m	2008/05/29 23:47:13	1.2
+++ MITgcm_contrib/jmc_script/read_StD.m	2008/10/09 01:09:48	1.4
@@ -1,5 +1,5 @@
-function [nIt,rList,tim,vvA,listV]=read_StD(namF,sufx,listV);
-% [nIt,rList,tim,vvA,listV]=read_StD(namF,sufx,listV);
+function [nIt,rList,tim,vvA,listV,kList]=read_StD(namF,sufx,listV);
+% [nIt,rList,tim,vvA,listV,kList]=read_StD(namF,sufx,listV);
 %
 % read ASCII stat-Diags output files (after splitted by script extract_StD)
 %
@@ -11,12 +11,16 @@
 %  rList    = list of region number
 %  tim(:,1) = iterations number ; tim(:,2) = time in simulation
 %  listV    = list of fields
+%  kList    = list of levels number
 %  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.2 2008/05/29 23:47:13 jmc Exp $
+% $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/jmc_script/read_StD.m,v 1.4 2008/10/09 01:09:48 jmc Exp $
 % $Name:  $
 
+%- Remove insignificant whitespace:
+%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,' :']);