quartz(width=18, height=7)

xx<-read.csv("8JulyTides.csv")
attach(xx)
day<- sort(as.real(Day))
junk<- as.POSIXlt(strptime(PMhi, "%H:%M"))
hrs<-junk$hour
hrs<- 12 + hrs - 12*(hrs==12)
mins<-junk$min
PMh<-hrs+mins/60
junk<- as.POSIXlt(strptime(PMlo, "%H:%M"))
hrs<-junk$hour
hrs<- 12 + hrs - 12*(hrs==12)
mins<-junk$min
PMl<-hrs+mins/60
Dayhi<- (day-1)+(PMh)/24
Daylo<- (day-1)+(PMl)/24
par(mar= c(5,5,4,2))
plot(c(0,62), c(0,24), typ='n', 
     main="Old Orchard Beach PM Tides July & August 2008, black = high, red = low",
     ylab="PM hi/lo black/red",
     xlab="day of the month",
     cex=1,
     cex.main=2,
     cex.axis=1.5,
     cex.lab=1.5)
points(Dayhi, PMh, cex=2, lwd=2)
points(Daylo,PMl, col='red', cex=2, lwd=2)
junk<- as.POSIXlt(strptime(AMhi, "%H:%M"))
hrs<-junk$hour
hrs<- hrs - 12*(hrs==12)
mins<-junk$min
AMh<-hrs+mins/60
junk<- as.POSIXlt(strptime(AMlo, "%H:%M"))
hrs<-junk$hour
hrs<- hrs - 12*(hrs==12)
mins<-junk$min
AMl<-hrs+mins/60
Dayhi<- (day-1)+(AMh)/24
Daylo<- (day-1)+(AMl)/24
points(Dayhi, AMh, cex=2, lwd=2)
points(Daylo,AMl, col='red', cex=2, lwd=2)
browser()
detach(xx)

xx<-read.csv("8AugTides.csv")
attach(xx)
day<- 31+ sort(as.real(Day))
junk<- as.POSIXlt(strptime(PMhi, "%H:%M"))
hrs<-junk$hour
hrs<- 12 + hrs - 12*(hrs==12)
mins<-junk$min
PMh<-hrs+mins/60
junk<- as.POSIXlt(strptime(PMlo, "%H:%M"))
hrs<-junk$hour
hrs<- 12 + hrs - 12*(hrs==12)
mins<-junk$min
PMl<-hrs+mins/60
Dayhi<- (day-1)+(PMh)/24
Daylo<- (day-1)+(PMl)/24
par(new=TRUE, ann=FALSE)
plot(c(0,62), c(0,24), typ='n', 
     main="Old Orchard Beach PM Tides Aug 2008, black = high, red = low",
     ylab="PM hi/lo black/red",
     xlab="day of the month",
     xaxt='n', yaxt='n',
     cex=1,
     cex.main=3)
points(Dayhi, PMh, cex=2, lwd=2)
points(Daylo,PMl, col='red', cex=2, lwd=2)
junk<- as.POSIXlt(strptime(AMhi, "%H:%M"))
hrs<-junk$hour
hrs<- hrs - 12*(hrs==12)
mins<-junk$min
AMh<-hrs+mins/60
junk<- as.POSIXlt(strptime(AMlo, "%H:%M"))
hrs<-junk$hour
hrs<- hrs - 12*(hrs==12)
mins<-junk$min
AMl<-hrs+mins/60
Dayhi<- (day-1)+(AMh)/24
Daylo<- (day-1)+(AMl)/24
points(Dayhi, AMh, cex=2, lwd=2)
points(Daylo,AMl, col='red', cex=2, lwd=2)
browser()
detach(xx)