lakemain()

subroutine clm_lake::lakemain ( 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)  prec, real(kind_lake), dimension(1), intent(in)  sabg, real(kind_lake), dimension(1), intent(in)  lat, real(kind_lake), dimension(1,nlevlake), intent(in)  z_lake, real(kind_lake), dimension(1,nlevlake), intent(in)  dz_lake, real(kind_lake), dimension(1), intent(in)  lakedepth, logical, dimension(1), intent(in)  do_capsnow, real(kind_lake), dimension(1), intent(inout)  h2osno, real(kind_lake), dimension(1), intent(inout)  snowdp, integer, dimension(1), intent(inout)  snl, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(inout)  z, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(inout)  dz, real(kind_lake), dimension(1,-nlevsnow+0:nlevsoil), intent(inout)  zi, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(inout)  h2osoi_vol, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(inout)  h2osoi_liq, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(inout)  h2osoi_ice, real(kind_lake), dimension(1), intent(inout)  t_grnd, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(inout)  t_soisno, real(kind_lake), dimension(1,nlevlake), intent(inout)  t_lake, real(kind_lake), dimension(1), intent(inout)  savedtke1, real(kind_lake), dimension(1,nlevlake), intent(inout)  lake_icefrac, real(kind_lake), dimension(1), intent(out)  eflx_lwrad_net, real(kind_lake), dimension(1), intent(out)  eflx_gnet, 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)  t_ref2m, real(kind_lake), dimension(1), intent(out)  q_ref2m, real(kind_lake), intent(in)  dtime, real(kind_lake), dimension(1,nlevsoil), intent(in)  watsat, real(kind_lake), dimension(1,nlevsoil), intent(in)  tksatu, real(kind_lake), dimension(1,nlevsoil), intent(in)  tkmg, real(kind_lake), dimension(1,nlevsoil), intent(in)  tkdry, real(kind_lake), dimension(1,nlevsoil), intent(in)  csol, 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)  z0mg, real(kind_lake), dimension(1), intent(out)  ustar_out, character(*), intent(inout)  errmsg, integer, intent(inout)  errflg, real(kind_lake), intent(in)  xlat_d, real(kind_lake), intent(in)  xlon_d  )

private

Parameters

[in]	dtime	timestep
[in]	xlon_d	grid location for debugging
[in]	forc_t	atmospheric temperature (Kelvin)
[in]	forc_pbot	atm bottom level 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]	prec	snow or rain rate [mm/s]
[in]	sabg	solar radiation absorbed by ground (W/m**2)
[in]	lat	latitude (radians)
[in]	z_lake	layer depth for lake (m)
[in]	dz_lake	layer thickness for lake (m)
[out]	ustar_out	friction velocity [m/s]
[in]	lakedepth	column lake depth (m) tep(in),hydro(in)
[in]	do_capsnow	true => do snow capping
[in]	watsat	volumetric soil water at saturation (porosity)
[in]	tksatu	thermal conductivity, saturated soil [W/m-K]
[in]	tkmg	thermal conductivity, soil minerals [W/m-K]
[in]	tkdry	thermal conductivity, dry soil (W/m/Kelvin)
[in]	csol	heat capacity, soil solids (J/m**3/Kelvin)
[in,out]	h2osoi_vol	volumetric soil water (0<=h2osoi_vol<=watsat)[m3/m3]
[in,out]	t_grnd	ground temperature (Kelvin)
[in,out]	h2osno	snow water (mm H2O)
[in,out]	snowdp	snow height (m)
[in,out]	z	layer depth for snow & soil (m)
[in,out]	dz	layer thickness for soil or snow (m)
[in,out]	t_soisno	soil (or snow) temperature (Kelvin)
[in,out]	t_lake	lake temperature (Kelvin)
[in,out]	snl	number of snow layers
[in,out]	h2osoi_liq	liquid water (kg/m2)
[in,out]	h2osoi_ice	ice lens (kg/m2)
[in,out]	savedtke1	top level eddy conductivity from previous timestep (W/m.K)
[in,out]	zi	interface level below a "z" level (m)
[in,out]	lake_icefrac	mass fraction of lake layer that is frozen
[out]	eflx_gnet	net heat flux into ground (W/m**2)
[out]	eflx_lwrad_net	net infrared (longwave) rad (W/m**2) [+ = to atm]
[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]	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) for calculation of decay of eddy diffusivity with depth Change the type variable to pass back to WRF.
[out]	z0mg	roughness length over ground, momentum (m(