include sim variants for CanESM2 model, cPouta

Main.R

@@ -20,23 +20,38 @@ regime <- read.csv(paste0(path, "params/regimes.csv"), sep = ",", stringsAsFacto
 
 
 ### Simulation variant that will be load
-# specify which one:  "CC45"        climate change with RCP scenraio 4.5 
-#                     "CC45_SA"     RCP 4.5 and set aside without deadwood extraction (no other management regimes!)
-#                     "CC85"        climate change with RCP scenario 8.5 
-#                     "CC85_SA"     RCP 8.5 and set aside without deadwood extraction (no other management regimes!)
-#                     "without"     no climate change
-#                     "without_SA"  no climate change and set aside without deadwood extraction (no other management regimes!)
-#
-# variants simulated with SIMO Linux version on cPouta virtual machine are indicated by "p" -> output comparison.
-#
-#                     "CC45_p"      RCP scenario 4.5 
-#                     "CC85_p"      RCP scenario 8.5
-#                     "CC26_p"      NEW: RCP scenario 2.6, was not simulated before
-#                     "CC0_p        NEW: RCP scenario 0, represent NO CLIMATE CHANGE, NEEDS TO BE validated!!
+### The corresponding db files need to be stored in folder "input_data"
+### specify which one should be read:
+
+
+# FOR MSC STUDENTS - total watershed have been already simulated:
+# variant simulated with WIN version of NON-climate sensitive SIMO version
+#                    "without"     no climate change
+#                    "without_SA"  no climate change and set aside without deadwood extraction (no other management regimes!)
+# 
+# variants simulated with WIN version of climate sensitive SIMO version, climate model HadGEM2
+#                    "CC45"        climate change with RCP scenraio 4.5 
+#                    "CC45_SA"     RCP 4.5 and set aside without deadwood extraction (no other management regimes!)
+#                    "CC85"        climate change with RCP scenario 8.5 
+#                    "CC85_SA"     RCP 8.5 and set aside without deadwood extraction (no other management regimes!)
+
+
+# FOR COMPARISON OF CLIMATE MODEL AND SIMO VERSIONS - ONLY 5 stands per watershed so far simulated
+# variants simulated with SIMO Linux version on cPouta, indicated by "p"
+# with climate model HadGEM2
+#                   "CC45_p"      RCP scenario 4.5
+#                   "CC85_p"      RCP scenario 8.5
+#                   "CC26_p"      NEW: RCP scenario 2.6
+#                   "CC0_p        NEW: RCP scenario 0, represents NO CLIMATE CHANGE ?? 
 #
+# with climate model CanESM2 (like in FutureBioEcon), indicated further by "canesm"
+#                   "CC45_p_canesm"      RCP sceanrio 4.5
+#                   "CC85_p_canesm"      RCP scenario 8.5
+#                   "CC26_p_canesm"      NEW: RCP scenario 2.6
+#                   "CC0_p_canesm"       NEW: RCP scenario 0
 
 
-sim_variant <- "CC0_p"  
+sim_variant <- "CC0_p_canesm"  
 
 
 ### Define the names of the databases (SIMO-output for 10 watersheds) that will be imported
@@ -67,7 +82,7 @@ first_load = FALSE
 if(first_load == TRUE){
 
   # If one of the follwing simulation variants is read ...
-  if(sim_variant %in% c("CC45", "CC85", "without", "CC45_p", "CC85_p", "CC26_p", "CC0_p")) {
+  if(sim_variant %in% c("CC45", "CC85", "without", "CC45_p", "CC85_p", "CC26_p", "CC0_p", "CC45_p_canesm", "CC85_p_canesm", "CC26_p_canesm", "CC0_p_canesm")) {
     
     # Run the script with the SQL query for all management regimes
     source(paste0(path, "structure_SIMO_rslDB_FBE.R"))  

loadDB.R

@@ -41,7 +41,7 @@ for (name in db_names){
   
   
   ### Rename set aside scenario if it considers deadwood extraction
-  if(sim_variant %in% c("CC45", "CC85", "without", "CC45_p", "CC85_p", "CC26_p", "CC0_p")) {
+  if(sim_variant %in% c("CC45", "CC85", "without", "CC45_p", "CC85_p", "CC26_p", "CC0_p", "CC45_p_canesm", "CC85_p_canesm", "CC26_p_canesm", "CC0_p_canesm")) {
     
    rsl <- rsl %>% mutate(regime = ifelse(regime %in% "SA", "SA_DWextract", regime))
     
@@ -76,7 +76,7 @@ rslt <- rslt %>%
 
 
 ### Rename set aside scenario if it considers deadwood extraction
-if(sim_variant %in% c("CC45", "CC85", "without", "CC45_p", "CC85_p", "CC26_p", "CC0_p")) {
+if(sim_variant %in% c("CC45", "CC85", "without", "CC45_p", "CC85_p", "CC26_p", "CC0_p", "CC45_p_canesm", "CC85_p_canesm", "CC26_p_canesm", "CC0_p_canesm")) {
   
   rslt <- rslt %>% mutate(regime = ifelse(regime %in% "SA", "SA_DWextract", regime))