BTd=fBTd(mode='serie')
lat=37.2
sol=calcSol(lat, BTd[100])
print(as.zooD(sol))
library(lattice)
xyplot(as.zooI(sol))
solStrous=calcSol(lat, BTd[100], method='strous')
print(as.zooD(solStrous))
solSpencer=calcSol(lat, BTd[100], method='spencer')
print(as.zooD(solSpencer))
solCooper=calcSol(lat, BTd[100], method='cooper')
print(as.zooD(solCooper))
Run the code above in your browser using DataLab