R Under development (unstable) (2024-01-28 r85838 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ## Test soil properties in APSIM-X > > require(apsimx) Loading required package: apsimx APSIM(X) not found. If APSIM(X) is installed in an alternative location, set paths manually using 'apsimx_options' or 'apsim_options'. You can still try as the package will look into the registry (under Windows) > require(ggplot2) Loading required package: ggplot2 > > apsimx_options(warn.versions = FALSE) > > run.soil <- get(".run.local.tests", envir = apsimx.options) > > tmp.dir <- tempdir() > > if(run.soil){ + + ex.dir <- auto_detect_apsimx_examples() + + wheat <- apsimx("Wheat.apsimx", src.dir = ex.dir) + + inspect_apsimx("Wheat.apsimx", src.dir = ex.dir, + node = "Soil", soil.child = "Physical") + + duls <- c(0.52, 0.49, 0.488, 0.48, 0.47, 0.457, 0.452) * 1.05 + + edit_apsimx("Wheat.apsimx", src.dir = ex.dir, + wrt.dir = tmp.dir, node = "Soil", + soil.child = "Physical", parm = "DUL", + value = duls) + + inspect_apsimx("Wheat-edited.apsimx", src.dir = tmp.dir, + node = "Soil", soil.child = "Physical") + + wheat2 <- apsimx("Wheat-edited.apsimx", src.dir = tmp.dir) + + ggplot() + + geom_line(data = wheat, aes(x = Date, y = Wheat.AboveGround.Wt, color = "Original DUL")) + + geom_line(data = wheat2, aes(x = Date, y = Wheat.AboveGround.Wt, color = "DUL times 1.05")) + + } > > if(run.soil){ + + ## For now + extd.dir <- system.file("extdata", package = "apsimx") + + inspect_apsimx("Wheat.apsimx", src.dir = extd.dir, + node = "Soil", soil.child = "Physical") + + inspect_apsimx("Maize.apsimx", src.dir = extd.dir, + node = "Soil", soil.child = "Physical") + + inspect_apsimx("Wheat_soil.apsimx", src.dir = extd.dir, + node = "Soil", soil.child = "Physical") + + ## CanI now just edit the Crop specific soil physical? + edit_apsimx("Wheat_soil.apsimx", src.dir = extd.dir, + wrt.dir = tmp.dir, + node = "Soil", soil.child = "Physical", + parm = "Maize KL", value = rep(0.07, 7)) + + edit_apsimx("Wheat_soil-edited.apsimx", + src.dir = tmp.dir, + overwrite = TRUE, + node = "Soil", soil.child = "Physical", + parm = "Wheat KL", value = rep(0.05, 7)) + + inspect_apsimx("Wheat_soil-edited.apsimx", src.dir = tmp.dir, + node = "Soil", soil.child = "Physical") + + } > > proc.time() user system elapsed 1.14 0.18 1.31