Unit issue with compartmentless models

[phantomas1234]

constructModel[{reaction["b1", {}, {metabolite["x1", None]}, {1}, False]}, Parameters -> {rateconst["b1", True] -> 0.1 Hour^-1}]

does not raise an error although the correct unit should be

rateconst["b1", True] -> 0.1 * Hour^-1 * Liter]

[avsastry]

Seems like this is a problem with the unit checking function in general. There is no error for either constructModel[{reaction["b1", {metabolite["x0", None]}, {metabolite["x1", None]}, {1, 1}, False]}, Parameters -> {rateconst["b1", True] -> 0.1*Hour^-1}, UnitChecking -> True] or constructModel[{reaction["b1", {metabolite["x0", None]}, {metabolite["x1", None]}, {1, 1}, False]}, Parameters -> {rateconst["b1", True] -> 0.1*Liter*Hour^-1}, UnitChecking -> True], even though only the first is correct.

[avsastry]

Also, shouldn't the units of the compartmentless rate constant be Liter/(Mole*Hour)?

[phantomas1234]

Its a tricky issue, although I've probably didn't intend this ambiguity,
both results are kinda ok (in the one case one ends up with a concentration
flux, in the other with a particle flux). Concentration fluxes are a
dangerous thing only in multi-compartment models but for no-compartment
models using both in the ODEs is fine I guess. Check out
http://www.copasi.org/tiki-download_file.php?fileId=73 for an explanation.
See also my notebook.

From what I can see, there should definitely be no quantity unit in the
rate constant's units because there is only one metabolite in the reaction,
right?

Best,

Niko

On Wed, Dec 3, 2014 at 9:53 AM, avsastry [email protected] wrote:

Also, shouldn't the units of the compartmentless rate constant be Liter *
Mole^-1 * Hour^-1?

—
Reply to this email directly or view it on GitHub
#21 (comment)
.