C
C There are a total of 13 entries in the algebraic variable array.
C There are a total of 8 entries in each of the rate and state variable arrays.
C There are a total of 25 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is Pyr in component Pyr (micromolar).
C CONSTS(25) is v1 in component v1 (micromolar_per_second).
C ALGBRC(1) is v2 in component v2 (micromolar_per_second).
C ALGBRC(6) is v7 in component v7 (micromolar_per_second).
C STATES(2) is AcCoA in component AcCoA (micromolar).
C ALGBRC(2) is v3 in component v3 (micromolar_per_second).
C STATES(3) is Cit in component Cit (micromolar).
C ALGBRC(3) is v4 in component v4 (micromolar_per_second).
C STATES(4) is KG in component KG (micromolar).
C ALGBRC(4) is v5 in component v5 (micromolar_per_second).
C ALGBRC(5) is v6 in component v6 (micromolar_per_second).
C STATES(5) is OAA in component OAA (micromolar).
C ALGBRC(7) is v8 in component v8 (micromolar_per_second).
C STATES(6) is NAD in component NAD (micromolar).
C ALGBRC(10) is vresp in component vresp (micromolar_per_second).
C STATES(7) is ATP in component ATP (micromolar).
C ALGBRC(13) is vATP in component vATP (micromolar_per_second).
C ALGBRC(8) is vANT in component vANT (micromolar_per_second).
C STATES(8) is delta_psi in component delta_psi (millivolt).
C CONSTS(1) is C in component delta_psi (millimolar_per_millivolt).
C ALGBRC(9) is vleak in component vleak (micromolar_per_second).
C CONSTS(2) is k1 in component v1 (micromolar_per_second).
C CONSTS(3) is k2 in component v2 (second_order_rate_constant).
C CONSTS(4) is k3 in component v3 (second_order_rate_constant).
C CONSTS(5) is k4 in component v4 (second_order_rate_constant).
C CONSTS(6) is k5 in component v5 (third_order_rate_constant).
C CONSTS(7) is At in component model_parameters (millimolar).
C CONSTS(8) is k6 in component v6 (first_order_rate_constant).
C CONSTS(9) is k7 in component v7 (second_order_rate_constant).
C CONSTS(10) is k8 in component v8 (first_order_rate_constant).
C CONSTS(11) is kANT in component vANT (first_order_rate_constant).
C CONSTS(12) is kleak in component vleak (molar_per_millivolt_per_second).
C CONSTS(13) is kresp in component vresp (millimolar_per_second).
C CONSTS(14) is K in component vresp (millimolar).
C CONSTS(15) is a in component vresp (per_millivolt).
C CONSTS(16) is delta_psi_m in component vresp (millivolt).
C CONSTS(17) is Nt in component model_parameters (millimolar).
C CONSTS(18) is kATP in component vATP (millimolar_per_second).
C CONSTS(19) is b in component vATP (per_micromolar).
C ALGBRC(12) is ATP_crit_delta_psi in component ATP_crit_delta_psi (micromolar).
C CONSTS(20) is R in component ATP_crit_delta_psi (joule_per_mole_kelvin).
C CONSTS(21) is T in component ATP_crit_delta_psi (kelvin).
C CONSTS(22) is F in component ATP_crit_delta_psi (coulomb_per_mole).
C CONSTS(23) is Kapp in component ATP_crit_delta_psi (per_millimolar).
C CONSTS(24) is Pi in component ATP_crit_delta_psi (millimolar).
C ALGBRC(11) is delta_G_transport in component ATP_crit_delta_psi (joule_per_mole).
C RATES(1) is d/dt Pyr in component Pyr (micromolar).
C RATES(2) is d/dt AcCoA in component AcCoA (micromolar).
C RATES(3) is d/dt Cit in component Cit (micromolar).
C RATES(4) is d/dt KG in component KG (micromolar).
C RATES(5) is d/dt OAA in component OAA (micromolar).
C RATES(6) is d/dt NAD in component NAD (micromolar).
C RATES(7) is d/dt ATP in component ATP (micromolar).
C RATES(8) is d/dt delta_psi in component delta_psi (millivolt).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.154
      STATES(2) = 0.063
      STATES(3) = 0.44
      STATES(4) = 0.225
      STATES(5) = 0.005
      STATES(6) = 0.856
      STATES(7) = 3.536
      STATES(8) = 150.0
      CONSTS(1) = 6.75e-06
      CONSTS(2) = 38.0
      CONSTS(3) = 152.0
      CONSTS(4) = 57142.0
      CONSTS(5) = 53.0
      CONSTS(6) = 82361.0
      CONSTS(7) = 4.160
      CONSTS(8) = 3.2e-3
      CONSTS(9) = 40.0
      CONSTS(10) = 3.6
      CONSTS(11) = 0.1
      CONSTS(12) = 0.426
      CONSTS(13) = 2.5
      CONSTS(14) = 2
      CONSTS(15) = 0.1
      CONSTS(16) = 150.0
      CONSTS(17) = 1.070
      CONSTS(18) = 131.9
      CONSTS(19) = 4
      CONSTS(20) = 8.314
      CONSTS(21) = 298
      CONSTS(22) = 96485
      CONSTS(23) = 4.4e-6
      CONSTS(24) = 2.440
      CONSTS(25) = CONSTS(2)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(3)*STATES(1)*STATES(6)
      ALGBRC(2) =  CONSTS(4)*STATES(5)*STATES(2)
      RATES(2) = ALGBRC(1) - ALGBRC(2)
      ALGBRC(3) =  CONSTS(5)*STATES(3)*STATES(6)
      RATES(3) = ALGBRC(2) - ALGBRC(3)
      ALGBRC(4) =  CONSTS(6)*STATES(4)*STATES(6)*(CONSTS(7) - STATES(7))
      ALGBRC(5) =  CONSTS(8)*(STATES(5) - STATES(4))
      RATES(4) = (ALGBRC(3)+ALGBRC(5)) - ALGBRC(4)
      ALGBRC(6) =  CONSTS(9)*STATES(1)*STATES(7)
      RATES(1) = CONSTS(25) - (ALGBRC(1)+ALGBRC(6))
      ALGBRC(7) =  CONSTS(10)*STATES(5)
      RATES(5) = (ALGBRC(4)+ALGBRC(6)) - (ALGBRC(2)+ALGBRC(7)+ALGBRC(5))
      ALGBRC(10) =  CONSTS(13)*((CONSTS(17) - STATES(6))/((CONSTS(14)+CONSTS(17)) - STATES(6)))*(1.00000/(1.00000+EXP( CONSTS(15)*(STATES(8) - CONSTS(16)))))
      RATES(6) = ALGBRC(10) - (ALGBRC(1)+ALGBRC(3)+ 2.00000*ALGBRC(4))
      ALGBRC(11) =  0.00120000*CONSTS(22)*STATES(8)
      ALGBRC(12) = CONSTS(7)/(1.00000+EXP(( -3.00000*ALGBRC(11))/( CONSTS(20)*CONSTS(21)))/( CONSTS(23)*CONSTS(24)))
      ALGBRC(13) =  CONSTS(18)*(2.00000/(1.00000+EXP( CONSTS(19)*(STATES(7) - ALGBRC(12)))) - 1.00000)
      ALGBRC(8) =  CONSTS(11)*STATES(7)
      RATES(7) = (ALGBRC(13)+ALGBRC(4)) - (ALGBRC(8)+ALGBRC(6))
      ALGBRC(9) =  CONSTS(12)*STATES(8)
      RATES(8) =  (1.00000/CONSTS(1))*( 10.0000*ALGBRC(10) - ( 3.00000*ALGBRC(13)+ALGBRC(9)+ALGBRC(8)))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(3)*STATES(1)*STATES(6)
      ALGBRC(2) =  CONSTS(4)*STATES(5)*STATES(2)
      ALGBRC(3) =  CONSTS(5)*STATES(3)*STATES(6)
      ALGBRC(4) =  CONSTS(6)*STATES(4)*STATES(6)*(CONSTS(7) - STATES(7))
      ALGBRC(5) =  CONSTS(8)*(STATES(5) - STATES(4))
      ALGBRC(6) =  CONSTS(9)*STATES(1)*STATES(7)
      ALGBRC(7) =  CONSTS(10)*STATES(5)
      ALGBRC(10) =  CONSTS(13)*((CONSTS(17) - STATES(6))/((CONSTS(14)+CONSTS(17)) - STATES(6)))*(1.00000/(1.00000+EXP( CONSTS(15)*(STATES(8) - CONSTS(16)))))
      ALGBRC(11) =  0.00120000*CONSTS(22)*STATES(8)
      ALGBRC(12) = CONSTS(7)/(1.00000+EXP(( -3.00000*ALGBRC(11))/( CONSTS(20)*CONSTS(21)))/( CONSTS(23)*CONSTS(24)))
      ALGBRC(13) =  CONSTS(18)*(2.00000/(1.00000+EXP( CONSTS(19)*(STATES(7) - ALGBRC(12)))) - 1.00000)
      ALGBRC(8) =  CONSTS(11)*STATES(7)
      ALGBRC(9) =  CONSTS(12)*STATES(8)
      RETURN
      END