This is not a question about R... it is either a bug in the sirad package or normal behaviour for that function. Either way you should be corresponding with the package maintainer. -- Sent from my phone. Please excuse my brevity.
On March 26, 2017 7:05:42 AM PDT, Frederic Ntirenganya <ntfr...@gmail.com> wrote: >Dear R users, > >I do have a challenge in a srcipt of computing Sora Radiation. One line >is >taking long time without giving the results. The line in which I have >problem is : >dat1$RS.daily<-ap(days=days,lat=lat,lon=lon,A=A,B=B,SSD=Sunshine,extraT=NULL) > >Could you please assist me? Thanks > >Here is the code: > >library("sirad", lib.loc="~/R/x86_64-pc-linux-gnu-library/3.0") >library("Evapotranspiration", >lib.loc="~/R/x86_64-pc-linux-gnu-library/3.0") >library("agridat", lib.loc="~/R/x86_64-pc-linux-gnu-library/3.0") >require(zoo) > >dat1<-read.csv("/home/fredo/Documents/Meteo Data/Meteo Rwa >data.csv",header=T,na.string="9999") >dat1$u.daily<-dat1$Wind.Speed*0.514444 >dat1$Date<-as.Date(paste(dat1$Year,dat1$Month,dat1$Day,sep="-")) >dat1$DOY<-dayOfYear(dat1$Date) >dat1$Longitude<-30.11 >dat1$Latitude<--1.95 >dat1$Elevation<-1490 >A <- 0.21 >B <- 0.57 >dat1$RHmax.daily<-dat1$RHmax.daily >dat1$RHmin.daily<-dat1$RHmin.daily >dat1$Tmax.daily<-dat1$Tmax.daily >dat1$Tmin.daily<-dat1$Tmin.daily >dat1<-subset(select=-Wind.direction.measured.in.Degrees) >Sunshine<-dat1$n.daily >lat<-as.numeric(dat1$Latitude) >lon<-as.numeric(dat1$Longitude) >days<-dat1$Date > >dat1$RS.daily<-ap(days=days,lat=lat,lon=lon,A=A,B=B,SSD=Sunshine,extraT=NULL) > >#View(dat1) > >DATA: > >dput(head(dat1)) >structure(list(Year = c(1984L, 1984L, 1984L, 1984L, 1984L, 1984L >), Month = c(1L, 1L, 1L, 1L, 1L, 1L), Day = 1:6, Wind.Speed = c(5L, >4L, 4L, 3L, 5L, 6L), n.daily = c(6.3, 4.8, 0.6, 8.2, 7.3, 1.7 >), Tmax.daily = c(27.4, 26.3, 22.9, 27.7, 28.5, 25.5), Tmin.daily = >c(14.5, >16, 14.4, 14.8, 16.6, 15.4), RHmax.daily = c(100L, 95L, 97L, >100L, 97L, 99L), RHmin.daily = c(45L, 62L, 72L, 55L, 54L, 63L >), Station.Name = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = "KIGALI >AERO", class = "factor"), >Elevation = c(1490, 1490, 1490, 1490, 1490, 1490), Longitude = c(30.11, > 30.11, 30.11, 30.11, 30.11, 30.11), Latitude = c(-1.95, -1.95, > -1.95, -1.95, -1.95, -1.95), u.daily = c(2.57222, 2.057776, > 2.057776, 1.543332, 2.57222, 3.086664), Date = structure(c(5113, > 5114, 5115, 5116, 5117, 5118), class = "Date"), DOY = c(1, > 2, 3, 4, 5, 6)), .Names = c("Year", "Month", "Day", "Wind.Speed", >"n.daily", "Tmax.daily", "Tmin.daily", "RHmax.daily", "RHmin.daily", >"Station.Name", "Elevation", "Longitude", "Latitude", "u.daily", >"Date", "DOY"), row.names = c(NA, 6L), class = "data.frame") > >Any help is appreciated! > >Frederic Ntirenganya >Maseno University, >African Maths Initiative, >Kenya. >Mobile:(+254)718492836 >Email: fr...@aims.ac.za >https://sites.google.com/a/aims.ac.za/fredo/ > > [[alternative HTML version deleted]] > >______________________________________________ >R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >https://stat.ethz.ch/mailman/listinfo/r-help >PLEASE do read the posting guide >http://www.R-project.org/posting-guide.html >and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.