library(SEF)
inventory <- read.table("C:/Users/elinl/Documents/UniBe/Instrumentelle data/ISTI/Inventorys/Inventory_Alps2.txt", colClasses="character")
colnames(inventory) <- c("Other_ID", "City", "Modern_Country", "Lat.degN", "Lon.degE", "Station_Elevation.m", "Start_Year", "End_Year" )
allfiles = list.files("C:/Users/elinl/Documents/Unibe/Instrumentelle data/ISTI/Alps", full.names=TRUE)
for (currentfile in inventory$Other_ID) {
currentfile innholder nå navnet på filen som skal leses inn
her fyller du på med kommandoene som skal kjøres for hver fil
index <- grep(currentfile, allfiles)
#read in data -the x
x <- read.table(allfiles[index], header = FALSE, fill=TRUE, sep="",stringsAsFactors=FALSE, na.strings="-9999")
colnames(x) <- c("City", "Lat", "Lon", "Station Elevation", "Date", "value_tmin","value_tmax","value_tavg","other")
year<- as.numeric(substr(x$Date,1,4))
months<- as.numeric(substr(x$Date,5,6))
value_tavg <- x$value_tavg/100
#inventory <- read.csv("Inventory_Canada.txt", sep=";", stringsAsFactors=F)
index <- which(inventory$Other_ID==currentfile)
lat <- inventory$Lat.degN[index]
lon <- inventory$Lon.degE[index]
alt <- inventory$Station_Elevation.m[index]
station_name <- inventory$City[index]
df <- data.frame(y=year, m=months, d="NA",
hh=rep("",nrow(x)), mm=rep("",nrow(x)), value_tavg, stringsAsFactors=FALSE)
write_sef(df, variable="ta", cod=inventory$Other_ID[index], nam=station_name, lat=lat,
lon=lon, alt=alt, sou="ISTI", units="C", stat="mean", period="month")
first_year <- min(df$y)
last_year <- max(df$y)
file.rename(list.files(pattern=as.character(inventory$Other_ID[index])),
paste("ISTI",inventory$Other_ID[index],station_name,first_year,last_year,"ta_monthly.tsv",
sep="_"))
}