lat=37.2;
Nm=1
SolD<-fSolD(lat,dn=100)
SolI<-fSolI(SolD,Nm=Nm)
G0d=5000
CompDPage<-fCompD(SolD,G0d, corr = "Page")
CompDCPR<-fCompD(SolD,G0d, corr = "CPR")
#define a function fKtd with the correlation of CPR
fKTd=function(x){(0.99*(x<=0.17))+(x>0.17)*(1.188-2.272*x+9.473*x^2-21.856*x^3+14.648*x^4)}
#CompDuser should be the same as CompDCPR
CompDuser<-fCompD(SolD,G0d, corr="user",f=fKTd)
lat=-37.2;
SolDs<-fSolD(lat,dn=283)
SolIs<-fSolI(SolDs,Nm=1)
CompDs<-fCompD(SolDs,G0d, corr = "CPR")
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')
Run the code above in your browser using DataLab