effect of temporal resolution

annual cycleの場合

multi-year simulationの場合

結論からいうとGOTMのサブブランチにあるWintonモデルはなんかおかしい。おそらく未完成なんだろう。

% dt = 60, 300, 600, 1200, 1800, 3600 seconds
% nsave = 1440, 288, 144, 72, 48, 24 timesteps
%
% out_ice: no eco. no T restoring.
% out_temp: no ice, no eco. no T restoring.
% out_eco: no ice. T restored every 60 sec.

directory = 'out_ice';
nfile = 6;
colid = 'kbgcrm';
figure;
for it = 1:nfile
    file=[directory '/' num2str(it) '.nc'];
subplot(221);
n = squeeze(ncread(file,'sst'));
plot(n,colid(it));
hold on;
title('SST');
subplot(222);
n = squeeze(ncread(file,'temp'));
n = squeeze(n(100,:));
plot(n,colid(it));
hold on;
title('Temperature at the first layer');
subplot(223);
n = squeeze(ncread(file,'uvic_eco_ph1'));
n = squeeze(n(100,:));
plot(n,colid(it));
hold on;
title('Small Phytoplankton at the first layer');
subplot(224);
n = squeeze(ncread(file,'uvic_eco_no3'));
n = squeeze(n(100,:));
plot(n,colid(it));
hold on;
title('Nitrate at the first layer');
end
legend('1 min','5 min','10 min','20 min','30 min','60 min');

figure;
for it = 1:nfile
file=[directory '/' num2str(it) '.nc'];
n=squeeze(ncread(file,'uvic_eco_no3'));
subplot(3,2,it);
[ax, h]=contourf(n,40);
shading flat
caxis([0 14]);
end
colorbar('location','west');