# Read header lines of 5 sbe csv files

cat(readLines("E0501_sbe16_pres_realtime_5m.csv", n=1),'\n\n')    
cat(readLines("E0501_sbe37_pres_realtime_10m.csv", n=1),'\n\n')
cat(readLines("E0501_sbe37_pres_realtime_15m.csv", n=1),'\n\n')
cat(readLines("E0501_sbe37_realtime_20m.csv", n=1),'\n\n')
cat(readLines("E0501_sbe37_realtime_2m.csv", n=1),'\n\n')

# header for 5 files manually renamed
Lsbe_5m=  c('T5','temp_5m','sal_5m','sig_5m','cond_5m','pres_5m')
Lsbe_10m= c('T10','sal_10m','sig_10m','cond_10m','pres_10m','temp_10m')
Lsbe_15m= c('T15','sal_15m','sig_15m','cond_15m','pres_15m','temp_15m')
Lsbe_20m= c('T20','sal_20m','sig_20m','cond_20m','temp_20m')
Lsbe_2m=  c('T2','sal_2m','sig_2m','cond_2m','temp_2m')

X5=read.csv("E0501_sbe16_pres_realtime_5m.csv")
X10=read.csv("E0501_sbe37_pres_realtime_10m.csv")
X15=read.csv("E0501_sbe37_pres_realtime_15m.csv")
X20=read.csv("E0501_sbe37_realtime_20m.csv")
X2=read.csv("E0501_sbe37_realtime_2m.csv")

cat('before\n')
cat(names(X5),'\n\n')
cat(names(X10),'\n\n')
cat(names(X15),'\n\n')
cat(names(X20),'\n\n')
cat(names(X2),'\n\n')

names(X5)=Lsbe_5m
names(X10)=Lsbe_10m
names(X15)=Lsbe_15m
names(X20)=Lsbe_20m
names(X2)=Lsbe_2m

cat('after\n')
cat(names(X5),'\n\n')
cat(names(X10),'\n\n')
cat(names(X15),'\n\n')
cat(names(X20),'\n\n')
cat(names(X2),'\n\n')

sec2day=24*60*60
attach(X2)
attach(X5)
attach(X10)
attach(X15)
attach(X20)
datenow=T2[length(T2)]
datefirst=T2[1]
time2_0=as.double(as.POSIXct(strptime(T2[1], "%Y-%m-%d")))/sec2day
time2=as.double(as.POSIXct(strptime(T2, "%Y-%m-%d %H:%M")))/sec2day - time2_0
time5=as.double(as.POSIXct(strptime(T5, "%Y-%m-%d %H:%M")))/sec2day - time2_0
time10=as.double(as.POSIXct(strptime(T10, "%Y-%m-%d %H:%M")))/sec2day - time2_0
time15=as.double(as.POSIXct(strptime(T15, "%Y-%m-%d %H:%M")))/sec2day - time2_0
time20=as.double(as.POSIXct(strptime(T20, "%Y-%m-%d %H:%M")))/sec2day - time2_0


print(X5[1,])
minl=min(length(temp_2m),length(temp_5m),length(temp_10m),length(temp_15m),length(temp_20m))
mini=1:minl
mint=min(temp_2m,temp_5m,temp_10m,temp_15m,temp_20m)
mint2=min(temp_2m)
maxt=max(temp_2m,temp_5m,temp_10m,temp_15m,temp_20m)
maxtim=max(time2,time5,time10,time15,time20)
plot(c(mint2,maxt),c(mint,maxt),typ='n', main='Temps 5m, 10m, 15m, 20m vs 2m')
points(temp_2m[mini],temp_5m[mini])
points(temp_2m[mini],temp_10m[mini], col='red')
points(temp_2m[mini],temp_15m[mini], col='green')
points(temp_2m[mini],temp_20m[mini], col='blue')

quartz(width=11,height=6)

plot(c(0,maxtim),
     c(mint,maxt),
     typ='n',
     xlab='days',
     ylab='Temp °C', 
     main='Temps 2m, 5m, 10m, 15m, 20m vs days')
points(time2,temp_2m,typ='l', col='black')
#points(time5,temp_5m,typ='l', col='red')
#points(time10,temp_10m,typ='l', col='green')
#points(time15,temp_15m,typ='l', col='blue')
#points(time20,temp_20m,typ='l', col='purple')

out=lowess(time5, temp_5m,f=0.01)
lines(out, col='red', lwd=4)
out=lowess(time10, temp_10m,f=0.01)
lines(out, col='green', lwd=4)
out=lowess(time15, temp_15m,f=0.01)
lines(out, col='blue', lwd=4)
out=lowess(time20, temp_20m,f=0.01)
lines(out, col='purple', lwd=4)

out=lowess(time5, temp_5m,f=0.1)
lines(out, col='red', lwd=4)
out=lowess(time10, temp_10m,f=0.1)
lines(out, col='green', lwd=4)
out=lowess(time15, temp_15m,f=0.1)
lines(out, col='blue', lwd=4)
out=lowess(time20, temp_20m,f=0.1)
lines(out, col='purple', lwd=4)

mtext(datenow, adj=0.95)
mtext(datefirst, adj=0.05)
mtext("Phenomena occuring at 12 and 17 days?")
mtext("Temp 2m", line= -1, adj=0.3)
mtext("Temp 5m", line= -1, adj=0.45, col='red')
mtext("Temp 10m", line= -1, adj=0.6, col='green')
mtext("Temp 15m", line= -1, adj=0.75, col='blue')
mtext("Temp 20m", line= -1, adj=0.9, col='purple')
mtext("2 levels of Smoothed Temperature profiles w thick colored lines", line= -2, adj=0.80, col='orange')
detach(X2)
detach(X5)
detach(X10)
detach(X15)
detach(X20)