-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathJC_Fig5_XFL_Section_Profiles.m
105 lines (80 loc) · 2.93 KB
/
JC_Fig5_XFL_Section_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
% clc;clear
run XFL_Section_Profiles_Version10.m
%Atlantic Ocean+Southern Ocean
load('lonlatDepth.mat','depth');
D=-depth(38:end);
%%
Lfont=2
Tempmin=-10
Tempmax=10
LineW=2
FontS=17
AX.F=10
%
fig=figure
set(gcf,'Units','inches')
set(gcf,'position',[0,0, 12 6])
%%
%
sp21=subplot(1,4,1)
errorbar(BetaTI05*100,D,StdTI05*100,'horizontal','Linewidth',Lfont)
hold on; grid on
title('I05')
ylim([-6000,-2000]);
xlim([Tempmin Tempmax])
text(min(xlim)+1, min(ylim), sprintf('(a)', 1), 'Horiz','left', 'Vert','bottom', 'FontSize',FontS,'Fontweight','bold')
% errorbar(BetaFI04*100,D,StdFI04*100,'horizontal','Linewidth',Lfont)
% errorbar(BetaPPI04*100,D,StdPPI04*100,'horizontal','Linewidth',Lfont)
errorbar(BetaPI05*100,D,StdPI05*100,'horizontal','Linewidth',Lfont)
set(gca,'FontSize',14);
%%legend
lgnd=legend('Section truth','Pesudo section','Location','Southeast');
legend boxoff
set(lgnd,'color','none')
ylabel('Depth (m)','FontSize',14)
sp22=subplot(1,4,2)
errorbar(BetaTP21*100,D,StdTP21*100,'horizontal','Linewidth',Lfont)
hold on; grid on
title('P21')
ylim([-6000,-2000]);
xlim([Tempmin Tempmax])
text(min(xlim)+1, min(ylim), sprintf('(b)', 1), 'Horiz','left', 'Vert','bottom', 'FontSize',FontS,'Fontweight','bold')
% errorbar(BetaFA02*100,D,StdFA02*100,'horizontal','Linewidth',Lfont)
% errorbar(BetaPPA02*100,D,StdPPA02*100,'horizontal','Linewidth',Lfont)
errorbar(BetaPP21*100,D,StdPP21*100,'horizontal','Linewidth',Lfont)
ylabel('Depth (m)','FontSize',14)
set(gca,'FontSize',14);
%
xlabel(' Potential Temperature Trend (10^-^3 ^oC year^-^1)','fontsize',16)
%
sp23=subplot(1,4,3)
errorbar(BetaTA16N*100,D,StdTA16N*100,'horizontal','Linewidth',Lfont)
hold on; grid on
title('A16N')
ylim([-6000,-2000]);
xlim([Tempmin Tempmax])
text(min(xlim)+1, min(ylim), sprintf('(c)', 1), 'Horiz','left', 'Vert','bottom', 'FontSize',FontS,'Fontweight','bold')
% errorbar(BetaFA01*100,D,StdFA01*100,'horizontal','Linewidth',Lfont)
% errorbar(BetaPPA01*100,D,StdPPA01*100,'horizontal','Linewidth',Lfont)
errorbar(BetaPA16N*100,D,StdPA16N*100,'horizontal','Linewidth',Lfont)
set(gca,'FontSize',14);
%
%
sp24=subplot(1,4,4)
errorbar(BetaTS04P*100,D,StdTS04P*100,'horizontal','Linewidth',Lfont)
hold on; grid on
title('S04P')
ylim([-6000,-2000]);
xlim([Tempmin Tempmax])
text(min(xlim)+1, min(ylim), sprintf('(d)', 1), 'Horiz','left', 'Vert','bottom', 'FontSize',FontS,'Fontweight','bold')
% errorbar(BetaFA12*100,D,StdFA12*100,'horizontal','Linewidth',Lfont)
% errorbar(BetaPPA12*100,D,StdPPA12*100,'horizontal','Linewidth',Lfont)
errorbar(BetaPS04P*100,D,StdPS04P*100,'horizontal','Linewidth',Lfont)
set(gca,'FontSize',14);
% ylabel('Depth (m)','FontSize',14)
%
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_Fig5_XFL_Section_ProfilesV11.png')