R/SatWatVap.R

## These functions justify the use of Arden Buck instead of
## Goff-Gratch
##
## References
##
## http://cires.colorado.edu/~voemel/vp.html


fnpsvp <- function(Tkelvin){
#   water boiling point = 373.16 oK

  u = Tkelvin / 373.16;
  v = 373.16 / Tkelvin;

  tmp <- -7.90298 * (v - 1) + 5.02808 * log10(v) 
          - 1.3816 * 10^-7 * (10^(11.344 * (1 - u)) - 1)   
          + 8.1328 * 10^-3 * (10^(-3.49149 * (v - 1)) - 1)
       + log10(1013.246);

  esat = 10^tmp; ## This is in hecto Pascals
  esat
}

## Arden Buck Equation

ardBuck <- function(Tcelsius){

  a <- (18.678 - Tcelsius/234.5)*Tcelsius
  b <- 257.14 + Tcelsius
  res <- (6.1121  * exp(a/b))
  res ## This is hecto Pascals
}

tetens <- function(Tcelsius){

    a = 0.611
    b = 17.502
    c = 240.97

    swvp <- a * exp((b * Tcelsius)/(Tcelsius + c))
    swvp
}

Try the BioCro package in your browser

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

BioCro documentation built on May 2, 2019, 6:15 p.m.