GFS Noah LSM Model

This is Noah LSM driver module, with the functionality of preparing variables to run Noah LSM gfssflx(), calling Noah LSM and post-processing variables for return to the parent model suite including unit conversion, as well as diagnotics calculation. More...

Detailed Description

Argument Table

local_name	standard_name	long_name	units	type	dimensions	kind	intent	optional
im	horizontal_loop_extent	horizontal loop extent	count	integer	()		in	False
km	soil_vertical_dimension	soil vertical layer dimension	count	integer	()		in	False
grav	gravitational_acceleration	gravitational acceleration	m s-2	real	()	kind_phys	in	False
cp	specific_heat_of_dry_air_at_constant_pressure	specific heat of dry air at constant pressure	J kg-1 K-1	real	()	kind_phys	in	False
hvap	latent_heat_of_vaporization_of_water_at_0C	latent heat of evaporation/sublimation	J kg-1	real	()	kind_phys	in	False
rd	gas_constant_dry_air	ideal gas constant for dry air	J kg-1 K-1	real	()	kind_phys	in	False
eps	ratio_of_dry_air_to_water_vapor_gas_constants	rd/rv	none	real	()	kind_phys	in	False
epsm1	ratio_of_dry_air_to_water_vapor_gas_constants_minus_one	(rd/rv) - 1	none	real	()	kind_phys	in	False
rvrdm1	ratio_of_vapor_to_dry_air_gas_constants_minus_one	(rv/rd) - 1 (rv = ideal gas constant for water vapor)	none	real	()	kind_phys	in	False
ps	surface_air_pressure	surface pressure	Pa	real	(horizontal_dimension)	kind_phys	in	False
t1	air_temperature_at_lowest_model_layer	1st model layer air temperature	K	real	(horizontal_dimension)	kind_phys	in	False
q1	water_vapor_specific_humidity_at_lowest_model_layer	1st model layer specific humidity	kg kg-1	real	(horizontal_dimension)	kind_phys	in	False
soiltyp	soil_type_classification	soil type at each grid cell	index	integer	(horizontal_dimension)		in	False
vegtype	vegetation_type_classification	vegetation type at each grid cell	index	integer	(horizontal_dimension)		in	False
sigmaf	bounded_vegetation_area_fraction	areal fractional cover of green vegetation bounded on the bottom	frac	real	(horizontal_dimension)	kind_phys	in	False
sfcemis	surface_longwave_emissivity_over_land_interstitial	surface lw emissivity in fraction over land (temporary use as interstitial)	frac	real	(horizontal_dimension)	kind_phys	in	False
dlwflx	surface_downwelling_longwave_flux_absorbed_by_ground_over_land	total sky surface downward longwave flux absorbed by the ground over land	W m-2	real	(horizontal_dimension)	kind_phys	in	False
dswsfc	surface_downwelling_shortwave_flux	total sky surface downward shortwave flux	W m-2	real	(horizontal_dimension)	kind_phys	in	False
snet	surface_net_downwelling_shortwave_flux	total sky surface net shortwave flux	W m-2	real	(horizontal_dimension)	kind_phys	in	False
delt	time_step_for_dynamics	dynamics time step	s	real	()	kind_phys	in	False
tg3	deep_soil_temperature	bottom soil temperature	K	real	(horizontal_dimension)	kind_phys	in	False
cm	surface_drag_coefficient_for_momentum_in_air_over_land	surface exchange coeff for momentum over land	none	real	(horizontal_dimension)	kind_phys	in	False
ch	surface_drag_coefficient_for_heat_and_moisture_in_air_over_land	surface exchange coeff heat & moisture over land	none	real	(horizontal_dimension)	kind_phys	in	False
prsl1	air_pressure_at_lowest_model_layer	Model layer 1 mean pressure	Pa	real	(horizontal_dimension)	kind_phys	in	False
prslki	ratio_of_exner_function_between_midlayer_and_interface_at_lowest_model_layer	Exner function ratio bt midlayer and interface at 1st layer	ratio	real	(horizontal_dimension)	kind_phys	in	False
zf	height_above_ground_at_lowest_model_layer	height above ground at 1st model layer	m	real	(horizontal_dimension)	kind_phys	in	False
land	flag_nonzero_land_surface_fraction	flag indicating presence of some land surface area fraction	flag	logical	(horizontal_dimension)		in	False
wind	wind_speed_at_lowest_model_layer	wind speed at lowest model level	m s-1	real	(horizontal_dimension)	kind_phys	in	False
slopetyp	surface_slope_classification	surface slope type at each grid cell	index	integer	(horizontal_dimension)		in	False
shdmin	minimum_vegetation_area_fraction	min fractional coverage of green veg	frac	real	(horizontal_dimension)	kind_phys	in	False
shdmax	maximum_vegetation_area_fraction	max fractnl cover of green veg (not used)	frac	real	(horizontal_dimension)	kind_phys	in	False
snoalb	upper_bound_on_max_albedo_over_deep_snow	upper bound on max albedo over deep snow	frac	real	(horizontal_dimension)	kind_phys	in	False
sfalb	surface_diffused_shortwave_albedo	mean surface diffused shortwave albedo	frac	real	(horizontal_dimension)	kind_phys	in	False
flag_iter	flag_for_iteration	flag for iteration	flag	logical	(horizontal_dimension)		in	False
flag_guess	flag_for_guess_run	flag for guess run	flag	logical	(horizontal_dimension)		in	False
lheatstrg	flag_for_canopy_heat_storage	flag for canopy heat storage parameterization	flag	logical	()		in	False
isot	soil_type_dataset_choice	soil type dataset choice	index	integer	()		in	False
ivegsrc	vegetation_type_dataset_choice	land use dataset choice	index	integer	()		in	False
bexppert	perturbation_of_soil_type_b_parameter	perturbation of soil type "b" parameter	frac	real	(horizontal_dimension)	kind_phys	in	False
xlaipert	perturbation_of_leaf_area_index	perturbation of leaf area index	frac	real	(horizontal_dimension)	kind_phys	in	False
vegfpert	perturbation_of_vegetation_fraction	perturbation of vegetation fraction	frac	real	(horizontal_dimension)	kind_phys	in	False
pertvegf	magnitude_of_perturbation_of_vegetation_fraction	magnitude of perturbation of vegetation fraction	frac	real	(5)	kind_phys	in	False
weasd	water_equivalent_accumulated_snow_depth_over_land	water equiv of acc snow depth over land	mm	real	(horizontal_dimension)	kind_phys	inout	False
snwdph	surface_snow_thickness_water_equivalent_over_land	water equivalent snow depth over land	mm	real	(horizontal_dimension)	kind_phys	inout	False
tskin	surface_skin_temperature_over_land_interstitial	surface skin temperature over land (temporary use as interstitial)	K	real	(horizontal_dimension)	kind_phys	inout	False
tprcp	nonnegative_lwe_thickness_of_precipitation_amount_on_dynamics_timestep_over_land	total precipitation amount in each time step over land	m	real	(horizontal_dimension)	kind_phys	inout	False
srflag	flag_for_precipitation_type	flag for snow or rain precipitation	flag	real	(horizontal_dimension)	kind_phys	inout	False
smc	volume_fraction_of_soil_moisture	volumetric fraction of soil moisture	frac	real	(horizontal_dimension, soil_vertical_dimension)	kind_phys	inout	False
stc	soil_temperature	soil temperature	K	real	(horizontal_dimension, soil_vertical_dimension)	kind_phys	inout	False
slc	volume_fraction_of_unfrozen_soil_moisture	volume fraction of unfrozen soil moisture	frac	real	(horizontal_dimension, soil_vertical_dimension)	kind_phys	inout	False
canopy	canopy_water_amount	canopy moisture content	kg m-2	real	(horizontal_dimension)	kind_phys	inout	False
trans	transpiration_flux	total plant transpiration rate	W m-2	real	(horizontal_dimension)	kind_phys	inout	False
tsurf	surface_skin_temperature_after_iteration_over_land	surface skin temperature after iteration over land	K	real	(horizontal_dimension)	kind_phys	inout	False
zorl	surface_roughness_length_over_land_interstitial	surface roughness length over land (temporary use as interstitial)	cm	real	(horizontal_dimension)	kind_phys	inout	False
sncovr1	surface_snow_area_fraction_over_land	surface snow area fraction	frac	real	(horizontal_dimension)	kind_phys	inout	False
qsurf	surface_specific_humidity_over_land	surface air saturation specific humidity over land	kg kg-1	real	(horizontal_dimension)	kind_phys	inout	False
gflux	upward_heat_flux_in_soil_over_land	soil heat flux over land	W m-2	real	(horizontal_dimension)	kind_phys	inout	False
drain	subsurface_runoff_flux	subsurface runoff flux	kg m-2 s-1	real	(horizontal_dimension)	kind_phys	inout	False
evap	kinematic_surface_upward_latent_heat_flux_over_land	kinematic surface upward latent heat flux over land	kg kg-1 m s-1	real	(horizontal_dimension)	kind_phys	inout	False
hflx	kinematic_surface_upward_sensible_heat_flux_over_land	kinematic surface upward sensible heat flux over land	K m s-1	real	(horizontal_dimension)	kind_phys	inout	False
ep	surface_upward_potential_latent_heat_flux_over_land	surface upward potential latent heat flux over land	W m-2	real	(horizontal_dimension)	kind_phys	inout	False
runoff	surface_runoff_flux	surface runoff flux	kg m-2 s-1	real	(horizontal_dimension)	kind_phys	inout	False
cmm	surface_drag_wind_speed_for_momentum_in_air_over_land	momentum exchange coefficient over land	m s-1	real	(horizontal_dimension)	kind_phys	inout	False
chh	surface_drag_mass_flux_for_heat_and_moisture_in_air_over_land	thermal exchange coefficient over land	kg m-2 s-1	real	(horizontal_dimension)	kind_phys	inout	False
evbs	soil_upward_latent_heat_flux	soil upward latent heat flux	W m-2	real	(horizontal_dimension)	kind_phys	inout	False
evcw	canopy_upward_latent_heat_flux	canopy upward latent heat flux	W m-2	real	(horizontal_dimension)	kind_phys	inout	False
sbsno	snow_deposition_sublimation_upward_latent_heat_flux	latent heat flux from snow depo/subl	W m-2	real	(horizontal_dimension)	kind_phys	inout	False
snowc	surface_snow_area_fraction	surface snow area fraction	frac	real	(horizontal_dimension)	kind_phys	inout	False
stm	soil_moisture_content	soil moisture content	kg m-2	real	(horizontal_dimension)	kind_phys	inout	False
snohf	snow_freezing_rain_upward_latent_heat_flux	latent heat flux due to snow and frz rain	W m-2	real	(horizontal_dimension)	kind_phys	inout	False
smcwlt2	volume_fraction_of_condensed_water_in_soil_at_wilting_point	soil water fraction at wilting point	frac	real	(horizontal_dimension)	kind_phys	inout	False
smcref2	threshold_volume_fraction_of_condensed_water_in_soil	soil moisture threshold	frac	real	(horizontal_dimension)	kind_phys	inout	False
wet1	normalized_soil_wetness	normalized soil wetness	frac	real	(horizontal_dimension)	kind_phys	inout	False
errmsg	ccpp_error_message	error message for error handling in CCPP	none	character	()	len=*	out	False
errflg	ccpp_error_flag	error flag for error handling in CCPP	flag	integer	()		out	False

GFS sfc_drv General Algorithm

Collaboration diagram for GFS Noah LSM Model:

Modules
	GFS sfcsub Module
	This module contains grib code for each parameter-used in subroutines sfccycle() and setrmsk().

Modules
module	namelist_soilveg

Functions/Subroutines
subroutine, public	lsm_noah::lsm_noah_init (me, isot, ivegsrc, nlunit,
	This subroutine contains the CCPP-compliant lsm_noah_init to initialize soil vegetation. More...
subroutine, public	lsm_noah::lsm_noah_run (im, km, grav, cp, hvap, rd, eps, epsm1, rvrdm1, ps, t1, q1, soiltyp, vegtype, sigmaf, sfcemis, dlwflx, dswsfc, snet, delt, tg3, cm, ch, prsl1, prslki, zf, land, wind, slopetyp, shdmin, shdmax, snoalb, sfalb, flag_iter, flag_guess, lheatstrg, isot, ivegsrc, bexppert, xlaipert, vegfpert, pertvegf,
subroutine	alcalc
	This subroutine calculates albedo including snow effect (0 -> 1). More...
subroutine	sfcdif
	This subroutine calculates surface layer exchange coefficients via iterative process(see Chen et al.(1997) [33]). More...
subroutine, public	set_soilveg_mod::set_soilveg (me, isot, ivet, nlunit)
	This subroutine initializes soil and vegetation. More...
subroutine	gfssflx (nsoil, couple, icein, ffrozp, dt, zlvl, sldpth, swdn, swnet, lwdn, sfcems, sfcprs, sfctmp, sfcspd, prcp, q2, q2sat, dqsdt2, th2, ivegsrc, vegtyp, soiltyp, slopetyp, shdmin, alb, snoalb, bexpp, xlaip, lheatstrg, tbot, cmc, t1, stc, smc, sh2o, sneqv, ch, cm, z0, nroot, shdfac, snowh, albedo, eta, sheat, ec, edir, et, ett, esnow, drip, dew, beta, etp, ssoil, flx1, flx2, flx3, runoff1, runoff2, runoff3, snomlt, sncovr, rc, pc, rsmin, xlai, rcs, rct, rcq, rcsoil, soilw, soilm, smcwlt, smcdry, smcref, smcmax)
	This is the entity of GFS Noah LSM model of physics subroutines. It is a soil/veg/snowpack land-surface model to update soil moisture, soil ice, soil temperature, skin temperature, snowpack water content, snowdepth, and all terms of the surface energy balance and surface water balance (excluding input atmospheric forcings of downward radiation and precipitation). More...
subroutine	canres
	This subroutine calculates canopy resistance which depends on incoming solar radiation, air temperature, atmospheric water vapor pressure deficit at the lowest model level, and soil moisture (preferably unfrozen soil moisture rather than total). More...
subroutine	csnow
	This subroutine calculates snow termal conductivity. More...
subroutine	nopac
	This subroutine calculates soil moisture and heat flux values and update soil moisture content and soil heat content values for the case when no snow pack is present. More...
subroutine	penman
	This subroutine calculates potential evaporation for the current point. various partial sums/products are also calculated and passed back to the calling routine for later use. More...
subroutine	redprm
	This subroutine internally sets default values or optionally read-in via namelist i/o, all soil and vegetation parateters requied for the execusion of the Noah LSM. More...
subroutine	snfrac
	This subroutine calculates snow fraction (0->1). More...
subroutine	snopac
	This subroutine calculates soil moisture and heat flux values and update soil moisture content and soil heat content values for the case when a snow pack is present. More...
subroutine	snow_new
	This subroutine calculates snow depth and densitity to account for the new snowfall. new values of snow depth & density returned. More...
subroutine	snowz0
	This subroutine calculates total roughness length over snow. More...
subroutine	tdfcnd
	This subroutine calculates thermal diffusivity and conductivity of the soil for a given point and time. More...
subroutine	evapo
	This subroutine calculates soil moisture flux. The soil moisture content (smc - a per unit volume measurement) is a dependent variable that is updated with prognostic equations. The canopy moisture content (cmc) is also updated. Frozen ground version: new states added: sh2o, and frozen ground correction factor, frzfact and paramter slope. More...
subroutine	shflx
	This subroutine updates the temperature state of the soil column based on the thermal diffusion equation and update the frozen soil moisture content based on the temperature. More...
subroutine	smflx
	This subroutine calculates soil moisture flux. The soil moisture content (smc - a per unit vulume measurement) is a dependent variable that is updated with prognostic equations. The canopy moisture content (cmc) is also updated. Frozen ground version: new states added: sh2o and frozen ground correction factor, frzx and parameter slope. More...
subroutine	snowpack
	This subroutine calculates compaction of a snowpack under conditions of increasing snow density, as obtained from an approximate solution of E. Anderson's differential equation (3.29),NOAA technical report NWS 19, by Victor Koren, 03/25/95. subroutine will return new values of snowh and sndens . More...
subroutine	devap
	This subrtouine calculates direct soil evaporation. More...
subroutine	frh2o
	This subroutine calculates amount of supercooled liquid soil water content if temperature is below 273.15K (t0). It requires Newton-type iteration to solve the nonlinear implicit equation given in eqn 17 of Koren et al.(1999) [100]. More...
subroutine	hrt
	This subroutine calculates the right hand side of the time tendency term of the soil thermal diffusion equation. Also to compute (prepare) the matrix coefficients for the tri-diagonal matrix of the implicit time scheme. More...
subroutine	hrtice
	This subroutine calculates the right hand side of the time tendency term of the soil thermal diffusion equation for sea-ice (ice = 1) or glacial-ice (ice). More...
subroutine	hstep
	This subroutine calculates/updates the soil temperature field. More...
subroutine	rosr12
	This subroutine inverts (solve) the tri-diagonal matrix problem. More...
subroutine	snksrc
	This subroutine calculates sink/source term of the termal diffusion equation. More...
subroutine	srt
	This subroutine calculates the right hand side of the time tendency term of the soil water diffusion equation. Also to compute (prepare) the matrix coefficients for the tri-diagonal matrix of the implicit time scheme. More...
subroutine	sstep
	This subroutine calculates/updates soil moisture content values and canopy moisture content values. More...
subroutine	tbnd
	This subroutine calculates temperature on the boundary of the layer by interpolation of the middle layer temperatures. More...
subroutine	tmpavg
	This subroutine calculates soil layer average temperature (tavg) in freezing/thawing layer using up, down, and middle layer temperature (tup, tdn, tm), where tup is at top boundary of layer, tdn is at bottom boundary of layer. tm is layer prognostic state temperature. More...
subroutine	transp
	This subroutine calculates transpiration for the veg class. More...
subroutine	wdfcnd
	This subroutine calculates soil water diffusivity and soil hydraulic conductivity. More...