Difference between revisions of "Electricity use in a house with geothermal heating"

# This code only works from a computer within THL network.
dat1 <- read.csv("//cesium/yhteiset/YMAL/Projects/Urgenche/WP4 Buildings/report2013b.csv")
dat1$Kulutus <- as.numeric(gsub(",", ".", dat1$Kulutus))
dat1 <- dat1[!is.na(dat1$Kulutus) , ]
dat1$Aika <- as.POSIXct("2013-01-19 0:00:00 EET") + (0:(nrow(dat1) - 1)) * 3600
dat1 <- dat1[c(4,2)]

write.csv(dat1 , file = "//cesium/yhteiset/YMAL/Projects/Urgenche/WP4 Buildings/report2013c.csv", row.names = FALSE)

# The function below does not work from own computer. The same was done using OpasnetBaseImport.
#opbase.upload(dat1, ident = "Op_en6315", name = "Electricity use in a house with geothermal heating", subset = "2013", 
#    obj_type = "variable", act_type = "replace", language = "eng", 
#    unit = "kWh", who = "Jouni", rescol = "Kulutus"
#)

library(fmi)

request <- FMIWFSRequest$new(apiKey=apiKey)
#request$setParameters(request="getFeature",
#                      storedquery_id="fmi::observations::weather::daily::timevaluepair",
#                      starttime="2014-01-01T00:00:00Z",
#                      endtime="2014-01-01T00:00:00Z",
#                      place="savilahti,Kuopio",
#                      parameters="rrday,snow,tday,tmin,tmax")

request$setParameters(request="getFeature",
                      storedquery_id="fmi::observations::weather::daily::timevaluepair",
                      starttime="2014-01-01T00:00:00Z",
                      endtime="2014-01-01T00:00:00Z",
                      bbox="19.09,59.3,31.59,70.13",
                      parameters="rrday,snow,tday,tmin,tmax")

client <- FMIWFSClient$new(request=request)
client
layers <- client$listLayers()
layers
#response <- client$getLayer(layer=layers[1], parameters=list(splitListFields=TRUE))
#response

request <- FMIWFSRequest$new(apiKey=apiKey)
client <- FMIWFSClient$new(request=request)
response <- client$getDailyWeather(startDateTime="2014-01-01", endDateTime="2014-01-01", bbox=getFinlandBBox(), parameters=list(explodeCollections=TRUE))
reponse