This MODULE contains the SUBROUTINE update_zetar needed to update
 the value of the stability parameter canopy%zetar=. 
This SUBROUTINE updates the value of the stability parameter 
 during the iteration loop of the Monin-Obukhov (MO) similarity theory in define_canopy.
 quantifies the role that the surface fluxes play in setting the 
 efficiency of turbulent transfer from the land to the atmosphere, and hence
 the aerodynamic component of the resistance network for those same surface 
 fluxes (an implicit problem which requires iteration to solve).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | mp | number of land points (-) | ||
| integer, | intent(out) | :: | iterplus | |||
| integer, | intent(in) | :: | NITER | number of MO-iterations (-) | ||
| real, | intent(out) | :: | canopy_zetar(mp,NITER) | stability parameter (-) | ||
| integer, | intent(in) | :: | iter | iteration counter (-) | ||
| integer, | intent(in) | :: | nrb | number of radiation bands (-) | ||
| real, | intent(in) | :: | CVONK | |||
| real, | intent(in) | :: | CGRAV | |||
| real, | intent(in) | :: | CCAPP | |||
| real, | intent(in) | :: | CLAI_THRESH | |||
| real, | intent(in) | :: | CZETmul | |||
| real, | intent(in) | :: | CZETPOS | |||
| real, | intent(in) | :: | CZETNEG | |||
| character(len=1), | intent(in) | :: | cable_user_soil_struc | name of soil model used | ||
| real, | intent(in) | :: | air_rho(mp) | air density (kg m) | ||
| real, | intent(in) | :: | met_tk(mp) | reference level air temperature (K) | ||
| real, | intent(in) | :: | met_fsd(mp,nrb) | downwelling shortwave (Wm) | ||
| real, | intent(in) | :: | rough_zref_tq(mp) | reference height for T and q (m) | ||
| real, | intent(in) | :: | rough_hruff(mp) | height of canopy (above snow) (m) | ||
| real, | intent(in) | :: | rough_term6a(mp) | term from ruff_resist (-) | ||
| real, | intent(in) | :: | rough_z0soilsn(mp) | roughness length of soil or snow (m) | ||
| real, | intent(in) | :: | canopy_vlaiw(mp) | canopy leaf area (mm) | ||
| real, | intent(out) | :: | canopy_zetash(mp,NITER) | stability parameter for soil (-) | ||
| real, | intent(in) | :: | canopy_us(mp) | friction velocity (ms) | ||
| real, | intent(in) | :: | canopy_fh(mp) | sensible heat flux (Wm) | ||
| real, | intent(in) | :: | canopy_fe(mp) | latent heat flux (Wm) | ||
| real, | intent(in) | :: | canopy_fhs(mp) | soil sensible heat flux (Wm) | ||
| real, | intent(in) | :: | canopy_fes(mp) | soil latent heat flux (Wm) |