-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathheat_transport_timedev_plots.m
74 lines (62 loc) · 4.33 KB
/
heat_transport_timedev_plots.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
% plot inter-model comparisons of time-deviation contributions to transport, and decomposition in POP
path(path,'~/POP/')
path(path,'~/plotting_scripts/')
cd('/indopac/adelman/Global_mesoscale/POP/')
load('N_Atlantic_tmean_contribs_multi_models.mat')
heat_transport_timedev_tmean = heat_transport_zonmean_lowfreq_tmean + heat_transport_zonmean_highfreq_tmean + heat_transport_largescale_lowfreq_tmean + heat_transport_largescale_highfreq_tmean + heat_transport_mesoscale_lowfreq_tmean + heat_transport_mesoscale_highfreq_tmean + heat_transport_undef_lowfreq_tmean + heat_transport_undef_highfreq_tmean;
heat_transport_eddy_Treguier2017_interp = interp2_fft(lat_Treguier2017,[0 1],repmat(heat_transport_eddy_Treguier2017,[1 2]),(20:0.1:60)',0.5*ones([401 1]));
heat_transport_eddy_Treguier2017_interpfilt = bandpass_err_fcn(heat_transport_eddy_Treguier2017_interp,1,0.1,1/(3*40),1/(5*0.2),5,1,1,0,0);
fig1 = figure(1);
h = plot(lat_POP,(1e-15)*heat_transport_timedev_tmean,lat_JM2002,(1e-15)*heat_transport_eddy_JM2002,(20:0.1:60)',(1e-15)*heat_transport_eddy_Treguier2017_interpfilt);
leg = legend('POP','Jayne and Marotzke (2002)','Treguier et al. (2017)','location','northeast');
set(leg,'FontSize',12)
set(gca,'FontSize',12,'xlim',[min(lat_POP) max(lat_POP)],'xtick',(-80):10:80,'xgrid','on','ylim',[-0.2 0.4],'ytick',(-2):0.1:2)
% daspect([(aspect_ratio_baseline + (0.1*(max(lat_transects_vec) - min(lat_transects_vec)))) 1 1])
set(h(1),'Color',[0 0 0],'LineWidth',2)
set(h(2),'Color',[0.8 0 0],'LineWidth',1.5)
set(h(3),'Color',[0.2 0.2 0.8],'LineWidth',1.5)
hold on
line([min(lat_POP) max(lat_POP)],[0 0],[1 1],'Color','k','LineStyle','-')
hold off
xlabel('Latitude')
ylabel('Time-deviation contribution to time-mean MHT (PW)')
title({'Total time mean heat transport contribution from time deviation components'; ''; ''; ''; ''; ''},'FontSize',10)
saveas(fig1,'Heat_transport_timedev_N_Atlantic_POP_JM2002_Treguier2017.pdf')
close(fig1)
fig2 = figure(2);
close(figure(1))
h = plot(lat_POP,(1e-15)*(heat_transport_zonmean_highfreq_tmean + heat_transport_largescale_highfreq_tmean + heat_transport_mesoscale_highfreq_tmean + heat_transport_undef_highfreq_tmean),lat_Volkov2008,(1e-15)*heat_transport_eddy_Volkov2008);
leg = legend('POP','Volkov et al. (2008)','location','northeast');
set(leg,'FontSize',12)
set(gca,'FontSize',12,'xlim',[min(lat_POP) max(lat_POP)],'xtick',(-80):10:80,'xgrid','on','ylim',[-0.2 0.4],'ytick',(-0.2):0.1:0.4)
% daspect([(aspect_ratio_baseline + (0.1*(max(lat_transects_vec) - min(lat_transects_vec)))) 1 1])
set(h(1),'Color',[0 0 0],'LineWidth',2)
set(h(2),'Color',[0.5 0 0.5],'LineWidth',1.5)
hold on
line([min(lat_POP) max(lat_POP)],[0 0],[1 1],'Color','k','LineStyle','-')
hold off
xlabel('Latitude')
ylabel('High frequency contribution to time-mean MHT (PW)')
title({'Total time mean heat transport contribution from high frequency components'; ''; ''; ''; ''; ''},'FontSize',10)
saveas(fig2,'Heat_transport_highfreq_N_Atlantic_POP_Volkov2008.pdf')
close(fig2)
fig3 = figure(3);
close(figure(1))
h = plot(lat_POP,(1e-15)*(heat_transport_zonmean_lowfreq_tmean + heat_transport_zonmean_highfreq_tmean),lat_POP,(1e-15)*(heat_transport_largescale_lowfreq_tmean + heat_transport_largescale_highfreq_tmean),lat_POP,(1e-15)*(heat_transport_mesoscale_lowfreq_tmean + heat_transport_mesoscale_highfreq_tmean),lat_POP,(1e-15)*(heat_transport_undef_lowfreq_tmean + heat_transport_undef_highfreq_tmean),lat_POP,(1e-15)*heat_transport_timedev_tmean);
leg = legend('Overturning','Large scale','Mesoscale','Residual','Total','location','northeast');
set(leg,'FontSize',12)
set(gca,'FontSize',12,'xlim',[min(lat_POP) max(lat_POP)],'xtick',(-80):10:80,'xgrid','on','ylim',[-0.2 0.4],'ytick',(-0.2):0.1:0.4)
% daspect([(aspect_ratio_baseline + (0.1*(max(lat_transects_vec) - min(lat_transects_vec)))) 1 1])
set(h(1),'Color',[0.55 0 0.55],'LineWidth',2)
set(h(2),'Color',[0 0 0.8],'LineWidth',2)
set(h(3),'Color',[0.8 0 0],'LineWidth',2)
set(h(4),'Color',[0.8 0.6 0.4],'LineWidth',2)
set(h(5),'Color',[0 0 0],'LineWidth',2)
hold on
line([min(lat_POP) max(lat_POP)],[0 0],[1 1],'Color','k','LineStyle','-')
hold off
xlabel('Latitude')
ylabel('Time-deviation contribution to time-mean MHT (PW)')
title({'Components of time mean heat transport contribution from time deviation components'; ''; ''; ''; ''; ''},'FontSize',10)
saveas(fig3,'Heat_transport_timedev_N_Atlantic_POP_components.pdf')
close(fig3)