Good night.
You should be able to start R & RStudio and copy and paste the play function (see below) into R. Try pasting it into the console window of RStudio hand hitting return.
Then do
play
If it is in your work area you should see it print out.
Another approach would be to save the function in an R file, call it play.r in your working directory and issue the command
source("play.r"
This should load and execute the function so that it is ready to use.
John
*Play Function
play <- function (Data, outpath, variable, cod, nam = "", lat = "", lon = "",
alt = "", sou = "", link = "", units, stat, metaHead = "",
meta = "", period = "", time_offset = 0, note = "", keep_na = FALSE,
outfile = NA)
{
for (i in 1:ncol(Data)) Data[, i] <- as.character(Data[,
i])
header <- array(dim = c(12, 2), data = "")
header[1, ] <- c("SEF", "1.0.0")
header[2, ] <- c("ID", trimws(as.character(cod)))
header[3, ] <- c("Name", trimws(as.character(nam)))
header[4, ] <- c("Lat", trimws(as.character(lat)))
header[5, ] <- c("Lon", trimws(as.character(lon)))
header[6, ] <- c("Alt", trimws(as.character(alt)))
header[7, ] <- c("Source", trimws(as.character(sou)))
header[8, ] <- c("Link", trimws(as.character(link)))
header[9, ] <- c("Vbl", trimws(as.character(variable)))
header[10, ] <- c("Stat", trimws(as.character(stat)))
header[11, ] <- c("Units", trimws(as.character(units)))
header[12, ] <- c("Meta", trimws(as.character(metaHead)))
if (stat == "point" & !all(as.character(period) == "0")) {
period <- "0"
warning("Period forced to 0 because of 'stat'")
}
if (!all(time_offset == 0) & !all(is.na(as.integer(Data[,
4]) + as.integer(Data[, 5])))) {
times <- ISOdate(Data[, 1], Data[, 2], Data[, 3], Data[,
4], Data[, 5])
times <- times - time_offset * 3600
Data[which(!is.na(times)), 1] <- as.integer(substr(times[which(!is.na(times))],
1, 4))
Data[which(!is.na(times)), 2] <- as.integer(substr(times[which(!is.na(times))],
6, 7))
Data[which(!is.na(times)), 3] <- as.integer(substr(times[which(!is.na(times))],
9, 10))
Data[which(!is.na(times)), 4] <- as.integer(substr(times[which(!is.na(times))],
12, 13))
Data[which(!is.na(times)), 5] <- as.integer(substr(times[which(!is.na(times))],
15, 16))
}
DataNew <- data.frame(Year = Data[, 1], Month = Data[, 2],
Day = Data[, 3], Hour = Data[, 4], Minute = Data[, 5],
Period = as.character(period), Value = Data[, 6], Meta = as.character(meta),
stringsAsFactors = FALSE)
if (!keep_na)
DataNew <- DataNew[which(!is.na(DataNew$Value)), ]
if (substr(outpath, nchar(outpath), nchar(outpath)) != "/") {
outpath <- paste0(outpath, "/")
}
if (is.na(outfile)) {
j <- 3
if (is.na(as.integer(DataNew[1, 3])))
j <- 2
if (is.na(as.integer(DataNew[1, 2])))
j <- 1
datemin <- paste(formatC(unlist(as.integer(DataNew[1,
1:j])), width = 2, flag = 0), collapse = "")
datemax <- paste(formatC(unlist(as.integer(DataNew[dim(DataNew)[1],
1:j])), width = 2, flag = 0), collapse = "")
dates <- paste(datemin, datemax, sep = "-")
filename <- paste(sou, cod, nam, dates, variable, sep = "_")
if (sou %in% c(NA, ""))
filename <- sub("_", "", filename)
if (note != "") {
note <- paste0("_", gsub(" ", "_", note))
}
filename <- gsub(" ", "", filename)
filename <- paste0(outpath, filename, note, ".tsv")
}
else {
filename <- paste0(outpath, outfile)
if (substr(filename, nchar(filename) - 3, nchar(filename)) !=
".tsv") {
filename <- paste0(filename, ".tsv")
}
}
write.table(header, file = filename, quote = FALSE, row.names = FALSE,
col.names = FALSE, sep = "\t", dec = ".", fileEncoding = "UTF-8")
write.table(t(names(DataNew)), file = filename, quote = FALSE,
row.names = FALSE, col.names = FALSE, sep = "\t", fileEncoding = "UTF-8",
append = TRUE)
write.table(DataNew, file = filename, quote = FALSE, row.names = FALSE,
col.names = FALSE, sep = "\t", dec = ".", fileEncoding = "UTF-8",
append = TRUE)
message(paste("Data written to file", filename))
}
