global_ocean.90x40x15/diags_matlab/mit_timesteps.m

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

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