Remove soilWaterInit

Author: Qianyu Li

models/sipnet/R/write.configs.SIPNET.R

@@ -557,7 +557,7 @@ write.config.SIPNET <- function(defaults, trait.values, settings, run.id, inputs
     }
     ## soilWater 
     if ("soilWater" %in% ic.names) {
-      param[which(param[, 1] == "soilWaterInit"), 2] <- IC$soilWater
+      param[which(param[, 1] == "soilWFracInit"), 2] <- IC$soilWater/param[which(param[, 1] == "soilWHC"), 2]
     }
     ## soilWFracInit fraction
     if ("soilWFrac" %in% ic.names) {

models/sipnet/inst/template.param

@@ -4,7 +4,6 @@ litterInit		280		0	130		1200		25
 soilInit		10000		1	3300		19000		3000.000000
 litterWFracInit		0.5		0	0.0		1.0		0.1
 soilWFracInit		0.5		1	0.0		1.0		0.1
-soilWaterInit			50		100	0		100		10
 snowInit		1.0		0	0.0		0.0		1.0
 microbeInit		0.50		0	0.020000	1.000000	0.001000
 fineRootFrac		0.200000	0	0.000000	1.000000	0.001000

modules/assim.sequential/inst/MultiSite-Exs/SDA/Bartlett.param

@@ -4,7 +4,6 @@ litterInit			200			0	130			1200				25
 soilInit			1600		0	1400	1700		1000
 litterWFracInit			0.5		1	0		1		0.1
 soilWFracInit			0.5		1	0		1		0.1
-soilWaterInit			50		100	0		100		10
 snowInit			0.0			0	0.0			0.0			1.0
 microbeInit		0.50		0	0.020000	1.000000	0.001000
 fineRootFrac		0.200000	0	0.000000	1.000000	0.001000

modules/assim.sequential/inst/MultiSite-Exs/SDA/Create_Multi_settings.R

@@ -122,6 +122,7 @@ template <- PEcAn.settings::Settings(list(
       #you could add more state variables here
       variable = structure(list(variable.name = "AbvGrndWood", unit = "MgC/ha", min_value = 0, max_value = 9999)),
       variable = structure(list(variable.name = "LAI", unit = "", min_value = 0, max_value = 9999)),
+      variable = structure(list(variable.name = "SoilMoist", unit = "cm", min_value = 0, max_value = 100)),
       variable = structure(list(variable.name = "SoilMoistFrac", unit = "", min_value = 0, max_value = 1)),#soilWFracInit
       variable = structure(list(variable.name = "TotSoilCarb", unit = "kg/m^2", min_value = 0, max_value = 9999))
     )),