// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// v1: DHAP => G3P 
// v2: G3P => Gly 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real compartment = 1 L;
  real F16BP = 6.01;
  real NADH = 1.87;
  real NAD = 1.45;
  real ATP = 2.37;
  real ADP = 2.17;
  real Phi = 1;
  real Vf1 = 47;
  real Keq1 = 1E4;
  real K1g3p = 1.2;
  real K1nadh = .023;
  real K1nad = .93;
  real K1dhap = .54;
  real K1f16bp = 4.8;
  real K1atp = .73;
  real K1adp = 2;
  real V2 = 53;
  real K2g3p = 3.5;
  real K2phi = 1;
  real DHAP = .59 M;
  real G3P(time) M;
  real Gly = 0 M;
  real v1(time) katal;
  real v2(time) katal;

  // equations
  when (time=time.min) G3P = 0;
  G3P:time = (v1 + -1*v2)/compartment;
  v1 = Vf1/(K1nadh*K1dhap)*(NADH*DHAP-NAD*G3P/Keq1)/((1+F16BP/K1f16bp+ATP/K1atp+ADP/K1adp)*(1+NADH/K1nadh+NAD/K1nad)*(1+DHAP/K1dhap+G3P/K1g3p));
  v2 = V2*G3P/K2g3p/((1+G3P/K2g3p)*(1+Phi/K2phi));

  // variable properties
  compartment.sbmlRole="compartment";
  F16BP.sbmlRole="parameter";
  NADH.sbmlRole="parameter";
  NAD.sbmlRole="parameter";
  ATP.sbmlRole="parameter";
  ADP.sbmlRole="parameter";
  Phi.sbmlRole="parameter";
  Vf1.sbmlRole="parameter";
  Keq1.sbmlRole="parameter";
  K1g3p.sbmlRole="parameter";
  K1nadh.sbmlRole="parameter";
  K1nad.sbmlRole="parameter";
  K1dhap.sbmlRole="parameter";
  K1f16bp.sbmlRole="parameter";
  K1atp.sbmlRole="parameter";
  K1adp.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  K2g3p.sbmlRole="parameter";
  K2phi.sbmlRole="parameter";
  DHAP.sbmlRole="species";
  DHAP.sbmlCompartment="compartment";
  G3P.sbmlRole="species";
  G3P.sbmlCompartment="compartment";
  Gly.sbmlRole="species";
  Gly.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
}