Plotting for poster.

Figure used for poster

%This is for FABM model output (not GOTM)!

clear all;

filename='resolute.nc';
dep=100;
avgdep=30;
day1=366;
colid='bcgmrk';

figure;

subplot(4,1,1);
n1=squeeze(ncread(filename,'sst'));
n2=squeeze(ncread(filename,'ice_hi'))*100;
xlen=day1:length(n1);
x=linspace(datenum('01-01-1990'),datenum('01-01-1994'),length(xlen));
[ax,h1,h2]=plotyy(x,n1(xlen),x,n2(xlen));
set(h1,'linewidth',2,'color','r');
set(h2,'linewidth',2,'color','b');
set(ax(1),'YLim',[-2 4],'ycolor','r')
set(ax(1),'YTick',-2:1:4)
set(ax(2),'YLim',[0 60],'ycolor','b')
set(ax(2),'YTick',0:10:60)
set(get(ax(1),'Ylabel'),'string','^oC')
set(get(ax(2),'Ylabel'),'string','cm')
axes(ax(2));
title('(a)');
axes(ax(1));
hold on;
datetick('x','mmmyy');
axes(ax(2));
hold on;
datetick('x','mmmyy');
%set(ax(1),'xtick',day1:365:length(xlen),'xlim',[day1 length(xlen)]);set(ax(2),'xtick',day1:365:length(xlen),'xlim',[day1 length(xlen)]);
grid on;
%xlabel('Days since January 1, 1990');
%legend('SST','Ice Thickness');

subplot(4,1,2);
parname={'uvic_eco_ph1' 'uvic_eco_ph2' 'uvic_eco_zo1' 'uvic_eco_zo2'};
n=squeeze(ncread(filename,'uvic_eco_ph1'));
n1=squeeze(mean(n(end-avgdep-1:end,xlen),1))*14/8.75;
n=squeeze(ncread(filename,'uvic_eco_zo1'));
n2=squeeze(mean(n(end-avgdep-1:end,xlen),1))*106/16*12;
[ax,h1,h2]=plotyy(x,n1,x,n2);
set(ax(1),'ycolor','r');
set(ax(2),'ycolor','b');
set(get(ax(1),'Ylabel'),'string','mg-Chl-\alpha/m^3');
set(get(ax(2),'Ylabel'),'string','mg-C/m^3');
set(h1,'linewidth',2,'color','r');
set(h2,'linewidth',2,'color','b');
set(ax(1),'YLim',[0 7],'ycolor','r')
set(ax(1),'YTick',0:1:7)
set(ax(2),'YLim',[0 112],'ycolor','b')
set(ax(2),'YTick',0:16:112)
%xlabel('Days since January 1, 1990');
n=squeeze(ncread(filename,'uvic_eco_ph2'));
n1=squeeze(mean(n(end-avgdep-1:end,xlen),1))*14/8.75;
axes(ax(1));
hold on;
hh=plot(x,n1,'--','color','r');
set(hh,'linewidth',2);
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
datetick('x','mmmyy');
n=squeeze(ncread(filename,'uvic_eco_zo2'));
n2=squeeze(mean(n(end-avgdep-1:end,xlen),1))*106/16*12;
axes(ax(2));
hold on;
hh=plot(x,n2,'--','color','b');
set(hh,'linewidth',2);
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
%set(ax(1),'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
%set(ax(2),'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
datetick('x','mmmyy');
grid on;
title('(b)');
%legend('Flagellates','Diatoms','Microzooplankton','Mesozooplankton','orientation','horizontal');

subplot(4,1,3);
n=squeeze(ncread(filename,'uvic_dms_dms'));
[ax, h]=contour(x,-dep:-1,n(1:dep,xlen));
set(h,'linewidth',2);
colorbar('location','west');
ylim([-50 -1]);
title('(c)');
%xlabel('Days since January 1, 1990');
datetick('x','mmmyy');
ylabel('Depth (m)');
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
grid on;

subplot(4,1,4);
n=squeeze(ncread(filename,'uvic_dms_dmspd'));
[ax, h]=contour(x,-dep:-1,n(1:dep,xlen));
set(h,'linewidth',2);
colorbar('location','west');
ylim([-50 -1]);
title('(d)');
%xlabel('Days since January 1, 1990');
datetick('x','mmmyy');
ylabel('Depth (m)');
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
grid on;