tutorial_global_oce_latlon/diags_matlab/mit_plotmeandrift.m

% m-file: mit_plotmeandrift.m
% mean temperature and salinity drift

% $Header:  $
% $Name:  $

if ~isempty(msg_spinup)
  % extracted from the monitor output with a script by jmc
  global_mt = load('spinup.t_A');
  global_ms = load('spinup.s_A');
  global_mt(itim) = [];
  global_ms(itim) = [];
else
  % or computed `on the fly'
  global_mt = zeros(size(kt));
  global_ms = zeros(size(kt));
  for k=kt;
    if meanfields
      % read the snapshots
      tk=rdmds('T',timesteps(k));
      sk=rdmds('S',timesteps(k));
      if strcmp(grd.buoyancy,'OCEANICP')
        % turn fields upside down
        tk=tk(:,:,end:-1:1).*grd.cmask;
        sk=sk(:,:,end:-1:1).*grd.cmask;
      end
    else
      % snapshots are already available
      if iscell(t)
        tk = t{k};
      else
        tk = t(:,:,:,k);
      end
      if iscell(s)
        sk = s{k};
      else
        sk = s(:,:,:,k);
      end
    end
    global_mt(k)= nansum(tk(:).*grd.volc(:))/nansum(grd.volc(:));
    global_ms(k)= nansum(sk(:).*grd.volc(:))/nansum(grd.volc(:));
  end
end
figure
subplot(2,1,1);
plot(tim,global_mt)
title('Drift of Mean Temperature'); xlabel(['Time [' timeunit ']']); ylabel('T [degC]')
subplot(2,1,2);
plot(tim,global_ms)
title('Drift of Mean Salinity'); xlabel(['Time [' timeunit ']']); ylabel('S [PSU]')
suptitle(['experiment ' dname])

1	% m-file: mit_plotmeandrift.m
2	% mean temperature and salinity drift
3
4	% $Header: $
5	% $Name: $
6
7	if ~isempty(msg_spinup)
8	% extracted from the monitor output with a script by jmc
9	global_mt = load('spinup.t_A');
10	global_ms = load('spinup.s_A');
11	global_mt(itim) = [];
12	global_ms(itim) = [];
13	else
14	% or computed `on the fly'
15	global_mt = zeros(size(kt));
16	global_ms = zeros(size(kt));
17	for k=kt;
18	if meanfields
19	% read the snapshots
20	tk=rdmds('T',timesteps(k));
21	sk=rdmds('S',timesteps(k));
22	if strcmp(grd.buoyancy,'OCEANICP')
23	% turn fields upside down
24	tk=tk(:,:,end:-1:1).*grd.cmask;
25	sk=sk(:,:,end:-1:1).*grd.cmask;
26	end
27	else
28	% snapshots are already available
29	if iscell(t)
30	tk = t{k};
31	else
32	tk = t(:,:,:,k);
33	end
34	if iscell(s)
35	sk = s{k};
36	else
37	sk = s(:,:,:,k);
38	end
39	end
40	global_mt(k)= nansum(tk(:).*grd.volc(:))/nansum(grd.volc(:));
41	global_ms(k)= nansum(sk(:).*grd.volc(:))/nansum(grd.volc(:));
42	end
43	end
44	figure
45	subplot(2,1,1);
46	plot(tim,global_mt)
47	title('Drift of Mean Temperature'); xlabel(['Time [' timeunit ']']); ylabel('T [degC]')
48	subplot(2,1,2);
49	plot(tim,global_ms)
50	title('Drift of Mean Salinity'); xlabel(['Time [' timeunit ']']); ylabel('S [PSU]')
51	suptitle(['experiment ' dname])
52