New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
albedo_cloud_correction.m in NEMO/trunk/doc/si3_doc/Figures/scripts – NEMO

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

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

Radiative transfer SI3 documentation

File size: 1.6 KB
Line 
1clear 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
10alpha_os = 0:0.001:1.;
11
12% SI3 correction
13f = ( -0.1010 * alpha_os.^2 + 0.1933 * alpha_os - 0.0148 );
14
15% clear sky sea ice albedo
16alpha_cs = alpha_os - f;
17
18% Grenfell and Perovich, Table 3 data
19alpha_GP04 = [ 0.114 0.151 0.141 0.266 0.367 0.309 0.332 0.492 0.678 0.928 ];
20f_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%--------------
25alpha_oce_cs = 0.05 / ( 1.1 * 0.40^1.4 + 0.15 ) % --- Briegleb and Ramanathan 1972
26alpha_oce_os = 0.06 %--- Payne 1972
27f_oce = alpha_oce_os - alpha_oce_cs
28
29%--------------
30% Plot
31%--------------
32plot(alpha_os,f,'k'); hold on
33zaddr = find( ( alpha_os >= alpha_GP04(3)  )   & ...
34              ( alpha_os <= alpha_GP04(10) ) );
35plot(alpha_os(zaddr),f(zaddr),'k', 'LineWidth', 3)
36plot(alpha_GP04(3:10),f_GP04(3:10),'ksq','MarkerFaceColor', 'k')
37plot(alpha_GP04(5),f_GP04(5),'rsq','MarkerFaceColor','r')
38text(alpha_GP04(5)+0.02,f_GP04(5),'dirty ice')
39plot(alpha_GP04(1:2),f_GP04(1:2),'rsq','MarkerFaceColor','r')
40text(alpha_GP04(1)+0.02,f_GP04(1),'wet tundra')
41text(alpha_GP04(2)+0.02,f_GP04(2),'dry tundra')
42plot(alpha_oce_os,f_oce,'bo','MarkerFaceColor', 'b')
43text(alpha_oce_os+0.02,f_oce,'ocean')
44plot( [ alpha_GP04(3) alpha_GP04(3) ], [ -0.06 0.08 ], 'k--' )
45plot( [ alpha_GP04(10) alpha_GP04(10) ], [ -0.06 0.08 ], 'k--' )
46
47xlabel('\alpha^{os}')
48ylabel('\alpha^{os}-\alpha^{cs}')
49
50set(gca,'fontsize', 16)
51set(gca, 'FontName', 'Helvetica LT Std')
Note: See TracBrowser for help on using the repository browser.