Line 35 appears incorrect:
interval<-interval*60*60
interval is already in hours (from line 18) but this appears to be converting from seconds to hours and needs to be checked.
