# Electricity use in a house with geothermal heating

## Question

What is the energy (i.e. electricity) use of geothermal heating in a detached house?

## Rationale

The answer is based on a single example house located near the Savilahti weather station, Kuopio, Finland.

### Unit

kWh (or kW, as the data is given in hourly slots)

### Calculations

 ```# 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" #) ```

 apiKey:```library(OpasnetUtils) library(fmi) library(rwfs) library(rgdal) request <- FMIWFSRequest\$new(apiKey=apiKey) if(FALSE) { #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 } request\$setParameters(request="getFeature", storedquery_id="fmi::observations::weather::daily::timevaluepair", parameters="rrday,snow,tday,tmin,tmax") client <- FMIWFSClient\$new(request=request) response <- client\$getDailyWeather(startDateTime="2014-01-01T00:00:00Z", endDateTime="2014-01-01T00:00:00Z", fmisid=101887) dat <- as.data.frame(response) oprint(dat) ```