phasechange_lake()

subroutine clm_lake::phasechange_lake ( integer, dimension(1), intent(in)  snl, real(kind_lake), dimension(1), intent(inout)  h2osno, 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(inout)  t_soisno, 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,nlevlake), intent(inout)  lake_icefrac, real(kind_lake), dimension(1,nlevlake), intent(inout)  t_lake, real(kind_lake), dimension(1), intent(inout)  snowdp, real(kind_lake), dimension(1), intent(out)  qflx_snomelt, real(kind_lake), dimension(1), intent(out)  eflx_snomelt, integer, dimension(1,-nlevsnow+1:nlevsoil), intent(out)  imelt, real(kind_lake), dimension(lbc:ubc,-nlevsnow+1:nlevsoil), intent(inout)  cv, real(kind_lake), dimension (lbc:ubc,1:nlevlake), intent(inout)  cv_lake, real(kind_lake), dimension(lbc:ubc), intent(out)  lhabs  )

private

The specific heats are assumed constant. Potential cycling errors resulting from this assumption will be trapped at the end of ShalLakeTemperature.

Parameters

[in]	snl	number of snow layers
[in,out]	h2osno	snow water (mm H2O)
[in]	dz	layer thickness (m)
[in]	dz_lake	lake layer thickness (m)
[in,out]	snowdp	snow height (m)
[in,out]	t_soisno	soil temperature (Kelvin)
[in,out]	h2osoi_liq	liquid water (kg/m2)
[in,out]	h2osoi_ice	ice lens (kg/m2)
[in,out]	lake_icefrac	mass fraction of lake layer that is frozen
[in,out]	t_lake	lake temperature (Kelvin)
[out]	qflx_snomelt	snow melt (mm H2O /s)
[out]	eflx_snomelt	snow melt heat flux (W/m**2)
[out]	imelt	flag for melting (=1), freezing (=2), Not=0 (new)
[in,out]	cv	heat capacity [J/(m2 K)]
[in,out]	cv_lake	heat capacity [J/(m2 K)]
[out]	lhabs	total per-column latent heat abs. (J/m^2)