C
C There are a total of 5 entries in the algebraic variable array.
C There are a total of 48 entries in each of the rate and state variable arrays.
C There are a total of 101 entries in the constant variable array.
C
C
C VOI is time in component environment (day).
C ALGBRC(1) is FAt in component FAt (litre_per_day).
C CONSTS(1) is FA1 in component FAt (litre_per_day).
C CONSTS(2) is tA in component FAt (day).
C CONSTS(3) is FA2 in component FAt (litre_per_day).
C ALGBRC(4) is FBt in component FBt (litre_per_day).
C CONSTS(4) is FB1 in component FBt (litre_per_day).
C CONSTS(5) is tB in component FBt (day).
C CONSTS(6) is FB2 in component FBt (litre_per_day).
C ALGBRC(5) is F0t in component F0t (litre_per_day).
C STATES(1) is c15 in component c15 (micromole_per_litre).
C CONSTS(7) is p0 in component model_parameters (litre_per_day).
C CONSTS(8) is V2 in component model_parameters (litre).
C STATES(2) is c015 in component c015 (micromole_per_litre).
C CONSTS(9) is p16 in component model_parameters (per_day).
C STATES(3) is c16 in component c16 (micromole_per_litre).
C CONSTS(10) is p18 in component model_parameters (per_day).
C STATES(4) is c17 in component c17 (micromole_per_litre).
C CONSTS(11) is p15 in component model_parameters (per_day_per_micromole_litre).
C STATES(5) is c19 in component c19 (micromole_per_litre).
C CONSTS(12) is p17 in component model_parameters (per_day).
C CONSTS(13) is p19 in component model_parameters (per_day_per_micromole_litre).
C CONSTS(14) is p20 in component model_parameters (micromole_per_litre).
C CONSTS(15) is p29 in component model_parameters (per_day).
C CONSTS(16) is pAA15 in component model_parameters (dimensionless).
C ALGBRC(2) is pp in component pp (dimensionless).
C STATES(6) is c016 in component c016 (micromole_per_litre).
C CONSTS(17) is pAA16 in component model_parameters (dimensionless).
C STATES(7) is c017 in component c017 (micromole_per_litre).
C CONSTS(18) is p22 in component model_parameters (per_day).
C STATES(8) is c18 in component c18 (micromole_per_litre).
C CONSTS(19) is p21 in component model_parameters (per_day).
C CONSTS(20) is p23 in component model_parameters (per_day).
C CONSTS(21) is p1 in component model_parameters (per_day).
C CONSTS(22) is s1 in component c17 (dimensionless).
C STATES(9) is c1 in component c1 (micromole_per_litre).
C CONSTS(23) is p2 in component model_parameters (per_day).
C CONSTS(24) is s2 in component c17 (dimensionless).
C STATES(10) is c2 in component c2 (micromole_per_litre).
C CONSTS(25) is p3 in component model_parameters (per_day).
C CONSTS(26) is s3 in component c17 (dimensionless).
C STATES(11) is c3 in component c3 (micromole_per_litre).
C CONSTS(27) is p4 in component model_parameters (per_day).
C CONSTS(28) is s4 in component c17 (dimensionless).
C STATES(12) is c4 in component c4 (micromole_per_litre).
C CONSTS(29) is p5 in component model_parameters (per_day).
C CONSTS(30) is s5 in component c17 (dimensionless).
C STATES(13) is c5 in component c5 (micromole_per_litre).
C CONSTS(31) is p6 in component model_parameters (per_day).
C CONSTS(32) is s6 in component c17 (dimensionless).
C STATES(14) is c6 in component c6 (micromole_per_litre).
C CONSTS(33) is p7 in component model_parameters (per_day).
C CONSTS(34) is s7 in component c17 (dimensionless).
C STATES(15) is c7 in component c7 (micromole_per_litre).
C CONSTS(35) is p8 in component model_parameters (per_day).
C CONSTS(36) is s8 in component c17 (dimensionless).
C STATES(16) is c8 in component c8 (micromole_per_litre).
C CONSTS(37) is p9 in component model_parameters (per_day).
C CONSTS(38) is s9 in component c17 (dimensionless).
C STATES(17) is c9 in component c9 (micromole_per_litre).
C CONSTS(39) is p14 in component model_parameters (per_day).
C CONSTS(40) is s14 in component c17 (dimensionless).
C STATES(18) is c14 in component c14 (micromole_per_litre).
C CONSTS(41) is p25 in component model_parameters (per_day_per_micromole_litre).
C CONSTS(42) is p26 in component model_parameters (per_day).
C CONSTS(43) is p28 in component model_parameters (per_day).
C CONSTS(44) is pAA17 in component model_parameters (dimensionless).
C STATES(19) is c018 in component c018 (micromole_per_litre).
C CONSTS(45) is pAA18 in component model_parameters (dimensionless).
C STATES(20) is c019 in component c019 (micromole_per_litre).
C ALGBRC(3) is gt in component gt (dimensionless).
C CONSTS(46) is p11 in component model_parameters (per_day).
C CONSTS(47) is s11 in component c19 (dimensionless).
C STATES(21) is c11 in component c11 (micromole_per_litre).
C CONSTS(48) is p12 in component model_parameters (per_day).
C CONSTS(49) is s12 in component c19 (dimensionless).
C STATES(22) is c12 in component c12 (micromole_per_litre).
C CONSTS(50) is p13 in component model_parameters (per_day).
C CONSTS(51) is s13 in component c19 (dimensionless).
C STATES(23) is c13 in component c13 (micromole_per_litre).
C CONSTS(52) is p10 in component model_parameters (per_day).
C CONSTS(53) is s10 in component c19 (dimensionless).
C STATES(24) is c10 in component c10 (micromole_per_litre).
C STATES(25) is c20 in component c20 (micromole_per_litre).
C STATES(26) is c020 in component c020 (micromole_per_litre).
C STATES(27) is c21 in component c21 (micromole_per_litre).
C STATES(28) is c021 in component c021 (micromole_per_litre).
C STATES(29) is c22 in component c22 (micromole_per_litre).
C STATES(30) is c022 in component c022 (micromole_per_litre).
C STATES(31) is c23 in component c23 (micromole_per_litre).
C STATES(32) is c023 in component c023 (micromole_per_litre).
C CONSTS(54) is p31 in component model_parameters (per_day).
C STATES(33) is c24 in component c24 (micromole_per_litre).
C STATES(34) is c024 in component c024 (micromole_per_litre).
C CONSTS(55) is p30 in component model_parameters (dimensionless).
C CONSTS(56) is p32 in component model_parameters (per_day).
C CONSTS(57) is p24 in component model_parameters (per_day).
C CONSTS(58) is p33 in component pp (micromole_per_litre_per_day).
C CONSTS(59) is t1 in component pp (day).
C CONSTS(60) is p27 in component model_parameters (dimensionless).
C STATES(35) is c01 in component c01 (micromole_per_litre).
C CONSTS(61) is pAA1 in component model_parameters (dimensionless).
C STATES(36) is c02 in component c02 (micromole_per_litre).
C CONSTS(62) is pAA2 in component model_parameters (dimensionless).
C STATES(37) is c03 in component c03 (micromole_per_litre).
C CONSTS(63) is pAA3 in component model_parameters (dimensionless).
C STATES(38) is c04 in component c04 (micromole_per_litre).
C CONSTS(64) is pAA4 in component model_parameters (dimensionless).
C STATES(39) is c05 in component c05 (micromole_per_litre).
C CONSTS(65) is pAA5 in component model_parameters (dimensionless).
C STATES(40) is c06 in component c06 (micromole_per_litre).
C CONSTS(66) is pAA6 in component model_parameters (dimensionless).
C STATES(41) is c07 in component c07 (micromole_per_litre).
C CONSTS(67) is pAA7 in component model_parameters (dimensionless).
C STATES(42) is c08 in component c08 (micromole_per_litre).
C CONSTS(68) is pAA8 in component model_parameters (dimensionless).
C STATES(43) is c09 in component c09 (micromole_per_litre).
C CONSTS(69) is pAA9 in component model_parameters (dimensionless).
C STATES(44) is c010 in component c010 (micromole_per_litre).
C CONSTS(70) is pAA10 in component model_parameters (dimensionless).
C STATES(45) is c011 in component c011 (micromole_per_litre).
C CONSTS(71) is pAA11 in component model_parameters (dimensionless).
C STATES(46) is c012 in component c012 (micromole_per_litre).
C CONSTS(72) is pAA12 in component model_parameters (dimensionless).
C STATES(47) is c013 in component c013 (micromole_per_litre).
C CONSTS(73) is pAA13 in component model_parameters (dimensionless).
C STATES(48) is c014 in component c014 (micromole_per_litre).
C CONSTS(74) is pAA14 in component model_parameters (dimensionless).
C CONSTS(75) is V1 in component model_parameters (litre).
C CONSTS(76) is CA1 in component c01 (micromole_per_litre).
C CONSTS(77) is CB1 in component model_parameters (micromole_per_litre).
C CONSTS(78) is CA2 in component c02 (micromole_per_litre).
C CONSTS(79) is CA3 in component c03 (micromole_per_litre).
C CONSTS(80) is CA4 in component c04 (micromole_per_litre).
C CONSTS(81) is CA5 in component c05 (micromole_per_litre).
C CONSTS(82) is CA6 in component c06 (micromole_per_litre).
C CONSTS(83) is CA7 in component c07 (micromole_per_litre).
C CONSTS(84) is CA8 in component c08 (micromole_per_litre).
C CONSTS(85) is CA9 in component c09 (micromole_per_litre).
C CONSTS(86) is CA10 in component c010 (micromole_per_litre).
C CONSTS(87) is CA11 in component c011 (micromole_per_litre).
C CONSTS(88) is CA12 in component c012 (micromole_per_litre).
C CONSTS(89) is CA13 in component c013 (micromole_per_litre).
C CONSTS(90) is CA14 in component c014 (micromole_per_litre).
C CONSTS(91) is CA15 in component c015 (micromole_per_litre).
C CONSTS(92) is CB15 in component model_parameters (micromole_per_litre).
C CONSTS(93) is CA16 in component c016 (micromole_per_litre).
C CONSTS(94) is CA17 in component c017 (micromole_per_litre).
C CONSTS(95) is CA18 in component c018 (micromole_per_litre).
C CONSTS(96) is CA19 in component c019 (micromole_per_litre).
C CONSTS(97) is CA20 in component c020 (micromole_per_litre).
C CONSTS(98) is CA21 in component c021 (micromole_per_litre).
C CONSTS(99) is CA22 in component c022 (micromole_per_litre).
C CONSTS(100) is CA23 in component c023 (micromole_per_litre).
C CONSTS(101) is CA24 in component c024 (micromole_per_litre).
C RATES(1) is d/dt c15 in component c15 (micromole_per_litre).
C RATES(3) is d/dt c16 in component c16 (micromole_per_litre).
C RATES(4) is d/dt c17 in component c17 (micromole_per_litre).
C RATES(8) is d/dt c18 in component c18 (micromole_per_litre).
C RATES(5) is d/dt c19 in component c19 (micromole_per_litre).
C RATES(25) is d/dt c20 in component c20 (micromole_per_litre).
C RATES(27) is d/dt c21 in component c21 (micromole_per_litre).
C RATES(29) is d/dt c22 in component c22 (micromole_per_litre).
C RATES(31) is d/dt c23 in component c23 (micromole_per_litre).
C RATES(33) is d/dt c24 in component c24 (micromole_per_litre).
C RATES(9) is d/dt c1 in component c1 (micromole_per_litre).
C RATES(10) is d/dt c2 in component c2 (micromole_per_litre).
C RATES(11) is d/dt c3 in component c3 (micromole_per_litre).
C RATES(12) is d/dt c4 in component c4 (micromole_per_litre).
C RATES(13) is d/dt c5 in component c5 (micromole_per_litre).
C RATES(14) is d/dt c6 in component c6 (micromole_per_litre).
C RATES(15) is d/dt c7 in component c7 (micromole_per_litre).
C RATES(16) is d/dt c8 in component c8 (micromole_per_litre).
C RATES(17) is d/dt c9 in component c9 (micromole_per_litre).
C RATES(24) is d/dt c10 in component c10 (micromole_per_litre).
C RATES(21) is d/dt c11 in component c11 (micromole_per_litre).
C RATES(22) is d/dt c12 in component c12 (micromole_per_litre).
C RATES(23) is d/dt c13 in component c13 (micromole_per_litre).
C RATES(18) is d/dt c14 in component c14 (micromole_per_litre).
C RATES(35) is d/dt c01 in component c01 (micromole_per_litre).
C RATES(36) is d/dt c02 in component c02 (micromole_per_litre).
C RATES(37) is d/dt c03 in component c03 (micromole_per_litre).
C RATES(38) is d/dt c04 in component c04 (micromole_per_litre).
C RATES(39) is d/dt c05 in component c05 (micromole_per_litre).
C RATES(40) is d/dt c06 in component c06 (micromole_per_litre).
C RATES(41) is d/dt c07 in component c07 (micromole_per_litre).
C RATES(42) is d/dt c08 in component c08 (micromole_per_litre).
C RATES(43) is d/dt c09 in component c09 (micromole_per_litre).
C RATES(44) is d/dt c010 in component c010 (micromole_per_litre).
C RATES(45) is d/dt c011 in component c011 (micromole_per_litre).
C RATES(46) is d/dt c012 in component c012 (micromole_per_litre).
C RATES(47) is d/dt c013 in component c013 (micromole_per_litre).
C RATES(48) is d/dt c014 in component c014 (micromole_per_litre).
C RATES(2) is d/dt c015 in component c015 (micromole_per_litre).
C RATES(6) is d/dt c016 in component c016 (micromole_per_litre).
C RATES(7) is d/dt c017 in component c017 (micromole_per_litre).
C RATES(19) is d/dt c018 in component c018 (micromole_per_litre).
C RATES(20) is d/dt c019 in component c019 (micromole_per_litre).
C RATES(26) is d/dt c020 in component c020 (micromole_per_litre).
C RATES(28) is d/dt c021 in component c021 (micromole_per_litre).
C RATES(30) is d/dt c022 in component c022 (micromole_per_litre).
C RATES(32) is d/dt c023 in component c023 (micromole_per_litre).
C RATES(34) is d/dt c024 in component c024 (micromole_per_litre).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 3.6
      CONSTS(2) = 1
      CONSTS(3) = 1.2
      CONSTS(4) = 0
      CONSTS(5) = 3
      CONSTS(6) = 0.024
      STATES(1) = 284
      CONSTS(7) = 1191.072
      CONSTS(8) = 0.9
      STATES(2) = 284
      CONSTS(9) = 3.593
      STATES(3) = 168
      CONSTS(10) = 56.73
      STATES(4) = 265
      CONSTS(11) = 0.0009
      STATES(5) = 41
      CONSTS(12) = 50.08
      CONSTS(13) = 225.76
      CONSTS(14) = 5
      CONSTS(15) = 2.8119
      CONSTS(16) = 0.31
      STATES(6) = 168
      CONSTS(17) = 0.25
      STATES(7) = 265
      CONSTS(18) = 0.0205
      STATES(8) = 687
      CONSTS(19) = 0.0033
      CONSTS(20) = 40.68
      CONSTS(21) = 1.816
      CONSTS(22) = 1
      STATES(9) = 2258
      CONSTS(23) = 16.29
      CONSTS(24) = 3
      STATES(10) = 940
      CONSTS(25) = 100
      CONSTS(26) = 4
      STATES(11) = 2604
      CONSTS(27) = 1.147
      CONSTS(28) = 1
      STATES(12) = 987
      CONSTS(29) = 15.99
      CONSTS(30) = 2
      STATES(13) = 464
      CONSTS(31) = 4.148
      CONSTS(32) = 1
      STATES(14) = 1300
      CONSTS(33) = 3.723
      CONSTS(34) = 1
      STATES(15) = 500
      CONSTS(35) = 1.573
      CONSTS(36) = 1
      STATES(16) = 1825
      CONSTS(37) = 1.135
      CONSTS(38) = 1
      STATES(17) = 2215
      CONSTS(39) = 9.326
      CONSTS(40) = 1
      STATES(18) = 859
      CONSTS(41) = 0
      CONSTS(42) = 0
      CONSTS(43) = 38.05
      CONSTS(44) = 1.27
      STATES(19) = 687
      CONSTS(45) = 0.06
      STATES(20) = 41
      CONSTS(46) = 34.53
      CONSTS(47) = 1
      STATES(21) = 259
      CONSTS(48) = 10.8
      CONSTS(49) = 1
      STATES(22) = 861
      CONSTS(50) = 10.76
      CONSTS(51) = 1
      STATES(23) = 1957
      CONSTS(52) = 0.192
      CONSTS(53) = 1
      STATES(24) = 327
      STATES(25) = 0
      STATES(26) = 0
      STATES(27) = 0
      STATES(28) = 0
      STATES(29) = 0
      STATES(30) = 0
      STATES(31) = 8159
      STATES(32) = 0
      CONSTS(54) = 1.31
      STATES(33) = 8159
      STATES(34) = 0
      CONSTS(55) = 0
      CONSTS(56) = 478.7
      CONSTS(57) = 0.5474
      CONSTS(58) = 0
      CONSTS(59) = 1
      CONSTS(60) = 0.6804
      STATES(35) = 2258
      CONSTS(61) = 1.4
      STATES(36) = 940
      CONSTS(62) = 0.18
      STATES(37) = 2604
      CONSTS(63) = 0
      STATES(38) = 987
      CONSTS(64) = 1.19
      STATES(39) = 464
      CONSTS(65) = 0.45
      STATES(40) = 1300
      CONSTS(66) = 1.05
      STATES(41) = 500
      CONSTS(67) = 0.62
      STATES(42) = 1825
      CONSTS(68) = 1.15
      STATES(43) = 2215
      CONSTS(69) = 1
      STATES(44) = 327
      CONSTS(70) = 0.15
      STATES(45) = 259
      CONSTS(71) = 0.77
      STATES(46) = 861
      CONSTS(72) = 1.55
      STATES(47) = 1957
      CONSTS(73) = 0.78
      STATES(48) = 859
      CONSTS(74) = 1.81
      CONSTS(75) = 0.6
      CONSTS(76) = 2258
      CONSTS(77) = 0
      CONSTS(78) = 940
      CONSTS(79) = 2604
      CONSTS(80) = 987
      CONSTS(81) = 464
      CONSTS(82) = 1300
      CONSTS(83) = 500
      CONSTS(84) = 1825
      CONSTS(85) = 2215
      CONSTS(86) = 327
      CONSTS(87) = 259
      CONSTS(88) = 861
      CONSTS(89) = 1957
      CONSTS(90) = 859
      CONSTS(91) = 284
      CONSTS(92) = 1500
      CONSTS(93) = 168
      CONSTS(94) = 265
      CONSTS(95) = 687
      CONSTS(96) = 41
      CONSTS(97) = 0
      CONSTS(98) = 0
      CONSTS(99) = 0
      CONSTS(100) = 0
      CONSTS(101) = 0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(25) =  (CONSTS(7)/CONSTS(8))*(STATES(26) - STATES(25))+ (( CONSTS(13)*STATES(1))/(CONSTS(14)+STATES(1)))*STATES(5)
      RATES(27) =  (CONSTS(7)/CONSTS(8))*(STATES(28) - STATES(27)) -  CONSTS(19)*STATES(27)
      RATES(29) =  (CONSTS(7)/CONSTS(8))*(STATES(30) - STATES(29)) -  CONSTS(18)*STATES(29)
      RATES(31) = (( (CONSTS(7)/CONSTS(8))*(STATES(32) - STATES(31))+ CONSTS(19)*STATES(27)+ CONSTS(18)*STATES(29)) -  CONSTS(20)*STATES(31))+ (( CONSTS(54)*STATES(1))/(CONSTS(14)+STATES(1)))*STATES(5)
      RATES(33) = ( (CONSTS(7)/CONSTS(8))*(STATES(34) - STATES(33))+ 2.00000*CONSTS(55)*CONSTS(20)*STATES(31)+ (( CONSTS(56)*STATES(1))/(CONSTS(14)+STATES(1)))*STATES(5)) -  CONSTS(57)*STATES(33)
      ALGBRC(2) = TERNRY(VOI.LT.CONSTS(59), CONSTS(58), 0.00000)
      RATES(1) = ( (CONSTS(7)/CONSTS(8))*(STATES(2) - STATES(1))+ CONSTS(9)*STATES(3)+ CONSTS(12)*STATES(4)+ CONSTS(16)*ALGBRC(2)) -  ( CONSTS(11)*STATES(5)+CONSTS(10)+( CONSTS(13)*STATES(5))/(CONSTS(14)+STATES(1))+CONSTS(15))*STATES(1)
      RATES(3) = (( (CONSTS(7)/CONSTS(8))*(STATES(6) - STATES(3))+ CONSTS(11)*STATES(1)*STATES(5)) -  CONSTS(9)*STATES(3))+ CONSTS(17)*ALGBRC(2)
      RATES(4) = (( (CONSTS(7)/CONSTS(8))*(STATES(7) - STATES(4))+ CONSTS(21)*CONSTS(22)*STATES(9)+ CONSTS(23)*CONSTS(24)*STATES(10)+ CONSTS(25)*CONSTS(26)*STATES(11)+ CONSTS(27)*CONSTS(28)*STATES(12)+ CONSTS(29)*CONSTS(30)*STATES(13)+ CONSTS(31)*CONSTS(32)*STATES(14)+ CONSTS(33)*CONSTS(34)*STATES(15)+ CONSTS(35)*CONSTS(36)*STATES(16)+ CONSTS(37)*CONSTS(38)*STATES(17)+ CONSTS(39)*CONSTS(40)*STATES(18)+ CONSTS(10)*STATES(1)+ CONSTS(42)*STATES(8)) -  (CONSTS(12)+ CONSTS(41)*STATES(5)+CONSTS(43))*STATES(4))+ CONSTS(44)*ALGBRC(2)
      RATES(8) = (( (CONSTS(7)/CONSTS(8))*(STATES(19) - STATES(8))+ CONSTS(41)*STATES(4)*STATES(5)) -  CONSTS(42)*STATES(8))+ CONSTS(45)*ALGBRC(2)
      ALGBRC(3) = TERNRY(VOI.LT.3.00000, 0.00000, CONSTS(60))
      RATES(5) = ( (CONSTS(7)/CONSTS(8))*(STATES(20) - STATES(5))+ CONSTS(46)*CONSTS(47)*STATES(21)+ CONSTS(48)*CONSTS(49)*STATES(22)+ CONSTS(50)*STATES(23)*CONSTS(51)+ CONSTS(52)*STATES(24)*CONSTS(53)+ CONSTS(20)*(1.00000 - ALGBRC(3))*STATES(4)+ CONSTS(9)*STATES(3)+ CONSTS(42)*STATES(8)) -  (( CONSTS(13)*STATES(1))/(CONSTS(14)+STATES(1))+ CONSTS(11)*STATES(1)+ CONSTS(41)*STATES(4))*STATES(5)
      RATES(9) = ( (CONSTS(7)/CONSTS(8))*(STATES(35) - STATES(9)) -  CONSTS(21)*STATES(9))+ CONSTS(61)*ALGBRC(2)
      RATES(10) = ( (CONSTS(7)/CONSTS(8))*(STATES(36) - STATES(10)) -  CONSTS(23)*STATES(10))+ CONSTS(62)*ALGBRC(2)
      RATES(11) = ( (CONSTS(7)/CONSTS(8))*(STATES(37) - STATES(11)) -  CONSTS(25)*STATES(11))+ CONSTS(63)*ALGBRC(2)
      RATES(12) = ( (CONSTS(7)/CONSTS(8))*(STATES(38) - STATES(12)) -  CONSTS(27)*STATES(12))+ CONSTS(64)*ALGBRC(2)
      RATES(13) = ( (CONSTS(7)/CONSTS(8))*(STATES(39) - STATES(13)) -  CONSTS(29)*STATES(13))+ CONSTS(65)*ALGBRC(2)
      RATES(14) = ( (CONSTS(7)/CONSTS(8))*(STATES(40) - STATES(14)) -  CONSTS(31)*STATES(14))+ CONSTS(66)*ALGBRC(2)
      RATES(15) = ( (CONSTS(7)/CONSTS(8))*(STATES(41) - STATES(15)) -  CONSTS(33)*STATES(15))+ CONSTS(67)*ALGBRC(2)
      RATES(16) = ( (CONSTS(7)/CONSTS(8))*(STATES(42) - STATES(16)) -  CONSTS(35)*STATES(16))+ CONSTS(68)*ALGBRC(2)
      RATES(17) = ( (CONSTS(7)/CONSTS(8))*(STATES(43) - STATES(17)) -  CONSTS(37)*STATES(17))+ CONSTS(69)*ALGBRC(2)
      RATES(24) = ( (CONSTS(7)/CONSTS(8))*(STATES(44) - STATES(24)) -  CONSTS(52)*STATES(24))+ CONSTS(70)*ALGBRC(2)
      RATES(21) = ( (CONSTS(7)/CONSTS(8))*(STATES(45) - STATES(21)) -  CONSTS(46)*STATES(21))+ CONSTS(71)*ALGBRC(2)
      RATES(22) = ( (CONSTS(7)/CONSTS(8))*(STATES(46) - STATES(22)) -  CONSTS(48)*STATES(22))+ CONSTS(72)*ALGBRC(2)
      RATES(23) = ( (CONSTS(7)/CONSTS(8))*(STATES(47) - STATES(23)) -  CONSTS(50)*STATES(23))+ CONSTS(73)*ALGBRC(2)
      RATES(18) = ( (CONSTS(7)/CONSTS(8))*(STATES(48) - STATES(18)) -  CONSTS(39)*STATES(18))+ CONSTS(74)*ALGBRC(2)
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(2), CONSTS(3), CONSTS(1))
      ALGBRC(4) = TERNRY(VOI.GE.CONSTS(5), CONSTS(6), CONSTS(4))
      ALGBRC(5) = ALGBRC(1)+ALGBRC(4)
      RATES(35) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(76)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(35)) -  (CONSTS(7)/CONSTS(75))*(STATES(35) - STATES(9))
      RATES(36) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(78)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(36)) -  (CONSTS(7)/CONSTS(75))*(STATES(36) - STATES(10))
      RATES(37) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(79)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(37)) -  (CONSTS(7)/CONSTS(75))*(STATES(37) - STATES(11))
      RATES(38) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(80)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(38)) -  (CONSTS(7)/CONSTS(75))*(STATES(38) - STATES(12))
      RATES(39) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(81)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(39)) -  (CONSTS(7)/CONSTS(75))*(STATES(39) - STATES(13))
      RATES(40) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(82)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(40)) -  (CONSTS(7)/CONSTS(75))*(STATES(40) - STATES(14))
      RATES(41) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(83)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(41)) -  (CONSTS(7)/CONSTS(75))*(STATES(41) - STATES(15))
      RATES(42) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(84)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(42)) -  (CONSTS(7)/CONSTS(75))*(STATES(42) - STATES(16))
      RATES(43) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(85)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(43)) -  (CONSTS(7)/CONSTS(75))*(STATES(43) - STATES(17))
      RATES(44) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(86)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(44)) -  (CONSTS(7)/CONSTS(75))*(STATES(44) - STATES(24))
      RATES(45) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(87)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(45)) -  (CONSTS(7)/CONSTS(75))*(STATES(45) - STATES(21))
      RATES(46) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(88)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(46)) -  (CONSTS(7)/CONSTS(75))*(STATES(46) - STATES(22))
      RATES(47) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(89)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(47)) -  (CONSTS(7)/CONSTS(75))*(STATES(47) - STATES(23))
      RATES(48) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(90)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(48)) -  (CONSTS(7)/CONSTS(75))*(STATES(48) - STATES(18))
      RATES(2) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(91)+ (ALGBRC(4)/CONSTS(75))*CONSTS(92)) -  (ALGBRC(5)/CONSTS(75))*STATES(2)) -  (CONSTS(7)/CONSTS(75))*(STATES(2) - STATES(1))
      RATES(6) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(93)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(6)) -  (CONSTS(7)/CONSTS(75))*(STATES(6) - STATES(3))
      RATES(7) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(94)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(7)) -  (CONSTS(7)/CONSTS(75))*(STATES(7) - STATES(4))
      RATES(19) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(95)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(19)) -  (CONSTS(7)/CONSTS(75))*(STATES(19) - STATES(8))
      RATES(20) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(96)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(20)) -  (CONSTS(7)/CONSTS(75))*(STATES(20) - STATES(5))
      RATES(26) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(97)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(26)) -  (CONSTS(7)/CONSTS(75))*(STATES(26) - STATES(25))
      RATES(28) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(98)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(28)) -  (CONSTS(7)/CONSTS(75))*(STATES(28) - STATES(27))
      RATES(30) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(99)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(30)) -  (CONSTS(7)/CONSTS(75))*(STATES(30) - STATES(29))
      RATES(32) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(100)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(32)) -  (CONSTS(7)/CONSTS(75))*(STATES(32) - STATES(31))
      RATES(34) = (( (ALGBRC(1)/CONSTS(75))*CONSTS(101)+ (ALGBRC(4)/CONSTS(75))*CONSTS(77)) -  (ALGBRC(5)/CONSTS(75))*STATES(34)) -  (CONSTS(7)/CONSTS(75))*(STATES(34) - STATES(33))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = TERNRY(VOI.LT.CONSTS(59), CONSTS(58), 0.00000)
      ALGBRC(3) = TERNRY(VOI.LT.3.00000, 0.00000, CONSTS(60))
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(2), CONSTS(3), CONSTS(1))
      ALGBRC(4) = TERNRY(VOI.GE.CONSTS(5), CONSTS(6), CONSTS(4))
      ALGBRC(5) = ALGBRC(1)+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