Hi,
I want to know the calculation of vegetation emissivity, ground emissivity, ground surface temperature and vegetation temperature in the code. However, when I use "grep -r emv_patch t_veg_patch t_grnd_col " in the ../src/biogeophys, I could not find the specific formula,do you have any suggestions?
All I found is :
biogeophys/SoilFluxesMod.F90:
! Calculate the skin temperature as a weighted sum of all the surface contributions (surface water table, snow, etc...)
! Note: This is the bare ground calculation of skin temperature
! The Urban and Vegetation are done in other place. Urban=Later in this function Veg=CanopyFluxMod
! t_skin_patch(p) = ((1._r8 - emv(p))*(1-frac_veg_nosno(p)) * sqrt(sqrt(lw_grnd))) + emv(p)*t_veg(p)
! if( frac_veg_nosno(p).eq.0 ) then
! t_skin_patch(p) = ((1._r8 - emv(p))*(1-frac_veg_nosno(p)) * sqrt(sqrt(lw_grnd))) + &
! emv(p) * frac_veg_nosno(p) * t_veg(p)
! end if
if(frac_veg_nosno(p).eq.0) t_skin_patch(p) = sqrt(sqrt(lw_grnd))
I want to know the calculation of vegetation emissivity, ground emissivity, ground surface temperature and vegetation temperature in the code. However, when I use "grep -r emv_patch t_veg_patch t_grnd_col " in the ../src/biogeophys, I could not find the specific formula,do you have any suggestions?
All I found is :
biogeophys/SoilFluxesMod.F90:
! Calculate the skin temperature as a weighted sum of all the surface contributions (surface water table, snow, etc...)
! Note: This is the bare ground calculation of skin temperature
! The Urban and Vegetation are done in other place. Urban=Later in this function Veg=CanopyFluxMod
! t_skin_patch(p) = ((1._r8 - emv(p))*(1-frac_veg_nosno(p)) * sqrt(sqrt(lw_grnd))) + emv(p)*t_veg(p)
! if( frac_veg_nosno(p).eq.0 ) then
! t_skin_patch(p) = ((1._r8 - emv(p))*(1-frac_veg_nosno(p)) * sqrt(sqrt(lw_grnd))) + &
! emv(p) * frac_veg_nosno(p) * t_veg(p)
! end if
if(frac_veg_nosno(p).eq.0) t_skin_patch(p) = sqrt(sqrt(lw_grnd))