C
C There are a total of 34 entries in the algebraic variable array.
C There are a total of 16 entries in each of the rate and state variable arrays.
C There are a total of 47 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C ALGBRC(1) is v1 in component v1 (flux).
C CONSTS(1) is Vmax in component v1 (flux).
C CONSTS(2) is Km in component v1 (umol_per_g).
C STATES(1) is GL in component GL (umol_per_g).
C ALGBRC(2) is v2 in component v2 (flux).
C CONSTS(3) is Vmax in component v2 (flux).
C CONSTS(4) is Km in component v2 (umol_per_g).
C STATES(2) is O2 in component O2 (umol_per_g).
C ALGBRC(16) is v3 in component v3 (flux).
C CONSTS(5) is Vmax in component v3 (first_order_rate_constant).
C ALGBRC(15) is k3 in component v3 (first_order_rate_constant).
C STATES(3) is FA in component FA (umol_per_g).
C ALGBRC(14) is PS in component PS (dimensionless).
C CONSTS(6) is PS0 in component PS0 (dimensionless).
C ALGBRC(31) is v4 in component v4 (flux).
C CONSTS(7) is Vmax in component v4 (first_order_rate_constant).
C ALGBRC(30) is k4 in component v4 (first_order_rate_constant).
C CONSTS(8) is epsilon in component v4 (dimensionless).
C ALGBRC(5) is RS in component RS (dimensionless).
C CONSTS(9) is RS0 in component RS0 (dimensionless).
C ALGBRC(29) is AF in component AF (dimensionless).
C CONSTS(10) is AF0 in component AF0 (dimensionless).
C ALGBRC(18) is v5 in component v5 (flux).
C CONSTS(11) is Vmax in component v5 (first_order_rate_constant).
C ALGBRC(17) is k5 in component v5 (first_order_rate_constant).
C CONSTS(12) is epsilon in component v5 (dimensionless).
C STATES(4) is GP in component GP (umol_per_g).
C ALGBRC(20) is v6 in component v6 (flux).
C CONSTS(13) is Vmax in component v6 (first_order_rate_constant).
C ALGBRC(19) is k6 in component v6 (first_order_rate_constant).
C CONSTS(14) is epsilon in component v6 (dimensionless).
C ALGBRC(4) is CS in component CS (dimensionless).
C CONSTS(15) is CS0 in component CS0 (dimensionless).
C ALGBRC(24) is v7 in component v7 (flux).
C CONSTS(16) is Vmax in component v7 (first_order_rate_constant).
C ALGBRC(22) is k7 in component v7 (first_order_rate_constant).
C CONSTS(17) is epsilon in component v7 (dimensionless).
C STATES(5) is GY in component GY (umol_per_g).
C ALGBRC(3) is v8 in component v8 (flux).
C CONSTS(18) is Vmax in component v8 (first_order_rate_constant).
C CONSTS(46) is k8 in component v8 (first_order_rate_constant).
C STATES(6) is TG in component TG (umol_per_g).
C ALGBRC(7) is v9 in component v9 (flux).
C CONSTS(19) is Vmax in component v9 (first_order_rate_constant).
C ALGBRC(6) is k9 in component v9 (first_order_rate_constant).
C STATES(7) is PY in component PY (umol_per_g).
C ALGBRC(33) is v10 in component v10 (flux).
C CONSTS(20) is Vmax in component v10 (first_order_rate_constant).
C ALGBRC(32) is k10 in component v10 (first_order_rate_constant).
C CONSTS(21) is epsilon in component v10 (dimensionless).
C ALGBRC(9) is v11 in component v11 (flux).
C CONSTS(22) is Vmax in component v11 (first_order_rate_constant).
C ALGBRC(8) is k11 in component v11 (first_order_rate_constant).
C STATES(8) is LA in component LA (umol_per_g).
C ALGBRC(23) is v12 in component v12 (flux).
C CONSTS(23) is Vmax in component v12 (first_order_rate_constant).
C ALGBRC(21) is k12 in component v12 (first_order_rate_constant).
C CONSTS(24) is epsilon in component v12 (dimensionless).
C STATES(9) is AC in component AC (umol_per_g).
C ALGBRC(26) is v13 in component v13 (flux).
C CONSTS(25) is Vmax in component v13 (first_order_rate_constant).
C ALGBRC(25) is k13 in component v13 (first_order_rate_constant).
C STATES(10) is CR in component CR (umol_per_g).
C ALGBRC(28) is v14 in component v14 (flux).
C CONSTS(26) is Vmax in component v14 (first_order_rate_constant).
C ALGBRC(27) is k14 in component v14 (first_order_rate_constant).
C STATES(11) is PC in component PC (umol_per_g).
C ALGBRC(11) is v15 in component v15 (flux).
C CONSTS(27) is Vmax in component v15 (first_order_rate_constant).
C ALGBRC(10) is k15 in component v15 (first_order_rate_constant).
C CONSTS(28) is epsilon in component v15 (dimensionless).
C STATES(12) is CoA_pool in component CoA_pool (umol_per_g).
C STATES(13) is FC in component FC (umol_per_g).
C CONSTS(29) is FC0 in component FC0 (umol_per_g).
C ALGBRC(13) is v16 in component v16 (flux).
C CONSTS(30) is Vmax in component v16 (first_order_rate_constant).
C ALGBRC(12) is k16 in component v16 (first_order_rate_constant).
C CONSTS(31) is epsilon in component v16 (dimensionless).
C CONSTS(47) is v17 in component v17 (flux).
C CONSTS(32) is Vmax in component v17 (first_order_rate_constant).
C CONSTS(45) is k17 in component v17 (first_order_rate_constant).
C CONSTS(33) is ATP in component ATP (umol_per_g).
C CONSTS(34) is aGL in component GL (umol_per_ml).
C CONSTS(35) is sigmaGL in component GL (g_per_ml).
C ALGBRC(34) is F in component model_parameters (ml_per_g_min).
C CONSTS(36) is aFA in component FA (umol_per_ml).
C CONSTS(37) is sigmaFA in component FA (g_per_ml).
C CONSTS(38) is aLA in component LA (umol_per_ml).
C CONSTS(39) is sigmaLA in component LA (g_per_ml).
C CONSTS(40) is aO2 in component O2 (umol_per_ml).
C CONSTS(41) is sigmaO2 in component O2 (g_per_ml).
C STATES(14) is CO2 in component CO2 (umol_per_g).
C CONSTS(42) is aCO2 in component CO2 (umol_per_ml).
C CONSTS(43) is sigmaCO2 in component CO2 (g_per_ml).
C STATES(15) is NAD in component NAD (umol_per_g).
C STATES(16) is ADP in component ADP (umol_per_g).
C CONSTS(44) is NADH in component NADH (umol_per_g).
C RATES(1) is d/dt GL in component GL (umol_per_g).
C RATES(3) is d/dt FA in component FA (umol_per_g).
C RATES(4) is d/dt GP in component GP (umol_per_g).
C RATES(5) is d/dt GY in component GY (umol_per_g).
C RATES(6) is d/dt TG in component TG (umol_per_g).
C RATES(7) is d/dt PY in component PY (umol_per_g).
C RATES(8) is d/dt LA in component LA (umol_per_g).
C RATES(9) is d/dt AC in component AC (umol_per_g).
C RATES(13) is d/dt FC in component FC (umol_per_g).
C RATES(12) is d/dt CoA_pool in component CoA_pool (umol_per_g).
C RATES(2) is d/dt O2 in component O2 (umol_per_g).
C RATES(14) is d/dt CO2 in component CO2 (umol_per_g).
C RATES(15) is d/dt NAD in component NAD (umol_per_g).
C RATES(16) is d/dt ADP in component ADP (umol_per_g).
C RATES(11) is d/dt PC in component PC (umol_per_g).
C RATES(10) is d/dt CR in component CR (umol_per_g).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 5.90
      CONSTS(2) = 0.05
      STATES(1) = 0.998
      CONSTS(3) = 67.6
      CONSTS(4) = 0.01
      STATES(2) = 0.963
      CONSTS(5) = 4.90
      STATES(3) = 0.021
      CONSTS(6) = 0.2
      CONSTS(7) = 21.3
      CONSTS(8) = 0.6
      CONSTS(9) = 0.111
      CONSTS(10) = 0.523
      CONSTS(11) = 2.82
      CONSTS(12) = 0.254
      STATES(4) = 0.171
      CONSTS(13) = 3.14
      CONSTS(14) = 0.5
      CONSTS(15) = 1.0
      CONSTS(16) = 0.0162
      CONSTS(17) = 0.5
      STATES(5) = 33.0
      CONSTS(18) = 0.005
      STATES(6) = 3.96
      CONSTS(19) = 1.8
      STATES(7) = 0.20
      CONSTS(20) = 12.6
      CONSTS(21) = 0.98
      CONSTS(22) = 0.96
      STATES(8) = 1.98
      CONSTS(23) = 695.7
      CONSTS(24) = 0.75
      STATES(9) = 0.0046
      CONSTS(25) = 0.455
      STATES(10) = 3.5
      CONSTS(26) = 0.455
      STATES(11) = 8.80
      CONSTS(27) = 626.1
      CONSTS(28) = 0.669
      STATES(12) = 0.043
      STATES(13) = 0.0088
      CONSTS(29) = 0.0088
      CONSTS(30) = 67.0
      CONSTS(31) = 0.775
      CONSTS(32) = 7.76
      CONSTS(33) = 4.5
      CONSTS(34) = 4.0
      CONSTS(35) = 3.76
      CONSTS(36) = 0.5
      CONSTS(37) = 13.2
      CONSTS(38) = 1.8
      CONSTS(39) = 0.51
      CONSTS(40) = 6.53
      CONSTS(41) = 1.0
      STATES(14) = 20.0
      CONSTS(42) = 15.5
      CONSTS(43) = 1.0
      STATES(15) = 1.81
      STATES(16) = 0.90
      CONSTS(44) = 0.19
      CONSTS(45) = CONSTS(32)
      CONSTS(46) = CONSTS(18)
      CONSTS(47) =  CONSTS(45)*CONSTS(33)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = CONSTS(44)/STATES(15)
      ALGBRC(10) =  CONSTS(27)*( CONSTS(28)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(28))*(STATES(13) ** -1.00000/(CONSTS(29) ** -1.00000+STATES(13) ** -1.00000)))
      ALGBRC(11) =  ALGBRC(10)*STATES(12)
      ALGBRC(12) =  CONSTS(30)*( CONSTS(31)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(31))*(STATES(13) ** -1.00000/(CONSTS(29) ** -1.00000+STATES(13) ** -1.00000)))
      ALGBRC(13) =  ALGBRC(12)*STATES(9)
      RATES(12) = ALGBRC(13) - ALGBRC(11)
      ALGBRC(14) = STATES(16)/CONSTS(33)
      ALGBRC(15) =  CONSTS(5)*(ALGBRC(14) ** -1.00000/(CONSTS(6) ** -1.00000+ALGBRC(14) ** -1.00000))
      ALGBRC(16) =  ALGBRC(15)*STATES(3)
      ALGBRC(3) =  CONSTS(46)*STATES(6)
      RATES(6) =  (1.00000/3.00000)*ALGBRC(16) - ALGBRC(3)
      ALGBRC(1) =  CONSTS(1)*(STATES(1)/(CONSTS(2)+STATES(1)))
      ALGBRC(17) =  CONSTS(11)*( CONSTS(12)*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14)))+ (1.00000 - CONSTS(12))*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000)))
      ALGBRC(18) =  ALGBRC(17)*STATES(4)
      ALGBRC(4) = STATES(10)/STATES(11)
      ALGBRC(19) =  CONSTS(13)*( CONSTS(14)*(ALGBRC(14) ** -1.00000/(CONSTS(6) ** -1.00000+ALGBRC(14) ** -1.00000))+ (1.00000 - CONSTS(14))*ALGBRC(4) ** -1.00000/(CONSTS(15) ** -1.00000+ALGBRC(4) ** -1.00000) ** 2.00000)
      ALGBRC(20) =  ALGBRC(19)*STATES(4)
      ALGBRC(22) =  CONSTS(16)*( CONSTS(17)*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14)))+ (1.00000 - CONSTS(17))*ALGBRC(4)/(CONSTS(15)+ALGBRC(4)) ** 2.00000)
      ALGBRC(24) =  ALGBRC(22)*STATES(5)
      RATES(4) = (ALGBRC(1)+ALGBRC(24)) - (ALGBRC(18)+ALGBRC(20))
      RATES(5) = ALGBRC(20) - ALGBRC(24)
      ALGBRC(25) =  CONSTS(25)*(ALGBRC(14) ** -1.00000/(CONSTS(6) ** -1.00000+ALGBRC(14) ** -1.00000))
      ALGBRC(26) =  ALGBRC(25)*STATES(10)
      ALGBRC(27) =  CONSTS(26)*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14)))
      ALGBRC(28) =  ALGBRC(27)*STATES(11)
      RATES(11) = ALGBRC(26) - ALGBRC(28)
      RATES(10) = ALGBRC(28) - ALGBRC(26)
      ALGBRC(2) =  CONSTS(3)*(STATES(2)/(CONSTS(4)+STATES(2)))
      ALGBRC(29) = STATES(9)/STATES(13)
      ALGBRC(30) =  CONSTS(7)*( CONSTS(8)*(ALGBRC(29) ** -1.00000/(CONSTS(10) ** -1.00000+ALGBRC(29) ** -1.00000))+ (1.00000 - CONSTS(8))*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000)))
      ALGBRC(31) =  ALGBRC(30)*STATES(3)
      ALGBRC(21) =  CONSTS(23)*( CONSTS(24)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(24))*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14))))
      ALGBRC(23) =  ALGBRC(21)*STATES(9)
      RATES(16) = (ALGBRC(1)+ALGBRC(20)+ 2.00000*ALGBRC(31)+ 2.00000*ALGBRC(16)+ALGBRC(26)+CONSTS(47)) - ( 3.00000*ALGBRC(18)+ALGBRC(23)+ 6.00000*ALGBRC(2)+ALGBRC(28))
      ALGBRC(6) =  CONSTS(19)*(ALGBRC(5)/(CONSTS(9)+ALGBRC(5)))
      ALGBRC(7) =  ALGBRC(6)*STATES(7)
      ALGBRC(32) =  CONSTS(20)*( CONSTS(21)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(21))*(ALGBRC(29) ** -1.00000/(CONSTS(10) ** -1.00000+ALGBRC(29) ** -1.00000)))
      ALGBRC(33) =  ALGBRC(32)*STATES(7)
      ALGBRC(8) =  CONSTS(22)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))
      ALGBRC(9) =  ALGBRC(8)*STATES(8)
      RATES(7) = ( 2.00000*ALGBRC(18)+ALGBRC(9)) - (ALGBRC(7)+ALGBRC(33))
      RATES(9) = (ALGBRC(33)+ALGBRC(31)) - (ALGBRC(23)+ALGBRC(13))
      RATES(13) = (ALGBRC(23)+ALGBRC(11)+ALGBRC(13)) - (ALGBRC(33)+ALGBRC(31))
      RATES(15) = (ALGBRC(7)+ 2.00000*ALGBRC(2)) - ( 2.00000*ALGBRC(18)+ALGBRC(33)+ALGBRC(9)+ (11.0000/3.00000)*ALGBRC(23)+ (35.0000/3.00000)*ALGBRC(31))
      ALGBRC(34) = TERNRY(VOI.GT.0.00000.AND.VOI.LT.5.00000, 1.00000, 0.400000)
      RATES(1) =  ALGBRC(34)*(CONSTS(34) -  CONSTS(35)*STATES(1)) - ALGBRC(1)
      RATES(3) = ( 3.00000*ALGBRC(3)+ ALGBRC(34)*(CONSTS(36) -  CONSTS(37)*STATES(3))) - (ALGBRC(16)+ALGBRC(31))
      RATES(8) = (ALGBRC(7)+ ALGBRC(34)*(CONSTS(38) -  CONSTS(39)*STATES(8))) - ALGBRC(9)
      RATES(2) =  ALGBRC(34)*(CONSTS(40) -  CONSTS(41)*STATES(2)) - ALGBRC(2)
      RATES(14) = ALGBRC(33)+ 2.00000*ALGBRC(23)+ ALGBRC(34)*(CONSTS(42) -  CONSTS(43)*STATES(14))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = CONSTS(44)/STATES(15)
      ALGBRC(10) =  CONSTS(27)*( CONSTS(28)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(28))*(STATES(13) ** -1.00000/(CONSTS(29) ** -1.00000+STATES(13) ** -1.00000)))
      ALGBRC(11) =  ALGBRC(10)*STATES(12)
      ALGBRC(12) =  CONSTS(30)*( CONSTS(31)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(31))*(STATES(13) ** -1.00000/(CONSTS(29) ** -1.00000+STATES(13) ** -1.00000)))
      ALGBRC(13) =  ALGBRC(12)*STATES(9)
      ALGBRC(14) = STATES(16)/CONSTS(33)
      ALGBRC(15) =  CONSTS(5)*(ALGBRC(14) ** -1.00000/(CONSTS(6) ** -1.00000+ALGBRC(14) ** -1.00000))
      ALGBRC(16) =  ALGBRC(15)*STATES(3)
      ALGBRC(3) =  CONSTS(46)*STATES(6)
      ALGBRC(1) =  CONSTS(1)*(STATES(1)/(CONSTS(2)+STATES(1)))
      ALGBRC(17) =  CONSTS(11)*( CONSTS(12)*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14)))+ (1.00000 - CONSTS(12))*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000)))
      ALGBRC(18) =  ALGBRC(17)*STATES(4)
      ALGBRC(4) = STATES(10)/STATES(11)
      ALGBRC(19) =  CONSTS(13)*( CONSTS(14)*(ALGBRC(14) ** -1.00000/(CONSTS(6) ** -1.00000+ALGBRC(14) ** -1.00000))+ (1.00000 - CONSTS(14))*ALGBRC(4) ** -1.00000/(CONSTS(15) ** -1.00000+ALGBRC(4) ** -1.00000) ** 2.00000)
      ALGBRC(20) =  ALGBRC(19)*STATES(4)
      ALGBRC(22) =  CONSTS(16)*( CONSTS(17)*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14)))+ (1.00000 - CONSTS(17))*ALGBRC(4)/(CONSTS(15)+ALGBRC(4)) ** 2.00000)
      ALGBRC(24) =  ALGBRC(22)*STATES(5)
      ALGBRC(25) =  CONSTS(25)*(ALGBRC(14) ** -1.00000/(CONSTS(6) ** -1.00000+ALGBRC(14) ** -1.00000))
      ALGBRC(26) =  ALGBRC(25)*STATES(10)
      ALGBRC(27) =  CONSTS(26)*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14)))
      ALGBRC(28) =  ALGBRC(27)*STATES(11)
      ALGBRC(2) =  CONSTS(3)*(STATES(2)/(CONSTS(4)+STATES(2)))
      ALGBRC(29) = STATES(9)/STATES(13)
      ALGBRC(30) =  CONSTS(7)*( CONSTS(8)*(ALGBRC(29) ** -1.00000/(CONSTS(10) ** -1.00000+ALGBRC(29) ** -1.00000))+ (1.00000 - CONSTS(8))*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000)))
      ALGBRC(31) =  ALGBRC(30)*STATES(3)
      ALGBRC(21) =  CONSTS(23)*( CONSTS(24)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(24))*(ALGBRC(14)/(CONSTS(6)+ALGBRC(14))))
      ALGBRC(23) =  ALGBRC(21)*STATES(9)
      ALGBRC(6) =  CONSTS(19)*(ALGBRC(5)/(CONSTS(9)+ALGBRC(5)))
      ALGBRC(7) =  ALGBRC(6)*STATES(7)
      ALGBRC(32) =  CONSTS(20)*( CONSTS(21)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))+ (1.00000 - CONSTS(21))*(ALGBRC(29) ** -1.00000/(CONSTS(10) ** -1.00000+ALGBRC(29) ** -1.00000)))
      ALGBRC(33) =  ALGBRC(32)*STATES(7)
      ALGBRC(8) =  CONSTS(22)*(ALGBRC(5) ** -1.00000/(CONSTS(9) ** -1.00000+ALGBRC(5) ** -1.00000))
      ALGBRC(9) =  ALGBRC(8)*STATES(8)
      ALGBRC(34) = TERNRY(VOI.GT.0.00000.AND.VOI.LT.5.00000, 1.00000, 0.400000)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END