tutorial_global_oce_latlon/diags_matlab/mit_timesteps.m

% m-file: mit_timesteps.m
% sets up the the variables tname and timesteps for the 4x4 global runs

% $Header:  $
% $Name:  $

deltat = mit_getparm('data','deltaTtracer');
if isempty(deltat)
  error('deltaTtracer is empty')
end

if strcmp(dname,'baseline.000')
  timesteps = [0:11]'*36e3; % one hundred years
  timesteps = [timesteps; [12:30]'*72e3]; % one hundred years (baseline experiment)
elseif strcmp(dname,'etopo5.000')
  timesteps = [0:30]'*72e3; % one hundred years (etopo5 experiments)
else
  if meanfields
    [dummy, timesteps] = rdmds('uVeltave',NaN);
    if isempty(timesteps)
      meanfields=0;
      error('There are no averaged fields, use meanfields=0')
    end
    timesteps = [0 timesteps];
  else
    [dummy, timesteps] = rdmds('U',NaN);
  end
  timesteps = timesteps';
  clear dummy
end

nt = length(timesteps);
kt = 1:nt;
tname = cell(nt,1);
for k = kt
  tname{k} = sprintf('%010i',timesteps(k));
end

%% alternative (unflexible way)
% $$$   timesteps = [0:11]'*36e3; % one hundred years
% $$$   timesteps = [timesteps; [12:30]'*72e3]; % one hundred years (baseline experiment)
% $$$ clear timesteps tname kt nt
% $$$ timesteps = [0:30]'*72e3; % one hundred years (etopo5 experiments)
% $$$ %tincr = 20; timesteps = [0:36]'*20; % one hundred years
% $$$ nt = length(timesteps);
% $$$ kt = 1:nt;
% $$$ tname = cell(nt,1);
% $$$ for k = kt
% $$$   tname{k} = sprintf('%010i',timesteps(k));
% $$$ end

% create a time axis
oneday = 3600*24;
onemonth =oneday*30;
oneyr=onemonth*12;
msg_spinup = dir('spinup.tim');
if isempty(msg_spinup)
  tim = deltat*timesteps';
else
  tim = load('spinup.tim');
  itim = find(diff(tim) == 0);
  tim(itim) = [];
end
% create reasonable unit;
if max(tim(:))/oneday <= 360
  tim = tim/oneday;
  timeunit = 'days';
  tuname = 'day';
elseif max(tim(:))/onemonth <= 120
  tim = tim/onemonth;
  timeunit = 'months';
  tuname = 'month';
else 
  tim = tim/oneyr;
  timeunit = 'yrs';
  tuname = 'year';
end

if ~exist('kmax','var')
  kmax=max(kt);
end
disp(['kmax = ' num2str(kmax) ', diplayed time = ' ...
      num2str(tim(kmax)) ' ' timeunit])

1	jmc	1.1	% m-file: mit_timesteps.m
2			% sets up the the variables tname and timesteps for the 4x4 global runs
3
4	jmc	1.3	% $Header: $
5	jmc	1.1	% $Name: $
6
7			deltat = mit_getparm('data','deltaTtracer');
8			if isempty(deltat)
9			error('deltaTtracer is empty')
10			end
11
12			if strcmp(dname,'baseline.000')
13			timesteps = [0:11]'*36e3; % one hundred years
14			timesteps = [timesteps; [12:30]'*72e3]; % one hundred years (baseline experiment)
15			elseif strcmp(dname,'etopo5.000')
16			timesteps = [0:30]'*72e3; % one hundred years (etopo5 experiments)
17			else
18			if meanfields
19			[dummy, timesteps] = rdmds('uVeltave',NaN);
20			if isempty(timesteps)
21			meanfields=0;
22			error('There are no averaged fields, use meanfields=0')
23			end
24			timesteps = [0 timesteps];
25			else
26			[dummy, timesteps] = rdmds('U',NaN);
27			end
28			timesteps = timesteps';
29			clear dummy
30			end
31
32			nt = length(timesteps);
33			kt = 1:nt;
34			tname = cell(nt,1);
35			for k = kt
36			tname{k} = sprintf('%010i',timesteps(k));
37			end
38
39			%% alternative (unflexible way)
40			% $$$ timesteps = [0:11]'*36e3; % one hundred years
41			% $$$ timesteps = [timesteps; [12:30]'*72e3]; % one hundred years (baseline experiment)
42			% $$$ clear timesteps tname kt nt
43			% $$$ timesteps = [0:30]'*72e3; % one hundred years (etopo5 experiments)
44			% $$$ %tincr = 20; timesteps = [0:36]'*20; % one hundred years
45			% $$$ nt = length(timesteps);
46			% $$$ kt = 1:nt;
47			% $$$ tname = cell(nt,1);
48			% $$$ for k = kt
49			% $$$ tname{k} = sprintf('%010i',timesteps(k));
50			% $$$ end
51
52			% create a time axis
53			oneday = 3600*24;
54			onemonth =oneday*30;
55			oneyr=onemonth*12;
56			msg_spinup = dir('spinup.tim');
57			if isempty(msg_spinup)
58			tim = deltat*timesteps';
59			else
60			tim = load('spinup.tim');
61			itim = find(diff(tim) == 0);
62			tim(itim) = [];
63			end
64			% create reasonable unit;
65			if max(tim(:))/oneday <= 360
66			tim = tim/oneday;
67			timeunit = 'days';
68			tuname = 'day';
69			elseif max(tim(:))/onemonth <= 120
70			tim = tim/onemonth;
71			timeunit = 'months';
72			tuname = 'month';
73			else
74			tim = tim/oneyr;
75			timeunit = 'yrs';
76			tuname = 'year';
77			end
78
79			if ~exist('kmax','var')
80			kmax=max(kt);
81			end
82			disp(['kmax = ' num2str(kmax) ', diplayed time = ' ...
83			num2str(tim(kmax)) ' ' timeunit])
84