Right now a step in runMemWithCohorts is to build the block of cohorts to start with before it runs the annual time steps of the sea-level rise scenario. It uses the initial conditions and the runToEquilibrium function to generate a set of cohorts given the assumption of current day equilibrium conditions.
After talking to @tilbud, it seems like we need to cover 2 new sets of potential starting points for the dynamic step before runMemWithCohorts resembles the latest versions Jim M has in the excel macros.
Supertidal peat: If the starting elevation is above the highest tide line, but below the maximum growing elevation of the plants, the sea-level rise scenario should start acting on a block of super-tidal peat. Jim M says that he does this by starting with a block of pure mineral matter, thinly sliced into regular arbitrary depth increments (1 cm cohorts?) with no real age. He then adds roots over and over until it's a block of peat. My questions are 1. is this an accurate description of how this initial condition gets made? 2. Can we use the existing addCohort or runToEquilibrium functions to do this @ktoddbrown? 3. At what point do you stop adding roots?
Upland soil: If the starting elevation is above the highest tide line and above the maximum growing elevation of wetland plants, but sea-level will raise one or both of those enough to affect the initial elevation, we need to initialize the MEM-scenario with a set of upland soil cohorts. After talking with Jim M I think this is a bit simpler. We have a generic upland set of updated 1-cm soil cohorts, 50% slow pool organic matter, 50% mineral matter by mass. We should also have an optional input so that users can input their own set up upland soil cohorts. Is this right Jim M?
Right now a step in
runMemWithCohortsis to build the block of cohorts to start with before it runs the annual time steps of the sea-level rise scenario. It uses the initial conditions and therunToEquilibriumfunction to generate a set of cohorts given the assumption of current day equilibrium conditions.After talking to @tilbud, it seems like we need to cover 2 new sets of potential starting points for the dynamic step before
runMemWithCohortsresembles the latest versions Jim M has in the excel macros.Supertidal peat: If the starting elevation is above the highest tide line, but below the maximum growing elevation of the plants, the sea-level rise scenario should start acting on a block of super-tidal peat. Jim M says that he does this by starting with a block of pure mineral matter, thinly sliced into regular arbitrary depth increments (1 cm cohorts?) with no real age. He then adds roots over and over until it's a block of peat. My questions are 1. is this an accurate description of how this initial condition gets made? 2. Can we use the existing
addCohortorrunToEquilibriumfunctions to do this @ktoddbrown? 3. At what point do you stop adding roots?Upland soil: If the starting elevation is above the highest tide line and above the maximum growing elevation of wetland plants, but sea-level will raise one or both of those enough to affect the initial elevation, we need to initialize the MEM-scenario with a set of upland soil cohorts. After talking with Jim M I think this is a bit simpler. We have a generic upland set of updated 1-cm soil cohorts, 50% slow pool organic matter, 50% mineral matter by mass. We should also have an optional input so that users can input their own set up upland soil cohorts. Is this right Jim M?