view redelacSticsInputGenerator.R @ 3:e344100f3354 draft

"planemo upload for repository https://forgemia.inra.fr/redelac commit ff64e62033c2485113100d60c38be73e59ef119c"
author siwaa
date Thu, 29 Jun 2023 10:50:58 +0000
parents c7d424481b03
children 4f8d87b9c246
line wrap: on
line source

library(SticsRFiles)
library(dplyr)

workspace <- paste0(getwd(), "/WS")
javastics_path <- getwd()

txt_path <- paste0(getwd(), "/WS/txt_files")
dir.create(txt_path,recursive = T)

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)) %>%
   select(!c(annee, date_min, date_max, Date_res_1, Date_res_2, Type_res_1, Type_res_2, Date_trav_1, Date_trav_2, Date_trav_3, Type_trav_1, Type_trav_2 ,Type_trav_3, plante_sem, date_sem, starts_with("nom_vari"), date_ferti_1, date_ferti_2, date_ferti_3, type_eng_1, type_eng_2, type_eng_3, date_rec, date_des, juleclair, nbinfloecl))
gen_tec_xml(param_df = tec_param, out_dir = workspace)
#   select(!c(annee, date_min, date_max, Date_res_1, Date_res_2, Type_res_1, Type_res_2, Date_trav_1, Date_trav_2, Date_trav_3, Type_trav_1, Type_trav_2 ,Type_trav_3, plante_sem, date_sem)

# 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)