lat=37.2;
Nm=1
SolD<-fSolD(lat,dn=100)
SolI<-fSolI(SolD,Nm=Nm)
G0d=5000
CompD<-fCompD(SolD,G0d)
CompI<-fCompI(CompD,SolI)
AngGen<-fTheta(SolI,beta=35);
Inclin<-fInclin(CompI, AngGen)
xyplot(Gef+Bef+Def~w,data=Inclin,type='l', auto.key=TRUE)
lat=37.2;
G0dm=c(2.766,3.491,4.494,5.912,6.989,7.742,7.919,7.027,5.369,3.562,2.814,2.179)*1000;
SolD<-fSolD(lat,BTd=fBTd(Modo='DiasProm'))
SolI<-fSolI(SolD,Nm=6)
CompD<-fCompD(SolD,G0d=G0dm, corr = 'Page')
CompD$Ta=25
CompI<-fCompI(CompD,SolI)
AngGen<-fTheta(SolI,beta=35);
Inclin<-fInclin(CompI, AngGen)
xyplot(Gef+Bef+Def~w|Mes,data=Inclin,type='l', auto.key=TRUE)
Run the code above in your browser using DataLab