MITgcm_contrib/mitgcm_tools/loadstate.m

function [STATE] = loadstate(varargin)
%loadstate()
%loadstate(DIRECTORY)
%loadstate(ITERATION)
%loadstate(DIRECTORY,ITERATION)
%
%Reads MITgcm parameter input files and output to create a STATE structure
%If ITERATION is not specified, the last model state is read by default.
%If DIRECTORY is not specified the current working directory is used.
%
%e.g.
% >> STATE=loadstate;
% >> STATE2=loadstate('/scratch/john/run2/');
% >> STATE3=loadstate(3000);
% >> STATE4=loadstate('/scratch/john/run2/',2460);
%
% Elements  (u,v,w,t,s,eta) correspond to snap-shots.
% Elements  (U,V,W,T,S,ETA) correspond to time averages.
%
%Written by adcroft@mit.edu, 2002
%$Header:

Dir='./';
iter=inf;
if nargin~=0
 for k = 1:nargin;
  if isstr(varargin{k})
   Dir=[varargin{k} '/'];
  else
   iter=varargin{k};
  end
 end
end

Its=[];
[STATE.t,its]=rdmds([Dir 'T'],iter); if isempty(Its); Its=its; end
[STATE.s,its]=rdmds([Dir 'S'],iter); if isempty(Its); Its=its; end
[STATE.u,its]=rdmds([Dir 'U'],iter); if isempty(Its); Its=its; end
[STATE.v,its]=rdmds([Dir 'V'],iter); if isempty(Its); Its=its; end
[STATE.w,its]=rdmds([Dir 'W'],iter); if isempty(Its); Its=its; end
[STATE.eta,its]=rdmds([Dir 'Eta'],iter); if isempty(Its); Its=its; end
[STATE.T,its]=rdmds([Dir 'Ttave'],iter); if isempty(Its); Its=its; end
[STATE.S,its]=rdmds([Dir 'Stave'],iter); if isempty(Its); Its=its; end
[STATE.U,its]=rdmds([Dir 'uVeltave'],iter); if isempty(Its); Its=its; end
[STATE.V,its]=rdmds([Dir 'vVeltave'],iter); if isempty(Its); Its=its; end
[STATE.W,its]=rdmds([Dir 'wVeltave'],iter); if isempty(Its); Its=its; end
[STATE.ETA,its]=rdmds([Dir 'ETAtave'],iter); if isempty(Its); Its=its; end
STATE.iters=Its;

datafile=[Dir 'data'];
dtClock = grepparameter(datafile,'deltatclock');
dtTracer = grepparameter(datafile,'deltattracer');
dtMom = grepparameter(datafile,'deltatmom');
dt = grepparameter(datafile,'deltat');
dt=[dt dtClock dtTracer dtMom]; dt=dt(1);
 
STATE.dt=dt;
STATE.time_secs=its*dt;
STATE.time_days=its*dt/86400;
STATE.time_years=its*dt/86400/360;
1	function [STATE] = loadstate(varargin)
2	%loadstate()
3	%loadstate(DIRECTORY)
4	%loadstate(ITERATION)
5	%loadstate(DIRECTORY,ITERATION)
6	%
7	%Reads MITgcm parameter input files and output to create a STATE structure
8	%If ITERATION is not specified, the last model state is read by default.
9	%If DIRECTORY is not specified the current working directory is used.
10	%
11	%e.g.
12	% >> STATE=loadstate;
13	% >> STATE2=loadstate('/scratch/john/run2/');
14	% >> STATE3=loadstate(3000);
15	% >> STATE4=loadstate('/scratch/john/run2/',2460);
16	%
17	% Elements (u,v,w,t,s,eta) correspond to snap-shots.
18	% Elements (U,V,W,T,S,ETA) correspond to time averages.
19	%
20	%Written by adcroft@mit.edu, 2002
21	%$Header:
22
23	Dir='./';
24	iter=inf;
25	if nargin~=0
26	for k = 1:nargin;
27	if isstr(varargin{k})
28	Dir=[varargin{k} '/'];
29	else
30	iter=varargin{k};
31	end
32	end
33	end
34
35	Its=[];
36	[STATE.t,its]=rdmds([Dir 'T'],iter); if isempty(Its); Its=its; end
37	[STATE.s,its]=rdmds([Dir 'S'],iter); if isempty(Its); Its=its; end
38	[STATE.u,its]=rdmds([Dir 'U'],iter); if isempty(Its); Its=its; end
39	[STATE.v,its]=rdmds([Dir 'V'],iter); if isempty(Its); Its=its; end
40	[STATE.w,its]=rdmds([Dir 'W'],iter); if isempty(Its); Its=its; end
41	[STATE.eta,its]=rdmds([Dir 'Eta'],iter); if isempty(Its); Its=its; end
42	[STATE.T,its]=rdmds([Dir 'Ttave'],iter); if isempty(Its); Its=its; end
43	[STATE.S,its]=rdmds([Dir 'Stave'],iter); if isempty(Its); Its=its; end
44	[STATE.U,its]=rdmds([Dir 'uVeltave'],iter); if isempty(Its); Its=its; end
45	[STATE.V,its]=rdmds([Dir 'vVeltave'],iter); if isempty(Its); Its=its; end
46	[STATE.W,its]=rdmds([Dir 'wVeltave'],iter); if isempty(Its); Its=its; end
47	[STATE.ETA,its]=rdmds([Dir 'ETAtave'],iter); if isempty(Its); Its=its; end
48	STATE.iters=Its;
49
50	datafile=[Dir 'data'];
51	dtClock = grepparameter(datafile,'deltatclock');
52	dtTracer = grepparameter(datafile,'deltattracer');
53	dtMom = grepparameter(datafile,'deltatmom');
54	dt = grepparameter(datafile,'deltat');
55	dt=[dt dtClock dtTracer dtMom]; dt=dt(1);
56
57	STATE.dt=dt;
58	STATE.time_secs=its*dt;
59	STATE.time_days=its*dt/86400;
60	STATE.time_years=its*dt/86400/360;