// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit substance = mole;
unit volume = litre;
unit area = metre^2;
unit length = metre;
unit time = second;

// 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 volume;
  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;
  private real DHAP.amt substance;
  real DHAP substance/volume;
  real DHAP.init substance/volume;
  private real G3P.amt(time) substance;
  real G3P(time) substance/volume;
  real G3P.init substance/volume;
  private real Gly.amt substance;
  real Gly substance/volume;
  real Gly.init substance/volume;
  real v1.rate(time) substance/time;
  real v2.rate(time) substance/time;

  // equations
  DHAP.amt = DHAP*compartment;
  DHAP = DHAP.init;
  DHAP.init = .59;
  when (time=time.min) G3P.amt = G3P.init*compartment;
  G3P.amt:time = v1.rate + -1*v2.rate;
  G3P = G3P.amt/compartment;
  G3P.init = 0;
  Gly.amt = Gly*compartment;
  Gly = Gly.init;
  Gly.init = 0;
  v1.rate = 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.rate = 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.amt.sbmlRole="speciesAmount";
  DHAP.sbmlRole="speciesConcentration";
  DHAP.sbmlCompartment="compartment";
  DHAP.init.sbmlRole="speciesInitialConcentration";
  G3P.amt.sbmlRole="speciesAmount";
  G3P.sbmlRole="speciesConcentration";
  G3P.sbmlCompartment="compartment";
  G3P.init.sbmlRole="speciesInitialConcentration";
  Gly.amt.sbmlRole="speciesAmount";
  Gly.sbmlRole="speciesConcentration";
  Gly.sbmlCompartment="compartment";
  Gly.init.sbmlRole="speciesInitialConcentration";
  v1.rate.sbmlRole="rate";
  v2.rate.sbmlRole="rate";
}