modification pour exclure non sol maisnon finalisé

dirsol/0_InitialisationsCLC.R

@@ -89,7 +89,7 @@ excludeWRB = 0
 excludeCLC = c(1,2,3,4,5,6,7,8,9,
                30,31,34,35,37,38,39,40,
                41,42,43,44,48,49,50
-)
+               )
 
 
 titreCarte = "Simulation du réseau de mesures de la SML avec CLC"

dirsol/0_InitialisationsLULCF.R

@@ -63,6 +63,9 @@ SIsample = 1000000
 
 excludeWRB = 0
 
+# the list of clc codes to exclude
+excludeCLC = 4:6
+
 titreCarte = "Simulation du réseau de mesures de la SML"
 
 ## Extraction pixels des stacks et clhs  --------

dirsol/Simul_DIRSOL_avec_tout_RMQS_LULUCF.R

@@ -34,7 +34,8 @@ sample3 <- sample3 %>%
          FR_WRB = WRBLEV1)
  
 #remove the non soil
-sample3 <- filter(sample3, !FR_WRB %in% excludeWRB )
+sample3 <- filter(sample3, !FR_WRB %in% excludeWRB ) %>%
+  filter(!FR_landuse %in% excludeWRB)
 
 
 # saveRDS( sample3, "~/smlsamplingeu_louis/data/octobre_avec_rmqs_sample3_lulucf.rds")
@@ -62,7 +63,7 @@ ListVarScenario <- c("FR_landuse",
 
 domains <- 1e6 * stackb[["FR_landuse"]] +
   1e4 * stackb[["code_supra"]]  +
-  1e2 *stackb[["WRBLEV1"]] +
+  1e2 * stackb[["WRBLEV1"]] +
   stackb[["climate_metzeger"]]