1 | clear all; close all; |
---|
2 | |
---|
3 | % Script to plot the cloud fraction dependency of albedo correction |
---|
4 | |
---|
5 | %------------------ |
---|
6 | % Sea Ice |
---|
7 | %------------------ |
---|
8 | |
---|
9 | % overcast albedo |
---|
10 | alpha_os = 0:0.001:1.; |
---|
11 | |
---|
12 | % SI3 correction |
---|
13 | f = ( -0.1010 * alpha_os.^2 + 0.1933 * alpha_os - 0.0148 ); |
---|
14 | |
---|
15 | % clear sky sea ice albedo |
---|
16 | alpha_cs = alpha_os - f; |
---|
17 | |
---|
18 | % Grenfell and Perovich, Table 3 data |
---|
19 | alpha_GP04 = [ 0.114 0.151 0.141 0.266 0.367 0.309 0.332 0.492 0.678 0.928 ]; |
---|
20 | f_GP04 = [ -0.010 -0.022 0.012 0.027 0.020 0.036 0.037 0.056 0.071 0.077 ]; |
---|
21 | |
---|
22 | %-------------- |
---|
23 | % Ocean |
---|
24 | %-------------- |
---|
25 | alpha_oce_cs = 0.05 / ( 1.1 * 0.40^1.4 + 0.15 ) % --- Briegleb and Ramanathan 1972 |
---|
26 | alpha_oce_os = 0.06 %--- Payne 1972 |
---|
27 | f_oce = alpha_oce_os - alpha_oce_cs |
---|
28 | |
---|
29 | %-------------- |
---|
30 | % Plot |
---|
31 | %-------------- |
---|
32 | plot(alpha_os,f,'k'); hold on |
---|
33 | zaddr = find( ( alpha_os >= alpha_GP04(3) ) & ... |
---|
34 | ( alpha_os <= alpha_GP04(10) ) ); |
---|
35 | plot(alpha_os(zaddr),f(zaddr),'k', 'LineWidth', 3) |
---|
36 | plot(alpha_GP04(3:10),f_GP04(3:10),'ksq','MarkerFaceColor', 'k') |
---|
37 | plot(alpha_GP04(5),f_GP04(5),'rsq','MarkerFaceColor','r') |
---|
38 | text(alpha_GP04(5)+0.02,f_GP04(5),'dirty ice') |
---|
39 | plot(alpha_GP04(1:2),f_GP04(1:2),'rsq','MarkerFaceColor','r') |
---|
40 | text(alpha_GP04(1)+0.02,f_GP04(1),'wet tundra') |
---|
41 | text(alpha_GP04(2)+0.02,f_GP04(2),'dry tundra') |
---|
42 | plot(alpha_oce_os,f_oce,'bo','MarkerFaceColor', 'b') |
---|
43 | text(alpha_oce_os+0.02,f_oce,'ocean') |
---|
44 | plot( [ alpha_GP04(3) alpha_GP04(3) ], [ -0.06 0.08 ], 'k--' ) |
---|
45 | plot( [ alpha_GP04(10) alpha_GP04(10) ], [ -0.06 0.08 ], 'k--' ) |
---|
46 | |
---|
47 | xlabel('\alpha^{os}') |
---|
48 | ylabel('\alpha^{os}-\alpha^{cs}') |
---|
49 | |
---|
50 | set(gca,'fontsize', 16) |
---|
51 | set(gca, 'FontName', 'Helvetica LT Std') |
---|