data preperation for rcp0 V7

Main_cPouta.R

@@ -23,10 +23,10 @@ path <- paste0(getwd())
 on_cPouta <- FALSE
 
 ### Do the SQL queries for the DB (only needed first time!)
-Process_SQL_Query <- FALSE
+Process_SQL_Query <- TRUE
 
 ### Export the defined columns from DB as CSV
-Process_CSV_export_from_DB <- FALSE
+Process_CSV_export_from_DB <- TRUE
 
 
 
@@ -34,7 +34,7 @@ Process_CSV_export_from_DB <- FALSE
 # Define name of genereated Optimization CSV file
 # ------------------------------------------------
 
-outp_name <- "rslt_RCP0_FIN_V4"
+outp_name <- "rslt_RCP0_FIN_V7"
 
 
 
@@ -59,7 +59,7 @@ if(on_cPouta == TRUE) {
 } else {
   # for Windows
   inputFolder  <- paste(path, "input_data/FIN", sep = "/")
-  outputFolder <- paste(path, "output/RCP0_FIN",     sep = "/")
+  outputFolder <- paste(path, "output/MF_FIN",     sep = "/")
 
 }
 

MultiForest/check_opt_data.Rmd

@@ -5,18 +5,23 @@ date: "7/28/2020"
 output: html_document
 ---
 
+
+# RCP0 whole Finland - version 7
+### (rslt_RCP0_FIN_V7.csv)
+
 ```{r setup, include=FALSE}
 knitr::opts_chunk$set(echo = TRUE)
 
 
 ### Running on cPouta or WIN local
-on_cPouta <- FALSE
+on_cPouta <- TRUE
 
 
 ### Set path to input files
 
 if(on_cPouta == TRUE) {
   # for pouta
+  datapath <- "/media/volume/outp_rcp0_test/"
   path <- "/home/ubuntu/workspace/MultiForest2/" 
   # path <- paste0(getwd(),"/")
   
@@ -43,11 +48,11 @@ library(egg)
 
 ```{r , echo=FALSE}
 
-# poutaxy <- read.csv(paste0(path, "rslt_RCP0_FIN_V3_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
-# pouta <- read.csv(paste0(path, "rslt_RCP0_FIN_V3.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
+poutaxy <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V7_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
+pouta <- read.csv(paste0(datapath, "rslt_RCP0_FIN_V7.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
 
-poutaxy <- read.csv(paste0(path, "rslt_test_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
-pouta <- read.csv(paste0(path, "rslt_test.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
+#poutaxy <- read.csv(paste0(datapath, "rslt_test_XY.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
+#pouta <- read.csv(paste0(datapath, "rslt_test.csv"), sep = ";", header = TRUE, stringsAsFactors = FALSE)
 
 ```
 
@@ -715,8 +720,6 @@ plotfkt("N_where_D_gt_40")
 
 ## Indicator: prc_V_deciduous
 
-## STILL WITH VALUES ABOVE 100% !!!
-
 ```{r , message=FALSE}
 
 plotfkt("prc_V_deciduous")
@@ -872,7 +875,18 @@ plotfkt_harvest("Harvested_V_pulp_under_bark")
 
 
 
+# Save as data set for the Optimization framework
 
+```{r , echo=TRUE}
+# For the Optimization (Jupyter Notebook) the data can be further ZIP compressed.
+
+# csv for optimization
+write.table(rslt.final, paste0(datapath, "rslt_RCP0_FIN_V7_no_outlier.csv" ),
+            sep = ";",
+            row.names = F,
+            col.names = TRUE)
+
+```