cbl_friction_vel_module

This MODULE contains the SUBROUTINE comp_friction_vel and two FUNCTIONS (psim and psis) needed to evaluate the friction velocity over each land point/tile given the wind speed and the current estimate of the Monin-Obukhov stability parameter $\xi$ . Outputs friction_vel, psim and psis are used in define_canopy.

Used by

Help

Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.
Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.

Functions

public function psim(zeta, mp, CPI_C) result(r)

Evaluates the integrated similarity function for momentum transfer, $\psi_m(\xi)$ .

Arguments

Type	Intent	Attributes		Name
real,	intent(in),	DIMENSION(mp)	::	zeta	current value of stability parameter $\xi$ (-)
integer,	intent(in)		::	mp	size of cable vector of land points (-)
real,	intent(in)		::	CPI_C	$\pi$ (-)

Return Value real, DIMENSION(mp)

public elemental function psis(zeta) result(r)

Evaluates the integrated similarity function for turbulent transfer of scalars, $\psi_s(\xi)$ .

Arguments

Type	Intent	Optional	Attributes		Name
real,	intent(in)			::	zeta	current value of stability parameter $\xi$ (-)

Return Value real

OUT $\psi_s (\xi)$

Subroutines

public subroutine comp_friction_vel(friction_vel, iter, mp, CVONK, CUMIN, CPI_C, zetar, zref_uv, zref_tq, z0m, ua)

This SUBROUTINE evaluates the value of the friction velocity $u_*$ as used during the iteration loop of the Monin-Obukhov (MO) similarity theory in define_canopy. $u_*$ is used when quantifying the stability parameter canopy%zetar in update_zetar, and when evaluating the resistance network components rt0 and rt1 in define_canopy. The friction velocity quantifies the magnitude of turbulent mixing that is occurring as well as the flux of momentum from the atmosphere to the land.

Arguments

Type	Intent		Name
real,	intent(out)	::	friction_vel(mp)	friction velocity (ms $^{-1}$ )
integer,	intent(in)	::	iter	MO iteration counter (-)
integer,	intent(in)	::	mp	size of cable vector of land points (-)
real,	intent(in)	::	CVONK	von Karman constant (-)
real,	intent(in)	::	CUMIN	minimum value of wind speed (ms $^{-1}$ )
real,	intent(in)	::	CPI_C	$\pi$ (-)
real,	intent(in)	::	zetar(mp,iter)	stability parameter - see update_zetar `canopy%zetar` (-)
real,	intent(in)	::	zref_uv(mp)	reference height for wind `rough%zref_uv` (m)
real,	intent(in)	::	zref_tq(mp)	reference height for temperature and humidity `rough%zref_tq` (m)
real,	intent(in)	::	z0m(mp)	roughness length `rough%z0m` (m)
real,	intent(in)	::	ua(mp)	wind speed `met%ua` (ms $^{-1}$ )

cbl_friction_vel_module Module

Contents

Functions

Subroutines

Used by

Functions

public function psim(zeta, mp, CPI_C) result(r)

Arguments

Return Value real, DIMENSION(mp)

public elemental function psis(zeta) result(r)

Arguments

Return Value real

Subroutines

public subroutine comp_friction_vel(friction_vel, iter, mp, CVONK, CUMIN, CPI_C, zetar, zref_uv, zref_tq, z0m, ua)

Arguments