-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathJC_Fig8_DeepOHCtrend.m
639 lines (528 loc) · 34.4 KB
/
JC_Fig8_DeepOHCtrend.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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
%Plot spatial map like Figure 5 in Garry et al., 2019
clc;clear;
load('XFL_HFspatialcolorMatrix_2K-6KallVersion10.mat')
load('Z_deepTruth_OHCtrend_ErrorV2.mat')
%load the small basin boundaries from PJ10;
load('bsn_boundaries_purkey.mat')
id=find(bsn_30(1,:)<=0);bsn_30(1,id)=bsn_30(1,id)+360;
y3=[-71:1:-55];
x1=[bsn_30(1,2:280)];%%start from 2:281
y1=[bsn_30(2,2:280)];
x2=[zeros(1,17),bsn_30(1,282:567)];
y2=[y3,bsn_30(2,282:567)];
bsn_30P1=[x1;y1];
bsn_30P2=[x2;y2];
%
bsn_33(1,:)=bsn_33(1,:) + 360;
%
id=find(bsn_28(1,:)<=0);bsn_28(1,id)=bsn_28(1,id)+360;
y3=[-32.5:1:5]; y3=sort(y3,'descend');
x1=[bsn_28(1,1:24),bsn_28(1,96:270)];
y1=[bsn_28(2,1:24),bsn_28(2,96:270)];
bsn_28P1=[x1;y1];
y3=[-32.5:1:5];
x2=[bsn_28(1,25:94),zeros(1,38)];
y2=[bsn_28(2,25:94),y3];
bsn_28P2=[x2;y2];
%
id=find(bsn_31(1,:)<=0);bsn_31(1,id)=bsn_31(1,id)+360;
y3=[-55:1:-33]; y3=sort(y3,'descend');
x1=[bsn_31(1,61:131)];
y1=[bsn_31(2,61:131)];
bsn_31P1=[x1;y1];
y3=[-55:1:-33];
x2=[bsn_31(1,1:59),zeros(1,23),bsn_31(1,133:156)];
y2=[bsn_31(2,1:59),y3,bsn_31(2,133:156)];
bsn_31P2=[x2;y2];
% load lon lat grid
load('lonlatDepth.mat')
lon(1:360)=lon(1:360)+360;
[xx,yy]=meshgrid(lon,lat);xx=xx';yy=yy';
% plot
% indian=[0.2 0.6 0.4]; so=[0.3 0.7 1]; pacific=[0.2 0.5 0.8]; alantic=[1 0.6 0.1];
level = 23; n = ceil(level/2);
cmap2 = [linspace(1, 1, n); linspace(1, 0, n); linspace(1, 0, n)]';
cmap1 = [linspace(0, 1, n); linspace(0, 1, n); linspace(1, 1, n)]';
cmap = [cmap1; cmap2(2:end, :)];
%%%%%%%%%%
%#Boudary layer thickness
BdyThick=0.5;
MarkerSize=5;
%%%%%%%%%%%%%%
figure
set(gcf,'Units','inches')
set(gcf,'position',[1,0, 8 21])
hh11=subplot(3,1,1)
m_proj('miller','lon',[0 360],'lat',[-80 75]);
hold on
m_pcolor(xx,yy,OHCL1);shading flat;
m_coast('patch',[.7 .7 .7],'edgecolor','none');
m_grid('tickdir','in','linewi',2,'fontsize',14,'xticklabels',[]);
colormap(gca,vivid(cmap, [.5 .5]));% colorbar
caxis([-1 1]);
% m_text(30,70,'Basin truth 2km-4km (W m^-^2)', 'fontweight','bold','fontsize',14)
m_text(10, 70, '(a)', 'fontweight','bold','fontsize',18)
title('Basin truth','fontsize',14)
% colormap('bluewhitered')
% h=colorbar('northoutside');
m_plot(bsn_1(1,:),bsn_1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_2(1,:),bsn_2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_3(1,:),bsn_3(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_4(1,:),bsn_4(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_5(1,:),bsn_5(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_6(1,:),bsn_6(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_20(1,:),bsn_20(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_27(1,:),bsn_27(2,:),'k','LineWidth',BdyThick)
%Southern Ocean
m_plot(bsn_7(1,:),bsn_7(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_16(1,:),bsn_16(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_30P1(1,:),bsn_30P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_30P2(1,:),bsn_30P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_33(1,:),bsn_33(2,:),'k','LineWidth',BdyThick)
%Pacific ocean
m_plot(bsn_8(1,:),bsn_8(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_9(1,:),bsn_9(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_10(1,:),bsn_10(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_11(1,:),bsn_11(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_12(1,:),bsn_12(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_13(1,:),bsn_13(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_14(1,:),bsn_14(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_17(1,:),bsn_17(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_18(1,:),bsn_18(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_19(1,:),bsn_19(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_21(1,:),bsn_21(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_22(1,:),bsn_22(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_32(1,:),bsn_32(2,:),'k','LineWidth',BdyThick)
%alantic
m_plot(bsn_15(1,:),bsn_15(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_23(1,:),bsn_23(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_24(1,:),bsn_24(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_25(1,:),bsn_25(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_26(1,:),bsn_26(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_28P1(1,:),bsn_28P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_28P2(1,:),bsn_28P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_29(1,:),bsn_29(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_31P1(1,:),bsn_31P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_31P2(1,:),bsn_31P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_1(1,20)-5,bsn_1(2,20)-5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_2(1,20)+10,bsn_2(2,20)-10,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_10(1,20)-8,bsn_10(2,20)-2,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_11(1,20)+7,bsn_11(2,20)+3.5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_32(1,20)-3,bsn_32(2,20)+1,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_21(1,20)+8,bsn_21(2,20)-5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_25(1,20)-4,bsn_25(2,20)+8,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
load('XFL_HF92_15_2K-6KallVersion11.mat')
HF_D=[HF2K_4Ktruth;HF2K_4KT;HF2K_4KP];
HF_D_SE=[HF2K_4Ktruth*0;HF2K_4KT_SE;HF2K_4KP_SE];
%$text sector
fontsize=10
%Indian
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(1)), char(177), sprintf('%.2f', ErrorL1(1)) )
m_text(30, 0, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(2)), char(177), sprintf('%.2f', ErrorL1(2)) )
m_text(42, 15, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(3)), char(177), sprintf('%.2f', ErrorL1(3)) )
m_text(40, -25, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(4)), char(177), sprintf('%.2f', ErrorL1(4)) )
m_text(60, -11, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(5)), char(177), sprintf('%.2f', ErrorL1(5)) )
m_text(90, -3, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(6)), char(177), sprintf('%.2f', ErrorL1(6)) )
m_text(50, -40, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(20)), char(177), sprintf('%.2f', ErrorL1(20)) )
m_text(98, -46, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(27)), char(177), sprintf('%.2f', ErrorL1(27)) )
m_text(10, -45, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$SO
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(30)), char(177), sprintf('%.2f', ErrorL1(30)) )
m_text(20, -65, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(30)), char(177), sprintf('%.2f', ErrorL1(30)) )
m_text(320, -70, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(7)), char(177), sprintf('%.2f', ErrorL1(7)) )
m_text(90, -63, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(16)), char(177), sprintf('%.2f', ErrorL1(16)) )
m_text(228, -65, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$pacific
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(18)), char(177), sprintf('%.2f', ErrorL1(18)) )
m_text(180, 25, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(14)), char(177), sprintf('%.2f', ErrorL1(14)) )
m_text(188, -48, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(8)), char(177), sprintf('%.2f', ErrorL1(8)) )
m_text(115, 25, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(32)), char(177), sprintf('%.2f', ErrorL1(32)) )
m_text(120, 40, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(21)), char(177), sprintf('%.2f', ErrorL1(21)) )
m_text(128, 52, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(22)), char(177), sprintf('%.2f', ErrorL1(22)) )
m_text(168, 57, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(10)), char(177), sprintf('%.2f', ErrorL1(10)) )
m_text(230, 10, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(11)), char(177), sprintf('%.2f', ErrorL1(11)) )
m_text(262, 2, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(12)), char(177), sprintf('%.2f', ErrorL1(12)) )
m_text(248, -15, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(13)), char(177), sprintf('%.2f', ErrorL1(13)) )
m_text(253, -30, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(19)), char(177), sprintf('%.2f', ErrorL1(19)) )
m_text(144, -45, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(17)), char(177), sprintf('%.2f', ErrorL1(17)) )
m_text(135, -15, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(9)), char(177), sprintf('%.2f', ErrorL1(9)) )
m_text(167, -23, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$Atlantic
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(26)), char(177), sprintf('%.2f', ErrorL1(26)) )
m_text(288, 37, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(29)), char(177), sprintf('%.2f', ErrorL1(29)) )
m_text(318, 28, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(24)), char(177), sprintf('%.2f', ErrorL1(24)) )
m_text(315, -24, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(28)))
m_text(347, -17, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( char(177), sprintf('%.2f', ErrorL1(28)) )
m_text(2, -10, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(23)), char(177), sprintf('%.2f', ErrorL1(23)) )
m_text(305, -45, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(33)), char(177), sprintf('%.2f', ErrorL1(33)) )
m_text(294, -56, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(31)) )
m_text(347, -43, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( char(177), sprintf('%.2f', ErrorL1(31)) )
m_text(2, -35, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(25)), char(177), sprintf('%.2f', ErrorL1(25)) )
m_text(255, 25, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTruth=strcat( sprintf('%.2f', HF2K_4Ktruth(15)), char(177), sprintf('%.2f', ErrorL1(15)) )
m_text(275, 15, StringTruth, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% saveas(gcf,'OHC2K_4K.png')
%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %$text sector
fontsize=8
COL='b';
% %Indian
% m_text(50, 0, 'S01','Color','g', 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% m_text(62, 15, 'S02', 'Color','g', 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(50, -15, 'S03','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(75, 0, 'S04', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(95, -20, 'S05','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(60, -45, 'S06','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(115, -40, 'S20', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(28, -32, 'S27','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% %$SO
m_text(50, -60, 'S30', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(340, -65, 'S30','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(115, -58, 'S07','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(245, -58, 'S16', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% %$pacific
m_text(190, 15, 'S18', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(200, -40, 'S14','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(127, 17, 'S08','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% m_text(129, 46, 'S32', 'Color','g', 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% m_text(145, 52, 'S21', 'Color','g', 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(181, 61, 'S22','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% m_text(260, 10, 'S10', 'Color','g', 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% m_text(272, 2, 'S11', 'Color','g', 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(262, -7, 'S12','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(272, -39, 'S13', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(147, -55, 'S19','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(153, -20, 'S17', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(169, -31, 'S09', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% %$Atlantic
m_text(302, 43, 'S26', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(320, 20, 'S29', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(328, -17, 'S24','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(347, -7, 'S28', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(0, -17, 'S28', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(330, -50, 'S23','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(295, -60, 'S33','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(348, -47, 'S31','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(0, -42, 'S31', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(260, 32, 'S25','Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
m_text(285, 8, 'S15', 'Color',COL, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%%
hh12=subplot(3,1,2)% truth for abyssal layer
m_proj('miller','lon',[0 360],'lat',[-80 75]);
hold on
m_pcolor(xx,yy,HFpL1);shading flat;
m_coast('patch',[.7 .7 .7],'edgecolor','none');
m_grid('tickdir','in','linewi',2,'fontsize',14,'xticklabels',[]);
colormap(gca,vivid(cmap, [.5 .5]));% colorbar
%colorbar
caxis([-1 1]);
% m_text(30,70,'Pseudo section 2km-4km (W m^-^2)', 'fontweight','bold','fontsize',14)
m_text(10, 70, '(b)', 'fontweight','bold','fontsize',18)
title('Pseudo section','fontsize',14)
m_plot(bsn_1(1,:),bsn_1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_2(1,:),bsn_2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_3(1,:),bsn_3(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_4(1,:),bsn_4(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_5(1,:),bsn_5(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_6(1,:),bsn_6(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_20(1,:),bsn_20(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_27(1,:),bsn_27(2,:),'k','LineWidth',BdyThick)
%Southern Ocean
m_plot(bsn_7(1,:),bsn_7(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_16(1,:),bsn_16(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_30P1(1,:),bsn_30P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_30P2(1,:),bsn_30P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_33(1,:),bsn_33(2,:),'k','LineWidth',BdyThick)
%Pacific ocean
m_plot(bsn_8(1,:),bsn_8(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_9(1,:),bsn_9(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_10(1,:),bsn_10(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_11(1,:),bsn_11(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_12(1,:),bsn_12(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_13(1,:),bsn_13(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_14(1,:),bsn_14(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_17(1,:),bsn_17(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_18(1,:),bsn_18(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_19(1,:),bsn_19(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_21(1,:),bsn_21(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_22(1,:),bsn_22(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_32(1,:),bsn_32(2,:),'k','LineWidth',BdyThick)
%alantic
m_plot(bsn_15(1,:),bsn_15(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_23(1,:),bsn_23(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_24(1,:),bsn_24(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_25(1,:),bsn_25(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_26(1,:),bsn_26(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_28P1(1,:),bsn_28P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_28P2(1,:),bsn_28P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_29(1,:),bsn_29(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_31P1(1,:),bsn_31P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_31P2(1,:),bsn_31P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_1(1,20)-5,bsn_1(2,20)-5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_2(1,20)+10,bsn_2(2,20)-10,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_10(1,20)-8,bsn_10(2,20)-2,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_11(1,20)+7,bsn_11(2,20)+3.5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_32(1,20)-3,bsn_32(2,20)+1,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_21(1,20)+8,bsn_21(2,20)-5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_25(1,20)-4,bsn_25(2,20)+8,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
%$text sector
fontsize=10
%Indian
StringP=strcat( sprintf('%.2f', HF2K_4KP(3)), char(177), sprintf('%.2f', HF2K_4KP_SE(3)) )
m_text(40, -25, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(4)), char(177), sprintf('%.2f', HF2K_4KP_SE(4)))
m_text(65, -2, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% StringP=strcat( char(177), sprintf('%.2f', HF2K_4KP_SE(4)) )
% m_text(73, -12, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(5)), char(177), sprintf('%.2f', HF2K_4KP_SE(5)) )
m_text(90, -18, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(6)), char(177), sprintf('%.2f', HF2K_4KP_SE(6)) )
m_text(47, -40, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(20)), char(177), sprintf('%.2f', HF2K_4KP_SE(20)) )
m_text(95, -40, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(27)), char(177), sprintf('%.2f', HF2K_4KP_SE(27)) )
m_text(10, -47, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$SO
StringP=strcat( sprintf('%.2f', HF2K_4KP(30)), char(177), sprintf('%.2f', HF2K_4KP_SE(30)) )
m_text(20, -60, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(30)), char(177), sprintf('%.2f', HF2K_4KP_SE(30)) )
m_text(320, -65, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(7)), char(177), sprintf('%.2f', HF2K_4KP_SE(7)) )
m_text(90, -58, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(16)), char(177), sprintf('%.2f', HF2K_4KP_SE(16)) )
m_text(220, -65, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$pacific
StringP=strcat( sprintf('%.2f', HF2K_4KP(18)), char(177), sprintf('%.2f', HF2K_4KP_SE(18)) )
m_text(170, 25, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(14)), char(177), sprintf('%.2f', HF2K_4KP_SE(14)) )
m_text(190, -40, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(8)), char(177), sprintf('%.2f', HF2K_4KP_SE(8)) )
m_text(115, 25, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(22)), char(177), sprintf('%.2f', HF2K_4KP_SE(22)) )
m_text(165, 58, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(12)), char(177), sprintf('%.2f', HF2K_4KP_SE(12)) )
m_text(252, -15, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(13)), char(177), sprintf('%.2f', HF2K_4KP_SE(13)) )
m_text(257, -33, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(19)), char(177), sprintf('%.2f', HF2K_4KP_SE(19)) )
m_text(145, -45, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(17)), char(177), sprintf('%.2f', HF2K_4KP_SE(17)) )
m_text(135, -15, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(9)), char(177), sprintf('%.2f', HF2K_4KP_SE(9)) )
m_text(165, -27, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$Atlantic
StringP=strcat( sprintf('%.2f', HF2K_4KP(26)), char(177), sprintf('%.2f', HF2K_4KP_SE(26)) )
m_text(282, 37, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(29)), char(177), sprintf('%.2f', HF2K_4KP_SE(29)) )
m_text(318, 28, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(24)), char(177), sprintf('%.2f', HF2K_4KP_SE(24)) )
m_text(315, -23, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(23)), char(177), sprintf('%.2f', HF2K_4KP_SE(23)) )
m_text(305, -45, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(33)), char(177), sprintf('%.2f', HF2K_4KP_SE(33)) )
m_text(294, -57, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(28)) )
m_text(345, -17, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( char(177), sprintf('%.2f', HF2K_4KP_SE(28)) )
m_text(2, -10, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(31)))
m_text(343, -43, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( char(177), sprintf('%.2f', HF2K_4KP_SE(31)) )
m_text(2, -35, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringP=strcat( sprintf('%.2f', HF2K_4KP(15)), char(177), sprintf('%.2f', HF2K_4KP_SE(15)) )
m_text(275, 15, StringP, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% saveas(gcf,'OHC2K_4K.png')
%%
hcolor1=colorbar('eastoutside','fontsize',14)
set(hcolor1,'position',[0.90 0.390 0.030 0.55])
clims = [-1 1];
hcolor1.Label.String = '(W m^-^2)';
hcolor1.Label.Position = [0.9, 1.15];
hcolor1.Label.Rotation = 360;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
load('XFL_HFspatialcolorMatrix_2K-6KallVersion10.mat')
% indian=[0.2 0.6 0.4]; so=[0.3 0.7 1]; pacific=[0.2 0.5 0.8]; alantic=[1 0.6 0.1];
mycolormap = customcolormap(linspace(0,1,11), {'#860454','#c51b7c','#dc75ab','#f0b7da','#ffdeef','#f8f7f7','#e5f4d9','#b9e084','#7fbc42','#4d921e','#276418'});
%%%%%%%%%%
%#Boudary layer thickness
BdyThick=0.5;
MarkerSize=5;
%%%%%%%%%%%%%%
hh13=subplot(3,1,3)
m_proj('miller','lon',[0 360],'lat',[-80 75]);
hold on
m_pcolor(xx,yy,HFpL1-OHCL1);shading flat;
m_coast('patch',[.7 .7 .7],'edgecolor','none');
m_grid('tickdir','in','linewi',2,'fontsize',14);
colormap(gca,mycolormap);% colorbar
caxis([-0.25 0.25]);
% m_text(20,70,'Total Bias (Pseudo section-Basin truth) 2km-4km (W m^-^2)', 'fontweight','bold','fontsize',13)
title('Total Bias','fontsize',14)
m_text(1, 70, '(c)', 'fontweight','bold','fontsize',18)
% colormap('bluewhitered')
% h=colorbar('northoutside');
m_plot(bsn_1(1,:),bsn_1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_2(1,:),bsn_2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_3(1,:),bsn_3(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_4(1,:),bsn_4(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_5(1,:),bsn_5(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_6(1,:),bsn_6(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_20(1,:),bsn_20(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_27(1,:),bsn_27(2,:),'k','LineWidth',BdyThick)
%Southern Ocean
m_plot(bsn_7(1,:),bsn_7(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_16(1,:),bsn_16(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_30P1(1,:),bsn_30P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_30P2(1,:),bsn_30P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_33(1,:),bsn_33(2,:),'k','LineWidth',BdyThick)
%Pacific ocean
m_plot(bsn_8(1,:),bsn_8(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_9(1,:),bsn_9(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_10(1,:),bsn_10(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_11(1,:),bsn_11(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_12(1,:),bsn_12(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_13(1,:),bsn_13(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_14(1,:),bsn_14(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_17(1,:),bsn_17(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_18(1,:),bsn_18(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_19(1,:),bsn_19(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_21(1,:),bsn_21(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_22(1,:),bsn_22(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_32(1,:),bsn_32(2,:),'k','LineWidth',BdyThick)
%alantic
m_plot(bsn_15(1,:),bsn_15(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_23(1,:),bsn_23(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_24(1,:),bsn_24(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_25(1,:),bsn_25(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_26(1,:),bsn_26(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_28P1(1,:),bsn_28P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_28P2(1,:),bsn_28P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_29(1,:),bsn_29(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_31P1(1,:),bsn_31P1(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_31P2(1,:),bsn_31P2(2,:),'k','LineWidth',BdyThick)
m_plot(bsn_1(1,20)-5,bsn_1(2,20)-5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_2(1,20)+10,bsn_2(2,20)-10,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_10(1,20)-8,bsn_10(2,20)-2,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_11(1,20)+7,bsn_11(2,20)+3.5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_32(1,20)-3,bsn_32(2,20)+1,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_21(1,20)+8,bsn_21(2,20)-5,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
m_plot(bsn_25(1,20)-4,bsn_25(2,20)+8,'-p','MarkerSize',MarkerSize,'MarkerFaceColor','r','MarkerEdgeColor','r');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
load('XFL_HF92_15_2K-6KallVersion11.mat')
HFbias=cat(1,HF2K_4KP-HF2K_4Ktruth,HF2K_4KT-HF2K_4Ktruth,HF2K_4KP-HF2K_4KT);
HF2K_4K_Btotal=HF2K_4KP-HF2K_4Ktruth;
BiasError=sqrt((HF2K_4KP_SE).^2+ErrorL1'.^2);
%$text sector
fontsize=10
%Indian
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(3)),char(177), sprintf('%.2f', BiasError(3)) )
m_text(40, -25, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(4)),char(177), sprintf('%.2f', BiasError(4)) )
m_text(60, -4, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(5)),char(177), sprintf('%.2f', BiasError(5)) )
m_text(90, -15, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(6)),char(177), sprintf('%.2f', BiasError(6)) )
m_text(47, -40, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(20)),char(177), sprintf('%.2f', BiasError(20)) )
m_text(100, -40, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(27)),char(177), sprintf('%.2f', BiasError(27)) )
m_text(10, -45, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$SO
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(30)),char(177), sprintf('%.2f', BiasError(30)) )
m_text(20, -60, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(30)),char(177), sprintf('%.2f', BiasError(30)) )
m_text(320, -65, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(7)),char(177), sprintf('%.2f', BiasError(7)) )
m_text(90, -58, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(16)),char(177), sprintf('%.2f', BiasError(16)) )
m_text(220, -65, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$pacific
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(18)),char(177), sprintf('%.2f', BiasError(18)) )
m_text(170, 25, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(14)),char(177), sprintf('%.2f', BiasError(14)) )
m_text(180, -40, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(8)),char(177), sprintf('%.2f', BiasError(8)) )
m_text(115, 25, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(22)),char(177), sprintf('%.2f', BiasError(22)) )
m_text(171, 58, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(12)),char(177), sprintf('%.2f', BiasError(12)) )
m_text(252, -15, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(13)),char(177), sprintf('%.2f', BiasError(13)) )
m_text(257, -33, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(19)),char(177), sprintf('%.2f', BiasError(19)) )
m_text(145, -45, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(17)),char(177), sprintf('%.2f', BiasError(17)) )
m_text(135, -15, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(9)),char(177), sprintf('%.2f', BiasError(9)) )
m_text(167, -23, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
%$Atlantic
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(26)),char(177), sprintf('%.2f', BiasError(26)) )
m_text(282, 33, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(29)),char(177), sprintf('%.2f', BiasError(29)) )
m_text(318, 28, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(24)),char(177), sprintf('%.2f', BiasError(24)) )
m_text(315, -23, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(28)) )
m_text(345, -17, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( char(177), sprintf('%.2f', BiasError(28)) )
m_text(2, -10, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(23)),char(177), sprintf('%.2f', BiasError(23)) )
m_text(305, -45, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(33)),char(177), sprintf('%.2f', BiasError(33)) )
m_text(294, -57, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(31)) )
m_text(345, -43, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat(char(177), sprintf('%.2f', BiasError(31)) )
m_text(2, -35, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
StringTotalB=strcat( sprintf('%.2f', HF2K_4K_Btotal(15)),char(177), sprintf('%.2f', BiasError(15)) )
m_text(275, 15, StringTotalB, 'fontweight','bold','fontsize',fontsize,'FontName', 'Times')
% saveas(gcf,'OHC2K_4K.png')
%%
%%
hcolor2=colorbar('eastoutside','fontsize',14)
set(hcolor2,'position',[0.90 0.040 0.030 0.25])
clims = [-0.25 0.25];
hcolor2.Label.String = '(W m^-^2)';
hcolor2.Label.Position = [0.9, 0.35];
hcolor2.Label.Rotation = 360;
%
set(hh11,'position',[0.090 0.6800 0.85 0.29])
set(hh12,'position',[0.0900 0.3600 0.85 0.29])
set(hh13,'position',[0.090 0.04 0.85 0.29])
% indian=[0.2 0.6
saveas(gcf,'JC_Fig7_OHC_spatialMapDeepVersion12_Revision.png')