shallakefluxes()

subroutine clm_lake::shallakefluxes ( real(kind_lake), dimension(1), intent(in)  forc_t, real(kind_lake), dimension(1), intent(in)  forc_pbot, real(kind_lake), dimension(1), intent(in)  forc_psrf, real(kind_lake), dimension(1), intent(in)  forc_hgt, real(kind_lake), dimension(1), intent(in)  forc_hgt_q, real(kind_lake), dimension(1), intent(in)  forc_hgt_t, real(kind_lake), dimension(1), intent(in)  forc_hgt_u, real(kind_lake), dimension(1), intent(in)  forc_q, real(kind_lake), dimension(1), intent(in)  forc_u, real(kind_lake), dimension(1), intent(in)  forc_v, real(kind_lake), dimension(1), intent(in)  forc_lwrad, real(kind_lake), dimension(1), intent(in)  forc_snow, real(kind_lake), dimension(1), intent(in)  forc_rain, real(kind_lake), dimension(1), intent(inout)  t_grnd, real(kind_lake), dimension(1), intent(in)  h2osno, real(kind_lake), dimension(1), intent(in)  snowdp, real(kind_lake), dimension(1), intent(in)  sabg, real(kind_lake), dimension(1), intent(in)  lat, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in)  dz, real(kind_lake), dimension(1,nlevlake), intent(in)  dz_lake, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in)  t_soisno, real(kind_lake), dimension(1,nlevlake), intent(in)  t_lake, integer, dimension(1), intent(in)  snl, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in)  h2osoi_liq, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in)  h2osoi_ice, real(kind_lake), dimension(1), intent(in)  savedtke1, real(kind_lake), dimension(1), intent(out)  qflx_prec_grnd, real(kind_lake), dimension(1), intent(out)  qflx_evap_soi, real(kind_lake), dimension(1), intent(out)  qflx_evap_tot, real(kind_lake), dimension(1), intent(out)  eflx_sh_grnd, real(kind_lake), dimension(1), intent(out)  eflx_lwrad_out, real(kind_lake), dimension(1), intent(out)  eflx_lwrad_net, real(kind_lake), dimension(1), intent(out)  eflx_soil_grnd, real(kind_lake), dimension(1), intent(out)  eflx_sh_tot, real(kind_lake), dimension(1), intent(out)  eflx_lh_tot, real(kind_lake), dimension(1), intent(out)  eflx_lh_grnd, real(kind_lake), dimension(1), intent(out)  t_veg, real(kind_lake), dimension(1), intent(out)  t_ref2m, real(kind_lake), dimension(1), intent(out)  q_ref2m, real(kind_lake), dimension(1), intent(out)  taux, real(kind_lake), dimension(1), intent(out)  tauy, real(kind_lake), dimension(1), intent(out)  ram1, real(kind_lake), dimension(1), intent(out)  ws, real(kind_lake), dimension(1), intent(out)  ks, real(kind_lake), dimension(1), intent(out)  eflx_gnet, real(kind_lake), dimension(1), intent(out)  z0mg, real(kind_lake), dimension(1), intent(out)  ustar_out, character(len=*), intent(inout)  errmsg, integer, intent(inout)  errflg, real(kind_lake), intent(in)  xlat_d, real(kind_lake), intent(in)  xlon_d  )

private

Shallow lakes have variable depth, possible snow layers above, freezing & thawing of lake water, and soil layers with active temperature and gas diffusion below.

WARNING: This subroutine assumes lake columns have one and only one pft.

Parameters

[in]	forc_t	atmospheric temperature (Kelvin)
[in]	forc_pbot	atmospheric pressure (Pa)
[in]	forc_psrf	atmospheric surface pressure (Pa)
[in]	forc_hgt	atmospheric reference height (m)
[in]	forc_hgt_q	observational height of humidity [m]
[in]	forc_hgt_t	observational height of temperature [m]
[in]	forc_hgt_u	observational height of wind [m]
[in]	forc_q	atmospheric specific humidity (kg/kg)
[in]	forc_u	atmospheric wind speed in east direction (m/s)
[in]	forc_v	atmospheric wind speed in north direction (m/s)
[in]	forc_lwrad	downward infrared (longwave) radiation (W/m**2)
[in]	forc_snow	snow rate [mm/s]
[in]	forc_rain	rain rate [mm/s]
[in]	h2osno	snow water (mm H2O)
[in]	snowdp	snow height (m)
[in]	sabg	solar radiation absorbed by ground (W/m**2)
[in]	lat	latitude (radians)
[in]	dz	layer thickness for soil or snow (m)
[in]	dz_lake	layer thickness for lake (m)
[in]	t_soisno	soil (or snow) temperature (Kelvin)
[in]	t_lake	lake temperature (Kelvin)
[in]	snl	number of snow layers
[in]	h2osoi_liq	liquid water (kg/m2)
[in]	h2osoi_ice	ice lens (kg/m2)
[in]	savedtke1	top level eddy conductivity from previous timestep (W/m.K)
[in,out]	t_grnd	ground temperature (Kelvin)
[out]	ustar_out	friction velocity [m/s]
[out]	qflx_prec_grnd	water onto ground including canopy runoff [kg/(m2 s)]
[out]	qflx_evap_soi	soil evaporation (mm H2O/s) (+ = to atm)
[out]	qflx_evap_tot	qflx_evap_soi + qflx_evap_veg + qflx_tran_veg
[out]	eflx_sh_grnd	sensible heat flux from ground (W/m**2) [+ to atm]
[out]	eflx_lwrad_out	emitted infrared (longwave) radiation (W/m**2)
[out]	eflx_lwrad_net	net infrared (longwave) rad (W/m**2) [+ = to atm]
[out]	eflx_soil_grnd	soil heat flux (W/m**2) [+ = into soil]
[out]	eflx_sh_tot	total sensible heat flux (W/m**2) [+ to atm]
[out]	eflx_lh_tot	total latent heat flux (W/m8*2) [+ to atm]
[out]	eflx_lh_grnd	ground evaporation heat flux (W/m**2) [+ to atm]
[out]	t_veg	vegetation temperature (Kelvin)
[out]	t_ref2m	2 m height surface air temperature (Kelvin)
[out]	q_ref2m	2 m height surface specific humidity (kg/kg)
[out]	taux	wind (shear) stress: e-w (kg/m/s**2)
[out]	tauy	wind (shear) stress: n-s (kg/m/s**2)
[out]	ram1	aerodynamical resistance (s/m)
[out]	ws	surface friction velocity (m/s)
[out]	ks	coefficient passed to ShalLakeTemperature for calculation of decay of eddy diffusivity with depth
[out]	eflx_gnet	net heat flux into ground (W/m**2) Change the type variable to pass back to WRF.
[out]	z0mg	roughness length over ground, momentum (m(