/[MITgcm]/MITgcm_contrib/jmc_script/read_StD.m
ViewVC logotype

Diff of /MITgcm_contrib/jmc_script/read_StD.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by jmc, Thu May 29 23:47:13 2008 UTC revision 1.5 by jmc, Wed May 14 16:53:04 2014 UTC
# Line 1  Line 1 
1  function [nIt,rList,tim,vvA,listV]=read_StD(namF,sufx,listV);  function [nIt,rList,tim,vvA,listV,kList]=read_StD(namF,sufx,listV);
2  % [nIt,rList,tim,vvA,listV]=read_StD(namF,sufx,listV);  % [nIt,rList,tim,vvA,listV,kList]=read_StD(namF,sufx,listV);
3  %  %
4  % read ASCII stat-Diags output files (after splitted by script extract_StD)  % read ASCII stat-Diags output files (after splitted by script extract_StD)
5  %  %
# Line 11  function [nIt,rList,tim,vvA,listV]=read_ Line 11  function [nIt,rList,tim,vvA,listV]=read_
11  %  rList    = list of region number  %  rList    = list of region number
12  %  tim(:,1) = iterations number ; tim(:,2) = time in simulation  %  tim(:,1) = iterations number ; tim(:,2) = time in simulation
13  %  listV    = list of fields  %  listV    = list of fields
14    %  kList    = list of levels number
15  %  vvA      = 5 dims output array:  %  vvA      = 5 dims output array:
16  %           ( kLev, time_rec, region_rec, [ave,std,min,max,vol], var_rec )  %           ( kLev, time_rec, region_rec, [ave,std,min,max,vol], var_rec )
17    
18  % $Header$  % $Header$
19  % $Name$  % $Name$
20    
21    %- Remove insignificant whitespace:
22    %sufx=strtrim(char(sufx)); % <-- only with matlab-7 or more recent
23    sufx=strrep(char(sufx),' ','');
24  namfhd=[namF,'_head','.',sufx];  namfhd=[namF,'_head','.',sufx];
25  namfil=[namF,'_Iter','.',sufx];  namfil=[namF,'_Iter','.',sufx];
26   fprintf(['read ',sufx,' :']);   fprintf(['read ',sufx,' :']);
# Line 125  if frq > 0, dIt=1./max(1,dIt); else dIt= Line 129  if frq > 0, dIt=1./max(1,dIt); else dIt=
129    
130  %- build time:  %- build time:
131  delT=0;  delT=0;
132  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
133  if delT > 0,  if delT > 0,
134    delT=abs(frq)/delT;    delT=abs(frq)/delT;
135    delta=delT-round(delT);    delta=delT-round(delT);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22