Reaction-Diffusion systems for reaction schemes. The schemes are the same as are used in RScalc, but concentrations are allowed to vary in space as well as time. The geometry may be either spherical or linear. The D-diff fields in the individual species give the diffusion constants in micron/s.
The system is computed on a grid of npoints points from rmin to rmax (micron) using a multi-dimensional Newton's method to iterate to the implicit solution at each time step. The points are spaced logarithmically in r for spherical diffusion and uniformly for linear diffusion.
Other features are the same as for RSCalc, with the exception that the conc-hist and injection options for individual pools apply only to the first point. In the first case, the concentration at the first point follows the specified concentration history, and all other points continue as before. In the second case, material is injected at the first point according to the given time profile, where the scaling is such that a Unit injection would raise the concentration in a cubic length Unit by one concentration Unit in one time Unit. So, if you choose to measure lengths in microns, concentrations in mM and time in ms, the Units for injections would be 10^-15 moles per second, which, as a flux of monovalent ions corresponds to about 96 pA.
The results of the calculation are the concentrations as a function of time and position for each of the species.
|    field||    type  ||    contents / default value||   |
|timestep||double||0.1||variable: timestep type: double|
|runtime||double||10.0||variable: runtime type: double|
|geometry||choice||variable: geometry type: choice|
|npoints||integer||0||variable: npoints type: integer|
|rmax||double||0.0||variable: rmax type: double|
|rmin||double||0.0||variable: rmin type: double|
|maxtimepts||integer||0||variable: maxtimepts type: integer|
|scheme||list selection; boolean; one only||ReactionScheme|