Skip to content

Commit 1a7481b

Browse files
authored
SIP111 Zero out envi values when relevant modes are off (#143)
* Zeros out microbeC and litter when modes are off * Updates after zeroing microbeC and litter
1 parent 3b9bd9f commit 1a7481b

File tree

5 files changed

+22796
-22786
lines changed

5 files changed

+22796
-22786
lines changed

src/sipnet/sipnet.c

Lines changed: 15 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -1898,17 +1898,27 @@ void setupModel(void) {
18981898
(1 - params.coarseRootFrac - params.fineRootFrac) * params.plantWoodInit;
18991899
envi.plantLeafC = params.laiInit * params.leafCSpWt;
19001900

1901-
envi.litter = params.litterInit;
1901+
if (ctx.litterPool) {
1902+
envi.litter = params.litterInit;
1903+
} else {
1904+
// Don't set a value if the litter pool is off
1905+
envi.litter = 0.0;
1906+
}
19021907
envi.soil = params.soilInit;
19031908

19041909
// change from per hour to per day rate
19051910
params.maxIngestionRate = params.maxIngestionRate * 24;
19061911

1907-
envi.microbeC = params.microbeInit * params.soilInit / 1000; // convert to gC
1908-
// m-2
1912+
if (ctx.microbes) {
1913+
// convert to gC m-2
1914+
envi.microbeC = params.microbeInit * params.soilInit / 1000;
1915+
} else {
1916+
// Don't set a value if microbes is off
1917+
envi.microbeC = 0.0;
1918+
}
19091919

1910-
params.totNitrogen = params.totNitrogen * params.soilInit; // convert to gC
1911-
// m-2
1920+
// convert to gC m-2
1921+
params.totNitrogen = params.totNitrogen * params.soilInit;
19121922

19131923
params.fineRootTurnoverRate /= 365.0;
19141924
params.coarseRootTurnoverRate /= 365.0;

0 commit comments

Comments
 (0)