fCalcPotRadiation: fCalcPotRadiation

In bgctw/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	by default corrects hour (given in local winter time) for latitude to solar time (where noon is exactly at 12:00). Set this to FALSE to directly use local winter time
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)

bgctw/REddyProc documentation built on March 26, 2024, 11:35 p.m.