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