Description Usage Arguments Value Author(s) References See Also Examples
Compute the angles which describe the intradaily apparent movement of the Sun from the Earth.
1 
solD 
A 
sample 
Increment of the intradaily sequence. It is a character
string, containing one of ‘"sec"’, ‘"min"’, ‘"hour"’. This can
optionally be preceded by a (positive or negative) integer and a
space, or followed by ‘"s"’. It is used by It is not considered when 
BTi 
Intradaily time base, a 
EoT 
logical, if 
keep.night 
logical, if 
method 

A zoo
object is returned with these components:
w 
numeric, solar hour angle (radians) 
aman 
logical, 
cosThzS 
numeric, cosine of the solar zenith angle 
AzS 
numeric, solar acimuth angle (radians) 
AlS 
numeric, solar elevation angle (radians) 
Bo0 
numeric, extraatmospheric irradiance (W/m2) 
rd, rg 
numeric, relation between irradiance and irradiation of diffuse and global values, respectively, following the correlations proposed by CollaresPereira and Rabl 
The latitude is stored as the attribute lat
of this object.
Oscar Perpiñán Lamigueiro.
Cooper, P.I., Solar Energy, 12, 3 (1969). "The Absorption of Solar Radiation in Solar Stills"
Spencer, Search 2 (5), 172, https://www.mailarchive.com/sundial@unikoeln.de/msg01050.html
Michalsky, J., 1988: The Astronomical Almanac's algorithm for approximate solar position (19502050), Solar Energy 40, 227235
CollaresPereira, M. y Rabl, A., The average distribution of solar radiation: correlations between diffuse and hemispherical and between daily and hourly insolation values. Solar Energy, 22:155–164, 1979.
Perpiñán, O, Energía Solar Fotovoltaica, 2015. (https://oscarperpinan.github.io/esf/)
Perpiñán, O. (2012), "solaR: Solar Radiation and Photovoltaic Systems with R", Journal of Statistical Software, 50(9), 132, doi: 10.18637/jss.v050.i09
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36  ###Angles for one day
BTd = fBTd(mode = 'serie')
#North hemisphere
lat = 37.2
solD < fSolD(lat,BTd[100])
solI < fSolI(solD, sample = 'hour')
print(solI)
#South hemisphere
lat = 37.2;
solDs < fSolD(lat,BTd[283])
solIs < fSolI(solDs, sample = 'hour')
print(solIs)
###Angles for the 12 average days
lat = 37.2;
solD < fSolD(lat,BTd = fBTd(mode = 'prom'))
solI < fSolI(solD, sample = '10 min', keep.night = FALSE)
library(lattice)
library(latticeExtra)
###Solar elevation angle vs. azimuth.
#This kind of graphics is useful for shadows calculations
mon = month.abb
p < xyplot(r2d(AlS)~r2d(AzS),
groups = month,
data = solI, type = 'l', col = 'black',
xlab = expression(psi[s]),ylab = expression(gamma[s]))
plab < p + glayer({
idx < round(length(x)/2+1)
panel.text(x[idx], y[idx], mon[group.value], pos = 3, offset = 0.2, cex = 0.8)})
print(plab)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.