-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathJC_Fig6_XFL_Sector_Profiles.m
126 lines (106 loc) · 4.13 KB
/
JC_Fig6_XFL_Sector_Profiles.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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
% clc;clear
run XFL_Section_Profiles_Version10.m
%Atlantic Ocean+Southern Ocean
load('lonlatDepth.mat','depth');
D=-depth(38:end);
fig=figure
set(gcf,'Units','inches')
set(gcf,'position',[0,0, 12 6])
%%
%Atlantic Ocean+Southern Ocean
load('lonlatDepth.mat'); clear lon lat
D=-depth(38:end);
load('TTrend_13Layers_33Sectors.mat')
%%
Lfont=2
Tempmin=-5
Tempmax=5
LineW=2
FontS=17
AX.F=10
%%
load('TrendSector06Version10.mat')
sp21=subplot(1,4,1)
plot(TrendTheta(:,6),D,'k', 'Linewidth',LineW)
hold on; grid on;box on
title('Ind-S06','FontSize',FontS)
ylim([-6000,-2000]);
xlim([-5 5])
text(min(xlim)+1, min(ylim), sprintf('(a)', 1), 'Horiz','left', 'Vert','bottom', 'FontSize',FontS,'FontWeight','bold')
% e1=errorbar(TrendTS06*100,D,SETS06*100,'horizontal','Linewidth',Lfont)
% e1.Color=[80/255,29/255,138/255];
% errorbar(TrendFS06*100,D,SEFS06*100,'horizontal','Linewidth',Lfont)
% errorbar(TrendPPS06*100,D,SEPPS06*100,'horizontal','Linewidth',Lfont)
e2=errorbar(TrendPS06*100,D,SEPS06*100,'horizontal','Linewidth',Lfont)
e2.Color=[229/255,87/255,9/255];
set(gca,'FontSize',14);
ylabel('Depth (m)','FontSize',14)
% xlabel('PT Trend (10^-^3 ^oC year^-^1)','fontsize',14)
lgnd=legend('ECCO truth','Pesudo section','Location','Southeast');
legend boxoff
set(lgnd,'color','none')
%%
load('TrendSector14Version10.mat')
sp22=subplot(1,4,2)
plot(TrendTheta(:,14),D,'k','Linewidth',LineW)
hold on; grid on;box on
title('Pac-S14','FontSize',FontS)
ylim([-6000,-2000]);
xlim([-2 2])
text(min(xlim)+0.5, min(ylim), sprintf('(b)', 1), 'Horiz','left', 'Vert','bottom', 'FontSize',FontS,'Fontweight','bold')
% e1=errorbar(TrendTS14*100,D,SETS14*100,'horizontal','Linewidth',Lfont)
% e1.Color=[80/255,29/255,138/255];
% errorbar(TrendFS14*100,D,SEFS14*100,'horizontal','Linewidth',Lfont)
% errorbar(TrendPPS14*100,D,SEPPS14*100,'horizontal','Linewidth',Lfont)
e2=errorbar(TrendPS14*100,D,SEPS14*100,'horizontal','Linewidth',Lfont)
e2.Color=[229/255,87/255,9/255];
set(gca,'FontSize',14);
xlabel(' Potential Temperature Trend (10^-^3 ^oC year^-^1)','fontsize',16)
%%
load('TrendSector26Version10.mat')
sp23=subplot(1,4,3)
plot(TrendTheta(:,26),D,'k', 'Linewidth',LineW)
hold on; grid on
title('Atl-S26','FontSize',FontS)
ylim([-6000,-2000]);
xlim([-5 5])
text(min(xlim)+1, min(ylim), sprintf('(c)', 5), 'Horiz','left', 'Vert','bottom','FontSize',FontS,'Fontweight','bold')
% e1=errorbar(TrendTS26*100,D,SETS26*100,'horizontal','Linewidth',Lfont)
% e1.Color=[80/255,29/255,138/255];
% errorbar(TrendFS24*100,D,SEFS24*100,'horizontal','Linewidth',Lfont)
% errorbar(TrendPPS24*100,D,SEPPS24*100,'horizontal','Linewidth',Lfont)
e2=errorbar(TrendPS26*100,D,SEPS26*100,'horizontal','Linewidth',Lfont)
e2.Color=[229/255,87/255,9/255];
set(gca,'FontSize',14);
% ylabel('Depth (m)','FontSize',14)
% xlabel('PT Trend (10^-^3 ^oC year^-^1)','fontsize',14)
%%
load('TrendSector30Version10.mat')
sp24=subplot(1,4,4)
plot(TrendTheta(:,30),D,'k', 'Linewidth',LineW)
hold on; grid on
title('SO-S30','FontSize',FontS)
ylim([-6000,-2000]);
xlim([0 10])
text(min(xlim)+1, min(ylim), sprintf('(d)', 1), 'Horiz','left', 'Vert','bottom', 'FontSize',FontS,'Fontweight','bold')
% e1=errorbar(TrendTS30*100,D,SETS30*100,'horizontal','Linewidth',Lfont)
% e1.Color=[80/255,29/255,138/255];
% errorbar(TrendFS16*100,D,SEFS16*100,'horizontal','Linewidth',Lfont)
% errorbar(TrendPPS16*100,D,SEPPS16*100,'horizontal','Linewidth',Lfont)
e2=errorbar(TrendPS30*100,D,SEPS30*100,'horizontal','Linewidth',Lfont)
e2.Color=[229/255,87/255,9/255];
set(gca,'FontSize',14);
% ylabel('Depth (m)','FontSize',14)
% xlabel('PT Trend (10^-^3 ^oC year^-^1)','fontsize',14)
%%legend
%%
% han=axes(fig,'visible','off');
% han.XLabel.Visible='on';
% han.YLabel.Visible='on';
% ylabel(han,'Depth (m)','FontName', 'Times','fontsize',15);
% xlabel(han,'Temperature trend (m ^oC year^-^1)','FontName', 'Times','fontsize',15);
set(sp21,'position', [0.08 0.11 0.19 0.8])
set(sp22,'position', [0.32 0.11 0.19 0.8])
set(sp23,'position', [0.56 0.11 0.19 0.8])
set(sp24,'position', [0.80 0.11 0.19 0.8])
saveas(gcf,'JC_Fig6_XFL_Sector_ProfilesV6.png')