PUBLIC ROUTINES
call compute_ocean_albedo (ocean, coszen, albedo, lat)
INPUT
ocean: logical flag; = true if ocean point
[logical, dimension(:,:)]
coszen: cosine of zenith angle (radians)
[real, dimension(:,:)]
lat latitude (radians)
[real, dimension(:,:)]
INPUT/OUTPUT
albedo albedo of ocean
[real, dimension(:,:)]
NOTES
All arguments dimensioned (longitude, latitude)
NAMELIST
&ocean_albedo_nml
logical :: ocean_albedo_option
ocean_albedo_option = 1 : Albedos used by GFDL Experimental Prediction Group
during 1980s and 1990s; source not currently documented.
(tabulated dependence on zenith angle)
ocean_albedo_option = 2 : Albedos used by GFDL Climate Dynamics Group during
1980s and 1990s; source not currently documented.
(tabulated dependence on latitude)
ocean_albedo_option = 3 : simple analytic dependence on zenith angle
used by J. E. Taylor, et. al.,
QJRMS, 1996, Vol. 122, 839-861
albedo = 0.037/[1.1*(cos(Z)**1.4) + 0.15]
ocean_albedo_option = 4 : constant uniform albedo.
Set by namelist variable const_alb
real :: const_alb (See ocean_albedo_option = 4)