R/PTpet.R

PTpet <-
function(Rn, T_C, PTconstant = 1.26){
    LatentHtEvap <- 2500	# kJ/kg
    DensityWater <- 1000	# kg/m3
    PsychConstant <- 0.066	# kPa/K
    potentialET <- PTconstant * SatVapPresSlope(T_C) * Rn/((SatVapPresSlope(T_C) + PsychConstant) * (LatentHtEvap * DensityWater))
    potentialET[which(potentialET < 0)] <- 0
    return(signif(potentialET,2))

}

Try the EcoHydRology package in your browser

Any scripts or data that you put into this service are public.

EcoHydRology documentation built on May 2, 2019, 4:47 p.m.