Changes between Initial Version and Version 1 of Ticket #1482
- Timestamp:
- 2016-06-07 02:14:54 (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Ticket #1482 – Description
initial v1 3 3 * Moved to an adaptive number of sub-timesteps per ecosim timestep (used to be fixed at 3). The number of timesteps is based on how far the concentration is away from the estimated equilibrium, and the fastest loss-rate parameter. The largest timestep is 0.1 times 1/(max loss rate). 4 4 * Changed the time-integration method to explicit Adams-Bashforth multistep. Previously, a semi-analytical solution was assumed over the timestep, but this was not great at conserving tracer mass. 5 * Bug fix: In SaveEcosimTimeStepData, removed the division by biomass. The tracer concentration was being divided by biomass to get the conc/biomass, but then being divided again in the graphing/output section.6 5 7 6 For Ecospace Ecotracer calcuations: 8 7 9 * Fixed initialization of spatial !CCell array for the environment (as a function of !RelCin and !CZero)8 * Fixed initialization of spatial CCell array for the environment (as a function of !RelCin and CZero) 10 9 * Same as Ecosim version, moved to an adaptive number of sub-timesteps per ecosapce timestep (used to be fixed at 1). Same criteria as ecosim to get the timestep. 11 10 * Changed integration method same as Ecosim. Previously this was a single-step backward Euler implicit integration. The old method performed very poorly if the loss rates were large compared to the timestep. 12 11 * Added a !SolveCellC routine that performs all the Ecotracer sections of the !SolveCell code since this needs to be done multiple times per ecospace timestep (while the rest of !SolveCell does not). 13 * Added multi-threading for the !SolveCellC calculations.12 * Added multi-threading for the SolveCellC calculations. 14 13 15 14 For *ALL* Ecospace calculations: 16 15 17 * In !setMovementParameters: set the Depth, width, !HabCap, Xvel and Yvel variables on the boundary to inherit the values just inside the boundary. Previously, open boundaries would sometimes appear to be closed under advection/diffusion since these values were not being set properly.16 * In setMovementParameters: set the Depth, width, !HabCap, Xvel and Yvel variables on the boundary to inherit the values just inside the boundary. Previously, open boundaries would sometimes appear to be closed under advection/diffusion since these values were not being set properly. 18 17 19 18 As far as interface stuff goes, I think we need one more column for the tracer input: excretion or metabolism rate. This should be the rate that the contaminant is excreted from itself into the environment. Currently there is a !CexcretionRate parameter, but the name should probably be changed to !CexcretionFraction or something, since it is the fraction of contaminant that is assimilated into tissue, and is not a rate at all. Maybe the !CexcretionRate should be used for excretion/metabolism, and a new variable made like !CassimFrac or something to represent the fraction assimilated by the organism.
