Merge pull request #108 from SciML/aj/askem

disable converting nonconstant parameters to states

Author: paulflang

src/SBMLToolkit.jl

@@ -4,6 +4,7 @@ using Catalyst
 using SBML
 using SymbolicUtils
 
+include("drafts.jl")
 include("systems.jl")
 include("reactions.jl")
 include("rules.jl")

src/drafts.jl

@@ -0,0 +1,10 @@
+function is_event_assignment(k, model)
+    for ev in values(model.events)
+        for as in ev.event_assignments
+            if as.variable == k
+                return true
+            end
+        end
+    end
+    return false
+end

src/systems.jl

@@ -82,7 +82,7 @@ function get_mappings(model::SBML.Model)
         end
     end
     for (k, v) in model.parameters
-        if v.constant == false && SBML.seemsdefined(k, model)
+        if v.constant == false && (SBML.seemsdefined(k, model) || is_event_assignment(k, model))
             var = create_var(k, IV; isbcspecies = true)
             push!(u0map, var => v.value)
         else