diff redelacSticsInputGenerator.R @ 0:c7d424481b03 draft

"planemo upload for repository https://forgemia.inra.fr/redelac"
author siwaa
date Wed, 07 Jun 2023 06:08:27 +0000
parents
children e344100f3354
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/redelacSticsInputGenerator.R	Wed Jun 07 06:08:27 2023 +0000
@@ -0,0 +1,52 @@
+library(SticsRFiles)
+library(dplyr)
+
+workspace <- paste0(getwd(), "/WS")
+javastics_path <- getwd()
+
+txt_path <- paste0(getwd(), "/WS/txt_files")
+dir.create(txt_path,recursive = T)
+
+gen_usms_xml2txt(javastics = javastics_path,
+                 workspace = workspace,
+		 out_dir = txt_path,
+		 check = FALSE)
+
+USMsFile <- "USMs.csv"
+TecFile <- "Tec.csv"
+IniFile <- "Ini.csv"
+StationFile <- "Station.csv"
+
+# Generate STICS xml files------------------------------------------------------------------------------
+
+# USMs file
+usms_param <- read_params_table(file.path(workspace, USMsFile)) %>%
+  select(usm_name,datedebut,datefin,finit,nomsol,fstation,fclim1,fclim2,culturean,nbplantes,codesimul,fplt_1,ftec_1,flai_1,fplt_2,ftec_2,flai_2)
+gen_usms_xml(file = file.path(workspace, "usms.xml"), param_df = usms_param)
+
+# Tec files
+tec_param <- read_params_table(file.path(workspace, TecFile))
+gen_tec_xml(param_table = tec_param, out_path = workspace)
+
+# Ini files
+ini_param <- read_params_table(file.path(workspace, IniFile))
+gen_ini_xml(param_df = ini_param, out_dir = workspace)
+
+# Sta files
+sta_param <- read_params_table(file.path(workspace, StationFile))
+gen_sta_xml(param_df = sta_param, out_dir = workspace)
+
+# Generate STICS txt files------------------------------------------------------------------------------
+# List of USMs
+list_usms <- SticsRFiles::get_usms_list(file = file.path(workspace,"usms.xml"))
+
+### Generate txt files into individual folders
+
+# generate text files
+gen_usms_xml2txt(javastics = javastics_path, 
+                 workspace = workspace, 
+                 out_dir = txt_path,
+                 verbose = TRUE,
+                 usm = list_usms,
+                 check=FALSE) #je mets l'option a FALSE car sinon j'ai soit disant des fichiers climatiques manquants du type 14050.198,14050.199,14050.2 (pas de lecture des z?ros)
+    
\ No newline at end of file