Potential bug when trying to run BNG model
I am having difficulties running an imported BNG model.
Two things I noticed:
-
So apparently CellBlender parses the imported BNG model and creates its own *.bngl file in .../BNGL_import_test_files. This internal BNG model is then used for CellBlender simulations and further parsing (not entirely sure what's going on under the hood).
-
Every kinetic rate parameter in an BNG model that is NOT defined explicitely but instead in a way like
begin parameters k 0.01 end parameters (...) begin reaction rules A() + B() -> C() k*10 end reaction rulesgets redefined in the internal BNG model as
begin parameters k 0.01 und_rateLaw1 k*0.1 end parameters (...) begin reaction rules A() + B() -> C() _rateLaw1 end reaction rules
However, notice that while in the parameters block the parameter is redefined explicitely as und_rateLaw1, the reaction rule uses _rateLaw1, which is naturally undefined.
I suspect this is unwanted behaviour and leads to errors further down the line.