def model{Glycolysis_model} Glycolysis_model as
    // define components 

    def unit millimolar as
        unit mole {pref: milli};
        unit litre {expo: -1};
    enddef;

    def unit minute as
        unit second {mult: 60.0};
    enddef;

    def unit flux as
        unit millimolar;
        unit minute {expo: -1};
    enddef;

    def comp environment as
        var time: minute {pub: out};
    enddef;

    def comp GLC as
        var GLC: millimolar {init: 10, pub: out};
        var v_GK: flux {pub: in};
        var time: minute {pub: in};

        ode(GLC, time) = -v_GK;
    enddef;

    def comp G6P as
        var G6P: millimolar {init: 0, pub: out};
        var v_GK: flux {pub: in};
        var v_PGI: flux {pub: in};
        var time: minute {pub: in};

        ode(G6P, time) = v_GK-v_PGI;
    enddef;

    def comp F6P as
        var F6P: millimolar {init: 0, pub: out};
        var v_PGI: flux {pub: in};
        var v_PFK: flux {pub: in};
        var time: minute {pub: in};

        ode(F6P, time) = v_PGI-v_PFK;
    enddef;

    def comp FBP as
        var FBP: millimolar {init: 0, pub: out};
        var v_PFK: flux {pub: in};
        var v_FBPA: flux {pub: in};
        var time: minute {pub: in};

        ode(FBP, time) = v_PFK-v_FBPA;
    enddef;

    def comp TRIO as
        var TRIO: millimolar {init: 0, pub: out};
        var v_FBPA: flux {pub: in};
        var v_GAPD: flux {pub: in};
        var time: minute {pub: in};

        ode(TRIO, time) = 2.0{dimensionless}*v_FBPA-v_GAPD;
    enddef;

    def comp BPG as
        var BPG: millimolar {init: 0, pub: out};
        var v_GAPD: flux {pub: in};
        var time: minute {pub: in};

        ode(BPG, time) = v_GAPD;
    enddef;

    def comp P as
        var P: millimolar {init: 100, pub: out};
        var time: minute {pub: in};
        var v_GAPD: flux {pub: in};

        ode(P, time) = -v_GAPD;
    enddef;

    def comp NADH as
        var NADH: millimolar {init: 0, pub: out};
        var v_GAPD: flux {pub: in};
        var time: minute {pub: in};

        ode(NADH, time) = v_GAPD;
    enddef;

    def comp NAD as
        var NAD: millimolar {init: 10, pub: out};
        var v_GAPD: flux {pub: in};
        var time: minute {pub: in};

        ode(NAD, time) = -v_GAPD;
    enddef;

    def comp DHAP as
        var DHAP: millimolar {pub: out};
        var Keq_TPI: dimensionless {init: 0.045, pub: out};
        var TRIO: millimolar {pub: in};

        DHAP = TRIO/(1.0{dimensionless}+Keq_TPI);
    enddef;

    def comp GAP as
        var GAP: millimolar {pub: out};
        var Keq_TPI: dimensionless {pub: in};
        var TRIO: millimolar {pub: in};

        GAP = TRIO*Keq_TPI/(1.0{dimensionless}+Keq_TPI);
    enddef;

    def comp ADP as
        var ADP: millimolar {init: 0, pub: out};
        var v_GK: flux {pub: in};
        var v_PFK: flux {pub: in};
        var time: minute {pub: in};

        ode(ADP, time) = v_GK+v_PFK;
    enddef;

    def comp ATP as
        var ATP: millimolar {init: 7, pub: out};
        var v_GK: flux {pub: in};
        var v_PFK: flux {pub: in};
        var time: minute {pub: in};

        ode(ATP, time) = -(v_GK+v_PFK);
    enddef;

    // define enzyme rates 

    def comp v_PGI as
        var v_PGI: flux {pub: out};
        var K_m_G6P: millimolar {init: 1.4};
        var V_max: flux {init: 50};
        var G6P: millimolar {pub: in};

        v_PGI = V_max*G6P/(K_m_G6P+G6P);
    enddef;

    def comp v_GK as
        var v_GK: flux {pub: out};
        var Km_GLC: millimolar {init: 0.11};
        var K_hlf_ATP: millimolar {init: 1.7};
        var h_ATP: dimensionless {init: 1.2};
        var v_max: flux {init: 0.28};
        var GLC: millimolar {pub: in};
        var ATP: millimolar {pub: in};

        v_GK = v_max*pow(ATP, h_ATP)*GLC/(pow(K_hlf_ATP, h_ATP)*Km_GLC+pow(ATP, h_ATP)*GLC);
    enddef;

    def comp v_GAPD as
        var v_GAPD: flux {pub: out};
        var V_max: flux {init: 1.3};
        var Km_GAP: millimolar {init: 2.4};
        var Km_NAD: millimolar {init: 0.9};
        var Km_P: millimolar {init: 13};
        var GAP: millimolar {pub: in};
        var NAD: millimolar {pub: in};
        var P: millimolar {pub: in};

        v_GAPD = V_max*GAP*NAD*P/(Km_GAP*Km_NAD*Km_P+GAP*NAD*P);
    enddef;

    def comp v_FBPA as
        var v_FBPA: flux {pub: out};
        var KM_FBP: millimolar {init: 4.4};
        var v_max: flux {init: 0.39};
        var K_i_FBP: millimolar {init: 10};
        var FBP: millimolar {pub: in};

        v_FBPA = v_max*FBP/(KM_FBP+FBP*(1.0{dimensionless}+FBP)/K_i_FBP);
    enddef;

    def comp v_PFK as
        var v_PFK: flux {pub: out};
        var v_max: flux {init: 0.035};
        var K_hlf_F6P: millimolar {init: 0.74};
        var h_F6P: dimensionless {init: 1.7};
        var Km_ATP: millimolar {init: 1.1};
        var F6P: millimolar {pub: in};
        var ATP: millimolar {pub: in};

        v_PFK = v_max*pow(F6P, h_F6P)*ATP/(pow(K_hlf_F6P, h_F6P)*Km_ATP+pow(F6P, h_F6P)*ATP);
    enddef;

    // Connections 

    def map between GLC and environment for
        vars time and time;
    enddef;

    def map between G6P and environment for
        vars time and time;
    enddef;

    def map between F6P and environment for
        vars time and time;
    enddef;

    def map between FBP and environment for
        vars time and time;
    enddef;

    def map between TRIO and environment for
        vars time and time;
    enddef;

    def map between BPG and environment for
        vars time and time;
    enddef;

    def map between P and environment for
        vars time and time;
    enddef;

    def map between NADH and environment for
        vars time and time;
    enddef;

    def map between NAD and environment for
        vars time and time;
    enddef;

    def map between ATP and environment for
        vars time and time;
    enddef;

    def map between ADP and environment for
        vars time and time;
    enddef;

    def map between GLC and v_GK for
        vars GLC and GLC;
        vars v_GK and v_GK;
    enddef;

    def map between G6P and v_GK for
        vars v_GK and v_GK;
    enddef;

    def map between G6P and v_PGI for
        vars G6P and G6P;
        vars v_PGI and v_PGI;
    enddef;

    def map between F6P and v_PGI for
        vars v_PGI and v_PGI;
    enddef;

    def map between F6P and v_PFK for
        vars F6P and F6P;
        vars v_PFK and v_PFK;
    enddef;

    def map between FBP and v_PFK for
        vars v_PFK and v_PFK;
    enddef;

    def map between FBP and v_FBPA for
        vars FBP and FBP;
        vars v_FBPA and v_FBPA;
    enddef;

    def map between TRIO and v_FBPA for
        vars v_FBPA and v_FBPA;
    enddef;

    def map between TRIO and v_GAPD for
        vars v_GAPD and v_GAPD;
    enddef;

    def map between BPG and v_GAPD for
        vars v_GAPD and v_GAPD;
    enddef;

    def map between ATP and v_GK for
        vars ATP and ATP;
        vars v_GK and v_GK;
    enddef;

    def map between ATP and v_PFK for
        vars ATP and ATP;
        vars v_PFK and v_PFK;
    enddef;

    def map between NADH and v_GAPD for
        vars v_GAPD and v_GAPD;
    enddef;

    def map between NAD and v_GAPD for
        vars NAD and NAD;
        vars v_GAPD and v_GAPD;
    enddef;

    def map between DHAP and TRIO for
        vars TRIO and TRIO;
    enddef;

    def map between GAP and TRIO for
        vars TRIO and TRIO;
    enddef;

    def map between DHAP and GAP for
        vars Keq_TPI and Keq_TPI;
    enddef;

    def map between P and v_GAPD for
        vars P and P;
    enddef;

    def map between GAP and v_GAPD for
        vars GAP and GAP;
    enddef;

    def map between ADP and v_GK for
        vars v_GK and v_GK;
    enddef;

    def map between ADP and v_PFK for
        vars v_PFK and v_PFK;
    enddef;

    def map between P and v_GAPD for
        vars v_GAPD and v_GAPD;
    enddef;
enddef;