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 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