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 24 entries in the algebraic variable array.
C There are a total of 33 entries in each of the rate and state variable arrays.
C There are a total of 97 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C STATES(1) is ATbl in component ATbl (nanomole).
C CONSTS(1) is k2T in component ATbl (nmol_hr).
C ALGBRC(18) is Qb in component model_parameters (l_hr).
C ALGBRC(16) is Qp in component Qp (l_hr).
C CONSTS(2) is Ql in component model_parameters (l_hr).
C CONSTS(3) is Qt in component model_parameters (l_hr).
C CONSTS(95) is CTbf in component model_parameters (nM).
C ALGBRC(1) is CTblf in component CTblf (nM).
C ALGBRC(11) is CTpf in component CTpf (nM).
C CONSTS(4) is CTbl in component model_parameters (nM).
C CONSTS(91) is CTlf in component CTlf (nM).
C CONSTS(94) is CTtv in component CTtv (nM).
C CONSTS(5) is ts in component model_parameters (dimensionless).
C STATES(2) is ADbl in component ADbl (nanomole).
C ALGBRC(19) is Qbt in component model_parameters (l_hr).
C CONSTS(96) is CDbf in component model_parameters (nM).
C ALGBRC(5) is CDblf in component CDblf (nM).
C ALGBRC(13) is CDpf in component CDpf (nM).
C CONSTS(93) is CDlf in component CDlf (nM).
C CONSTS(6) is CDbl in component model_parameters (nM).
C STATES(3) is CT_A in component CT_A (nM).
C CONSTS(7) is kTA_on in component CT_A (second_order_rate_constant).
C CONSTS(8) is kTA_off in component CT_A (first_order_rate_constant).
C ALGBRC(9) is CAf in component CAf (nM).
C STATES(4) is CD_A in component CD_A (nM).
C CONSTS(9) is kDA_on in component CD_A (second_order_rate_constant).
C CONSTS(10) is kDA_off in component CD_A (first_order_rate_constant).
C CONSTS(11) is CA in component model_parameters (nM).
C STATES(5) is ATl in component ATl (nanomole).
C CONSTS(12) is k1T in component model_parameters (first_order_rate_constant).
C CONSTS(13) is km5alpha in component model_parameters (nM).
C CONSTS(14) is Vmax1 in component model_parameters (nmol_hr).
C CONSTS(15) is PTl in component model_parameters (dimensionless).
C STATES(6) is ADl in component ADl (nanomole).
C CONSTS(16) is k1D in component model_parameters (first_order_rate_constant).
C CONSTS(17) is PDl in component model_parameters (dimensionless).
C CONSTS(18) is CTl in component model_parameters (nM).
C CONSTS(19) is CDl in component model_parameters (nM).
C STATES(7) is ATif in component ATif (nanomole).
C CONSTS(20) is mu in component model_parameters (l_hr).
C CONSTS(21) is PTif in component model_parameters (dimensionless).
C CONSTS(22) is CTif in component model_parameters (nM).
C CONSTS(23) is CTst in component model_parameters (nM).
C CONSTS(24) is CLH in component model_parameters (nanomole).
C CONSTS(25) is PTst in component model_parameters (dimensionless).
C STATES(8) is ATst in component ATst (nanomole).
C STATES(9) is ATp in component ATp (nanomole).
C CONSTS(26) is Vmax0 in component model_parameters (nmol_hr).
C STATES(10) is ADp in component ADp (nanomole).
C STATES(11) is CT_R in component CT_R (nM).
C CONSTS(27) is kTR_on in component CT_R (second_order_rate_constant).
C CONSTS(28) is kTR_off in component CT_R (first_order_rate_constant).
C ALGBRC(24) is CRf in component CRf (nM).
C STATES(12) is CD_R in component CD_R (nM).
C CONSTS(29) is kDR_on in component CD_R (second_order_rate_constant).
C CONSTS(30) is kDR_off in component CD_R (first_order_rate_constant).
C STATES(13) is CDT in component CDT (nM).
C CONSTS(31) is kDT_on in component CDT (second_order_rate_constant).
C CONSTS(32) is kDT_off in component CDT (first_order_rate_constant).
C STATES(14) is CDD in component CDD (nM).
C CONSTS(33) is kDD_on in component CDD (second_order_rate_constant).
C CONSTS(34) is kDD_off in component CDD (first_order_rate_constant).
C STATES(15) is CTT in component CTT (nM).
C CONSTS(35) is kTT_on in component CTT (second_order_rate_constant).
C CONSTS(36) is kTT_off in component CTT (first_order_rate_constant).
C STATES(16) is CDD5a in component CDD5a (nM).
C CONSTS(37) is kDNA_onDD5a in component CDD5a (second_order_rate_constant).
C CONSTS(38) is kDNA_offDD5a in component CDD5a (first_order_rate_constant).
C ALGBRC(20) is CDNA5af in component CDNA5af (nM).
C STATES(17) is CDDcd in component CDDcd (nM).
C CONSTS(39) is kDNA_onDDcd in component CDDcd (second_order_rate_constant).
C CONSTS(40) is kDNA_offDDcd in component CDDcd (first_order_rate_constant).
C ALGBRC(15) is CDNAcdf in component CDNAcdf (nM).
C STATES(18) is CDDsec in component CDDsec (nM).
C CONSTS(41) is kDNA_onDDsec in component CDDsec (second_order_rate_constant).
C CONSTS(42) is kDNA_offDDsec in component CDDsec (first_order_rate_constant).
C ALGBRC(17) is CDNAsecf in component CDNAsecf (nM).
C STATES(19) is CDDcp in component CDDcp (nM).
C CONSTS(43) is kDNA_onDDcp in component CDDcp (second_order_rate_constant).
C CONSTS(44) is kDNA_offDDcp in component CDDcp (first_order_rate_constant).
C ALGBRC(21) is CDNAcpf in component CDNAcpf (nM).
C STATES(20) is CDT5a in component CDT5a (nM).
C CONSTS(45) is kDNA_onDT5a in component CDT5a (second_order_rate_constant).
C CONSTS(46) is kDNA_offDT5a in component CDT5a (first_order_rate_constant).
C STATES(21) is CDTcd in component CDTcd (nM).
C CONSTS(47) is kDNA_onDTcd in component CDTcd (second_order_rate_constant).
C CONSTS(48) is kDNA_offDTcd in component CDTcd (first_order_rate_constant).
C STATES(22) is CDTsec in component CDTsec (nM).
C CONSTS(49) is kDNA_onDTsec in component CDTsec (second_order_rate_constant).
C CONSTS(50) is kDNA_offDTsec in component CDTsec (first_order_rate_constant).
C STATES(23) is CDTcp in component CDTcp (nM).
C CONSTS(51) is kDNA_onDTcp in component CDTcp (second_order_rate_constant).
C CONSTS(52) is kDNA_offDTcp in component CDTcp (first_order_rate_constant).
C STATES(24) is CTT5a in component CTT5a (nM).
C CONSTS(53) is kDNA_onTT5a in component CTT5a (second_order_rate_constant).
C CONSTS(54) is kDNA_offTT5a in component CTT5a (first_order_rate_constant).
C STATES(25) is CTTcd in component CTTcd (nM).
C CONSTS(55) is kDNA_onTTcd in component CTTcd (second_order_rate_constant).
C CONSTS(56) is kDNA_offTTcd in component CTTcd (first_order_rate_constant).
C STATES(26) is CTTsec in component CTTsec (nM).
C CONSTS(57) is kDNA_onTTsec in component CTTsec (second_order_rate_constant).
C CONSTS(58) is kDNA_offTTsec in component CTTsec (first_order_rate_constant).
C STATES(27) is CTTcp in component CTTcp (nM).
C CONSTS(59) is kDNA_onTTcp in component CTTcp (second_order_rate_constant).
C CONSTS(60) is kDNA_offTTcp in component CTTcp (first_order_rate_constant).
C ALGBRC(8) is T in component CTpf (nM).
C CONSTS(61) is NTp in component CTpf (dimensionless).
C CONSTS(62) is CTp in component model_parameters (nM).
C ALGBRC(12) is D in component CDpf (nM).
C CONSTS(63) is NDp in component CDpf (dimensionless).
C CONSTS(64) is CDp in component model_parameters (nM).
C ALGBRC(23) is CR in component CR (nM).
C CONSTS(65) is CDNAcd in component model_parameters (nM).
C CONSTS(66) is CDNAsec in component model_parameters (nM).
C CONSTS(67) is CDNA5a in component model_parameters (nM).
C CONSTS(68) is CDNAcp in component model_parameters (nM).
C ALGBRC(2) is DNAocp in component DNAocp (dimensionless).
C ALGBRC(3) is DNAosec in component DNAosec (dimensionless).
C ALGBRC(6) is DNAocd in component DNAocd (dimensionless).
C ALGBRC(4) is DNAo5a in component DNAo5a (dimensionless).
C STATES(28) is ATb in component ATb (nanomole).
C CONSTS(69) is CTb in component model_parameters (nM).
C CONSTS(70) is PTb in component model_parameters (dimensionless).
C STATES(29) is ADb in component ADb (nanomole).
C CONSTS(71) is CDb in component model_parameters (nM).
C CONSTS(72) is PDb in component model_parameters (dimensionless).
C STATES(30) is ALH in component ALH (nanomole).
C CONSTS(73) is kLH1 in component ALH (l_hr).
C CONSTS(74) is kLH2 in component ALH (nmol_hr).
C CONSTS(75) is kLH3 in component ALH (nmol_hr).
C CONSTS(76) is kLH4 in component ALH (first_order_rate_constant).
C CONSTS(77) is alpha_T in component ALH (dimensionless).
C CONSTS(78) is alpha_D in component ALH (dimensionless).
C ALGBRC(10) is A in component ALH (nmol_hr).
C CONSTS(92) is B in component ALH (nmol_hr).
C STATES(31) is AR in component AR (nanomole).
C CONSTS(79) is k1R in component AR (nmol_hr).
C CONSTS(80) is keR in component AR (l_hr).
C ALGBRC(22) is Vp in component Vp (litre).
C CONSTS(81) is kQp in component Qp (first_order_rate_constant).
C ALGBRC(14) is Vpc in component Qp (litre).
C STATES(32) is VPC1 in component VPC1 (litre).
C CONSTS(82) is VPC2 in component model_parameters (litre).
C ALGBRC(7) is Vmax in component Vmax (nmol_hr).
C CONSTS(83) is k5a in component Vmax (nmol_hr).
C CONSTS(84) is kcpl in component VPC1 (kg_hr).
C CONSTS(85) is kcd1 in component VPC1 (first_order_rate_constant).
C CONSTS(86) is VPC1b in component model_parameters (kg).
C STATES(33) is VPL1 in component VPL1 (litre).
C CONSTS(87) is ksec in component VPL1 (kg_hr).
C CONSTS(88) is kflo in component VPL1 (first_order_rate_constant).
C CONSTS(89) is VPL2 in component model_parameters (litre).
C CONSTS(90) is Qc in component model_parameters (l_hr).
C RATES(1) is d/dt ATbl in component ATbl (nanomole).
C RATES(2) is d/dt ADbl in component ADbl (nanomole).
C RATES(3) is d/dt CT_A in component CT_A (nM).
C RATES(4) is d/dt CD_A in component CD_A (nM).
C RATES(5) is d/dt ATl in component ATl (nanomole).
C RATES(6) is d/dt ADl in component ADl (nanomole).
C RATES(7) is d/dt ATif in component ATif (nanomole).
C RATES(8) is d/dt ATst in component ATst (nanomole).
C RATES(9) is d/dt ATp in component ATp (nanomole).
C RATES(10) is d/dt ADp in component ADp (nanomole).
C RATES(11) is d/dt CT_R in component CT_R (nM).
C RATES(12) is d/dt CD_R in component CD_R (nM).
C RATES(13) is d/dt CDT in component CDT (nM).
C RATES(14) is d/dt CDD in component CDD (nM).
C RATES(15) is d/dt CTT in component CTT (nM).
C RATES(16) is d/dt CDD5a in component CDD5a (nM).
C RATES(17) is d/dt CDDcd in component CDDcd (nM).
C RATES(18) is d/dt CDDsec in component CDDsec (nM).
C RATES(19) is d/dt CDDcp in component CDDcp (nM).
C RATES(20) is d/dt CDT5a in component CDT5a (nM).
C RATES(21) is d/dt CDTcd in component CDTcd (nM).
C RATES(22) is d/dt CDTsec in component CDTsec (nM).
C RATES(23) is d/dt CDTcp in component CDTcp (nM).
C RATES(24) is d/dt CTT5a in component CTT5a (nM).
C RATES(25) is d/dt CTTcd in component CTTcd (nM).
C RATES(26) is d/dt CTTsec in component CTTsec (nM).
C RATES(27) is d/dt CTTcp in component CTTcp (nM).
C RATES(28) is d/dt ATb in component ATb (nanomole).
C RATES(29) is d/dt ADb in component ADb (nanomole).
C RATES(30) is d/dt ALH in component ALH (nanomole).
C RATES(31) is d/dt AR in component AR (nanomole).
C RATES(32) is d/dt VPC1 in component VPC1 (litre).
C RATES(33) is d/dt VPL1 in component VPL1 (litre).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0.0
CONSTS(1) = 0.17
CONSTS(2) = 1.06
CONSTS(3) = 0.061
CONSTS(4) = 7.60
CONSTS(5) = 0.56
STATES(2) = 0.0
CONSTS(6) = 0.59
STATES(3) = 0.0
CONSTS(7) = 0.18
CONSTS(8) = 8131.0
STATES(4) = 0.0
CONSTS(9) = 0.15
CONSTS(10) = 5407.0
CONSTS(11) = 5.00E5
STATES(5) = 0.0
CONSTS(12) = 87.93
CONSTS(13) = 40.00
CONSTS(14) = 3.65
CONSTS(15) = 2.75
STATES(6) = 0.0
CONSTS(16) = 77.20
CONSTS(17) = 2.00
CONSTS(18) = 9.92
CONSTS(19) = 0.92
STATES(7) = 0.0
CONSTS(20) = 1000.0
CONSTS(21) = 1.96
CONSTS(22) = 307.29
CONSTS(23) = 307.29
CONSTS(24) = 0.104
CONSTS(25) = 1.00
STATES(8) = 0.0
STATES(9) = 0.0
CONSTS(26) = 2.89
STATES(10) = 0.0
STATES(11) = 0.0
CONSTS(27) = 0.14
CONSTS(28) = 0.069
STATES(12) = 0.0
CONSTS(29) = 0.053
CONSTS(30) = 0.018
STATES(13) = 0.0
CONSTS(31) = 0.14
CONSTS(32) = 3.13
STATES(14) = 0.0
CONSTS(33) = 0.14
CONSTS(34) = 3.13
STATES(15) = 0.0
CONSTS(35) = 0.14
CONSTS(36) = 3.13
STATES(16) = 0.0
CONSTS(37) = 0.14
CONSTS(38) = 0.07
STATES(17) = 0.0
CONSTS(39) = 0.14
CONSTS(40) = 0.1
STATES(18) = 0.0
CONSTS(41) = 0.14
CONSTS(42) = 0.2
STATES(19) = 0.0
CONSTS(43) = 0.14
CONSTS(44) = 0.7
STATES(20) = 0.0
CONSTS(45) = 0.14
CONSTS(46) = 0.56
STATES(21) = 0.0
CONSTS(47) = 0.14
CONSTS(48) = 0.8
STATES(22) = 0.0
CONSTS(49) = 0.14
CONSTS(50) = 1.6
STATES(23) = 0.0
CONSTS(51) = 0.14
CONSTS(52) = 5.6
STATES(24) = 0.0
CONSTS(53) = 0.14
CONSTS(54) = 0.84
STATES(25) = 0.0
CONSTS(55) = 0.14
CONSTS(56) = 1.2
STATES(26) = 0.0
CONSTS(57) = 0.14
CONSTS(58) = 2.4
STATES(27) = 0.0
CONSTS(59) = 0.14
CONSTS(60) = 8.4
CONSTS(61) = 0.79
CONSTS(62) = 11.38
CONSTS(63) = 1.00
CONSTS(64) = 52.46
CONSTS(65) = 0.075
CONSTS(66) = 0.075
CONSTS(67) = 0.075
CONSTS(68) = 0.075
STATES(28) = 0.0
CONSTS(69) = 0.6275
CONSTS(70) = 1.00
STATES(29) = 0.0
CONSTS(71) = 0.019884
CONSTS(72) = 0.62
STATES(30) = 0.0
CONSTS(73) = 0.13
CONSTS(74) = 0.026
CONSTS(75) = 1.68E-4
CONSTS(76) = 0.80
CONSTS(77) = 0.25
CONSTS(78) = 0.75
STATES(31) = 0.0
CONSTS(79) = 68.15
CONSTS(80) = 71.90
CONSTS(81) = 73.84
STATES(32) = 0.0
CONSTS(82) = 5.08E-5
CONSTS(83) = 3.02
CONSTS(84) = 1.17E-7
CONSTS(85) = 0.014
CONSTS(86) = 1E-4
STATES(33) = 0.0
CONSTS(87) = 8.44E-6
CONSTS(88) = 0.043
CONSTS(89) = 3.08E-6
CONSTS(90) = 6.08
CONSTS(91) = CONSTS(18)/CONSTS(15)
CONSTS(92) = CONSTS(75)
CONSTS(96) = CONSTS(20)*(CONSTS(22) - CONSTS(23)/CONSTS(25))
CONSTS(93) = CONSTS(19)/CONSTS(17)
CONSTS(94) = CONSTS(22)/CONSTS(21)
CONSTS(95) = CONSTS(69)/CONSTS(70)
CONSTS(96) = CONSTS(71)/CONSTS(72)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(8) = CONSTS(96)
RATES(5) = CONSTS(2)*(CONSTS(4) - CONSTS(91)) - (( CONSTS(14)*CONSTS(91))/(CONSTS(13)+CONSTS(91))+ CONSTS(12)*(STATES(5)/CONSTS(15)))
RATES(6) = CONSTS(2)*(CONSTS(6) - CONSTS(93)) - (( CONSTS(14)*CONSTS(91))/(CONSTS(13)+CONSTS(91))+ CONSTS(16)*(STATES(6)/CONSTS(17)))
RATES(13) = CONSTS(31)*STATES(12)*STATES(11) - CONSTS(32)*STATES(13)
RATES(14) = CONSTS(33)*STATES(12) ** 2.00000 - CONSTS(34)*STATES(14)
RATES(15) = CONSTS(35)*STATES(11) ** 2.00000 - CONSTS(36)*STATES(15)
ALGBRC(1) = CONSTS(4) - STATES(3)
RATES(7) = CONSTS(3)*(1.00000 - CONSTS(5))*(ALGBRC(1) - CONSTS(22)/CONSTS(21))+ CONSTS(20)*(CONSTS(23)/CONSTS(25) - CONSTS(22))+ CONSTS(12)*CONSTS(24)
ALGBRC(3) = (STATES(26)+STATES(18)+STATES(22))/CONSTS(66)
RATES(33) = (STATES(32)/CONSTS(86))*CONSTS(87)*ALGBRC(3) - CONSTS(88)*STATES(33)
ALGBRC(2) = (STATES(27)+STATES(19)+STATES(23))/CONSTS(68)
ALGBRC(6) = (STATES(25)+STATES(17)+STATES(21))/CONSTS(65)
RATES(32) = (STATES(32)/CONSTS(86))*CONSTS(84)*ALGBRC(2) - CONSTS(85)*(1.00000 - ALGBRC(6))*STATES(32)
ALGBRC(9) = CONSTS(11) - (STATES(3)+STATES(4))
RATES(3) = CONSTS(7)*ALGBRC(1)*ALGBRC(9) - CONSTS(8)*STATES(3)
ALGBRC(5) = CONSTS(6) - STATES(4)
RATES(4) = CONSTS(9)*ALGBRC(5)*ALGBRC(9) - CONSTS(10)*STATES(4)
ALGBRC(10) = - CONSTS(73)*( CONSTS(77)*ALGBRC(1)+ CONSTS(78)*ALGBRC(5))+CONSTS(74)
RATES(30) = TERNRY(ALGBRC(10).GT.CONSTS(92), ALGBRC(10) - CONSTS(76)*STATES(30), CONSTS(92) - CONSTS(76)*STATES(30))
ALGBRC(15) = CONSTS(65) - (STATES(25)+STATES(17)+STATES(21))
RATES(17) = CONSTS(39)*STATES(14)*ALGBRC(15) - CONSTS(40)*STATES(17)
RATES(21) = CONSTS(47)*STATES(13)*ALGBRC(15) - CONSTS(48)*STATES(21)
RATES(25) = CONSTS(55)*STATES(15)*ALGBRC(15) - CONSTS(56)*STATES(25)
ALGBRC(14) = STATES(32)+CONSTS(82)
ALGBRC(16) = CONSTS(81)*ALGBRC(14)
ALGBRC(8) = STATES(26)+STATES(22)+STATES(24)+STATES(20)+STATES(25)+STATES(21)+STATES(27)+STATES(23)
ALGBRC(11) = (CONSTS(62) - (STATES(11)+ 2.00000*STATES(15)+STATES(13)+ALGBRC(8)))/(1.00000+CONSTS(61))
RATES(9) = ALGBRC(16)*(CONSTS(4) - ALGBRC(11)) - ( CONSTS(26)*ALGBRC(11))/(CONSTS(13)+ALGBRC(11))
ALGBRC(12) = STATES(18)+STATES(22)+STATES(16)+STATES(20)+STATES(17)+STATES(21)+STATES(19)+STATES(23)
ALGBRC(13) = (CONSTS(64) - (STATES(12)+ 2.00000*STATES(14)+STATES(13)+ALGBRC(12)))/(1.00000+CONSTS(63))
RATES(10) = ALGBRC(16)*(CONSTS(6) - ALGBRC(13)) - ( CONSTS(26)*ALGBRC(11))/(CONSTS(13)+ALGBRC(11))
ALGBRC(17) = CONSTS(66) - (STATES(26)+STATES(18)+STATES(22))
RATES(18) = CONSTS(41)*STATES(14)*ALGBRC(17) - CONSTS(42)*STATES(18)
RATES(22) = CONSTS(49)*STATES(13)*ALGBRC(17) - CONSTS(50)*STATES(22)
RATES(26) = CONSTS(57)*STATES(15)*ALGBRC(17) - CONSTS(58)*STATES(26)
ALGBRC(18) = CONSTS(90) - (ALGBRC(16)+CONSTS(3)+CONSTS(2))
RATES(1) = ALGBRC(18)*(CONSTS(95) - ALGBRC(1))+ ALGBRC(16)*(ALGBRC(11) - CONSTS(4))+ CONSTS(2)*(CONSTS(91) - CONSTS(4))+ CONSTS(3)*(1.00000 - CONSTS(5))*(CONSTS(94) - ALGBRC(1))+CONSTS(1)
ALGBRC(19) = CONSTS(90) - (ALGBRC(16)+CONSTS(2))
RATES(2) = ALGBRC(19)*(CONSTS(96) - ALGBRC(5))+ ALGBRC(16)*(ALGBRC(13) - ALGBRC(5))+ CONSTS(2)*(CONSTS(93) - CONSTS(6))
ALGBRC(20) = CONSTS(67) - (STATES(24)+STATES(16)+STATES(20))
RATES(16) = CONSTS(37)*STATES(14)*ALGBRC(20) - CONSTS(38)*STATES(16)
RATES(20) = CONSTS(45)*STATES(13)*ALGBRC(20) - CONSTS(46)*STATES(20)
RATES(24) = CONSTS(53)*STATES(15)*ALGBRC(20) - CONSTS(54)*STATES(24)
RATES(28) = ALGBRC(18)*(ALGBRC(1) - CONSTS(69)/CONSTS(70))
RATES(29) = ALGBRC(19)*(ALGBRC(5) - CONSTS(71)/CONSTS(72))
ALGBRC(21) = CONSTS(68) - (STATES(27)+STATES(19)+STATES(23))
RATES(19) = CONSTS(43)*STATES(14)*ALGBRC(21) - CONSTS(44)*STATES(19)
RATES(23) = CONSTS(51)*STATES(13)*ALGBRC(21) - CONSTS(52)*STATES(23)
RATES(27) = CONSTS(59)*STATES(15)*ALGBRC(21) - CONSTS(60)*STATES(27)
ALGBRC(22) = STATES(32)+CONSTS(82)+STATES(33)+CONSTS(89)
ALGBRC(23) = STATES(31)/ALGBRC(22)
ALGBRC(24) = ALGBRC(23) - (STATES(11)+STATES(12)+ 2.00000*(STATES(15)+STATES(14)+STATES(13))+ 2.00000*(ALGBRC(15)+ALGBRC(17)+ALGBRC(21)+ALGBRC(20)))
RATES(11) = CONSTS(27)*ALGBRC(11)*ALGBRC(24) - CONSTS(28)*STATES(11)
RATES(12) = CONSTS(29)*ALGBRC(13)*ALGBRC(24) - CONSTS(30)*STATES(12)
RATES(31) = CONSTS(79) - CONSTS(80)*ALGBRC(24)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = CONSTS(4) - STATES(3)
ALGBRC(3) = (STATES(26)+STATES(18)+STATES(22))/CONSTS(66)
ALGBRC(2) = (STATES(27)+STATES(19)+STATES(23))/CONSTS(68)
ALGBRC(6) = (STATES(25)+STATES(17)+STATES(21))/CONSTS(65)
ALGBRC(9) = CONSTS(11) - (STATES(3)+STATES(4))
ALGBRC(5) = CONSTS(6) - STATES(4)
ALGBRC(10) = - CONSTS(73)*( CONSTS(77)*ALGBRC(1)+ CONSTS(78)*ALGBRC(5))+CONSTS(74)
ALGBRC(15) = CONSTS(65) - (STATES(25)+STATES(17)+STATES(21))
ALGBRC(14) = STATES(32)+CONSTS(82)
ALGBRC(16) = CONSTS(81)*ALGBRC(14)
ALGBRC(8) = STATES(26)+STATES(22)+STATES(24)+STATES(20)+STATES(25)+STATES(21)+STATES(27)+STATES(23)
ALGBRC(11) = (CONSTS(62) - (STATES(11)+ 2.00000*STATES(15)+STATES(13)+ALGBRC(8)))/(1.00000+CONSTS(61))
ALGBRC(12) = STATES(18)+STATES(22)+STATES(16)+STATES(20)+STATES(17)+STATES(21)+STATES(19)+STATES(23)
ALGBRC(13) = (CONSTS(64) - (STATES(12)+ 2.00000*STATES(14)+STATES(13)+ALGBRC(12)))/(1.00000+CONSTS(63))
ALGBRC(17) = CONSTS(66) - (STATES(26)+STATES(18)+STATES(22))
ALGBRC(18) = CONSTS(90) - (ALGBRC(16)+CONSTS(3)+CONSTS(2))
ALGBRC(19) = CONSTS(90) - (ALGBRC(16)+CONSTS(2))
ALGBRC(20) = CONSTS(67) - (STATES(24)+STATES(16)+STATES(20))
ALGBRC(21) = CONSTS(68) - (STATES(27)+STATES(19)+STATES(23))
ALGBRC(22) = STATES(32)+CONSTS(82)+STATES(33)+CONSTS(89)
ALGBRC(23) = STATES(31)/ALGBRC(22)
ALGBRC(24) = ALGBRC(23) - (STATES(11)+STATES(12)+ 2.00000*(STATES(15)+STATES(14)+STATES(13))+ 2.00000*(ALGBRC(15)+ALGBRC(17)+ALGBRC(21)+ALGBRC(20)))
ALGBRC(4) = (STATES(24)+STATES(16)+STATES(20))/CONSTS(67)
ALGBRC(7) = CONSTS(83)*ALGBRC(4)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END