Work Shop 2
Natasha Langton
5 March 2019
Load Appropriate libraries
library(readr)
library(maps)Set Wd
setwd("~/College/Msc Climate Change/GY667")Check directory for file name - Arklow
dir()## [1] "~$rkshop 2 instructions.docx"
## [2] "arklow.txt"
## [3] "arklowmonthlydata.csv"
## [4] "Assignment 1.docx"
## [5] "global_temps_monthly.csv"
## [6] "HadCRUT.4.6.0.0.median.nc"
## [7] "marking-from-WS1.xlsx"
## [8] "New folder"
## [9] "Tide Gauge.htm"
## [10] "Workshop 2 instructions.docx"
## [11] "Workshop 2.Rmd"
## [12] "workshop_1_instructions_script.R"
## [13] "workshop_1_solutions_script.R"
## [14] "Workshop_2.html"
## [15] "Workshop_2.Rmd"
## [16] "Workshop1-Guidelines_and_Instructions.docx"
## [17] "workshop2-skeleton-script.R"To open up the Arklow data use read.table function. Need to skip the first 6 lines as it only contains metadata.
Arklow = read.table(file.path(getwd(),'arklow.txt'),skip = 6,header=TRUE, sep = "\t")Investigate
head(Arklow)## Date Value Quality
## 1 2003/08/26 00:00:00 0.067 31
## 2 2003/08/26 00:15:00 0.003 31
## 3 2003/08/26 00:30:00 -0.050 31
## 4 2003/08/26 00:45:00 -0.102 31
## 5 2003/08/26 01:00:00 -0.161 31
## 6 2003/08/26 01:15:00 -0.214 31colnames(Arklow)## [1] "Date" "Value" "Quality"Map the location of Arklow Tide Gauge, as latitude and longitude data is not available in the data set, get the latitude and longitude information from https://erddap.marine.ie/erddap/tabledap/IrishNationalTideGaugeNetwork.subset?station_id,longitude,latitude&station_id=%22OPW%20Arklow%20Harbour%22&.viewDistinctMap=true&.viewDistinctData=1000&.viewRelatedData=0&distinct()&.last=station_id
map("world",c("ireland","uk"),fill=TRUE,xlim=c(-12,-4),ylim=c(51,56))
map.axes(cex.axis=1)
title(main="Location of Arklow Tide Gauge",xlab="Longitude",ylab="Latitude")
points(-6.145231,52.792046
,pch=21,col="gray",bg="red")
The red point indicates Arklows position on the map.
Need to Convert your data to POSIXlt date format.
Arklow.time <- as.POSIXlt(Arklow$Date,format="%Y/%m/%d %H:%M:%S",tz='UTC')Investigate
head(Arklow.time)## [1] "2003-08-26 00:00:00 UTC" "2003-08-26 00:15:00 UTC"
## [3] "2003-08-26 00:30:00 UTC" "2003-08-26 00:45:00 UTC"
## [5] "2003-08-26 01:00:00 UTC" "2003-08-26 01:15:00 UTC"Make a new data frame with variables year, month, day, hour, minute, second, sea level (h), and flag
ark <- data.frame("time"=Arklow.time,
"year"=Arklow.time$year+1900,
"month"=Arklow.time$mon+1,
"day"=Arklow.time$mday,
"hour"=Arklow.time$hour+1,
"min"=Arklow.time$min,
"sec"=Arklow.time$sec,
"h"=Arklow$Value,
"flag"=Arklow$Quality)Make an hourly data frame before aggregating by hour, chop off un even bits at start and end
head(ark, 10)## time year month day hour min sec h flag
## 1 2003-08-26 00:00:00 2003 8 26 1 0 0 0.067 31
## 2 2003-08-26 00:15:00 2003 8 26 1 15 0 0.003 31
## 3 2003-08-26 00:30:00 2003 8 26 1 30 0 -0.050 31
## 4 2003-08-26 00:45:00 2003 8 26 1 45 0 -0.102 31
## 5 2003-08-26 01:00:00 2003 8 26 2 0 0 -0.161 31
## 6 2003-08-26 01:15:00 2003 8 26 2 15 0 -0.214 31
## 7 2003-08-26 01:30:00 2003 8 26 2 30 0 -0.281 31
## 8 2003-08-26 01:45:00 2003 8 26 2 45 0 -0.355 31
## 9 2003-08-26 02:00:00 2003 8 26 3 0 0 -0.388 31
## 10 2003-08-26 02:15:00 2003 8 26 3 15 0 -0.435 31tail(ark, 10)## time year month day hour min sec h flag
## 520708 2019-01-31 21:45:00 2019 1 31 22 45 0 0.516 *
## 520709 2019-01-31 22:00:00 2019 1 31 23 0 0 0.511 *
## 520710 2019-01-31 22:15:00 2019 1 31 23 15 0 0.546 *
## 520711 2019-01-31 22:30:00 2019 1 31 23 30 0 0.522 *
## 520712 2019-01-31 22:45:00 2019 1 31 23 45 0 0.489 *
## 520713 2019-01-31 23:00:00 2019 1 31 24 0 0 0.416 *
## 520714 2019-01-31 23:15:00 2019 1 31 24 15 0 0.352 *
## 520715 2019-01-31 23:30:00 2019 1 31 24 30 0 0.253 *
## 520716 2019-01-31 23:45:00 2019 1 31 24 45 0 0.247 *
## 520717 2019-02-01 00:00:00 2019 2 1 1 0 0 0.139 *From the output of head and tail of data set ark, it is clear that I need to start at 1 and end at 520716.
ark.sub <- ark[1:520716,]
ark.hour <- aggregate(h~hour+day+month+year,ark.sub,mean)Create a new empty column for a decimal year
ark.hour <- cbind(ark.hour, NA)
names(ark.hour) <- c("hour","day","month", "year", "h", "time")Add a decimal year variable
ark.hour$time <- as.POSIXlt(sprintf("%s/%s/%s %s",
ark.hour$year, ark.hour$month,
ark.hour$day, ark.hour$hour),
format="%Y/%m/%d %H",tz='UTC')Investigate the data
summary(ark.hour)## hour day month year
## Min. : 1.0 Min. : 1.00 Min. : 1.000 Min. :2003
## 1st Qu.: 6.0 1st Qu.: 8.00 1st Qu.: 4.000 1st Qu.:2007
## Median :12.0 Median :16.00 Median : 7.000 Median :2011
## Mean :12.5 Mean :15.74 Mean : 6.565 Mean :2011
## 3rd Qu.:18.0 3rd Qu.:23.00 3rd Qu.:10.000 3rd Qu.:2015
## Max. :24.0 Max. :31.00 Max. :12.000 Max. :2019
## h time
## Min. :-1.4748 Min. :2003-08-26 01:00:00
## 1st Qu.:-0.2760 1st Qu.:2007-06-28 01:15:00
## Median :-0.0410 Median :2011-04-17 07:30:00
## Mean :-0.0657 Mean :2011-05-25 20:07:37
## 3rd Qu.: 0.1550 3rd Qu.:2015-05-10 17:45:00
## Max. : 1.2517 Max. :2019-02-01 00:00:00head(ark.hour)## hour day month year h time
## 1 1 26 8 2003 -0.02050 2003-08-26 01:00:00
## 2 2 26 8 2003 -0.25275 2003-08-26 02:00:00
## 3 3 26 8 2003 -0.44150 2003-08-26 03:00:00
## 4 4 26 8 2003 -0.48425 2003-08-26 04:00:00
## 5 5 26 8 2003 -0.37725 2003-08-26 05:00:00
## 6 6 26 8 2003 -0.16650 2003-08-26 06:00:00tail(ark.hour)## hour day month year h time
## 130189 19 31 1 2019 0.24550 2019-01-31 19:00:00
## 130190 20 31 1 2019 0.25225 2019-01-31 20:00:00
## 130191 21 31 1 2019 0.39450 2019-01-31 21:00:00
## 130192 22 31 1 2019 0.50650 2019-01-31 22:00:00
## 130193 23 31 1 2019 0.51700 2019-01-31 23:00:00
## 130194 24 31 1 2019 0.31700 2019-02-01 00:00:00Repeat the same steps for monthly data.
head(ark.hour, 200)## hour day month year h time
## 1 1 26 8 2003 -0.02050 2003-08-26 01:00:00
## 2 2 26 8 2003 -0.25275 2003-08-26 02:00:00
## 3 3 26 8 2003 -0.44150 2003-08-26 03:00:00
## 4 4 26 8 2003 -0.48425 2003-08-26 04:00:00
## 5 5 26 8 2003 -0.37725 2003-08-26 05:00:00
## 6 6 26 8 2003 -0.16650 2003-08-26 06:00:00
## 7 7 26 8 2003 -0.00825 2003-08-26 07:00:00
## 8 8 26 8 2003 0.03350 2003-08-26 08:00:00
## 9 9 26 8 2003 -0.00100 2003-08-26 09:00:00
## 10 10 26 8 2003 -0.00900 2003-08-26 10:00:00
## 11 11 26 8 2003 0.05300 2003-08-26 11:00:00
## 12 12 26 8 2003 0.05675 2003-08-26 12:00:00
## 13 13 26 8 2003 -0.04025 2003-08-26 13:00:00
## 14 14 26 8 2003 -0.21275 2003-08-26 14:00:00
## 15 15 26 8 2003 -0.36475 2003-08-26 15:00:00
## 16 16 26 8 2003 -0.41900 2003-08-26 16:00:00
## 17 17 26 8 2003 -0.29300 2003-08-26 17:00:00
## 18 18 26 8 2003 -0.04925 2003-08-26 18:00:00
## 19 19 26 8 2003 0.17625 2003-08-26 19:00:00
## 20 20 26 8 2003 0.27725 2003-08-26 20:00:00
## 21 21 26 8 2003 0.26850 2003-08-26 21:00:00
## 22 22 26 8 2003 0.24225 2003-08-26 22:00:00
## 23 23 26 8 2003 0.27375 2003-08-26 23:00:00
## 24 24 26 8 2003 0.25950 2003-08-27 00:00:00
## 25 1 27 8 2003 0.09125 2003-08-27 01:00:00
## 26 2 27 8 2003 -0.15000 2003-08-27 02:00:00
## 27 3 27 8 2003 -0.42950 2003-08-27 03:00:00
## 28 4 27 8 2003 -0.55750 2003-08-27 04:00:00
## 29 5 27 8 2003 -0.52875 2003-08-27 05:00:00
## 30 6 27 8 2003 -0.30050 2003-08-27 06:00:00
## 31 7 27 8 2003 -0.02350 2003-08-27 07:00:00
## 32 8 27 8 2003 0.13525 2003-08-27 08:00:00
## 33 9 27 8 2003 0.12975 2003-08-27 09:00:00
## 34 10 27 8 2003 0.07175 2003-08-27 10:00:00
## 35 11 27 8 2003 0.08750 2003-08-27 11:00:00
## 36 12 27 8 2003 0.09050 2003-08-27 12:00:00
## 37 13 27 8 2003 0.00600 2003-08-27 13:00:00
## 38 14 27 8 2003 -0.15375 2003-08-27 14:00:00
## 39 15 27 8 2003 -0.37325 2003-08-27 15:00:00
## 40 16 27 8 2003 -0.46075 2003-08-27 16:00:00
## 41 17 27 8 2003 -0.38425 2003-08-27 17:00:00
## 42 18 27 8 2003 -0.13050 2003-08-27 18:00:00
## 43 19 27 8 2003 0.20125 2003-08-27 19:00:00
## 44 20 27 8 2003 0.41500 2003-08-27 20:00:00
## 45 21 27 8 2003 0.45050 2003-08-27 21:00:00
## 46 22 27 8 2003 0.37900 2003-08-27 22:00:00
## 47 23 27 8 2003 0.34100 2003-08-27 23:00:00
## 48 24 27 8 2003 0.33850 2003-08-28 00:00:00
## 49 1 28 8 2003 0.20700 2003-08-28 01:00:00
## 50 2 28 8 2003 -0.04425 2003-08-28 02:00:00
## 51 3 28 8 2003 -0.38575 2003-08-28 03:00:00
## 52 4 28 8 2003 -0.64675 2003-08-28 04:00:00
## 53 5 28 8 2003 -0.70475 2003-08-28 05:00:00
## 54 6 28 8 2003 -0.52675 2003-08-28 06:00:00
## 55 7 28 8 2003 -0.18150 2003-08-28 07:00:00
## 56 8 28 8 2003 0.11175 2003-08-28 08:00:00
## 57 9 28 8 2003 0.21900 2003-08-28 09:00:00
## 58 10 28 8 2003 0.17600 2003-08-28 10:00:00
## 59 11 28 8 2003 0.10500 2003-08-28 11:00:00
## 60 12 28 8 2003 0.12500 2003-08-28 12:00:00
## 61 13 28 8 2003 0.08550 2003-08-28 13:00:00
## 62 14 28 8 2003 -0.06525 2003-08-28 14:00:00
## 63 15 28 8 2003 -0.30575 2003-08-28 15:00:00
## 64 16 28 8 2003 -0.52825 2003-08-28 16:00:00
## 65 17 28 8 2003 -0.56975 2003-08-28 17:00:00
## 66 18 28 8 2003 -0.38550 2003-08-28 18:00:00
## 67 19 28 8 2003 -0.01575 2003-08-28 19:00:00
## 68 20 28 8 2003 0.35900 2003-08-28 20:00:00
## 69 21 28 8 2003 0.53150 2003-08-28 21:00:00
## 70 22 28 8 2003 0.52600 2003-08-28 22:00:00
## 71 23 28 8 2003 0.39975 2003-08-28 23:00:00
## 72 24 28 8 2003 0.36375 2003-08-29 00:00:00
## 73 1 29 8 2003 0.28650 2003-08-29 01:00:00
## 74 2 29 8 2003 0.06250 2003-08-29 02:00:00
## 75 3 29 8 2003 -0.25925 2003-08-29 03:00:00
## 76 4 29 8 2003 -0.61275 2003-08-29 04:00:00
## 77 5 29 8 2003 -0.81650 2003-08-29 05:00:00
## 78 6 29 8 2003 -0.77100 2003-08-29 06:00:00
## 79 7 29 8 2003 -0.48725 2003-08-29 07:00:00
## 80 8 29 8 2003 -0.12425 2003-08-29 08:00:00
## 81 9 29 8 2003 0.15325 2003-08-29 09:00:00
## 82 10 29 8 2003 0.21125 2003-08-29 10:00:00
## 83 11 29 8 2003 0.13725 2003-08-29 11:00:00
## 84 12 29 8 2003 0.05700 2003-08-29 12:00:00
## 85 13 29 8 2003 0.05025 2003-08-29 13:00:00
## 86 14 29 8 2003 -0.08300 2003-08-29 14:00:00
## 87 15 29 8 2003 -0.29675 2003-08-29 15:00:00
## 88 16 29 8 2003 -0.55375 2003-08-29 16:00:00
## 89 17 29 8 2003 -0.71250 2003-08-29 17:00:00
## 90 18 29 8 2003 -0.65225 2003-08-29 18:00:00
## 91 19 29 8 2003 -0.34300 2003-08-29 19:00:00
## 92 20 29 8 2003 0.11175 2003-08-29 20:00:00
## 93 21 29 8 2003 0.44925 2003-08-29 21:00:00
## 94 22 29 8 2003 0.55550 2003-08-29 22:00:00
## 95 23 29 8 2003 0.47950 2003-08-29 23:00:00
## 96 24 29 8 2003 0.36225 2003-08-30 00:00:00
## 97 1 30 8 2003 0.31000 2003-08-30 01:00:00
## 98 2 30 8 2003 0.15625 2003-08-30 02:00:00
## 99 3 30 8 2003 -0.10900 2003-08-30 03:00:00
## 100 4 30 8 2003 -0.46750 2003-08-30 04:00:00
## 101 5 30 8 2003 -0.79525 2003-08-30 05:00:00
## 102 6 30 8 2003 -0.89450 2003-08-30 06:00:00
## 103 7 30 8 2003 -0.72900 2003-08-30 07:00:00
## 104 8 30 8 2003 -0.34500 2003-08-30 08:00:00
## 105 9 30 8 2003 0.03725 2003-08-30 09:00:00
## 106 10 30 8 2003 0.24650 2003-08-30 10:00:00
## 107 11 30 8 2003 0.22150 2003-08-30 11:00:00
## 108 12 30 8 2003 0.10225 2003-08-30 12:00:00
## 109 13 30 8 2003 0.07950 2003-08-30 13:00:00
## 110 14 30 8 2003 0.03925 2003-08-30 14:00:00
## 111 15 30 8 2003 -0.15525 2003-08-30 15:00:00
## 112 16 30 8 2003 -0.40150 2003-08-30 16:00:00
## 113 17 30 8 2003 -0.66800 2003-08-30 17:00:00
## 114 18 30 8 2003 -0.76725 2003-08-30 18:00:00
## 115 19 30 8 2003 -0.62975 2003-08-30 19:00:00
## 116 20 30 8 2003 -0.23975 2003-08-30 20:00:00
## 117 21 30 8 2003 0.20300 2003-08-30 21:00:00
## 118 22 30 8 2003 0.46875 2003-08-30 22:00:00
## 119 23 30 8 2003 0.49875 2003-08-30 23:00:00
## 120 24 30 8 2003 0.38225 2003-08-31 00:00:00
## 121 1 31 8 2003 0.29850 2003-08-31 01:00:00
## 122 2 31 8 2003 0.25250 2003-08-31 02:00:00
## 123 3 31 8 2003 0.07225 2003-08-31 03:00:00
## 124 4 31 8 2003 -0.21975 2003-08-31 04:00:00
## 125 5 31 8 2003 -0.57500 2003-08-31 05:00:00
## 126 6 31 8 2003 -0.83550 2003-08-31 06:00:00
## 127 7 31 8 2003 -0.85300 2003-08-31 07:00:00
## 128 8 31 8 2003 -0.62125 2003-08-31 08:00:00
## 129 9 31 8 2003 -0.22475 2003-08-31 09:00:00
## 130 10 31 8 2003 0.09450 2003-08-31 10:00:00
## 131 11 31 8 2003 0.19950 2003-08-31 11:00:00
## 132 12 31 8 2003 0.11675 2003-08-31 12:00:00
## 133 13 31 8 2003 0.01500 2003-08-31 13:00:00
## 134 14 31 8 2003 0.02975 2003-08-31 14:00:00
## 135 15 31 8 2003 -0.06600 2003-08-31 15:00:00
## 136 16 31 8 2003 -0.23975 2003-08-31 16:00:00
## 137 17 31 8 2003 -0.48300 2003-08-31 17:00:00
## 138 18 31 8 2003 -0.71575 2003-08-31 18:00:00
## 139 19 31 8 2003 -0.74100 2003-08-31 19:00:00
## 140 20 31 8 2003 -0.52175 2003-08-31 20:00:00
## 141 21 31 8 2003 -0.11725 2003-08-31 21:00:00
## 142 22 31 8 2003 0.25225 2003-08-31 22:00:00
## 143 23 31 8 2003 0.42950 2003-08-31 23:00:00
## 144 24 31 8 2003 0.39950 2003-09-01 00:00:00
## 145 1 1 9 2003 0.27425 2003-09-01 01:00:00
## 146 2 1 9 2003 0.24625 2003-09-01 02:00:00
## 147 3 1 9 2003 0.18275 2003-09-01 03:00:00
## 148 4 1 9 2003 0.00600 2003-09-01 04:00:00
## 149 5 1 9 2003 -0.28025 2003-09-01 05:00:00
## 150 6 1 9 2003 -0.61050 2003-09-01 06:00:00
## 151 7 1 9 2003 -0.80750 2003-09-01 07:00:00
## 152 8 1 9 2003 -0.78025 2003-09-01 08:00:00
## 153 9 1 9 2003 -0.52450 2003-09-01 09:00:00
## 154 10 1 9 2003 -0.16900 2003-09-01 10:00:00
## 155 11 1 9 2003 0.05750 2003-09-01 11:00:00
## 156 12 1 9 2003 0.08550 2003-09-01 12:00:00
## 157 13 1 9 2003 -0.00225 2003-09-01 13:00:00
## 158 14 1 9 2003 -0.03850 2003-09-01 14:00:00
## 159 15 1 9 2003 -0.00550 2003-09-01 15:00:00
## 160 16 1 9 2003 -0.08125 2003-09-01 16:00:00
## 161 17 1 9 2003 -0.24400 2003-09-01 17:00:00
## 162 18 1 9 2003 -0.47600 2003-09-01 18:00:00
## 163 19 1 9 2003 -0.66750 2003-09-01 19:00:00
## 164 20 1 9 2003 -0.67150 2003-09-01 20:00:00
## 165 21 1 9 2003 -0.45925 2003-09-01 21:00:00
## 166 22 1 9 2003 -0.10900 2003-09-01 22:00:00
## 167 23 1 9 2003 0.16700 2003-09-01 23:00:00
## 168 24 1 9 2003 0.26325 2003-09-02 00:00:00
## 169 1 2 9 2003 0.20300 2003-09-02 01:00:00
## 170 2 2 9 2003 0.12525 2003-09-02 02:00:00
## 171 3 2 9 2003 0.16425 2003-09-02 03:00:00
## 172 4 2 9 2003 0.13375 2003-09-02 04:00:00
## 173 5 2 9 2003 0.00600 2003-09-02 05:00:00
## 174 6 2 9 2003 -0.23625 2003-09-02 06:00:00
## 175 7 2 9 2003 -0.51050 2003-09-02 07:00:00
## 176 8 2 9 2003 -0.67425 2003-09-02 08:00:00
## 177 9 2 9 2003 -0.66000 2003-09-02 09:00:00
## 178 10 2 9 2003 -0.44600 2003-09-02 10:00:00
## 179 11 2 9 2003 -0.19350 2003-09-02 11:00:00
## 180 12 2 9 2003 -0.04150 2003-09-02 12:00:00
## 181 13 2 9 2003 -0.04825 2003-09-02 13:00:00
## 182 14 2 9 2003 -0.10975 2003-09-02 14:00:00
## 183 15 2 9 2003 -0.08800 2003-09-02 15:00:00
## 184 16 2 9 2003 -0.01250 2003-09-02 16:00:00
## 185 17 2 9 2003 -0.06525 2003-09-02 17:00:00
## 186 18 2 9 2003 -0.17950 2003-09-02 18:00:00
## 187 19 2 9 2003 -0.38825 2003-09-02 19:00:00
## 188 20 2 9 2003 -0.55575 2003-09-02 20:00:00
## 189 21 2 9 2003 -0.57525 2003-09-02 21:00:00
## 190 22 2 9 2003 -0.40875 2003-09-02 22:00:00
## 191 23 2 9 2003 -0.14700 2003-09-02 23:00:00
## 192 24 2 9 2003 0.04525 2003-09-03 00:00:00
## 193 1 3 9 2003 0.08500 2003-09-03 01:00:00
## 194 2 3 9 2003 0.03000 2003-09-03 02:00:00
## 195 3 3 9 2003 0.00675 2003-09-03 03:00:00
## 196 4 3 9 2003 0.10550 2003-09-03 04:00:00
## 197 5 3 9 2003 0.12925 2003-09-03 05:00:00
## 198 6 3 9 2003 0.06075 2003-09-03 06:00:00
## 199 7 3 9 2003 -0.12450 2003-09-03 07:00:00
## 200 8 3 9 2003 -0.37575 2003-09-03 08:00:00From this it is clear that the start of a month starts at line 145, which is the 1st of September 2003.
tail(ark.hour,200)## hour day month year h time
## 129995 17 23 1 2019 -0.77450 2019-01-23 17:00:00
## 129996 18 23 1 2019 -1.00400 2019-01-23 18:00:00
## 129997 19 23 1 2019 -0.90150 2019-01-23 19:00:00
## 129998 20 23 1 2019 -0.53525 2019-01-23 20:00:00
## 129999 21 23 1 2019 -0.08200 2019-01-23 21:00:00
## 130000 22 23 1 2019 0.24100 2019-01-23 22:00:00
## 130001 23 23 1 2019 0.29225 2019-01-23 23:00:00
## 130002 24 23 1 2019 0.18500 2019-01-24 00:00:00
## 130003 1 24 1 2019 0.12100 2019-01-24 01:00:00
## 130004 2 24 1 2019 0.12875 2019-01-24 02:00:00
## 130005 3 24 1 2019 -0.00025 2019-01-24 03:00:00
## 130006 4 24 1 2019 -0.24075 2019-01-24 04:00:00
## 130007 5 24 1 2019 -0.54575 2019-01-24 05:00:00
## 130008 6 24 1 2019 -0.76700 2019-01-24 06:00:00
## 130009 7 24 1 2019 -0.73400 2019-01-24 07:00:00
## 130010 8 24 1 2019 -0.42925 2019-01-24 08:00:00
## 130011 9 24 1 2019 0.03750 2019-01-24 09:00:00
## 130012 10 24 1 2019 0.42125 2019-01-24 10:00:00
## 130013 11 24 1 2019 0.55350 2019-01-24 11:00:00
## 130014 12 24 1 2019 0.47525 2019-01-24 12:00:00
## 130015 13 24 1 2019 0.33675 2019-01-24 13:00:00
## 130016 14 24 1 2019 0.32275 2019-01-24 14:00:00
## 130017 15 24 1 2019 0.22325 2019-01-24 15:00:00
## 130018 16 24 1 2019 -0.02875 2019-01-24 16:00:00
## 130019 17 24 1 2019 -0.41900 2019-01-24 17:00:00
## 130020 18 24 1 2019 -0.78525 2019-01-24 18:00:00
## 130021 19 24 1 2019 -0.93875 2019-01-24 19:00:00
## 130022 20 24 1 2019 -0.85200 2019-01-24 20:00:00
## 130023 21 24 1 2019 -0.48750 2019-01-24 21:00:00
## 130024 22 24 1 2019 -0.06350 2019-01-24 22:00:00
## 130025 23 24 1 2019 0.16475 2019-01-24 23:00:00
## 130026 24 24 1 2019 0.16200 2019-01-25 00:00:00
## 130027 1 25 1 2019 0.05300 2019-01-25 01:00:00
## 130028 2 25 1 2019 0.00775 2019-01-25 02:00:00
## 130029 3 25 1 2019 0.03175 2019-01-25 03:00:00
## 130030 4 25 1 2019 -0.06950 2019-01-25 04:00:00
## 130031 5 25 1 2019 -0.26300 2019-01-25 05:00:00
## 130032 6 25 1 2019 -0.51625 2019-01-25 06:00:00
## 130033 7 25 1 2019 -0.67450 2019-01-25 07:00:00
## 130034 8 25 1 2019 -0.61025 2019-01-25 08:00:00
## 130035 9 25 1 2019 -0.27775 2019-01-25 09:00:00
## 130036 10 25 1 2019 0.12625 2019-01-25 10:00:00
## 130037 11 25 1 2019 0.45425 2019-01-25 11:00:00
## 130038 12 25 1 2019 0.52900 2019-01-25 12:00:00
## 130039 13 25 1 2019 0.45200 2019-01-25 13:00:00
## 130040 14 25 1 2019 0.37200 2019-01-25 14:00:00
## 130041 15 25 1 2019 0.37850 2019-01-25 15:00:00
## 130042 16 25 1 2019 0.25450 2019-01-25 16:00:00
## 130043 17 25 1 2019 0.03450 2019-01-25 17:00:00
## 130044 18 25 1 2019 -0.32125 2019-01-25 18:00:00
## 130045 19 25 1 2019 -0.65675 2019-01-25 19:00:00
## 130046 20 25 1 2019 -0.79800 2019-01-25 20:00:00
## 130047 21 25 1 2019 -0.71050 2019-01-25 21:00:00
## 130048 22 25 1 2019 -0.40150 2019-01-25 22:00:00
## 130049 23 25 1 2019 -0.07225 2019-01-25 23:00:00
## 130050 24 25 1 2019 0.08500 2019-01-26 00:00:00
## 130051 1 26 1 2019 0.08850 2019-01-26 01:00:00
## 130052 2 26 1 2019 0.01575 2019-01-26 02:00:00
## 130053 3 26 1 2019 0.05025 2019-01-26 03:00:00
## 130054 4 26 1 2019 0.09925 2019-01-26 04:00:00
## 130055 5 26 1 2019 0.02900 2019-01-26 05:00:00
## 130056 6 26 1 2019 -0.13650 2019-01-26 06:00:00
## 130057 7 26 1 2019 -0.34800 2019-01-26 07:00:00
## 130058 8 26 1 2019 -0.46075 2019-01-26 08:00:00
## 130059 9 26 1 2019 -0.37050 2019-01-26 09:00:00
## 130060 10 26 1 2019 -0.06350 2019-01-26 10:00:00
## 130061 11 26 1 2019 0.30325 2019-01-26 11:00:00
## 130062 12 26 1 2019 0.55350 2019-01-26 12:00:00
## 130063 13 26 1 2019 0.62725 2019-01-26 13:00:00
## 130064 14 26 1 2019 0.56875 2019-01-26 14:00:00
## 130065 15 26 1 2019 0.57675 2019-01-26 15:00:00
## 130066 16 26 1 2019 0.67025 2019-01-26 16:00:00
## 130067 17 26 1 2019 0.66225 2019-01-26 17:00:00
## 130068 18 26 1 2019 0.50600 2019-01-26 18:00:00
## 130069 19 26 1 2019 0.19400 2019-01-26 19:00:00
## 130070 20 26 1 2019 -0.14825 2019-01-26 20:00:00
## 130071 21 26 1 2019 -0.31825 2019-01-26 21:00:00
## 130072 22 26 1 2019 -0.30725 2019-01-26 22:00:00
## 130073 23 26 1 2019 -0.13100 2019-01-26 23:00:00
## 130074 24 26 1 2019 0.07225 2019-01-27 00:00:00
## 130075 1 27 1 2019 0.14975 2019-01-27 01:00:00
## 130076 2 27 1 2019 0.09450 2019-01-27 02:00:00
## 130077 3 27 1 2019 -0.01300 2019-01-27 03:00:00
## 130078 4 27 1 2019 -0.07975 2019-01-27 04:00:00
## 130079 5 27 1 2019 -0.16925 2019-01-27 05:00:00
## 130080 6 27 1 2019 -0.35750 2019-01-27 06:00:00
## 130081 7 27 1 2019 -0.52150 2019-01-27 07:00:00
## 130082 8 27 1 2019 -0.73475 2019-01-27 08:00:00
## 130083 9 27 1 2019 -0.85850 2019-01-27 09:00:00
## 130084 10 27 1 2019 -0.76550 2019-01-27 10:00:00
## 130085 11 27 1 2019 -0.55125 2019-01-27 11:00:00
## 130086 12 27 1 2019 -0.24200 2019-01-27 12:00:00
## 130087 13 27 1 2019 -0.01550 2019-01-27 13:00:00
## 130088 14 27 1 2019 0.05250 2019-01-27 14:00:00
## 130089 15 27 1 2019 0.04550 2019-01-27 15:00:00
## 130090 16 27 1 2019 0.00775 2019-01-27 16:00:00
## 130091 17 27 1 2019 0.01975 2019-01-27 17:00:00
## 130092 18 27 1 2019 0.00075 2019-01-27 18:00:00
## 130093 19 27 1 2019 -0.11225 2019-01-27 19:00:00
## 130094 20 27 1 2019 -0.33800 2019-01-27 20:00:00
## 130095 21 27 1 2019 -0.55500 2019-01-27 21:00:00
## 130096 22 27 1 2019 -0.66150 2019-01-27 22:00:00
## 130097 23 27 1 2019 -0.63250 2019-01-27 23:00:00
## 130098 24 27 1 2019 -0.45000 2019-01-28 00:00:00
## 130099 1 28 1 2019 -0.28100 2019-01-28 01:00:00
## 130100 2 28 1 2019 -0.21325 2019-01-28 02:00:00
## 130101 3 28 1 2019 -0.27875 2019-01-28 03:00:00
## 130102 4 28 1 2019 -0.28300 2019-01-28 04:00:00
## 130103 5 28 1 2019 -0.22500 2019-01-28 05:00:00
## 130104 6 28 1 2019 -0.12875 2019-01-28 06:00:00
## 130105 7 28 1 2019 -0.08450 2019-01-28 07:00:00
## 130106 8 28 1 2019 -0.18575 2019-01-28 08:00:00
## 130107 9 28 1 2019 -0.33250 2019-01-28 09:00:00
## 130108 10 28 1 2019 -0.45950 2019-01-28 10:00:00
## 130109 11 28 1 2019 -0.46175 2019-01-28 11:00:00
## 130110 12 28 1 2019 -0.36675 2019-01-28 12:00:00
## 130111 13 28 1 2019 -0.16175 2019-01-28 13:00:00
## 130112 14 28 1 2019 -0.08500 2019-01-28 14:00:00
## 130113 15 28 1 2019 -0.10175 2019-01-28 15:00:00
## 130114 16 28 1 2019 -0.17025 2019-01-28 16:00:00
## 130115 17 28 1 2019 -0.11200 2019-01-28 17:00:00
## 130116 18 28 1 2019 0.01600 2019-01-28 18:00:00
## 130117 19 28 1 2019 0.08575 2019-01-28 19:00:00
## 130118 20 28 1 2019 0.04075 2019-01-28 20:00:00
## 130119 21 28 1 2019 -0.13775 2019-01-28 21:00:00
## 130120 22 28 1 2019 -0.34400 2019-01-28 22:00:00
## 130121 23 28 1 2019 -0.50175 2019-01-28 23:00:00
## 130122 24 28 1 2019 -0.53550 2019-01-29 00:00:00
## 130123 1 29 1 2019 -0.42250 2019-01-29 01:00:00
## 130124 2 29 1 2019 -0.28450 2019-01-29 02:00:00
## 130125 3 29 1 2019 -0.21750 2019-01-29 03:00:00
## 130126 4 29 1 2019 -0.21300 2019-01-29 04:00:00
## 130127 5 29 1 2019 -0.16625 2019-01-29 05:00:00
## 130128 6 29 1 2019 -0.05075 2019-01-29 06:00:00
## 130129 7 29 1 2019 0.11000 2019-01-29 07:00:00
## 130130 8 29 1 2019 0.18975 2019-01-29 08:00:00
## 130131 9 29 1 2019 0.14625 2019-01-29 09:00:00
## 130132 10 29 1 2019 -0.01050 2019-01-29 10:00:00
## 130133 11 29 1 2019 -0.13100 2019-01-29 11:00:00
## 130134 12 29 1 2019 -0.17250 2019-01-29 12:00:00
## 130135 13 29 1 2019 -0.11650 2019-01-29 13:00:00
## 130136 14 29 1 2019 -0.02775 2019-01-29 14:00:00
## 130137 15 29 1 2019 0.02025 2019-01-29 15:00:00
## 130138 16 29 1 2019 -0.02175 2019-01-29 16:00:00
## 130139 17 29 1 2019 -0.03475 2019-01-29 17:00:00
## 130140 18 29 1 2019 0.02325 2019-01-29 18:00:00
## 130141 19 29 1 2019 0.16275 2019-01-29 19:00:00
## 130142 20 29 1 2019 0.23725 2019-01-29 20:00:00
## 130143 21 29 1 2019 0.20250 2019-01-29 21:00:00
## 130144 22 29 1 2019 0.00950 2019-01-29 22:00:00
## 130145 23 29 1 2019 -0.23625 2019-01-29 23:00:00
## 130146 24 29 1 2019 -0.42050 2019-01-30 00:00:00
## 130147 1 30 1 2019 -0.50575 2019-01-30 01:00:00
## 130148 2 30 1 2019 -0.44575 2019-01-30 02:00:00
## 130149 3 30 1 2019 -0.33050 2019-01-30 03:00:00
## 130150 4 30 1 2019 -0.24950 2019-01-30 04:00:00
## 130151 5 30 1 2019 -0.19975 2019-01-30 05:00:00
## 130152 6 30 1 2019 -0.14225 2019-01-30 06:00:00
## 130153 7 30 1 2019 -0.04625 2019-01-30 07:00:00
## 130154 8 30 1 2019 0.14275 2019-01-30 08:00:00
## 130155 9 30 1 2019 0.16925 2019-01-30 09:00:00
## 130156 10 30 1 2019 0.08225 2019-01-30 10:00:00
## 130157 11 30 1 2019 -0.10475 2019-01-30 11:00:00
## 130158 12 30 1 2019 -0.25675 2019-01-30 12:00:00
## 130159 13 30 1 2019 -0.36450 2019-01-30 13:00:00
## 130160 14 30 1 2019 -0.37250 2019-01-30 14:00:00
## 130161 15 30 1 2019 -0.26725 2019-01-30 15:00:00
## 130162 16 30 1 2019 -0.18550 2019-01-30 16:00:00
## 130163 17 30 1 2019 -0.13875 2019-01-30 17:00:00
## 130164 18 30 1 2019 -0.12775 2019-01-30 18:00:00
## 130165 19 30 1 2019 -0.06625 2019-01-30 19:00:00
## 130166 20 30 1 2019 0.09925 2019-01-30 20:00:00
## 130167 21 30 1 2019 0.23975 2019-01-30 21:00:00
## 130168 22 30 1 2019 0.24050 2019-01-30 22:00:00
## 130169 23 30 1 2019 0.11875 2019-01-30 23:00:00
## 130170 24 30 1 2019 -0.07500 2019-01-31 00:00:00
## 130171 1 31 1 2019 -0.21325 2019-01-31 01:00:00
## 130172 2 31 1 2019 -0.26250 2019-01-31 02:00:00
## 130173 3 31 1 2019 -0.21625 2019-01-31 03:00:00
## 130174 4 31 1 2019 -0.08425 2019-01-31 04:00:00
## 130175 5 31 1 2019 0.02425 2019-01-31 05:00:00
## 130176 6 31 1 2019 0.06575 2019-01-31 06:00:00
## 130177 7 31 1 2019 0.10275 2019-01-31 07:00:00
## 130178 8 31 1 2019 0.22000 2019-01-31 08:00:00
## 130179 9 31 1 2019 0.39025 2019-01-31 09:00:00
## 130180 10 31 1 2019 0.46875 2019-01-31 10:00:00
## 130181 11 31 1 2019 0.47200 2019-01-31 11:00:00
## 130182 12 31 1 2019 0.30350 2019-01-31 12:00:00
## 130183 13 31 1 2019 0.18025 2019-01-31 13:00:00
## 130184 14 31 1 2019 0.03800 2019-01-31 14:00:00
## 130185 15 31 1 2019 0.08975 2019-01-31 15:00:00
## 130186 16 31 1 2019 0.14975 2019-01-31 16:00:00
## 130187 17 31 1 2019 0.25300 2019-01-31 17:00:00
## 130188 18 31 1 2019 0.29025 2019-01-31 18:00:00
## 130189 19 31 1 2019 0.24550 2019-01-31 19:00:00
## 130190 20 31 1 2019 0.25225 2019-01-31 20:00:00
## 130191 21 31 1 2019 0.39450 2019-01-31 21:00:00
## 130192 22 31 1 2019 0.50650 2019-01-31 22:00:00
## 130193 23 31 1 2019 0.51700 2019-01-31 23:00:00
## 130194 24 31 1 2019 0.31700 2019-02-01 00:00:00From this it is clear that it ends on line 130194, which is the 31st of January 2019.
min(which(ark.hour$year==2003&ark.hour$month==9))## [1] 145max(which(ark.hour$year==2019&ark.hour$month==1))## [1] 130194From the output cut the data from 145 and 130194
ark.month <- ark.hour[145:130194,]Create a new empty column for time
ark.month <- cbind(ark.month, NA)
names(ark.month) <- c("hour", "day", "month", "year", "h", "time")
ark.month$time <- as.POSIXlt(sprintf("%s/%s/%s %s", ark.month$year, ark.month$month, ark.month$day, ark.month$h), format="%Y/%m/%d %H",tz='UTC')ark.month <- aggregate(h~month+year,ark.month,mean)ark.month <- cbind(ark.month, NA)
names(ark.month) <- c("month", "year", "h", "time")Add a POSIXlt variable back in
ark.month$time <- as.POSIXlt(sprintf("%s/%s/15", ark.month$year, ark.month$month),tz="UTC")Plot the resultant data, add legends etc.
plot(ark.month$time,ark.month$h,col='blue',main = 'Monthly sea level',xlab = 'Years',ylab = 'Sea level [m]',type = 'o')
legend("topleft",
legend=c("Arklow tide gauge"),
col= c("blue"), lty=1,lwd=1, bty="n")
From the graph it is clear that there is distinct high and low peaks in tidal measurements.
Save Arklow monthly data as a CSV file
write.csv(ark.month, file = "arklowmonthlydata.csv")