Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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