fCalcPotRadiation: fCalcPotRadiation

In REddyProc: Data processing and plotting utilities of (half-)hourly eddy-covariance measurements

Description

Calculate the potential radiation

Usage

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

Arguments

DoY.V.n	Data vector with day of year (DoY), same length as Hour or length 1
Hour.V.n	Data vector with time as decimal hour
Lat_deg.n	Latitude in (decimal) degrees
Long_deg.n	Longitude in (decimal) degrees
TimeZone_h.n	Time zone (in hours)
useSolartime.b	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

Value

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

Author(s)

AMM (Department for Biogeochemical Integration at MPI-BGC, Jena, Germany)

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.b = 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)

Example output

Warning message:
In fCalcPotRadiation(160, hour, 39.94, -5.77, TimeZone = +1, useSolartime.b = FALSE) :
  Argument names useSolartime.b have been deprecated. Please, use instead useSolartime

REddyProc documentation built on May 2, 2019, 5:19 p.m.