fix cal_rH2

CUG-hydro · Aug 21, 2024 · 728bc4a · 728bc4a
1 parent a1e94be
commit 728bc4a
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/R/ET0_helper.R b/R/ET0_helper.R
@@ -153,9 +153,10 @@ cal_rH <- function(U2, h = 0.12) {
 cal_rH2 <- function(U2, Tair, Pa = atm) {
   # `f(U2) = 2.6 * (1 + 0.54U2)` is equivalent to Shuttleworth1993
   # rho_a * Cp / rH = f(U2)
-  f_U2 <- 2.6 * (1 + 0.54 * U2)
+  f_U2 <- 6.43 * (1 + 0.536 * U2)
   rho_a <- 3.486 * Pa / cal_TvK(Tair) # FAO56, Eq. 3-5, kg m-3
-  rH <- rho_a * Cp / f_U2 * 86400
+  gamma <- cal_gamma(Tair, Pa)
+  rH <- rho_a * Cp / f_U2 / gamma * 86400
   rH
 }