global_ocean.128x64x15/diags_matlab/mit_plotmeandrift.m

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