###Angles for one day
#North hemisphere
Nm=1
lat=37.2
SolD<-fSolD(lat,dn=100)
(SolI<-fSolI(SolD,Nm=Nm))
#South hemisphere
lat=-37.2;
SolDs<-fSolD(lat,dn=283)
(SolIs<-fSolI(SolDs,Nm=1))
###Angles for the 12 average days
lat=37.2;
SolD<-fSolD(lat,BTd=fBTd(Modo='DiasProm'))
SolI<-fSolI(SolD,Nm=6)
library(lattice)
library(latticedl)
###Solar elevation angle vs. azimuth.
#This kind of graphics is useful for shadows calculations
NombreMes=format(ISOdate(2000, 1:12, 1), "%b")
p<-xyplot(AlS*180/pi~AzS*180/pi,
groups=factor(NombreMes[Mes]),data=SolI,
col='black', type='l',
xlab=expression(psi[s]),ylab=expression(gamma[s]))
direct.label(p,method='top.points')
Run the code above in your browser using DataLab