# HG changeset patch # User siwaa # Date 1725024588 0 # Node ID 5fdc9d0c223bb9ae34f112ef25aa494694b92d4a # Parent 366449cab681b86f1d9e6991028134489dca8124 "planemo upload for repository https://forgemia.inra.fr/carboseq/record-projet-carboseq/-/tree/43e7f4ebfe4e552bf684936df0a050cb6def5aab/toos/CarboSeqSimulator commit 43e7f4ebfe4e552bf684936df0a050cb6def5aab-dirty" diff -r 366449cab681 -r 5fdc9d0c223b carboseqSimulator.xml --- a/carboseqSimulator.xml Thu Aug 29 14:58:19 2024 +0000 +++ b/carboseqSimulator.xml Fri Aug 30 13:29:48 2024 +0000 @@ -1,4 +1,4 @@ - + To run the CarboSeq simulator docker://registry.forgemia.inra.fr/csopra/csopralibs:siwaa diff -r 366449cab681 -r 5fdc9d0c223b carboseqUnsplit.R --- a/carboseqUnsplit.R Thu Aug 29 14:58:19 2024 +0000 +++ b/carboseqUnsplit.R Fri Aug 30 13:29:48 2024 +0000 @@ -1,13 +1,26 @@ -resU <- list(SOC=data.frame(), forcings=data.frame()) +resU <- data.frame() + +listFiles <- list.files('./inputs') + +# restoring the order first + +orderVec <- c() for (i in list.files('./inputs')) { load(file.path('.', 'inputs',i)) - resU$SOC <- rbind(resU$SOC, res$SOC); - resU$forcings <- rbind(resU$forcings, res$forcings) + orderVec <- c(orderVec, res$ID[1]); + rm(res) + } + +listFilesReOrder <- listFiles[order(orderVec)] + +for (i in listFilesReOrder) { + load(file.path('.', 'inputs',i)) + resU<- rbind(resU, res); rm(res) } res <-resU -save(res, file = 'csopraRes.Rdata') \ No newline at end of file +save(res, file = 'csopraRes.Rdata') \ No newline at end of file diff -r 366449cab681 -r 5fdc9d0c223b carboseqUnsplit.xml --- a/carboseqUnsplit.xml Thu Aug 29 14:58:19 2024 +0000 +++ b/carboseqUnsplit.xml Fri Aug 30 13:29:48 2024 +0000 @@ -1,4 +1,4 @@ - + To unsplit CarboSeq results docker://registry.forgemia.inra.fr/csopra/csopralibs:siwa diff -r 366449cab681 -r 5fdc9d0c223b run.r --- a/run.r Thu Aug 29 14:58:19 2024 +0000 +++ b/run.r Fri Aug 30 13:29:48 2024 +0000 @@ -118,8 +118,8 @@ preprocessorInputs = ppi) print("between load an sim") modeltoolbox::runSims(maxCores = NBCores) - # res <- modeltoolbox::getYearlySOCResults() - res <- modeltoolbox::readSocForcingsResults() + res <- modeltoolbox::getYearlySOCResults() + # res <- modeltoolbox::readSocForcingsResults() fileName <- "resRun.RData" filePath <- file.path(srcOut, fileName) save(res, file = filePath) @@ -148,8 +148,8 @@ modeltoolbox::loadCSVs(maxCores = NBCores, cinputCalculator = "cinest", preprocessorInputs = ppi) modeltoolbox::runSims(maxCores = NBCores) - # res <- modeltoolbox::getYearlySOCResults() - res <- modeltoolbox::readSocForcingsResults() + res <- modeltoolbox::getYearlySOCResults() + # res <- modeltoolbox::readSocForcingsResults() fileName <- "resRun.RData" filePath <- file.path(srcOut, fileName) save(res, file = filePath)