1 | clear all, close all |
---|
2 | |
---|
3 | % plots SI3 albedo dependence versus ice thickness |
---|
4 | |
---|
5 | z1_href_pnd = 1. / 0.05 |
---|
6 | z1_c1 = 1. / ( log(1.5) - log(0.05) ); |
---|
7 | z1_c2 = 1. / 0.05; |
---|
8 | z1_c3 = 1. / 0.02 |
---|
9 | z1_c4 = 1. / 0.03 |
---|
10 | |
---|
11 | rn_alb_idry = 0.60; |
---|
12 | rn_alb_imlt = 0.50; |
---|
13 | rn_alb_sdry = 0.85 |
---|
14 | rn_alb_smlt = 0.75 |
---|
15 | rn_alb_oce = 0.06; |
---|
16 | rn_alb_dpnd = 0.27; |
---|
17 | |
---|
18 | h_i = 0:0.001:3; |
---|
19 | |
---|
20 | %--- Bare ice albedo (for hi > 150cm) |
---|
21 | zalb_ice = h_i; |
---|
22 | zalb_ice(:) = rn_alb_idry; |
---|
23 | |
---|
24 | %--- Bare ice albedo (for hi < 150cm) |
---|
25 | zaddr = find( ( 0.05 < h_i ) & ( h_i <= 1.5 ) ); % 5cm < hi < 150cm |
---|
26 | zalb_ice(zaddr) = zalb_ice(zaddr) + ( 0.18 - zalb_ice(zaddr) ) * z1_c1 .* ( log(1.5) - log(h_i(zaddr)) ); |
---|
27 | zaddr = find( h_i <= 0.05 ); |
---|
28 | zalb_ice(zaddr) = rn_alb_oce + ( 0.18 - rn_alb_oce ) * z1_c2 * h_i(zaddr); |
---|
29 | |
---|
30 | figure; |
---|
31 | |
---|
32 | |
---|
33 | subplot(1,3,1) |
---|
34 | |
---|
35 | plot(h_i,zalb_ice,'k','LineWidth', 3) |
---|
36 | xlabel('h_i(m)') |
---|
37 | ylabel('\alpha') |
---|
38 | ylim([0.1 0.9]) |
---|
39 | set(gca,'fontsize', 16) |
---|
40 | set(gca, 'FontName', 'Helvetica LT Std') |
---|
41 | |
---|
42 | |
---|
43 | |
---|
44 | %--- Snow covered albedo |
---|
45 | |
---|
46 | h_s = 0:0.001:0.6; |
---|
47 | zalb_ice = 0.6 |
---|
48 | zalb_snw = rn_alb_sdry - ( rn_alb_sdry - zalb_ice ) * exp( -h_s * z1_c3 ) |
---|
49 | |
---|
50 | subplot(1,3,2) |
---|
51 | plot(h_s,zalb_snw,'k','LineWidth', 3) |
---|
52 | xlabel('h_s(m)') |
---|
53 | ylabel('\alpha') |
---|
54 | ylim([0.1 0.9]) |
---|
55 | set(gca,'fontsize', 16) |
---|
56 | set(gca, 'FontName', 'Helvetica LT Std') |
---|
57 | |
---|
58 | %--- Ponded ice albedo |
---|
59 | zalb_ice = 0.5 |
---|
60 | h_pnd = 0:0.001:0.6; |
---|
61 | zalb_pnd = rn_alb_dpnd - ( rn_alb_dpnd - zalb_ice ) * exp( -h_pnd * z1_href_pnd ) |
---|
62 | |
---|
63 | subplot(1,3,3) |
---|
64 | plot(h_pnd,zalb_pnd,'k','LineWidth', 3) |
---|
65 | xlabel('h_{pnd}(m)') |
---|
66 | ylabel('\alpha') |
---|
67 | ylim([0.1 0.9]) |
---|
68 | set(gca,'fontsize', 16) |
---|
69 | set(gca, 'FontName', 'Helvetica LT Std') |
---|
70 | |
---|
71 | |
---|