Remove V_Aspen - counted twice next to V_popolus

Main_cPouta.R

@@ -158,7 +158,6 @@ columns <-  paste0("id,
                    V,
                    V_penbirch,
                    V_pubbirch,
-                   V_aspen,
                    V_populus,
                    V_Alnus_incana,
                    V_Alnus_glutinosa,

prepare_opt_data_KE.R

@@ -214,7 +214,6 @@ rslt <- rslt %>%
   # sum up the volumes of deciduous
   mutate(V_deciduous = V_penbirch
          + V_pubbirch
-         + V_aspen
          + V_populus
          + V_Alnus_incana
          + V_Alnus_glutinosa
@@ -227,7 +226,7 @@ rslt <- rslt %>%
 mutate(prc_V_deciduous = ifelse(prc_V_deciduous %in% "NaN", 0, prc_V_deciduous))  # %>%
 
 # # remove unnecessary columns
-# select(!c("V_deciduous", "V_penbirch","V_pubbirch","V_aspen","V_populus","V_Alnus_incana","V_Alnus_glutinosa","V_o_decidious"))
+# select(!c("V_deciduous", "V_penbirch","V_pubbirch","V_populus","V_Alnus_incana","V_Alnus_glutinosa","V_o_decidious"))
 
 
 print("*** %deciduous trees calculated ***")
@@ -308,7 +307,6 @@ rslt <- rslt %>%
          V_deciduous =  round( V_deciduous , digits = 2), 
          V_penbirch =  round( V_penbirch , digits = 2), 
          V_pubbirch =  round( V_pubbirch , digits = 2), 
-         V_aspen =  round( V_aspen , digits = 2), 
          V_populus =  round( V_populus , digits = 2), 
          V_Alnus_incana =  round( V_Alnus_incana , digits = 2), 
          V_Alnus_glutinosa =  round( V_Alnus_glutinosa , digits = 2), 
@@ -557,9 +555,9 @@ unique(rslt$protection)
 
 
 print("*** status of protection added ***")
-cat("Number of stands status     strict:" , length(unique(rslt[rslt$protection %in% "strict",]$unique_id)))
-cat("Number of stands status  landscape:" , length(unique(rslt[rslt$protection %in% "landscape",]$unique_id)))
-cat("Number of stands status commercial:" , length(unique(rslt[rslt$protection %in% "commercial",]$unique_id)))
+print(paste("Number of stands status     strict:" , length(unique(rslt[rslt$protection %in% "strict",]$unique_id))))
+print(paste("Number of stands status  landscape:" , length(unique(rslt[rslt$protection %in% "landscape",]$unique_id))))
+print(paste("Number of stands status commercial:" , length(unique(rslt[rslt$protection %in% "commercial",]$unique_id))))
 
 
 
@@ -614,7 +612,7 @@ rslt <- rslt %>%
 rslt.xy_no <- rslt  %>%
   select(id, year, regime, X, Y, standid, ES_GAME, CARBON_STORAGE, scenario, region,
          # used for calculating % share broadleave
-         V_deciduous, V_penbirch, V_pubbirch, V_aspen, V_populus, V_Alnus_incana, V_Alnus_glutinosa, V_o_decidious,
+         V_deciduous, V_penbirch, V_pubbirch, V_populus, V_Alnus_incana, V_Alnus_glutinosa, V_o_decidious,
          # used for calculating meand HSI27
          HSI_RL_S1, HSI_RL_S2, HSI_RL_S3, HSI_RL_S4, HSI_RL_S5, HSI_RL_S6, HSI_RL_S7, HSI_RL_S8, HSI_RL_S9, HSI_RL_S10, HSI_RL_S11,
          HSI_RL_S12, HSI_RL_S13, HSI_RL_S14, HSI_RL_S15, HSI_RL_S16, HSI_RL_S17, HSI_RL_S18, HSI_RL_S19, HSI_RL_S20, HSI_RL_S21, HSI_RL_S22,
@@ -645,7 +643,7 @@ rslt.final <- rslt %>%
 
 
 print("*** rslt.final created ***")
-cat("Output contains all 39445 stand:", 39445 == length(unique(rslt.final$id)))
+print(paste("Output contains all 39445 stands:", 39445 == length(unique(rslt.final$id))))