R/UtilityFunctions.R
In trenchproject/TrenchR: Tools for Microclimate and Biophysical Ecology
Defines functions
daylength
Documented in
daylength
#' @title Calendar Day from Date
#'
#' @description The function converts a date (day, month, year) to Calendar Day (day of year).
#'
#' @param day \code{character} numerical date in standard format (e.g. \code{"2017-01-02"}, \code{"01-02"}, \code{"01/02/2017"} etc).
#'
#' @param format \code{character} date format following \code{\link[base]{as.POSIXlt}} conventions. Default value = \code{"\%Y-\%m-\%d"}.
#'
#' @return \code{numeric} Julian day number, 1-366 (e.g. 1 for January 1st).
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' day_of_year(day = "2017-04-22",
#' format = "%Y-%m-%d")
#' day_of_year(day = "2017-04-22")
#' day_of_year(day = "04/22/2017",
#' format = "%m/%d/%Y")
#'
day_of_year <- function (day,
format = "%Y-%m-%d") {
day <- as.POSIXlt(day, format = format)
as.numeric(strftime(day, format = "%j"))
}
#' @title Solar Declination in Radians
#'
#' @description The function calculates solar declination, which is the angular distance of the sun north or south of the earth’s equator, based on the day of year \insertCite{Campbell1998}{TrenchR}.
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @return \code{numeric} declination angle (radians).
#'
#' @family utility functions
#'
#' @references
#' \insertAllCited{}
#'
#' @export
#'
#' @examples
#' dec_angle(doy = 112)
#' dec_angle(doy = 360)
#'
dec_angle <- function (doy) {
stopifnot(doy > 0,
doy < 367)
RevAng <- 0.21631 + 2 * atan(0.967 * tan (0.0086 * (-186 + doy)))
asin(0.39795 * cos (RevAng))
}
#' @title Day Length
#'
#' @description The function calculates daylength in hours as a function of latitude and day of year. Uses the CMB model \insertCite{Campbell1998}{TrenchR}.
#'
#' @param lat \code{numeric} latitude (decimal degrees).
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @return \code{numeric} day length (hours).
#'
#' @references
#' \insertAllCited{}
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' daylength(lat = 47.61,
#' doy = 112)
#'
daylength <- function(lat,
doy){
stopifnot(lat >= -90,
lat <= 90,
doy > 0,
doy < 367)
lat_rad <- degrees_to_radians(lat)
DecAng <- dec_angle(doy)
subset <- (sin (6 * pi / 180) + sin (lat_rad) * sin (DecAng)) / (cos (lat_rad) * cos (DecAng))
subset[which(subset>1)] <- 1
subset[which(subset< -1)] <- -1
24 - (24 / pi) * acos(subset)
}
#' @title Time of Solar Noon
#'
#' @description The function calculates the time of solar noon in hours as a function of the day of year and longitude \insertCite{Campbell1998}{TrenchR}.
#'
#' @param lon \code{numeric} longitude (decimal degrees).
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from a standard date via \code{\link{day_of_year}}.
#'
#' @param offset \code{numeric} number of hours to add to UTC (Coordinated Universal Time) to get local time (improves accuracy but not always necessary). Defaults to NA.
#'
#' @return \code{numeric} time of solar noon (hours).
#'
#' @references
#' \insertAllCited{}
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' solar_noon(lon = -122.335,
#' doy = 112)
#'
solar_noon <- function (lon,
doy,
offset = NA){
stopifnot(lon >= -180,
lon <= 180,
doy > 0,
doy < 367)
# Calculates the time of solar noon for each day using longitude correction (LC), equation of time (ET), and a conversion (f)
f <- (279.575 + 0.9856 * doy) # f in degrees as a function of day of year, p.169 Campbell & Norman 1998
f[f > 360] = f[f > 360] - 360 # ensure 0 to 360 degrees
f <- f * pi / 180 # convert f in degrees to radians
# (11.4) Equation of time: ET is a 15-20 minute correction which depends on calendar day
ET <- (-104.7 * sin(f) + 596.2 * sin(2*f) + 4.3 * sin(3*f) - 12.7 * sin(4*f) - 429.3 * cos(f) - 2.0 * cos(2*f) + 19.3 * cos(3*f)) / 3600
lon[lon < 0] <- 360 + lon[lon < 0] # convert to 0 to 360
LC <- 1 / 15 * (lon %% 15) # longitude correction, 1/15h for each degree of standard meridian
LC[LC > 0.5] <- LC[LC > 0.5] - 1
t_0 <- 12 - LC - ET # solar noon, p168, 11.3, Campbell &Norman 1998
# Check if offset is as expected. (Is the timezone of the location the same as that of the meridian
# that's within 7.5 degrees from that location?)
lon[lon > 180] <- lon[lon > 180] - 360
if (!is.na(offset)) {
offset_theory <- as.integer(lon / 15) + lon / abs(lon) * as.integer(abs(lon) %% 15 / 7.5)
t_0 = t_0 - offset_theory + offset
}
t_0
}
#' @title Zenith Angle
#'
#' @description The function calculates the zenith angle, the location of the sun as an angle (in degrees) measured from vertical \insertCite{Campbell1998}{TrenchR}.
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @param lat \code{numeric} latitude (decimal degrees).
#'
#' @param lon \code{numeric} longitude (decimal degrees).
#'
#' @param hour \code{numeric} hour of the day (0-24).
#'
#' @param offset \code{numeric} the number of hours to add to UTC (Coordinated Universal Time) to get local time (improves accuracy but not always necessary). Optional. Defaults to NA.
#'
#' @return \code{numeric} zenith angle (degrees)
#'
#' @family utility functions
#'
#' @export
#'
#' @references
#' \insertAllCited{}
#'
#' @examples
#' zenith_angle(doy = 112,
#' lat = 47.61,
#' lon = -122.33,
#' hour = 12)
#'
zenith_angle <- function (doy,
lat,
lon,
hour,
offset = NA) {
stopifnot(doy > 0,
doy < 367,
lat >= -90,
lat <= 90,
lon >= -180,
lon <= 180,
hour >= 0,
hour <= 24)
lat <- lat * pi / 180 # to radians
RevAng <- 0.21631 + 2 * atan(0.967 * tan(0.0086 * (-186 + doy))) # Revolution angle in radians
DecAng <- asin(0.39795 * cos(RevAng)) # Declination angle in radians
f <- 279.575 + 0.9856 * doy # f in degrees as a function of day of year, p.169 Campbell & Norman 1998
f <- f * pi / 180 # convert f in degrees to radians
# (11.4) Equation of time: ET is a 15-20 minute correction which depends on calendar day
ET <- (-104.7 * sin(f) + 596.2 * sin(2 * f) + 4.3 * sin(3 * f) - 12.7 * sin(4 * f) - 429.3 * cos(f) - 2.0 * cos(2 * f) + 19.3 * cos(3 * f)) / 3600
lon[lon < 0] <- 360 + lon[lon < 0] #convert to 0 to 360
LC <- 1 / 15 * (lon %% 15) # longitude correction, 1/15h for each degree of standard meridian
LC[LC > 0.5] <- LC[LC > 0.5] - 1
t_0 <- 12 - LC - ET # solar noon
# Check if offset is as expected. (Is the timezone of the location the same as that of the meridian
# that's within 7.5 degrees from that location?)
lon[lon > 180] <- lon[lon > 180] - 360
if (!is.na(offset)) {
offset_theory <- as.integer(lon / 15) + lon / abs(lon) * as.integer(abs(lon) %% 15 / 7.5)
t_0 <- t_0 - offset_theory + offset
}
cos.zenith <- sin(DecAng) * sin(lat) + cos(DecAng) * cos(lat) * cos(pi / 12 * (hour - t_0)) # cos of zenith angle in radians
zenith <- acos(cos.zenith) * 180 / pi # zenith angle in degrees
zenith[zenith > 90] <- 90 # if measured from the vertical psi can't be greater than pi/2 (90 degrees)
zenith
}
#' @title Azimuth Angle
#'
#' @description The function calculates the azimuth angle, the angle (degrees) from which the sunlight is coming measured from true north or south measured in the horizontal plane. The azimuth angle is measured with respect to due south, increasing in the counter clockwise direction so 90 degrees is east \insertCite{Campbell1998}{TrenchR}.
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @param lat \code{numeric} latitude (decimal degrees).
#'
#' @param lon \code{numeric} longitude (decimal degrees).
#'
#' @param hour \code{numeric} hour of the day.
#'
#' @param offset \code{numeric} number of hours to add to UTC (Coordinated Universal Time) to get local time (to improve accuracy but not always necessary). Optional. Defaults to NA.
#'
#' @return \code{numeric} azimuth angle (degrees).
#'
#' @references
#' \insertAllCited{}
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' azimuth_angle(doy = 112,
#' lat = 47.61,
#' lon = -122.33,
#' hour = 12,
#' offset = -8)
#'
azimuth_angle <- function (doy,
lat,
lon,
hour,
offset = NA) {
stopifnot(doy > 0,
doy < 367,
lat >= -90,
lat <= 90,
lon >= -180,
lon <= 180,
hour >= 0,
hour <= 24)
lat <- lat * pi / 180 #to radians
RevAng <- 0.21631 + 2 * atan(0.967 * tan(0.0086 * (-186 + doy))) # Revolution angle in radians
DecAng <- asin(0.39795 * cos(RevAng)) # Declination angle in radians
f <- 279.575 + 0.9856 * doy # f in degrees as a function of day of year, p.169 Campbell & Norman 2000
f <- f * pi / 180 # convert f in degrees to radians
# (11.4) Equation of time: ET is a 15-20 minute correction which depends on calendar day
ET <- (-104.7 * sin(f) + 596.2 * sin(2 * f) + 4.3 * sin(3 * f) - 12.7 * sin(4 * f) - 429.3 * cos(f) - 2.0 * cos(2 * f) + 19.3 * cos(3 * f)) / 3600
lon[lon < 0] <- 360 + lon[lon < 0] # convert to 0 to 360
# Set up two Booleans indicating whether we need to apply a correction on azimuth angle at the end
azi_corr1 <- FALSE
azi_corr2 <- TRUE
LC <- 1 / 15 * (lon %% 15) # longitude correction, 1/15h for each degree of standard meridian
if (LC > 0.5) {
LC = LC - 1
} else {
azi_corr1 = TRUE
}
t_0 <- 12 - LC - ET # solar noon
# Check if offset is as expected. (Is the timezone of the location the same as that of the meridian
# that's within 7.5 degrees from that location?)
lon[lon > 180] <- lon[lon > 180] - 360
if (!is.na(offset)) {
offset_theory <- as.integer(lon / 15) + lon / abs(lon) * as.integer(abs(lon) %% 15 / 7.5)
if (offset_theory != offset) {
t_0 <- t_0 - offset_theory + offset
azi_corr2 <- FALSE
}
}
cos.zenith <- sin(DecAng) * sin(lat) + cos(DecAng) * cos(lat) * cos(pi / 12 * (hour - t_0)) #cos of zenith angle in radians
zenith <- acos(cos.zenith) # zenith angle in radians
if (zenith > pi / 2){
zenith <- pi / 2 # if measured from the vertical psi can't be greater than pi/2 (90 degrees)
}
cos.azimuth <- -(sin(DecAng) - cos(zenith) * sin(lat)) / (cos(lat) * sin(zenith))
azimuth <- acos(cos.azimuth) * 180 / pi # azimuth angle in degrees
if (azi_corr1 && azi_corr2) {
azimuth = 360 - azimuth
}
azimuth
}
#' @title Air Pressure from Elevation
#'
#' @description The function estimates air pressure (kPa) as a function of elevation \insertCite{engingeeringtoolbox}{TrenchR}.
#'
#' @param elev \code{numeric} elevation (meters).
#'
#' @family utility functions
#'
#' @return \code{numeric} air pressure (kPa).
#'
#' @export
#'
#' @references
#' \insertAllCited{}
#'
#' @examples
#' airpressure_from_elev(elev = 1500)
#'
airpressure_from_elev <- function (elev) {
stopifnot(elev >= 0)
101.3 * exp (-elev / 8200)
}
#' @title Convert Among Temperature Scales
#'
#' @description The function converts temperatures among Celsius, Fahrenheit, and Kelvin \insertCite{Blischak}{TrenchR}.
#'
#' @param temperature \code{numeric} temperature (Celsius, Fahrenheit, or Kelvin).
#'
#' @family utility functions
#'
#' @return \code{numeric} temperature (Celsius, Fahrenheit, or Kelvin).
#'
#' @export
#'
#' @references
#' \insertAllCited{}
#'
#' @name temperature conversions
#'
#' @examples
#' kelvin_to_celsius(temperature = 270)
#' fahrenheit_to_kelvin(temperature = 85)
#' fahrenheit_to_celsius(temperature = 85)
#' celsius_to_kelvin(temperature = -10)
#'
fahrenheit_to_kelvin <- function (temperature) {
((temperature - 32) * (5 / 9)) + 273.15
}
#' @name temperature conversions
#'
#' @export
#'
kelvin_to_celsius <- function (temperature) {
temperature - 273.15
}
#' @name temperature conversions
#'
#' @export
#'
celsius_to_kelvin <- function (temperature) {
temperature + 273.15
}
#' @name temperature conversions
#'
#' @export
#'
fahrenheit_to_celsius <- function (temperature) {
temp_k <- fahrenheit_to_kelvin(temperature)
kelvin_to_celsius(temp_k)
}
#' @title Convert Angles Between Radians and Degrees
#'
#' @description The function converts angles in radians to degrees or degrees to radians.
#'
#' @param rad \code{numeric} angle (radians).
#'
#' @param deg \code{numeric} angle (degrees).
#'
#' @return \code{numeric} angle (degrees or radians).
#'
#' @name angle conversions
#'
#' @export
#'
#' @examples
#' radians_to_degrees(0.831)
#' degrees_to_radians(47.608)
#'
radians_to_degrees <- function (rad) {
rad * 180 / pi
}
#' @name angle conversions
#'
#' @export
#'
degrees_to_radians <- function (deg) {
deg * pi / 180
}
trenchproject/TrenchR documentation built on Oct. 10, 2023, 10:12 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions daylength
Documented in daylength
#' @title Calendar Day from Date
#'
#' @description The function converts a date (day, month, year) to Calendar Day (day of year).
#'
#' @param day \code{character} numerical date in standard format (e.g. \code{"2017-01-02"}, \code{"01-02"}, \code{"01/02/2017"} etc).
#'
#' @param format \code{character} date format following \code{\link[base]{as.POSIXlt}} conventions. Default value = \code{"\%Y-\%m-\%d"}.
#'
#' @return \code{numeric} Julian day number, 1-366 (e.g. 1 for January 1st).
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' day_of_year(day = "2017-04-22",
#' format = "%Y-%m-%d")
#' day_of_year(day = "2017-04-22")
#' day_of_year(day = "04/22/2017",
#' format = "%m/%d/%Y")
#'
day_of_year <- function (day,
format = "%Y-%m-%d") {
day <- as.POSIXlt(day, format = format)
as.numeric(strftime(day, format = "%j"))
}
#' @title Solar Declination in Radians
#'
#' @description The function calculates solar declination, which is the angular distance of the sun north or south of the earth’s equator, based on the day of year \insertCite{Campbell1998}{TrenchR}.
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @return \code{numeric} declination angle (radians).
#'
#' @family utility functions
#'
#' @references
#' \insertAllCited{}
#'
#' @export
#'
#' @examples
#' dec_angle(doy = 112)
#' dec_angle(doy = 360)
#'
dec_angle <- function (doy) {
stopifnot(doy > 0,
doy < 367)
RevAng <- 0.21631 + 2 * atan(0.967 * tan (0.0086 * (-186 + doy)))
asin(0.39795 * cos (RevAng))
}
#' @title Day Length
#'
#' @description The function calculates daylength in hours as a function of latitude and day of year. Uses the CMB model \insertCite{Campbell1998}{TrenchR}.
#'
#' @param lat \code{numeric} latitude (decimal degrees).
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @return \code{numeric} day length (hours).
#'
#' @references
#' \insertAllCited{}
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' daylength(lat = 47.61,
#' doy = 112)
#'
daylength <- function(lat,
doy){
stopifnot(lat >= -90,
lat <= 90,
doy > 0,
doy < 367)
lat_rad <- degrees_to_radians(lat)
DecAng <- dec_angle(doy)
subset <- (sin (6 * pi / 180) + sin (lat_rad) * sin (DecAng)) / (cos (lat_rad) * cos (DecAng))
subset[which(subset>1)] <- 1
subset[which(subset< -1)] <- -1
24 - (24 / pi) * acos(subset)
}
#' @title Time of Solar Noon
#'
#' @description The function calculates the time of solar noon in hours as a function of the day of year and longitude \insertCite{Campbell1998}{TrenchR}.
#'
#' @param lon \code{numeric} longitude (decimal degrees).
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from a standard date via \code{\link{day_of_year}}.
#'
#' @param offset \code{numeric} number of hours to add to UTC (Coordinated Universal Time) to get local time (improves accuracy but not always necessary). Defaults to NA.
#'
#' @return \code{numeric} time of solar noon (hours).
#'
#' @references
#' \insertAllCited{}
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' solar_noon(lon = -122.335,
#' doy = 112)
#'
solar_noon <- function (lon,
doy,
offset = NA){
stopifnot(lon >= -180,
lon <= 180,
doy > 0,
doy < 367)
# Calculates the time of solar noon for each day using longitude correction (LC), equation of time (ET), and a conversion (f)
f <- (279.575 + 0.9856 * doy) # f in degrees as a function of day of year, p.169 Campbell & Norman 1998
f[f > 360] = f[f > 360] - 360 # ensure 0 to 360 degrees
f <- f * pi / 180 # convert f in degrees to radians
# (11.4) Equation of time: ET is a 15-20 minute correction which depends on calendar day
ET <- (-104.7 * sin(f) + 596.2 * sin(2*f) + 4.3 * sin(3*f) - 12.7 * sin(4*f) - 429.3 * cos(f) - 2.0 * cos(2*f) + 19.3 * cos(3*f)) / 3600
lon[lon < 0] <- 360 + lon[lon < 0] # convert to 0 to 360
LC <- 1 / 15 * (lon %% 15) # longitude correction, 1/15h for each degree of standard meridian
LC[LC > 0.5] <- LC[LC > 0.5] - 1
t_0 <- 12 - LC - ET # solar noon, p168, 11.3, Campbell &Norman 1998
# Check if offset is as expected. (Is the timezone of the location the same as that of the meridian
# that's within 7.5 degrees from that location?)
lon[lon > 180] <- lon[lon > 180] - 360
if (!is.na(offset)) {
offset_theory <- as.integer(lon / 15) + lon / abs(lon) * as.integer(abs(lon) %% 15 / 7.5)
t_0 = t_0 - offset_theory + offset
}
t_0
}
#' @title Zenith Angle
#'
#' @description The function calculates the zenith angle, the location of the sun as an angle (in degrees) measured from vertical \insertCite{Campbell1998}{TrenchR}.
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @param lat \code{numeric} latitude (decimal degrees).
#'
#' @param lon \code{numeric} longitude (decimal degrees).
#'
#' @param hour \code{numeric} hour of the day (0-24).
#'
#' @param offset \code{numeric} the number of hours to add to UTC (Coordinated Universal Time) to get local time (improves accuracy but not always necessary). Optional. Defaults to NA.
#'
#' @return \code{numeric} zenith angle (degrees)
#'
#' @family utility functions
#'
#' @export
#'
#' @references
#' \insertAllCited{}
#'
#' @examples
#' zenith_angle(doy = 112,
#' lat = 47.61,
#' lon = -122.33,
#' hour = 12)
#'
zenith_angle <- function (doy,
lat,
lon,
hour,
offset = NA) {
stopifnot(doy > 0,
doy < 367,
lat >= -90,
lat <= 90,
lon >= -180,
lon <= 180,
hour >= 0,
hour <= 24)
lat <- lat * pi / 180 # to radians
RevAng <- 0.21631 + 2 * atan(0.967 * tan(0.0086 * (-186 + doy))) # Revolution angle in radians
DecAng <- asin(0.39795 * cos(RevAng)) # Declination angle in radians
f <- 279.575 + 0.9856 * doy # f in degrees as a function of day of year, p.169 Campbell & Norman 1998
f <- f * pi / 180 # convert f in degrees to radians
# (11.4) Equation of time: ET is a 15-20 minute correction which depends on calendar day
ET <- (-104.7 * sin(f) + 596.2 * sin(2 * f) + 4.3 * sin(3 * f) - 12.7 * sin(4 * f) - 429.3 * cos(f) - 2.0 * cos(2 * f) + 19.3 * cos(3 * f)) / 3600
lon[lon < 0] <- 360 + lon[lon < 0] #convert to 0 to 360
LC <- 1 / 15 * (lon %% 15) # longitude correction, 1/15h for each degree of standard meridian
LC[LC > 0.5] <- LC[LC > 0.5] - 1
t_0 <- 12 - LC - ET # solar noon
# Check if offset is as expected. (Is the timezone of the location the same as that of the meridian
# that's within 7.5 degrees from that location?)
lon[lon > 180] <- lon[lon > 180] - 360
if (!is.na(offset)) {
offset_theory <- as.integer(lon / 15) + lon / abs(lon) * as.integer(abs(lon) %% 15 / 7.5)
t_0 <- t_0 - offset_theory + offset
}
cos.zenith <- sin(DecAng) * sin(lat) + cos(DecAng) * cos(lat) * cos(pi / 12 * (hour - t_0)) # cos of zenith angle in radians
zenith <- acos(cos.zenith) * 180 / pi # zenith angle in degrees
zenith[zenith > 90] <- 90 # if measured from the vertical psi can't be greater than pi/2 (90 degrees)
zenith
}
#' @title Azimuth Angle
#'
#' @description The function calculates the azimuth angle, the angle (degrees) from which the sunlight is coming measured from true north or south measured in the horizontal plane. The azimuth angle is measured with respect to due south, increasing in the counter clockwise direction so 90 degrees is east \insertCite{Campbell1998}{TrenchR}.
#'
#' @param doy \code{numeric} day of year (1-366). This can be obtained from standard date via \code{\link{day_of_year}}.
#'
#' @param lat \code{numeric} latitude (decimal degrees).
#'
#' @param lon \code{numeric} longitude (decimal degrees).
#'
#' @param hour \code{numeric} hour of the day.
#'
#' @param offset \code{numeric} number of hours to add to UTC (Coordinated Universal Time) to get local time (to improve accuracy but not always necessary). Optional. Defaults to NA.
#'
#' @return \code{numeric} azimuth angle (degrees).
#'
#' @references
#' \insertAllCited{}
#'
#' @family utility functions
#'
#' @export
#'
#' @examples
#' azimuth_angle(doy = 112,
#' lat = 47.61,
#' lon = -122.33,
#' hour = 12,
#' offset = -8)
#'
azimuth_angle <- function (doy,
lat,
lon,
hour,
offset = NA) {
stopifnot(doy > 0,
doy < 367,
lat >= -90,
lat <= 90,
lon >= -180,
lon <= 180,
hour >= 0,
hour <= 24)
lat <- lat * pi / 180 #to radians
RevAng <- 0.21631 + 2 * atan(0.967 * tan(0.0086 * (-186 + doy))) # Revolution angle in radians
DecAng <- asin(0.39795 * cos(RevAng)) # Declination angle in radians
f <- 279.575 + 0.9856 * doy # f in degrees as a function of day of year, p.169 Campbell & Norman 2000
f <- f * pi / 180 # convert f in degrees to radians
# (11.4) Equation of time: ET is a 15-20 minute correction which depends on calendar day
ET <- (-104.7 * sin(f) + 596.2 * sin(2 * f) + 4.3 * sin(3 * f) - 12.7 * sin(4 * f) - 429.3 * cos(f) - 2.0 * cos(2 * f) + 19.3 * cos(3 * f)) / 3600
lon[lon < 0] <- 360 + lon[lon < 0] # convert to 0 to 360
# Set up two Booleans indicating whether we need to apply a correction on azimuth angle at the end
azi_corr1 <- FALSE
azi_corr2 <- TRUE
LC <- 1 / 15 * (lon %% 15) # longitude correction, 1/15h for each degree of standard meridian
if (LC > 0.5) {
LC = LC - 1
} else {
azi_corr1 = TRUE
}
t_0 <- 12 - LC - ET # solar noon
# Check if offset is as expected. (Is the timezone of the location the same as that of the meridian
# that's within 7.5 degrees from that location?)
lon[lon > 180] <- lon[lon > 180] - 360
if (!is.na(offset)) {
offset_theory <- as.integer(lon / 15) + lon / abs(lon) * as.integer(abs(lon) %% 15 / 7.5)
if (offset_theory != offset) {
t_0 <- t_0 - offset_theory + offset
azi_corr2 <- FALSE
}
}
cos.zenith <- sin(DecAng) * sin(lat) + cos(DecAng) * cos(lat) * cos(pi / 12 * (hour - t_0)) #cos of zenith angle in radians
zenith <- acos(cos.zenith) # zenith angle in radians
if (zenith > pi / 2){
zenith <- pi / 2 # if measured from the vertical psi can't be greater than pi/2 (90 degrees)
}
cos.azimuth <- -(sin(DecAng) - cos(zenith) * sin(lat)) / (cos(lat) * sin(zenith))
azimuth <- acos(cos.azimuth) * 180 / pi # azimuth angle in degrees
if (azi_corr1 && azi_corr2) {
azimuth = 360 - azimuth
}
azimuth
}
#' @title Air Pressure from Elevation
#'
#' @description The function estimates air pressure (kPa) as a function of elevation \insertCite{engingeeringtoolbox}{TrenchR}.
#'
#' @param elev \code{numeric} elevation (meters).
#'
#' @family utility functions
#'
#' @return \code{numeric} air pressure (kPa).
#'
#' @export
#'
#' @references
#' \insertAllCited{}
#'
#' @examples
#' airpressure_from_elev(elev = 1500)
#'
airpressure_from_elev <- function (elev) {
stopifnot(elev >= 0)
101.3 * exp (-elev / 8200)
}
#' @title Convert Among Temperature Scales
#'
#' @description The function converts temperatures among Celsius, Fahrenheit, and Kelvin \insertCite{Blischak}{TrenchR}.
#'
#' @param temperature \code{numeric} temperature (Celsius, Fahrenheit, or Kelvin).
#'
#' @family utility functions
#'
#' @return \code{numeric} temperature (Celsius, Fahrenheit, or Kelvin).
#'
#' @export
#'
#' @references
#' \insertAllCited{}
#'
#' @name temperature conversions
#'
#' @examples
#' kelvin_to_celsius(temperature = 270)
#' fahrenheit_to_kelvin(temperature = 85)
#' fahrenheit_to_celsius(temperature = 85)
#' celsius_to_kelvin(temperature = -10)
#'
fahrenheit_to_kelvin <- function (temperature) {
((temperature - 32) * (5 / 9)) + 273.15
}
#' @name temperature conversions
#'
#' @export
#'
kelvin_to_celsius <- function (temperature) {
temperature - 273.15
}
#' @name temperature conversions
#'
#' @export
#'
celsius_to_kelvin <- function (temperature) {
temperature + 273.15
}
#' @name temperature conversions
#'
#' @export
#'
fahrenheit_to_celsius <- function (temperature) {
temp_k <- fahrenheit_to_kelvin(temperature)
kelvin_to_celsius(temp_k)
}
#' @title Convert Angles Between Radians and Degrees
#'
#' @description The function converts angles in radians to degrees or degrees to radians.
#'
#' @param rad \code{numeric} angle (radians).
#'
#' @param deg \code{numeric} angle (degrees).
#'
#' @return \code{numeric} angle (degrees or radians).
#'
#' @name angle conversions
#'
#' @export
#'
#' @examples
#' radians_to_degrees(0.831)
#' degrees_to_radians(47.608)
#'
radians_to_degrees <- function (rad) {
rad * 180 / pi
}
#' @name angle conversions
#'
#' @export
#'
degrees_to_radians <- function (deg) {
deg * pi / 180
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.