Need for species selection in findRxnFlux

[maringos]

Hi,

It is more a request than a bug in the findRxnFlux

BacArena/R/Arena.R

Line 3038 in d07b125

for(i in 1:length(sim@specs)){

There are cases (see below an example from a gapseq model) where the nomenclature of the exchange reactions does not follow the nomenclature of the metabolites and this breaks the code. As it is a pity that the whole function does not work owning to one model in multi-model simulations, I would like to ask you for adding an argument, so that we choose each time the species (one or many)

Regards,
Georgios

printReaction(model,react="EX_h2o(e)")
abbreviation	equation
EX_h2o(e)	(1) cpd00001[e0] <==>