source:
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/doc/latex/SI3/Figures/scripts/albedo_hi.m
@
10419
Last change on this file since 10419 was 9996, checked in by vancop, 6 years ago | |
---|---|
File size: 1.5 KB |
Line | |
---|---|
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 |
Note: See TracBrowser
for help on using the repository browser.