Generated Code

C
C There are a total of 29 entries in the algebraic variable array.
C There are a total of 29 entries in each of the rate and state variable arrays.
C There are a total of 67 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is Rs in component Rs (dimensionless).
C ALGBRC(1) is v1 in component v1 (first_order_rate_constant).
C ALGBRC(8) is v3 in component v3 (first_order_rate_constant).
C STATES(2) is RL in component RL (dimensionless).
C ALGBRC(5) is v2 in component v2 (first_order_rate_constant).
C ALGBRC(11) is v4 in component v4 (first_order_rate_constant).
C STATES(3) is Ri in component Ri (dimensionless).
C ALGBRC(16) is v6 in component v6 (first_order_rate_constant).
C STATES(4) is L in component L (dimensionless).
C STATES(5) is R2L2 in component R2L2 (dimensionless).
C ALGBRC(14) is v5 in component v5 (first_order_rate_constant).
C ALGBRC(19) is v7 in component v7 (first_order_rate_constant).
C STATES(6) is R2_CPP in component R2_CPP (dimensionless).
C ALGBRC(23) is v8 in component v8 (first_order_rate_constant).
C STATES(7) is Li in component Li (dimensionless).
C STATES(8) is R2i in component R2i (dimensionless).
C STATES(9) is Shc in component Shc (dimensionless).
C ALGBRC(2) is v9 in component v9 (first_order_rate_constant).
C ALGBRC(6) is v10 in component v10 (first_order_rate_constant).
C STATES(10) is ShcP in component ShcP (dimensionless).
C ALGBRC(9) is v11 in component v11 (first_order_rate_constant).
C ALGBRC(17) is v27 in component v27 (first_order_rate_constant).
C STATES(11) is ShcGS in component ShcGS (dimensionless).
C ALGBRC(15) is v13 in component v13 (first_order_rate_constant).
C ALGBRC(12) is v12 in component v12 (first_order_rate_constant).
C STATES(12) is GS in component GS (dimensionless).
C ALGBRC(20) is v28 in component v28 (first_order_rate_constant).
C STATES(13) is GSP in component GSP (dimensionless).
C STATES(14) is RasGDP in component RasGDP (dimensionless).
C ALGBRC(21) is v15 in component v15 (first_order_rate_constant).
C STATES(15) is Ras_ShcGS in component Ras_ShcGS (dimensionless).
C STATES(16) is RasGTP in component RasGTP (dimensionless).
C ALGBRC(24) is v17 in component v17 (first_order_rate_constant).
C ALGBRC(18) is v14 in component v14 (first_order_rate_constant).
C ALGBRC(22) is v16 in component v16 (first_order_rate_constant).
C STATES(17) is GAP in component GAP (dimensionless).
C STATES(18) is Ras_GAP in component Ras_GAP (dimensionless).
C STATES(19) is Raf in component Raf (dimensionless).
C ALGBRC(25) is v18 in component v18 (first_order_rate_constant).
C STATES(20) is Ras_Raf in component Ras_Raf (dimensionless).
C STATES(21) is Rafa in component Rafa (dimensionless).
C ALGBRC(26) is v19 in component v19 (first_order_rate_constant).
C ALGBRC(28) is v21 in component v21 (first_order_rate_constant).
C STATES(22) is MEK in component MEK (dimensionless).
C ALGBRC(27) is v20 in component v20 (first_order_rate_constant).
C STATES(23) is MEKP in component MEKP (dimensionless).
C ALGBRC(29) is v22 in component v22 (first_order_rate_constant).
C STATES(24) is MEKPP in component MEKPP (dimensionless).
C STATES(25) is ERK in component ERK (dimensionless).
C ALGBRC(7) is v24 in component v24 (first_order_rate_constant).
C ALGBRC(3) is v23 in component v23 (first_order_rate_constant).
C STATES(26) is ERKP in component ERKP (dimensionless).
C ALGBRC(13) is v26 in component v26 (first_order_rate_constant).
C ALGBRC(10) is v25 in component v25 (first_order_rate_constant).
C STATES(27) is ERKPP in component ERKPP (dimensionless).
C STATES(28) is t in component t (dimensionless).
C CONSTS(66) is v29 in component v29 (first_order_rate_constant).
C STATES(29) is X in component X (dimensionless).
C ALGBRC(4) is v30 in component v30 (first_order_rate_constant).
C CONSTS(1) is k1 in component v1 (first_order_rate_constant).
C CONSTS(2) is kn1 in component v1 (first_order_rate_constant).
C CONSTS(3) is DT in component v2 (dimensionless).
C CONSTS(4) is E in component v2 (dimensionless).
C CONSTS(5) is k2 in component v2 (first_order_rate_constant).
C CONSTS(6) is f in component v2 (dimensionless).
C CONSTS(7) is DT in component v3 (dimensionless).
C CONSTS(8) is E in component v3 (dimensionless).
C CONSTS(9) is kn3 in component v3 (first_order_rate_constant).
C CONSTS(10) is f in component v3 (dimensionless).
C CONSTS(11) is k3 in component v3 (first_order_rate_constant).
C CONSTS(12) is k2_4 in component v4 (first_order_rate_constant).
C CONSTS(13) is k4 in component v4 (first_order_rate_constant).
C CONSTS(14) is DT in component v5 (dimensionless).
C CONSTS(15) is E in component v5 (dimensionless).
C CONSTS(16) is k5 in component v5 (first_order_rate_constant).
C CONSTS(17) is f in component v5 (dimensionless).
C CONSTS(18) is DT in component v6 (dimensionless).
C CONSTS(19) is E in component v6 (dimensionless).
C CONSTS(20) is k6 in component v6 (first_order_rate_constant).
C CONSTS(21) is k7 in component v7 (first_order_rate_constant).
C CONSTS(22) is f in component v7 (dimensionless).
C CONSTS(23) is kn7 in component v7 (first_order_rate_constant).
C CONSTS(24) is DT in component v8 (dimensionless).
C CONSTS(25) is E in component v8 (dimensionless).
C CONSTS(26) is k8 in component v8 (first_order_rate_constant).
C CONSTS(27) is k9 in component v9 (first_order_rate_constant).
C CONSTS(28) is K_9 in component v9 (dimensionless).
C CONSTS(29) is V_10 in component v10 (first_order_rate_constant).
C CONSTS(30) is K_10 in component v10 (dimensionless).
C CONSTS(31) is k11 in component v11 (first_order_rate_constant).
C CONSTS(32) is kn11 in component v11 (first_order_rate_constant).
C CONSTS(33) is k12 in component v12 (first_order_rate_constant).
C CONSTS(34) is kn12 in component v12 (first_order_rate_constant).
C CONSTS(35) is k_13 in component v13 (first_order_rate_constant).
C CONSTS(36) is k14 in component v14 (first_order_rate_constant).
C CONSTS(37) is kn14 in component v14 (first_order_rate_constant).
C CONSTS(38) is k15 in component v15 (first_order_rate_constant).
C CONSTS(39) is k16 in component v16 (first_order_rate_constant).
C CONSTS(40) is kn16 in component v16 (first_order_rate_constant).
C CONSTS(41) is k17 in component v17 (first_order_rate_constant).
C CONSTS(42) is V_18 in component v18 (first_order_rate_constant).
C CONSTS(43) is K_18 in component v18 (dimensionless).
C CONSTS(44) is k19 in component v19 (first_order_rate_constant).
C CONSTS(45) is K_19 in component v19 (dimensionless).
C CONSTS(46) is V_20 in component v20 (first_order_rate_constant).
C CONSTS(47) is K_20 in component v20 (dimensionless).
C CONSTS(48) is k21 in component v21 (first_order_rate_constant).
C CONSTS(49) is K_21 in component v21 (dimensionless).
C CONSTS(50) is V_22 in component v22 (first_order_rate_constant).
C CONSTS(51) is K_22 in component v22 (dimensionless).
C CONSTS(52) is k23 in component v23 (first_order_rate_constant).
C CONSTS(53) is K_23 in component v23 (dimensionless).
C CONSTS(54) is V_24 in component v24 (first_order_rate_constant).
C CONSTS(55) is K_24 in component v24 (dimensionless).
C CONSTS(56) is k25 in component v25 (first_order_rate_constant).
C CONSTS(57) is K_25 in component v25 (dimensionless).
C CONSTS(58) is V_26 in component v26 (first_order_rate_constant).
C CONSTS(59) is K_26 in component v26 (dimensionless).
C CONSTS(60) is k27 in component v27 (first_order_rate_constant).
C CONSTS(61) is K_27 in component v27 (dimensionless).
C CONSTS(62) is V_28 in component v28 (first_order_rate_constant).
C CONSTS(63) is K_28 in component v28 (dimensionless).
C CONSTS(64) is v_1 in component v29 (first_order_rate_constant).
C CONSTS(65) is k_11 in component v30 (first_order_rate_constant).
C RATES(1) is d/dt Rs in component Rs (dimensionless).
C RATES(2) is d/dt RL in component RL (dimensionless).
C RATES(3) is d/dt Ri in component Ri (dimensionless).
C RATES(4) is d/dt L in component L (dimensionless).
C RATES(5) is d/dt R2L2 in component R2L2 (dimensionless).
C RATES(6) is d/dt R2_CPP in component R2_CPP (dimensionless).
C RATES(7) is d/dt Li in component Li (dimensionless).
C RATES(8) is d/dt R2i in component R2i (dimensionless).
C RATES(9) is d/dt Shc in component Shc (dimensionless).
C RATES(10) is d/dt ShcP in component ShcP (dimensionless).
C RATES(11) is d/dt ShcGS in component ShcGS (dimensionless).
C RATES(12) is d/dt GS in component GS (dimensionless).
C RATES(13) is d/dt GSP in component GSP (dimensionless).
C RATES(14) is d/dt RasGDP in component RasGDP (dimensionless).
C RATES(15) is d/dt Ras_ShcGS in component Ras_ShcGS (dimensionless).
C RATES(16) is d/dt RasGTP in component RasGTP (dimensionless).
C RATES(17) is d/dt GAP in component GAP (dimensionless).
C RATES(18) is d/dt Ras_GAP in component Ras_GAP (dimensionless).
C RATES(19) is d/dt Raf in component Raf (dimensionless).
C RATES(20) is d/dt Ras_Raf in component Ras_Raf (dimensionless).
C RATES(21) is d/dt Rafa in component Rafa (dimensionless).
C RATES(22) is d/dt MEK in component MEK (dimensionless).
C RATES(23) is d/dt MEKP in component MEKP (dimensionless).
C RATES(24) is d/dt MEKPP in component MEKPP (dimensionless).
C RATES(25) is d/dt ERK in component ERK (dimensionless).
C RATES(26) is d/dt ERKP in component ERKP (dimensionless).
C RATES(27) is d/dt ERKPP in component ERKPP (dimensionless).
C RATES(28) is d/dt t in component t (dimensionless).
C RATES(29) is d/dt X in component X (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 11100
      STATES(2) = 0
      STATES(3) = 3900
      STATES(4) = 0.0000001
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      STATES(9) = 30000
      STATES(10) = 0
      STATES(11) = 0
      STATES(12) = 20000
      STATES(13) = 0
      STATES(14) = 19800
      STATES(15) = 0
      STATES(16) = 200
      STATES(17) = 15000
      STATES(18) = 0
      STATES(19) = 10000
      STATES(20) = 0
      STATES(21) = 0
      STATES(22) = 360000
      STATES(23) = 0
      STATES(24) = 0
      STATES(25) = 750000
      STATES(26) = 0
      STATES(27) = 0
      STATES(28) = 0
      STATES(29) = 0
      CONSTS(1) = 384210000
      CONSTS(2) = 0.73
      CONSTS(3) = 6.5
      CONSTS(4) = 0.12
      CONSTS(5) = 0.7
      CONSTS(6) = 0.2
      CONSTS(7) = 6.5
      CONSTS(8) = 0.12
      CONSTS(9) = 0.7
      CONSTS(10) = 0.2
      CONSTS(11) = 0.0484
      CONSTS(12) = 0.0000001
      CONSTS(13) = 0.001383
      CONSTS(14) = 6.5
      CONSTS(15) = 0.12
      CONSTS(16) = 0.35
      CONSTS(17) = 0.2
      CONSTS(18) = 6.5
      CONSTS(19) = 0.12
      CONSTS(20) = 0.35
      CONSTS(21) = 1
      CONSTS(22) = 0.2
      CONSTS(23) = 0.000347
      CONSTS(24) = 6.5
      CONSTS(25) = 0.12
      CONSTS(26) = 0.35
      CONSTS(27) = 12
      CONSTS(28) = 6000
      CONSTS(29) = 300000
      CONSTS(30) = 6000
      CONSTS(31) = 0.002
      CONSTS(32) = 3.81
      CONSTS(33) = 0.0163
      CONSTS(34) = 10
      CONSTS(35) = 15
      CONSTS(36) = 0.005
      CONSTS(37) = 60
      CONSTS(38) = 720
      CONSTS(39) = 0.0012
      CONSTS(40) = 3
      CONSTS(41) = 27
      CONSTS(42) = 97000
      CONSTS(43) = 6000
      CONSTS(44) = 50
      CONSTS(45) = 9000
      CONSTS(46) = 920000
      CONSTS(47) = 600000
      CONSTS(48) = 50
      CONSTS(49) = 9000
      CONSTS(50) = 920000
      CONSTS(51) = 600000
      CONSTS(52) = 8.3
      CONSTS(53) = 90000
      CONSTS(54) = 200000
      CONSTS(55) = 600000
      CONSTS(56) = 8.3
      CONSTS(57) = 90000
      CONSTS(58) = 400000
      CONSTS(59) = 600000
      CONSTS(60) = 1.6
      CONSTS(61) = 600000
      CONSTS(62) = 75
      CONSTS(63) = 20000
      CONSTS(64) = 1
      CONSTS(65) = 0
      CONSTS(66) = CONSTS(64)
      CONSTS(66) = CONSTS(66)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(28) = CONSTS(66)
      ALGBRC(1) =  CONSTS(1)*STATES(1)*STATES(4) -  CONSTS(2)*STATES(2)
      RATES(4) = - ALGBRC(1)
      ALGBRC(4) =  CONSTS(65)*STATES(28)
      RATES(29) = ALGBRC(4)
      ALGBRC(2) = ( CONSTS(27)*2.00000*(STATES(5)+STATES(8)+STATES(6))*STATES(9))/(CONSTS(28)+STATES(9))
      ALGBRC(6) = ( CONSTS(29)*STATES(10))/(CONSTS(30)+STATES(10))
      RATES(9) = ALGBRC(6) - ALGBRC(2)
      ALGBRC(7) = ( CONSTS(54)*STATES(26))/(CONSTS(55)+STATES(26))
      ALGBRC(3) = ( CONSTS(52)*STATES(25)*(STATES(23)+STATES(24)))/(CONSTS(53)+STATES(25))
      RATES(25) = ALGBRC(7) - ALGBRC(3)
      ALGBRC(8) =  CONSTS(11)*STATES(3) -  CONSTS(10)*CONSTS(9)*(CONSTS(8)+ (1.00000 - CONSTS(8))*(1.00000 - EXP(- STATES(28)/CONSTS(7) ** 3.00000)))*STATES(1)
      RATES(1) = ALGBRC(8) - ALGBRC(1)
      ALGBRC(5) =  CONSTS(6)*CONSTS(5)*(CONSTS(4)+ (1.00000 - CONSTS(4))*(1.00000 - EXP(- STATES(28)/CONSTS(3) ** 3.00000)))*STATES(2)
      ALGBRC(11) =  CONSTS(13)*STATES(2)*STATES(2) -  CONSTS(12)*STATES(5)
      RATES(2) = (ALGBRC(1) - ALGBRC(5)) - ALGBRC(11)
      ALGBRC(13) = ( CONSTS(58)*STATES(27))/(CONSTS(59)+STATES(27))
      ALGBRC(10) = ( CONSTS(56)*STATES(26)*(STATES(23)+STATES(24)))/(CONSTS(57)+STATES(26))
      RATES(26) = ((ALGBRC(3)+ALGBRC(13)) - ALGBRC(7)) - ALGBRC(10)
      RATES(27) = ALGBRC(10) - ALGBRC(13)
      ALGBRC(15) =  CONSTS(35)*STATES(15)
      ALGBRC(12) =  CONSTS(33)*STATES(14)*STATES(11) -  CONSTS(34)*STATES(15)
      RATES(15) = ALGBRC(12) - ALGBRC(15)
      ALGBRC(16) =  CONSTS(20)*(CONSTS(19)+ (1.00000 - CONSTS(19))*(1.00000 - EXP(- STATES(28)/CONSTS(18) ** 3.00000)))*STATES(8)
      RATES(3) = ALGBRC(5)+ALGBRC(8)+ALGBRC(16)
      RATES(7) = ALGBRC(5)+ALGBRC(16)
      ALGBRC(9) =  CONSTS(31)*STATES(10)*STATES(12) -  CONSTS(32)*STATES(11)
      ALGBRC(17) = ( STATES(27)*STATES(11)*CONSTS(60))/(CONSTS(61)+STATES(11))
      RATES(10) = (ALGBRC(17) - ALGBRC(6)) - ALGBRC(9)
      RATES(11) = ((ALGBRC(15)+ALGBRC(9)) - ALGBRC(17)) - ALGBRC(12)
      ALGBRC(14) =  CONSTS(17)*CONSTS(16)*(CONSTS(15)+ (1.00000 - CONSTS(15))*(1.00000 - EXP(- STATES(28)/CONSTS(14) ** 3.00000)))*STATES(5)
      ALGBRC(19) =  CONSTS(21)*CONSTS(22)*STATES(5) -  CONSTS(23)*STATES(6)
      RATES(5) = (ALGBRC(11) - ALGBRC(14)) - ALGBRC(19)
      ALGBRC(20) = ( CONSTS(62)*STATES(13))/(CONSTS(63)+STATES(13))
      RATES(12) = ALGBRC(20) - ALGBRC(9)
      RATES(13) = ALGBRC(17) - ALGBRC(20)
      ALGBRC(21) =  CONSTS(38)*STATES(18)
      RATES(14) = ALGBRC(21) - ALGBRC(12)
      ALGBRC(18) =  CONSTS(36)*STATES(16)*STATES(17) -  CONSTS(37)*STATES(18)
      RATES(17) = ALGBRC(21) - ALGBRC(18)
      RATES(18) = ALGBRC(18) - ALGBRC(21)
      ALGBRC(23) =  CONSTS(26)*(CONSTS(25)+ (1.00000 - CONSTS(25))*(1.00000 - EXP(- STATES(28)/CONSTS(24) ** 3.00000)))*STATES(6)
      RATES(6) = ALGBRC(19) - ALGBRC(23)
      RATES(8) = (ALGBRC(14)+ALGBRC(23)) - ALGBRC(16)
      ALGBRC(24) =  CONSTS(41)*STATES(20)
      ALGBRC(22) =  CONSTS(39)*STATES(16)*STATES(19) -  CONSTS(40)*STATES(20)
      RATES(16) = ((ALGBRC(15)+ALGBRC(24)) - ALGBRC(18)) - ALGBRC(22)
      RATES(20) = ALGBRC(22) - ALGBRC(24)
      ALGBRC(25) = ( CONSTS(42)*STATES(21))/(CONSTS(43)+STATES(21))
      RATES(19) = ALGBRC(25) - ALGBRC(22)
      ALGBRC(26) = ( STATES(21)*STATES(22)*CONSTS(44))/(CONSTS(45)+STATES(22))
      ALGBRC(27) = ( CONSTS(46)*STATES(23))/(CONSTS(47)+STATES(23))
      RATES(22) = ALGBRC(27) - ALGBRC(26)
      ALGBRC(28) = ( STATES(21)*STATES(23)*CONSTS(48))/(CONSTS(49)+STATES(23))
      RATES(21) = ((ALGBRC(24) - ALGBRC(26)) - ALGBRC(28)) - ALGBRC(25)
      ALGBRC(29) = ( CONSTS(50)*STATES(24))/(CONSTS(51)+STATES(24))
      RATES(23) = ((ALGBRC(26)+ALGBRC(29)) - ALGBRC(27)) - ALGBRC(28)
      RATES(24) = ALGBRC(28) - ALGBRC(29)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*STATES(1)*STATES(4) -  CONSTS(2)*STATES(2)
      ALGBRC(4) =  CONSTS(65)*STATES(28)
      ALGBRC(2) = ( CONSTS(27)*2.00000*(STATES(5)+STATES(8)+STATES(6))*STATES(9))/(CONSTS(28)+STATES(9))
      ALGBRC(6) = ( CONSTS(29)*STATES(10))/(CONSTS(30)+STATES(10))
      ALGBRC(7) = ( CONSTS(54)*STATES(26))/(CONSTS(55)+STATES(26))
      ALGBRC(3) = ( CONSTS(52)*STATES(25)*(STATES(23)+STATES(24)))/(CONSTS(53)+STATES(25))
      ALGBRC(8) =  CONSTS(11)*STATES(3) -  CONSTS(10)*CONSTS(9)*(CONSTS(8)+ (1.00000 - CONSTS(8))*(1.00000 - EXP(- STATES(28)/CONSTS(7) ** 3.00000)))*STATES(1)
      ALGBRC(5) =  CONSTS(6)*CONSTS(5)*(CONSTS(4)+ (1.00000 - CONSTS(4))*(1.00000 - EXP(- STATES(28)/CONSTS(3) ** 3.00000)))*STATES(2)
      ALGBRC(11) =  CONSTS(13)*STATES(2)*STATES(2) -  CONSTS(12)*STATES(5)
      ALGBRC(13) = ( CONSTS(58)*STATES(27))/(CONSTS(59)+STATES(27))
      ALGBRC(10) = ( CONSTS(56)*STATES(26)*(STATES(23)+STATES(24)))/(CONSTS(57)+STATES(26))
      ALGBRC(15) =  CONSTS(35)*STATES(15)
      ALGBRC(12) =  CONSTS(33)*STATES(14)*STATES(11) -  CONSTS(34)*STATES(15)
      ALGBRC(16) =  CONSTS(20)*(CONSTS(19)+ (1.00000 - CONSTS(19))*(1.00000 - EXP(- STATES(28)/CONSTS(18) ** 3.00000)))*STATES(8)
      ALGBRC(9) =  CONSTS(31)*STATES(10)*STATES(12) -  CONSTS(32)*STATES(11)
      ALGBRC(17) = ( STATES(27)*STATES(11)*CONSTS(60))/(CONSTS(61)+STATES(11))
      ALGBRC(14) =  CONSTS(17)*CONSTS(16)*(CONSTS(15)+ (1.00000 - CONSTS(15))*(1.00000 - EXP(- STATES(28)/CONSTS(14) ** 3.00000)))*STATES(5)
      ALGBRC(19) =  CONSTS(21)*CONSTS(22)*STATES(5) -  CONSTS(23)*STATES(6)
      ALGBRC(20) = ( CONSTS(62)*STATES(13))/(CONSTS(63)+STATES(13))
      ALGBRC(21) =  CONSTS(38)*STATES(18)
      ALGBRC(18) =  CONSTS(36)*STATES(16)*STATES(17) -  CONSTS(37)*STATES(18)
      ALGBRC(23) =  CONSTS(26)*(CONSTS(25)+ (1.00000 - CONSTS(25))*(1.00000 - EXP(- STATES(28)/CONSTS(24) ** 3.00000)))*STATES(6)
      ALGBRC(24) =  CONSTS(41)*STATES(20)
      ALGBRC(22) =  CONSTS(39)*STATES(16)*STATES(19) -  CONSTS(40)*STATES(20)
      ALGBRC(25) = ( CONSTS(42)*STATES(21))/(CONSTS(43)+STATES(21))
      ALGBRC(26) = ( STATES(21)*STATES(22)*CONSTS(44))/(CONSTS(45)+STATES(22))
      ALGBRC(27) = ( CONSTS(46)*STATES(23))/(CONSTS(47)+STATES(23))
      ALGBRC(28) = ( STATES(21)*STATES(23)*CONSTS(48))/(CONSTS(49)+STATES(23))
      ALGBRC(29) = ( CONSTS(50)*STATES(24))/(CONSTS(51)+STATES(24))
      RETURN
      END