calc_PotRadiation_CosineResponsePower: Model potential solar radiation at the surface with power law...

In laubblatt/cleaRskyQuantileRegression: Calculate monthly fractional clear-sky shortwave transmission from global radiation records

Description

Attenuation of the solar beam in a curved, refractive atmosphere modelled with a simple power law with an exponent to the cosine of the solar zenit angle The function calles other functions to calculate the solar zenith angle, the extra terrestrial radiation copied from the R packages REddyProc::fCalcPotRadiation and solartime. Since their arguments where changing I copied part of these packages.

Usage

calc_PotRadiation_CosineResponsePower(doy, hour, latDeg, longDeg, timeZone,
  isCorrectSolartime = TRUE, cosineResponsePower = 1.2)

Arguments

doy	value or vector of day of year
hour	value or vector of hour of day also with decimal fraction
latDeg	latitude in degrees
timeZone	numeric, if local time is used then 0 otherwise the time must be adapted
cosineResponsePower	numeric defaults to 1.2, may be lower at high latitude sites, see Long and Ackerman 2000
latDeg	longitude in degrees

Value

vector of potential solar radiation at the surface

Author(s)

Maik Renner, mrenner [at] bgc-jena.mpg.de

References

Long and Ackerman 2000 JGR Renner et al., 2019 ESS

Examples

data(LIN2003)
# calculate potential surface solar radiation
LIN2003[ , IncomingShortwavePotential := calc_PotRadiation_CosineResponsePower(doy = yday(Date),hour = Time/3600 + 0.25, latDeg = 52.21, longDeg = 14.122, timeZone = 0) ]
LIN2003
plot(IncomingShortwave ~ IncomingShortwavePotential, data = LIN2003[month(Date) == 6, ], type = "p")

laubblatt/cleaRskyQuantileRegression documentation built on Nov. 27, 2019, 12:26 p.m.