Mercurial > repos > siwaa > redelac_stics_t

comparison redelacSticsTool.R @ 23:3b44cb0a1e76 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://forgemia.inra.fr/redelac commit 4cad1cb9719e5d2ca13cf8a5c8b98b2409303b59"
author siwaa
date Wed, 19 Aug 2026 14:35:50 +0000
parents 15317ef4816b
children
comparison
equal deleted inserted replaced
22:15317ef4816b 23:3b44cb0a1e76
249 force = TRUE, 249 force = TRUE,
250 successive = list_successive_usms 250 successive = list_successive_usms
251 ) 251 )
252 252
253 list_usms <- vec 253 list_usms <- vec
254 254
255 cat(paste0("--> succession : ", successive_usms_tab_split[[j]][i, 1], " launched\n")) 255 cat(paste0("--> succession : ", successive_usms_tab_split[[j]][i, 1], " launched\n"))
256 256
257 return(stics_wrapper(model_options = sim_options, situation = list_usms)) 257 return(stics_wrapper(model_options = sim_options, situation = list_usms))
258 258
259 } 259 }
265 265
266 # 266 #
267 # Output part 267 # Output part
268 # 268 #
269 269
270 beforeOutputsTime <- Sys.time()
271
272 simPCPlist <- list()
273
274 for (j in 1:length(casesUpcSol)) { 270 for (j in 1:length(casesUpcSol)) {
275
276 successive_usms_tab_split[[j]]$Treatment <- NULL 271 successive_usms_tab_split[[j]]$Treatment <- NULL
277 successive_usms_tab_split[[j]]$Number_USM <- NULL 272 successive_usms_tab_split[[j]]$Number_USM <- NULL
273 }
274
275 beforeOutputsTime <- Sys.time()
276
277 foreach(
278 j = 1:length(casesUpcSol),
279 .options.future = list(packages =c("foreach", "tidyr", "arrow"))
280 ) %dordopar% {
281
278 splited <- strsplit(casesUpcSol[j], split = "-")[[1]] 282 splited <- strsplit(casesUpcSol[j], split = "-")[[1]]
279 283
280 upc <- splited[1] 284 upc <- splited[1]
281 285
282 sol <- splited[2] 286 sol <- splited[2]
286 if (i == 1) { 290 if (i == 1) {
287 simPCP <- result[[j]][[i]] 291 simPCP <- result[[j]][[i]]
288 } else { 292 } else {
289 simPCP$sim_list <- c(simPCP$sim_list, result[[j]][[i]]$sim_list) 293 simPCP$sim_list <- c(simPCP$sim_list, result[[j]][[i]]$sim_list)
290 } 294 }
291 } 295
292 296 }
293 simPCPlist[[j]] <- simPCP
294 }
295
296 rm(casesUpcSol)
297
298 foreach(
299 j = 1:length(simPCPlist),
300 .options.future = list(packages =c("foreach", "tidyr", "arrow"))
301 ) %dordopar% {
302 297
303 #agregation des sorties dans un dataframe 298 #agregation des sorties dans un dataframe
304 sim_df <- 299 sim_df <-
305 rbindlist(simPCPlist[[j]]$sim_list, use.names = TRUE, idcol = "usm") 300 rbindlist(simPCP$sim_list, use.names = TRUE, idcol = "usm")
301
306 302
307 #ajout des colonnes explicatives + de l'horizon temporel 303 #ajout des colonnes explicatives + de l'horizon temporel
308 sim_df <- sim_df %>% 304 sim_df <- sim_df %>%
309 separate( 305 separate(
310 col = usm, 306 col = usm,
319 "culture" 315 "culture"
320 ), 316 ),
321 sep = "-", 317 sep = "-",
322 remove = FALSE 318 remove = FALSE
323 ) 319 )
320
321 cat(paste0("--> rotation upc sol : ", " ", rotation, " ", upc, " " , sol, " to be gathered\n"))
324 322
325 #sauvegarde sous la forme d'un fichier rds 323 #sauvegarde sous la forme d'un fichier rds
326 324
327 if (baseOut) { 325 if (baseOut) {
328 326