module_radiation_surface

This module sets up surface albedo for sw radiation and surface emissivity for lw radiation. More...

Detailed Description

This module sets up surface albedo for sw radiation and surface emissivity for lw radiation.

Version

NCEP-Radiation_surface v5.1 Nov 2012

Collaboration diagram for module_radiation_surface:

Variables
integer, parameter, public	module_radiation_surface::imxems = 360
	num of longitude points in global emis-type map
integer, parameter, public	module_radiation_surface::jmxems = 180
	num of latitude points in global emis-type map
real(kind=kind_phys), parameter	module_radiation_surface::f_zero = 0.0
real(kind=kind_phys), parameter	module_radiation_surface::f_one = 1.0
real(kind=kind_phys), parameter	module_radiation_surface::rad2dg = 180.0 / con_pi
integer, dimension(:,:), allocatable	module_radiation_surface::idxems
	global surface emissivity index array
integer	module_radiation_surface::iemslw = 0
	global surface emissivity contrl flag set up in 'sfc_init'
subroutine, public	module_radiation_surface::sfc_init
	This subroutine is the initialization program for surface radiation related quantities (albedo, emissivity, etc.) More...
subroutine, public	module_radiation_surface::setalb
	This subroutine computes four components of surface albedos (i.e., vis-nir, direct-diffused) according to control flag ialbflg. 1) climatological surface albedo scheme (Briegleb 1992 [5]) 2) MODIS retrieval based scheme from Boston univ. More...
subroutine, public	module_radiation_surface::setemis
	This subroutine computes surface emissivity for LW radiation. More...

Function/Subroutine Documentation

subroutine, public module_radiation_surface::setalb

(

)

This subroutine computes four components of surface albedos (i.e., vis-nir, direct-diffused) according to control flag ialbflg.
1) climatological surface albedo scheme (Briegleb 1992 [5])
2) MODIS retrieval based scheme from Boston univ.

Parameters

slmsk	(IMAX), sea(0),land(1),ice(2) mask on fcst model grid
snowf	(IMAX), snow depth water equivalent in mm
sncovr	(IMAX), snow cover over land
snoalb	(IMAX), maximum snow albedo over land (for deep snow)
zorlf	(IMAX), surface roughness in cm
coszf	(IMAX), cosin of solar zenith angle
tsknf	(IMAX), ground surface temperature in K
tairf	(IMAX), lowest model layer air temperature in K
hprif	(IMAX), topographic sdv in m — for ialbflg=0 climtological albedo scheme —
alvsf	(IMAX), 60 degree vis albedo with strong cosz dependency
alnsf	(IMAX), 60 degree nir albedo with strong cosz dependency
alvwf	(IMAX), 60 degree vis albedo with weak cosz dependency
alnwf	(IMAX), 60 degree nir albedo with weak cosz dependency — for ialbflg=1 MODIS based land albedo scheme —
alvsf	(IMAX), visible black sky albedo at zenith 60 degree
alnsf	(IMAX), near-ir black sky albedo at zenith 60 degree
alvwf	(IMAX), visible white sky albedo
alnwf	(IMAX), near-ir white sky albedo
facsf	(IMAX), fractional coverage with strong cosz dependency
facwf	(IMAX), fractional coverage with weak cosz dependency
fice	(IMAX), sea-ice fraction
tisfc	(IMAX), sea-ice surface temperature
IMAX	array horizontal dimension
sfcalb	(IMAX,NF_ALBD), mean sfc albedo ( :, 1) - near ir direct beam albedo ( :, 2) - near ir diffused albedo ( :, 3) - uv+vis direct beam albedo ( :, 4) - uv+vis diffused albedo

General Algorithm

1. If use climatological albedo scheme:
   • Modified snow albedo scheme - units convert to m (originally snowf in mm; zorlf in cm)
   • Calculate diffused sea surface albedo
   • Calculate diffused snow albedo.
   • Calculate direct snow albedo.
   • Calculate direct sea surface albedo.
2. If use modis based albedo for land area:
   • Calculate snow cover input directly for land model, no conversion needed.
   • Calculate diffused sea surface albedo.
   • Calculate diffused snow albedo, land area use input max snow albedo.
   • Calculate direct snow albedo.
   • Calculate direct sea surface albedo, use fanglin's zenith angle treatment.

Definition at line 309 of file radiation_surface.f.

References physcons::con_t0c, physcons::con_tice, physcons::con_ttp, and physparam::ialbflg.

Referenced by module_radiation_driver::grrad().

Here is the caller graph for this function:

subroutine, public module_radiation_surface::setemis

(

)

This subroutine computes surface emissivity for LW radiation.

Parameters

xlon	(IMAX), longitude in radiance, ok for both 0->2pi or -pi -> +pi ranges
xlat	(IMAX), latitude in radiance, default to pi/2 -> -pi/2 range, otherwise see in-line comment
slmsk	(IMAX), sea(0),land(1),ice(2) mask on fcst model grid
snowf	(IMAX), snow depth water equivalent in mm
sncovr	(IMAX), snow cover over land
zorlf	(IMAX), surface roughness in cm
tsknf	(IMAX), ground surface temperature in K
tairf	(IMAX), lowest model layer air temperature in K
hprif	(IMAX), topographic standard deviation in m
IMAX	array horizontal dimension
sfcemis	(IMAX), surface emissivity

General Algorithm

1. Set sfcemis default to 1.0 or by surface type and condition.
2. Check for snow covered area.

Definition at line 641 of file radiation_surface.f.

References physparam::ialbflg, idxems, iemslw, imxems, and jmxems.

Referenced by module_radiation_driver::grrad().

Here is the caller graph for this function:

subroutine, public module_radiation_surface::sfc_init

(

)

This subroutine is the initialization program for surface radiation related quantities (albedo, emissivity, etc.)

Parameters

me	print control flag

General Algorithm

• Initialization of surface albedo section
  physparam::ialbflg
  • = 0: using climatology surface albedo scheme for SW
  • = 1: using MODIS based land surface albedo for SW
• Initialization of surface emissivity section
  physparam::iemsflg
  • = 0: fixed SFC emissivity at 1.0
  • = 1: input SFC emissivity type map from "semis_file"

Definition at line 134 of file radiation_surface.f.

References physparam::ialbflg, idxems, physparam::iemsflg, iemslw, imxems, jmxems, and physparam::semis_file.

Referenced by module_radiation_driver::radinit().

Here is the caller graph for this function: