fCalcPotRadiation: fCalcPotRadiation

In REddyProc: Post Processing of (Half-)Hourly Eddy-Covariance Measurements

fCalcPotRadiation

Description

Calculate the potential radiation

Usage

fCalcPotRadiation(DoY = DoY.V.n, Hour = Hour.V.n, 
    LatDeg = Lat_deg.n, LongDeg = Long_deg.n, 
    TimeZone = TimeZone_h.n, useSolartime = TRUE, 
    DoY.V.n, Hour.V.n, Lat_deg.n, Long_deg.n, 
    TimeZone_h.n, useSolartime.b = TRUE)

Arguments

DoY	Data vector with day of year (DoY), same length as Hour or length 1
Hour	Data vector with time as decimal hour of local time zone
LatDeg	Latitude in (decimal) degrees
LongDeg	Longitude in (decimal) degrees
TimeZone	Time zone (in hours)
useSolartime
DoY.V.n	deprecated
Hour.V.n	deprecated
Lat_deg.n	deprecated
Long_deg.n	deprecated
TimeZone_h.n	deprecated
useSolartime.b	deprecated

Value

Data vector of potential radiation (PotRad, W_m-2)

Author(s)

AMM Department for Biogeochemical Integration at MPI-BGC, Jena, Germany <[e-mail]> [cph], Thomas Wutzler <[e-mail]> [aut, cre], Markus Reichstein <[e-mail]> [aut], Antje Maria Moffat <[e-mail]> [aut, trl], Olaf Menzer <[e-mail]> [ctb], Mirco Migliavacca <[e-mail]> [aut], Kerstin Sickel <[e-mail]> [ctb, trl], Ladislav <U+0160>igut <[e-mail]> [ctb]

Examples

hour <- seq(8, 16, by = 0.1)
potRadSolar <- fCalcPotRadiation(160, hour, 39.94, -5.77, TimeZone = +1)
potRadLocal <- fCalcPotRadiation(160, hour, 39.94, -5.77, TimeZone = +1
	  , useSolartime = FALSE)
plot(potRadSolar ~ hour, type = 'l')
abline(v = 13, lty = "dotted")
lines(potRadLocal ~  hour, col = "blue")
abline(v = 12, col = "blue", lty = "dotted")
legend("bottomright", legend = c("solar time", "local winter time")
	, col = c("black", "blue"), inset = 0.05, lty = 1)

REddyProc documentation built on March 18, 2022, 5:41 p.m.