source: NEMO/trunk/doc/si3_doc/Figures/scripts/albedo_hi.m @ 9996

Last change on this file since 9996 was 9996, checked in by vancop, 2 years ago

Radiative transfer SI3 documentation

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