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 124 entries in the algebraic variable array. C There are a total of 19 entries in each of the rate and state variable arrays. C There are a total of 58 entries in the constant variable array. C C C VOI is t in component environment (second). C CONSTS(1) is R in component environment (J_per_K_per_mol). C CONSTS(2) is T in component environment (kelvin). C CONSTS(3) is F in component environment (C_per_mol). C CONSTS(4) is C_m in component environment (fF). C STATES(1) is q_Na_o in component environment (fmol). C STATES(2) is q_Na_i in component environment (fmol). C STATES(3) is q_S000_Na in component environment (fmol). C STATES(4) is q_S010_Na in component environment (fmol). C STATES(5) is q_S100_Na in component environment (fmol). C STATES(6) is q_S110_Na in component environment (fmol). C STATES(7) is q_S200_Na in component environment (fmol). C STATES(8) is q_S210_Na in component environment (fmol). C STATES(9) is q_S300_Na in component environment (fmol). C STATES(10) is q_S310_Na in component environment (fmol). C STATES(11) is q_S001_Na in component environment (fmol). C STATES(12) is q_S011_Na in component environment (fmol). C STATES(13) is q_S101_Na in component environment (fmol). C STATES(14) is q_S111_Na in component environment (fmol). C STATES(15) is q_S201_Na in component environment (fmol). C STATES(16) is q_S211_Na in component environment (fmol). C STATES(17) is q_S301_Na in component environment (fmol). C STATES(18) is q_S311_Na in component environment (fmol). C ALGBRC(121) is v_Na in component fast_Na (fmol_per_sec). C STATES(19) is q_mem in component environment (fC). C ALGBRC(123) is I_mem_Na in component fast_Na (fA). C CONSTS(5) is K_000_Na in component Na_parameters (per_fmol). C CONSTS(6) is K_001_Na in component Na_parameters (per_fmol). C CONSTS(7) is K_010_Na in component Na_parameters (per_fmol). C CONSTS(8) is K_011_Na in component Na_parameters (per_fmol). C CONSTS(9) is K_100_Na in component Na_parameters (per_fmol). C CONSTS(10) is K_101_Na in component Na_parameters (per_fmol). C CONSTS(11) is K_110_Na in component Na_parameters (per_fmol). C CONSTS(12) is K_111_Na in component Na_parameters (per_fmol). C CONSTS(13) is K_200_Na in component Na_parameters (per_fmol). C CONSTS(14) is K_201_Na in component Na_parameters (per_fmol). C CONSTS(15) is K_210_Na in component Na_parameters (per_fmol). C CONSTS(16) is K_211_Na in component Na_parameters (per_fmol). C CONSTS(17) is K_300_Na in component Na_parameters (per_fmol). C CONSTS(18) is K_301_Na in component Na_parameters (per_fmol). C CONSTS(19) is K_310_Na in component Na_parameters (per_fmol). C CONSTS(20) is K_311_Na in component Na_parameters (per_fmol). C CONSTS(21) is K_Na_o in component Na_parameters (per_fmol). C CONSTS(22) is K_Na_i in component Na_parameters (per_fmol). C CONSTS(23) is kappa_Na in component Na_parameters (fmol_per_sec). C CONSTS(24) is kappa_h000 in component Na_parameters (fmol_per_sec). C CONSTS(25) is kappa_h001 in component Na_parameters (fmol_per_sec). C CONSTS(26) is kappa_h100 in component Na_parameters (fmol_per_sec). C CONSTS(27) is kappa_h101 in component Na_parameters (fmol_per_sec). C CONSTS(28) is kappa_h200 in component Na_parameters (fmol_per_sec). C CONSTS(29) is kappa_h201 in component Na_parameters (fmol_per_sec). C CONSTS(30) is kappa_h300 in component Na_parameters (fmol_per_sec). C CONSTS(31) is kappa_h301 in component Na_parameters (fmol_per_sec). C CONSTS(32) is kappa_j000 in component Na_parameters (fmol_per_sec). C CONSTS(33) is kappa_j010 in component Na_parameters (fmol_per_sec). C CONSTS(34) is kappa_j100 in component Na_parameters (fmol_per_sec). C CONSTS(35) is kappa_j110 in component Na_parameters (fmol_per_sec). C CONSTS(36) is kappa_j200 in component Na_parameters (fmol_per_sec). C CONSTS(37) is kappa_j210 in component Na_parameters (fmol_per_sec). C CONSTS(38) is kappa_j300 in component Na_parameters (fmol_per_sec). C CONSTS(39) is kappa_j310 in component Na_parameters (fmol_per_sec). C CONSTS(40) is kappa_m000 in component Na_parameters (fmol_per_sec). C CONSTS(41) is kappa_m001 in component Na_parameters (fmol_per_sec). C CONSTS(42) is kappa_m010 in component Na_parameters (fmol_per_sec). C CONSTS(43) is kappa_m011 in component Na_parameters (fmol_per_sec). C CONSTS(44) is kappa_m100 in component Na_parameters (fmol_per_sec). C CONSTS(45) is kappa_m101 in component Na_parameters (fmol_per_sec). C CONSTS(46) is kappa_m110 in component Na_parameters (fmol_per_sec). C CONSTS(47) is kappa_m111 in component Na_parameters (fmol_per_sec). C CONSTS(48) is kappa_m200 in component Na_parameters (fmol_per_sec). C CONSTS(49) is kappa_m201 in component Na_parameters (fmol_per_sec). C CONSTS(50) is kappa_m210 in component Na_parameters (fmol_per_sec). C CONSTS(51) is kappa_m211 in component Na_parameters (fmol_per_sec). C CONSTS(52) is zNa in component Na_parameters (dimensionless). C CONSTS(53) is z_fh in component Na_parameters (dimensionless). C CONSTS(54) is z_fj in component Na_parameters (dimensionless). C CONSTS(55) is z_fm in component Na_parameters (dimensionless). C CONSTS(56) is z_rh in component Na_parameters (dimensionless). C CONSTS(57) is z_rj in component Na_parameters (dimensionless). C CONSTS(58) is z_rm in component Na_parameters (dimensionless). C ALGBRC(17) is mu_Na_o in component fast_Na (J_per_mol). C ALGBRC(18) is mu_Na_i in component fast_Na (J_per_mol). C ALGBRC(2) is mu_S000_Na in component fast_Na (J_per_mol). C ALGBRC(3) is mu_S010_Na in component fast_Na (J_per_mol). C ALGBRC(4) is mu_S100_Na in component fast_Na (J_per_mol). C ALGBRC(5) is mu_S110_Na in component fast_Na (J_per_mol). C ALGBRC(6) is mu_S200_Na in component fast_Na (J_per_mol). C ALGBRC(7) is mu_S210_Na in component fast_Na (J_per_mol). C ALGBRC(8) is mu_S300_Na in component fast_Na (J_per_mol). C ALGBRC(9) is mu_S310_Na in component fast_Na (J_per_mol). C ALGBRC(10) is mu_S001_Na in component fast_Na (J_per_mol). C ALGBRC(11) is mu_S011_Na in component fast_Na (J_per_mol). C ALGBRC(12) is mu_S101_Na in component fast_Na (J_per_mol). C ALGBRC(13) is mu_S111_Na in component fast_Na (J_per_mol). C ALGBRC(14) is mu_S201_Na in component fast_Na (J_per_mol). C ALGBRC(15) is mu_S211_Na in component fast_Na (J_per_mol). C ALGBRC(16) is mu_S301_Na in component fast_Na (J_per_mol). C ALGBRC(19) is mu_S311_Na in component fast_Na (J_per_mol). C ALGBRC(49) is Am_Na in component fast_Na (J_per_mol). C ALGBRC(20) is Af_Na in component fast_Na (J_per_mol). C ALGBRC(50) is Ar_Na in component fast_Na (J_per_mol). C ALGBRC(21) is Af_h000 in component fast_Na (J_per_mol). C ALGBRC(51) is Ar_h000 in component fast_Na (J_per_mol). C ALGBRC(52) is v_h000 in component fast_Na (fmol_per_sec). C ALGBRC(22) is Af_h001 in component fast_Na (J_per_mol). C ALGBRC(53) is Ar_h001 in component fast_Na (J_per_mol). C ALGBRC(54) is v_h001 in component fast_Na (fmol_per_sec). C ALGBRC(23) is Af_h100 in component fast_Na (J_per_mol). C ALGBRC(55) is Ar_h100 in component fast_Na (J_per_mol). C ALGBRC(56) is v_h100 in component fast_Na (fmol_per_sec). C ALGBRC(24) is Af_h101 in component fast_Na (J_per_mol). C ALGBRC(57) is Ar_h101 in component fast_Na (J_per_mol). C ALGBRC(58) is v_h101 in component fast_Na (fmol_per_sec). C ALGBRC(25) is Af_h200 in component fast_Na (J_per_mol). C ALGBRC(59) is Ar_h200 in component fast_Na (J_per_mol). C ALGBRC(60) is v_h200 in component fast_Na (fmol_per_sec). C ALGBRC(26) is Af_h201 in component fast_Na (J_per_mol). C ALGBRC(61) is Ar_h201 in component fast_Na (J_per_mol). C ALGBRC(62) is v_h201 in component fast_Na (fmol_per_sec). C ALGBRC(27) is Af_h300 in component fast_Na (J_per_mol). C ALGBRC(63) is Ar_h300 in component fast_Na (J_per_mol). C ALGBRC(64) is v_h300 in component fast_Na (fmol_per_sec). C ALGBRC(28) is Af_h301 in component fast_Na (J_per_mol). C ALGBRC(65) is Ar_h301 in component fast_Na (J_per_mol). C ALGBRC(66) is v_h301 in component fast_Na (fmol_per_sec). C ALGBRC(29) is Af_j000 in component fast_Na (J_per_mol). C ALGBRC(67) is Ar_j000 in component fast_Na (J_per_mol). C ALGBRC(68) is v_j000 in component fast_Na (fmol_per_sec). C ALGBRC(30) is Af_j010 in component fast_Na (J_per_mol). C ALGBRC(69) is Ar_j010 in component fast_Na (J_per_mol). C ALGBRC(70) is v_j010 in component fast_Na (fmol_per_sec). C ALGBRC(31) is Af_j100 in component fast_Na (J_per_mol). C ALGBRC(71) is Ar_j100 in component fast_Na (J_per_mol). C ALGBRC(72) is v_j100 in component fast_Na (fmol_per_sec). C ALGBRC(32) is Af_j110 in component fast_Na (J_per_mol). C ALGBRC(73) is Ar_j110 in component fast_Na (J_per_mol). C ALGBRC(74) is v_j110 in component fast_Na (fmol_per_sec). C ALGBRC(33) is Af_j200 in component fast_Na (J_per_mol). C ALGBRC(75) is Ar_j200 in component fast_Na (J_per_mol). C ALGBRC(76) is v_j200 in component fast_Na (fmol_per_sec). C ALGBRC(34) is Af_j210 in component fast_Na (J_per_mol). C ALGBRC(77) is Ar_j210 in component fast_Na (J_per_mol). C ALGBRC(78) is v_j210 in component fast_Na (fmol_per_sec). C ALGBRC(35) is Af_j300 in component fast_Na (J_per_mol). C ALGBRC(79) is Ar_j300 in component fast_Na (J_per_mol). C ALGBRC(80) is v_j300 in component fast_Na (fmol_per_sec). C ALGBRC(36) is Af_j310 in component fast_Na (J_per_mol). C ALGBRC(81) is Ar_j310 in component fast_Na (J_per_mol). C ALGBRC(82) is v_j310 in component fast_Na (fmol_per_sec). C ALGBRC(37) is Af_m000 in component fast_Na (J_per_mol). C ALGBRC(83) is Ar_m000 in component fast_Na (J_per_mol). C ALGBRC(84) is v_m000 in component fast_Na (fmol_per_sec). C ALGBRC(38) is Af_m001 in component fast_Na (J_per_mol). C ALGBRC(85) is Ar_m001 in component fast_Na (J_per_mol). C ALGBRC(87) is v_m001 in component fast_Na (fmol_per_sec). C ALGBRC(39) is Af_m010 in component fast_Na (J_per_mol). C ALGBRC(88) is Ar_m010 in component fast_Na (J_per_mol). C ALGBRC(90) is v_m010 in component fast_Na (fmol_per_sec). C ALGBRC(40) is Af_m011 in component fast_Na (J_per_mol). C ALGBRC(91) is Ar_m011 in component fast_Na (J_per_mol). C ALGBRC(93) is v_m011 in component fast_Na (fmol_per_sec). C ALGBRC(41) is Af_m100 in component fast_Na (J_per_mol). C ALGBRC(94) is Ar_m100 in component fast_Na (J_per_mol). C ALGBRC(96) is v_m100 in component fast_Na (fmol_per_sec). C ALGBRC(42) is Af_m101 in component fast_Na (J_per_mol). C ALGBRC(97) is Ar_m101 in component fast_Na (J_per_mol). C ALGBRC(99) is v_m101 in component fast_Na (fmol_per_sec). C ALGBRC(43) is Af_m110 in component fast_Na (J_per_mol). C ALGBRC(100) is Ar_m110 in component fast_Na (J_per_mol). C ALGBRC(102) is v_m110 in component fast_Na (fmol_per_sec). C ALGBRC(44) is Af_m111 in component fast_Na (J_per_mol). C ALGBRC(103) is Ar_m111 in component fast_Na (J_per_mol). C ALGBRC(105) is v_m111 in component fast_Na (fmol_per_sec). C ALGBRC(45) is Af_m200 in component fast_Na (J_per_mol). C ALGBRC(106) is Ar_m200 in component fast_Na (J_per_mol). C ALGBRC(108) is v_m200 in component fast_Na (fmol_per_sec). C ALGBRC(46) is Af_m201 in component fast_Na (J_per_mol). C ALGBRC(109) is Ar_m201 in component fast_Na (J_per_mol). C ALGBRC(112) is v_m201 in component fast_Na (fmol_per_sec). C ALGBRC(47) is Af_m210 in component fast_Na (J_per_mol). C ALGBRC(113) is Ar_m210 in component fast_Na (J_per_mol). C ALGBRC(116) is v_m210 in component fast_Na (fmol_per_sec). C ALGBRC(48) is Af_m211 in component fast_Na (J_per_mol). C ALGBRC(117) is Ar_m211 in component fast_Na (J_per_mol). C ALGBRC(120) is v_m211 in component fast_Na (fmol_per_sec). C ALGBRC(86) is v_S000_Na in component fast_Na (fmol_per_sec). C ALGBRC(92) is v_S010_Na in component fast_Na (fmol_per_sec). C ALGBRC(98) is v_S100_Na in component fast_Na (fmol_per_sec). C ALGBRC(104) is v_S110_Na in component fast_Na (fmol_per_sec). C ALGBRC(110) is v_S200_Na in component fast_Na (fmol_per_sec). C ALGBRC(118) is v_S210_Na in component fast_Na (fmol_per_sec). C ALGBRC(111) is v_S300_Na in component fast_Na (fmol_per_sec). C ALGBRC(119) is v_S310_Na in component fast_Na (fmol_per_sec). C ALGBRC(89) is v_S001_Na in component fast_Na (fmol_per_sec). C ALGBRC(95) is v_S011_Na in component fast_Na (fmol_per_sec). C ALGBRC(101) is v_S101_Na in component fast_Na (fmol_per_sec). C ALGBRC(107) is v_S111_Na in component fast_Na (fmol_per_sec). C ALGBRC(114) is v_S201_Na in component fast_Na (fmol_per_sec). C ALGBRC(122) is v_S211_Na in component fast_Na (fmol_per_sec). C ALGBRC(115) is v_S301_Na in component fast_Na (fmol_per_sec). C ALGBRC(124) is v_S311_Na in component fast_Na (fmol_per_sec). C ALGBRC(1) is V_mem in component fast_Na (volt). C RATES(1) is d/dt q_Na_o in component environment (fmol). C RATES(2) is d/dt q_Na_i in component environment (fmol). C RATES(19) is d/dt q_mem in component environment (fC). C RATES(3) is d/dt q_S000_Na in component environment (fmol). C RATES(4) is d/dt q_S010_Na in component environment (fmol). C RATES(5) is d/dt q_S100_Na in component environment (fmol). C RATES(6) is d/dt q_S110_Na in component environment (fmol). C RATES(7) is d/dt q_S200_Na in component environment (fmol). C RATES(8) is d/dt q_S210_Na in component environment (fmol). C RATES(9) is d/dt q_S300_Na in component environment (fmol). C RATES(10) is d/dt q_S310_Na in component environment (fmol). C RATES(11) is d/dt q_S001_Na in component environment (fmol). C RATES(12) is d/dt q_S011_Na in component environment (fmol). C RATES(13) is d/dt q_S101_Na in component environment (fmol). C RATES(14) is d/dt q_S111_Na in component environment (fmol). C RATES(15) is d/dt q_S201_Na in component environment (fmol). C RATES(16) is d/dt q_S211_Na in component environment (fmol). C RATES(17) is d/dt q_S301_Na in component environment (fmol). C RATES(18) is d/dt q_S311_Na in component environment (fmol). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = 8.314 CONSTS(2) = 310 CONSTS(3) = 96485 CONSTS(4) = 153400 STATES(1) = 725.48 STATES(2) = 380 STATES(3) = 1.72572E-05 STATES(4) = 1.91747E-06 STATES(5) = 5.75256E-06 STATES(6) = 6.39176E-07 STATES(7) = 6.39176E-07 STATES(8) = 7.10196E-08 STATES(9) = 2.36737E-08 STATES(10) = 2.63039E-09 STATES(11) = 1.91747E-06 STATES(12) = 2.13059E-07 STATES(13) = 6.39176E-07 STATES(14) = 7.10196E-08 STATES(15) = 7.10196E-08 STATES(16) = 7.89103E-09 STATES(17) = 2.63039E-09 STATES(18) = 2.92267E-10 STATES(19) = -13039 CONSTS(5) = 0.0073376 CONSTS(6) = 23718.0222 CONSTS(7) = 17499.8767 CONSTS(8) = 56566767727.4036 CONSTS(9) = 1.5633e-05 CONSTS(10) = 50.5317 CONSTS(11) = 37.2839 CONSTS(12) = 120516694.0958 CONSTS(13) = 9.9918e-08 CONSTS(14) = 0.32298 CONSTS(15) = 0.2383 CONSTS(16) = 770290.0911 CONSTS(17) = 1.9159e-09 CONSTS(18) = 0.006193 CONSTS(19) = 0.0045694 CONSTS(20) = 14770.0739 CONSTS(21) = 0.017678 CONSTS(22) = 0.0024108 CONSTS(23) = 32.858 CONSTS(24) = 0.045709 CONSTS(25) = 1.4141e-08 CONSTS(26) = 21.4544 CONSTS(27) = 6.6373e-06 CONSTS(28) = 3356.675 CONSTS(29) = 0.0010384 CONSTS(30) = 175057.5863 CONSTS(31) = 0.054157 CONSTS(32) = 0.017824 CONSTS(33) = 7.4736e-09 CONSTS(34) = 8.3662 CONSTS(35) = 3.5079e-06 CONSTS(36) = 1308.9427 CONSTS(37) = 0.00054883 CONSTS(38) = 68264.0876 CONSTS(39) = 0.028623 CONSTS(40) = 5117405.6443 CONSTS(41) = 1.5832 CONSTS(42) = 2.1457 CONSTS(43) = 6.6381e-07 CONSTS(44) = 1601300127.6001 CONSTS(45) = 495.389 CONSTS(46) = 671.4132 CONSTS(47) = 0.00020771 CONSTS(48) = 125266701374.7363 CONSTS(49) = 38753.3515 CONSTS(50) = 52523.3901 CONSTS(51) = 0.016249 CONSTS(52) = 1 CONSTS(53) = -4.1892 CONSTS(54) = -4.0381 CONSTS(55) = 0.49541 CONSTS(56) = 1.2995 CONSTS(57) = 1.4281 CONSTS(58) = -2.4284 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = STATES(19)/CONSTS(4) ALGBRC(21) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(51) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(52) = CONSTS(24)*(EXP(ALGBRC(21)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(51)/( CONSTS(1)*CONSTS(2)))) ALGBRC(29) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(67) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(68) = CONSTS(32)*(EXP(ALGBRC(29)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(67)/( CONSTS(1)*CONSTS(2)))) ALGBRC(37) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(83) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(84) = CONSTS(40)*(EXP(ALGBRC(37)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(83)/( CONSTS(1)*CONSTS(2)))) ALGBRC(86) = (- ALGBRC(84) - ALGBRC(52)) - ALGBRC(68) RATES(3) = ALGBRC(86) ALGBRC(22) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(53) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(54) = CONSTS(25)*(EXP(ALGBRC(22)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(53)/( CONSTS(1)*CONSTS(2)))) ALGBRC(38) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(85) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(87) = CONSTS(41)*(EXP(ALGBRC(38)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(85)/( CONSTS(1)*CONSTS(2)))) ALGBRC(89) = (- ALGBRC(87) - ALGBRC(54)) - - ALGBRC(68) RATES(11) = ALGBRC(89) ALGBRC(30) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(69) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(70) = CONSTS(33)*(EXP(ALGBRC(30)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(69)/( CONSTS(1)*CONSTS(2)))) ALGBRC(39) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(88) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(90) = CONSTS(42)*(EXP(ALGBRC(39)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(88)/( CONSTS(1)*CONSTS(2)))) ALGBRC(92) = (- ALGBRC(90)+ALGBRC(52)) - ALGBRC(70) RATES(4) = ALGBRC(92) ALGBRC(40) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(91) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(93) = CONSTS(43)*(EXP(ALGBRC(40)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(91)/( CONSTS(1)*CONSTS(2)))) ALGBRC(95) = (- ALGBRC(93)+ALGBRC(54)) - - ALGBRC(70) RATES(12) = ALGBRC(95) ALGBRC(23) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(55) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(56) = CONSTS(26)*(EXP(ALGBRC(23)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(55)/( CONSTS(1)*CONSTS(2)))) ALGBRC(31) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(71) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(72) = CONSTS(34)*(EXP(ALGBRC(31)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(71)/( CONSTS(1)*CONSTS(2)))) ALGBRC(41) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(94) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(96) = CONSTS(44)*(EXP(ALGBRC(41)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(94)/( CONSTS(1)*CONSTS(2)))) ALGBRC(98) = ((ALGBRC(84) - ALGBRC(96)) - ALGBRC(56)) - ALGBRC(72) RATES(5) = ALGBRC(98) ALGBRC(24) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(57) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(58) = CONSTS(27)*(EXP(ALGBRC(24)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(57)/( CONSTS(1)*CONSTS(2)))) ALGBRC(42) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(97) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(99) = CONSTS(45)*(EXP(ALGBRC(42)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(97)/( CONSTS(1)*CONSTS(2)))) ALGBRC(101) = ((ALGBRC(87) - ALGBRC(99)) - ALGBRC(58)) - - ALGBRC(72) RATES(13) = ALGBRC(101) ALGBRC(32) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(73) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(74) = CONSTS(35)*(EXP(ALGBRC(32)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(73)/( CONSTS(1)*CONSTS(2)))) ALGBRC(43) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(100) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(102) = CONSTS(46)*(EXP(ALGBRC(43)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(100)/( CONSTS(1)*CONSTS(2)))) ALGBRC(104) = ((ALGBRC(90) - ALGBRC(102))+ALGBRC(56)) - ALGBRC(74) RATES(6) = ALGBRC(104) ALGBRC(44) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(103) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(105) = CONSTS(47)*(EXP(ALGBRC(44)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(103)/( CONSTS(1)*CONSTS(2)))) ALGBRC(107) = ((ALGBRC(93) - ALGBRC(105))+ALGBRC(58)) - - ALGBRC(74) RATES(14) = ALGBRC(107) ALGBRC(25) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(59) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(60) = CONSTS(28)*(EXP(ALGBRC(25)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(59)/( CONSTS(1)*CONSTS(2)))) ALGBRC(33) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(75) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(76) = CONSTS(36)*(EXP(ALGBRC(33)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(75)/( CONSTS(1)*CONSTS(2)))) ALGBRC(45) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(106) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(108) = CONSTS(48)*(EXP(ALGBRC(45)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(106)/( CONSTS(1)*CONSTS(2)))) ALGBRC(110) = ((ALGBRC(96) - ALGBRC(108)) - ALGBRC(60)) - ALGBRC(76) RATES(7) = ALGBRC(110) ALGBRC(27) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(63) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(64) = CONSTS(30)*(EXP(ALGBRC(27)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(63)/( CONSTS(1)*CONSTS(2)))) ALGBRC(35) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(79) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(80) = CONSTS(38)*(EXP(ALGBRC(35)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(79)/( CONSTS(1)*CONSTS(2)))) ALGBRC(111) = (ALGBRC(108) - ALGBRC(64)) - ALGBRC(80) RATES(9) = ALGBRC(111) ALGBRC(26) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(61) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(62) = CONSTS(29)*(EXP(ALGBRC(26)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(61)/( CONSTS(1)*CONSTS(2)))) ALGBRC(46) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(109) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(112) = CONSTS(49)*(EXP(ALGBRC(46)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(109)/( CONSTS(1)*CONSTS(2)))) ALGBRC(114) = ((ALGBRC(99) - ALGBRC(112)) - ALGBRC(62)) - - ALGBRC(76) RATES(15) = ALGBRC(114) ALGBRC(28) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(65) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(66) = CONSTS(31)*(EXP(ALGBRC(28)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(65)/( CONSTS(1)*CONSTS(2)))) ALGBRC(115) = (ALGBRC(112) - ALGBRC(66)) - - ALGBRC(80) RATES(17) = ALGBRC(115) ALGBRC(34) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(77) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(78) = CONSTS(37)*(EXP(ALGBRC(34)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(77)/( CONSTS(1)*CONSTS(2)))) ALGBRC(47) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(113) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(116) = CONSTS(50)*(EXP(ALGBRC(47)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(113)/( CONSTS(1)*CONSTS(2)))) ALGBRC(118) = ((ALGBRC(102) - ALGBRC(116))+ALGBRC(60)) - ALGBRC(78) RATES(8) = ALGBRC(118) ALGBRC(36) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(81) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(82) = CONSTS(39)*(EXP(ALGBRC(36)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(81)/( CONSTS(1)*CONSTS(2)))) ALGBRC(119) = (ALGBRC(116)+ALGBRC(64)) - ALGBRC(82) RATES(10) = ALGBRC(119) ALGBRC(49) = CONSTS(52)*CONSTS(3)*ALGBRC(1) ALGBRC(18) = CONSTS(1)*CONSTS(2)*log( CONSTS(22)*STATES(2)) ALGBRC(19) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18)) ALGBRC(20) = ALGBRC(18)+ CONSTS(52)*CONSTS(3)*ALGBRC(1)+ALGBRC(19) ALGBRC(17) = CONSTS(1)*CONSTS(2)*log( CONSTS(21)*STATES(1)) ALGBRC(50) = ALGBRC(17)+ALGBRC(19) ALGBRC(121) = TERNRY(ALGBRC(49).EQ.0.00000, CONSTS(23)*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2)))), ((( CONSTS(23)*ALGBRC(49))/( CONSTS(1)*CONSTS(2)))/(EXP(ALGBRC(49)/( CONSTS(1)*CONSTS(2))) - 1.00000))*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2))))) RATES(1) = ALGBRC(121) RATES(2) = - ALGBRC(121) ALGBRC(48) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(117) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(120) = CONSTS(51)*(EXP(ALGBRC(48)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(117)/( CONSTS(1)*CONSTS(2)))) ALGBRC(122) = ((ALGBRC(105) - ALGBRC(120))+ALGBRC(62)) - - ALGBRC(78) RATES(16) = ALGBRC(122) ALGBRC(123) = CONSTS(3)*(((((((((((((((((((((((((((((((((((((((((((((((((((((((( - CONSTS(52)*ALGBRC(121)+ CONSTS(58)*ALGBRC(84)) - CONSTS(55)*ALGBRC(84))+ CONSTS(58)*ALGBRC(96)) - CONSTS(55)*ALGBRC(96))+ CONSTS(58)*ALGBRC(90)) - CONSTS(55)*ALGBRC(90))+ CONSTS(58)*ALGBRC(102)) - CONSTS(55)*ALGBRC(102))+ CONSTS(58)*ALGBRC(116)) - CONSTS(55)*ALGBRC(116))+ CONSTS(56)*ALGBRC(52)) - CONSTS(53)*ALGBRC(52))+ CONSTS(56)*ALGBRC(60)) - CONSTS(53)*ALGBRC(60))+ CONSTS(57)*ALGBRC(72)) - CONSTS(54)*ALGBRC(72))+ CONSTS(57)*ALGBRC(76)) - CONSTS(54)*ALGBRC(76))+ CONSTS(57)*ALGBRC(80)) - CONSTS(54)*ALGBRC(80))+ CONSTS(58)*ALGBRC(108)) - CONSTS(55)*ALGBRC(108))+ CONSTS(56)*ALGBRC(56)) - CONSTS(53)*ALGBRC(56))+ CONSTS(56)*ALGBRC(64)) - CONSTS(53)*ALGBRC(64))+ CONSTS(57)*ALGBRC(68)) - CONSTS(54)*ALGBRC(68))+ CONSTS(58)*ALGBRC(87)) - CONSTS(55)*ALGBRC(87))+ CONSTS(58)*ALGBRC(99)) - CONSTS(55)*ALGBRC(99))+ CONSTS(58)*ALGBRC(93)) - CONSTS(55)*ALGBRC(93))+ CONSTS(58)*ALGBRC(105)) - CONSTS(55)*ALGBRC(105))+ CONSTS(58)*ALGBRC(120)) - CONSTS(55)*ALGBRC(120))+ CONSTS(56)*ALGBRC(54)) - CONSTS(53)*ALGBRC(54))+ CONSTS(56)*ALGBRC(62)) - CONSTS(53)*ALGBRC(62))+ CONSTS(57)*ALGBRC(70)) - CONSTS(54)*ALGBRC(70))+ CONSTS(57)*ALGBRC(74)) - CONSTS(54)*ALGBRC(74))+ CONSTS(57)*ALGBRC(78)) - CONSTS(54)*ALGBRC(78))+ CONSTS(57)*ALGBRC(82)) - CONSTS(54)*ALGBRC(82))+ CONSTS(58)*ALGBRC(112)) - CONSTS(55)*ALGBRC(112))+ CONSTS(56)*ALGBRC(58)) - CONSTS(53)*ALGBRC(58))+ CONSTS(56)*ALGBRC(66)) - CONSTS(53)*ALGBRC(66)) RATES(19) = ALGBRC(123) ALGBRC(124) = (ALGBRC(120)+ALGBRC(66)) - ((- ALGBRC(82)+ALGBRC(121)) - ALGBRC(121)) RATES(18) = ALGBRC(124) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = STATES(19)/CONSTS(4) ALGBRC(21) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(51) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(52) = CONSTS(24)*(EXP(ALGBRC(21)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(51)/( CONSTS(1)*CONSTS(2)))) ALGBRC(29) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(67) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(68) = CONSTS(32)*(EXP(ALGBRC(29)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(67)/( CONSTS(1)*CONSTS(2)))) ALGBRC(37) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(83) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(84) = CONSTS(40)*(EXP(ALGBRC(37)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(83)/( CONSTS(1)*CONSTS(2)))) ALGBRC(86) = (- ALGBRC(84) - ALGBRC(52)) - ALGBRC(68) ALGBRC(22) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(53) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(54) = CONSTS(25)*(EXP(ALGBRC(22)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(53)/( CONSTS(1)*CONSTS(2)))) ALGBRC(38) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(85) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(87) = CONSTS(41)*(EXP(ALGBRC(38)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(85)/( CONSTS(1)*CONSTS(2)))) ALGBRC(89) = (- ALGBRC(87) - ALGBRC(54)) - - ALGBRC(68) ALGBRC(30) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(69) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(70) = CONSTS(33)*(EXP(ALGBRC(30)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(69)/( CONSTS(1)*CONSTS(2)))) ALGBRC(39) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(88) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(90) = CONSTS(42)*(EXP(ALGBRC(39)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(88)/( CONSTS(1)*CONSTS(2)))) ALGBRC(92) = (- ALGBRC(90)+ALGBRC(52)) - ALGBRC(70) ALGBRC(40) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(91) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(93) = CONSTS(43)*(EXP(ALGBRC(40)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(91)/( CONSTS(1)*CONSTS(2)))) ALGBRC(95) = (- ALGBRC(93)+ALGBRC(54)) - - ALGBRC(70) ALGBRC(23) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(55) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(56) = CONSTS(26)*(EXP(ALGBRC(23)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(55)/( CONSTS(1)*CONSTS(2)))) ALGBRC(31) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(71) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(72) = CONSTS(34)*(EXP(ALGBRC(31)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(71)/( CONSTS(1)*CONSTS(2)))) ALGBRC(41) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(94) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(96) = CONSTS(44)*(EXP(ALGBRC(41)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(94)/( CONSTS(1)*CONSTS(2)))) ALGBRC(98) = ((ALGBRC(84) - ALGBRC(96)) - ALGBRC(56)) - ALGBRC(72) ALGBRC(24) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(57) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(58) = CONSTS(27)*(EXP(ALGBRC(24)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(57)/( CONSTS(1)*CONSTS(2)))) ALGBRC(42) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(97) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(99) = CONSTS(45)*(EXP(ALGBRC(42)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(97)/( CONSTS(1)*CONSTS(2)))) ALGBRC(101) = ((ALGBRC(87) - ALGBRC(99)) - ALGBRC(58)) - - ALGBRC(72) ALGBRC(32) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(73) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(74) = CONSTS(35)*(EXP(ALGBRC(32)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(73)/( CONSTS(1)*CONSTS(2)))) ALGBRC(43) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(100) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(102) = CONSTS(46)*(EXP(ALGBRC(43)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(100)/( CONSTS(1)*CONSTS(2)))) ALGBRC(104) = ((ALGBRC(90) - ALGBRC(102))+ALGBRC(56)) - ALGBRC(74) ALGBRC(44) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(103) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(105) = CONSTS(47)*(EXP(ALGBRC(44)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(103)/( CONSTS(1)*CONSTS(2)))) ALGBRC(107) = ((ALGBRC(93) - ALGBRC(105))+ALGBRC(58)) - - ALGBRC(74) ALGBRC(25) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(59) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(60) = CONSTS(28)*(EXP(ALGBRC(25)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(59)/( CONSTS(1)*CONSTS(2)))) ALGBRC(33) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(75) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(76) = CONSTS(36)*(EXP(ALGBRC(33)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(75)/( CONSTS(1)*CONSTS(2)))) ALGBRC(45) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(106) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(108) = CONSTS(48)*(EXP(ALGBRC(45)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(106)/( CONSTS(1)*CONSTS(2)))) ALGBRC(110) = ((ALGBRC(96) - ALGBRC(108)) - ALGBRC(60)) - ALGBRC(76) ALGBRC(27) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(63) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(64) = CONSTS(30)*(EXP(ALGBRC(27)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(63)/( CONSTS(1)*CONSTS(2)))) ALGBRC(35) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(79) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(80) = CONSTS(38)*(EXP(ALGBRC(35)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(79)/( CONSTS(1)*CONSTS(2)))) ALGBRC(111) = (ALGBRC(108) - ALGBRC(64)) - ALGBRC(80) ALGBRC(26) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(61) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(62) = CONSTS(29)*(EXP(ALGBRC(26)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(61)/( CONSTS(1)*CONSTS(2)))) ALGBRC(46) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(109) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(112) = CONSTS(49)*(EXP(ALGBRC(46)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(109)/( CONSTS(1)*CONSTS(2)))) ALGBRC(114) = ((ALGBRC(99) - ALGBRC(112)) - ALGBRC(62)) - - ALGBRC(76) ALGBRC(28) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17))+ CONSTS(53)*CONSTS(3)*ALGBRC(1) ALGBRC(65) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(56)*CONSTS(3)*ALGBRC(1) ALGBRC(66) = CONSTS(31)*(EXP(ALGBRC(28)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(65)/( CONSTS(1)*CONSTS(2)))) ALGBRC(115) = (ALGBRC(112) - ALGBRC(66)) - - ALGBRC(80) ALGBRC(34) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(77) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(78) = CONSTS(37)*(EXP(ALGBRC(34)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(77)/( CONSTS(1)*CONSTS(2)))) ALGBRC(47) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(113) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(116) = CONSTS(50)*(EXP(ALGBRC(47)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(113)/( CONSTS(1)*CONSTS(2)))) ALGBRC(118) = ((ALGBRC(102) - ALGBRC(116))+ALGBRC(60)) - ALGBRC(78) ALGBRC(36) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10))+ CONSTS(54)*CONSTS(3)*ALGBRC(1) ALGBRC(81) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(57)*CONSTS(3)*ALGBRC(1) ALGBRC(82) = CONSTS(39)*(EXP(ALGBRC(36)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(81)/( CONSTS(1)*CONSTS(2)))) ALGBRC(119) = (ALGBRC(116)+ALGBRC(64)) - ALGBRC(82) ALGBRC(49) = CONSTS(52)*CONSTS(3)*ALGBRC(1) ALGBRC(18) = CONSTS(1)*CONSTS(2)*log( CONSTS(22)*STATES(2)) ALGBRC(19) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18)) ALGBRC(20) = ALGBRC(18)+ CONSTS(52)*CONSTS(3)*ALGBRC(1)+ALGBRC(19) ALGBRC(17) = CONSTS(1)*CONSTS(2)*log( CONSTS(21)*STATES(1)) ALGBRC(50) = ALGBRC(17)+ALGBRC(19) ALGBRC(121) = TERNRY(ALGBRC(49).EQ.0.00000, CONSTS(23)*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2)))), ((( CONSTS(23)*ALGBRC(49))/( CONSTS(1)*CONSTS(2)))/(EXP(ALGBRC(49)/( CONSTS(1)*CONSTS(2))) - 1.00000))*(EXP(ALGBRC(20)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(50)/( CONSTS(1)*CONSTS(2))))) ALGBRC(48) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16))+ CONSTS(55)*CONSTS(3)*ALGBRC(1) ALGBRC(117) = CONSTS(1)*CONSTS(2)*log( CONSTS(20)*STATES(18))+ CONSTS(58)*CONSTS(3)*ALGBRC(1) ALGBRC(120) = CONSTS(51)*(EXP(ALGBRC(48)/( CONSTS(1)*CONSTS(2))) - EXP(ALGBRC(117)/( CONSTS(1)*CONSTS(2)))) ALGBRC(122) = ((ALGBRC(105) - ALGBRC(120))+ALGBRC(62)) - - ALGBRC(78) ALGBRC(123) = CONSTS(3)*(((((((((((((((((((((((((((((((((((((((((((((((((((((((( - CONSTS(52)*ALGBRC(121)+ CONSTS(58)*ALGBRC(84)) - CONSTS(55)*ALGBRC(84))+ CONSTS(58)*ALGBRC(96)) - CONSTS(55)*ALGBRC(96))+ CONSTS(58)*ALGBRC(90)) - CONSTS(55)*ALGBRC(90))+ CONSTS(58)*ALGBRC(102)) - CONSTS(55)*ALGBRC(102))+ CONSTS(58)*ALGBRC(116)) - CONSTS(55)*ALGBRC(116))+ CONSTS(56)*ALGBRC(52)) - CONSTS(53)*ALGBRC(52))+ CONSTS(56)*ALGBRC(60)) - CONSTS(53)*ALGBRC(60))+ CONSTS(57)*ALGBRC(72)) - CONSTS(54)*ALGBRC(72))+ CONSTS(57)*ALGBRC(76)) - CONSTS(54)*ALGBRC(76))+ CONSTS(57)*ALGBRC(80)) - CONSTS(54)*ALGBRC(80))+ CONSTS(58)*ALGBRC(108)) - CONSTS(55)*ALGBRC(108))+ CONSTS(56)*ALGBRC(56)) - CONSTS(53)*ALGBRC(56))+ CONSTS(56)*ALGBRC(64)) - CONSTS(53)*ALGBRC(64))+ CONSTS(57)*ALGBRC(68)) - CONSTS(54)*ALGBRC(68))+ CONSTS(58)*ALGBRC(87)) - CONSTS(55)*ALGBRC(87))+ CONSTS(58)*ALGBRC(99)) - CONSTS(55)*ALGBRC(99))+ CONSTS(58)*ALGBRC(93)) - CONSTS(55)*ALGBRC(93))+ CONSTS(58)*ALGBRC(105)) - CONSTS(55)*ALGBRC(105))+ CONSTS(58)*ALGBRC(120)) - CONSTS(55)*ALGBRC(120))+ CONSTS(56)*ALGBRC(54)) - CONSTS(53)*ALGBRC(54))+ CONSTS(56)*ALGBRC(62)) - CONSTS(53)*ALGBRC(62))+ CONSTS(57)*ALGBRC(70)) - CONSTS(54)*ALGBRC(70))+ CONSTS(57)*ALGBRC(74)) - CONSTS(54)*ALGBRC(74))+ CONSTS(57)*ALGBRC(78)) - CONSTS(54)*ALGBRC(78))+ CONSTS(57)*ALGBRC(82)) - CONSTS(54)*ALGBRC(82))+ CONSTS(58)*ALGBRC(112)) - CONSTS(55)*ALGBRC(112))+ CONSTS(56)*ALGBRC(58)) - CONSTS(53)*ALGBRC(58))+ CONSTS(56)*ALGBRC(66)) - CONSTS(53)*ALGBRC(66)) ALGBRC(124) = (ALGBRC(120)+ALGBRC(66)) - ((- ALGBRC(82)+ALGBRC(121)) - ALGBRC(121)) ALGBRC(2) = CONSTS(1)*CONSTS(2)*log( CONSTS(5)*STATES(3)) ALGBRC(3) = CONSTS(1)*CONSTS(2)*log( CONSTS(7)*STATES(4)) ALGBRC(4) = CONSTS(1)*CONSTS(2)*log( CONSTS(9)*STATES(5)) ALGBRC(5) = CONSTS(1)*CONSTS(2)*log( CONSTS(11)*STATES(6)) ALGBRC(6) = CONSTS(1)*CONSTS(2)*log( CONSTS(13)*STATES(7)) ALGBRC(7) = CONSTS(1)*CONSTS(2)*log( CONSTS(15)*STATES(8)) ALGBRC(8) = CONSTS(1)*CONSTS(2)*log( CONSTS(17)*STATES(9)) ALGBRC(9) = CONSTS(1)*CONSTS(2)*log( CONSTS(19)*STATES(10)) ALGBRC(10) = CONSTS(1)*CONSTS(2)*log( CONSTS(6)*STATES(11)) ALGBRC(11) = CONSTS(1)*CONSTS(2)*log( CONSTS(8)*STATES(12)) ALGBRC(12) = CONSTS(1)*CONSTS(2)*log( CONSTS(10)*STATES(13)) ALGBRC(13) = CONSTS(1)*CONSTS(2)*log( CONSTS(12)*STATES(14)) ALGBRC(14) = CONSTS(1)*CONSTS(2)*log( CONSTS(14)*STATES(15)) ALGBRC(15) = CONSTS(1)*CONSTS(2)*log( CONSTS(16)*STATES(16)) ALGBRC(16) = CONSTS(1)*CONSTS(2)*log( CONSTS(18)*STATES(17)) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END