library(maptools) RM215 <- matrix(c(33.5959109, -82.1461363), nrow=1) RM215.sp <- SpatialPoints(RM215, proj4string=CRS("+proj=longlat +datum=WGS84")) d060101 <- as.POSIXct("2006-01-01") study_seq <- seq(from=d060101, length.out=761, by="days") up <- sunriset(RM215, study_seq, direction="sunrise", POSIXct.out=TRUE) #why am i getting NA. I am sure that i am missing something. #Stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis [[alternative HTML version deleted]]