C
C There are a total of 98 entries in the algebraic variable array.
C There are a total of 30 entries in each of the rate and state variable arrays.
C There are a total of 108 entries in the constant variable array.
C
C
C VOI is t in component lpc (s).
C CONSTS(1) is Vlvrd in component lpc (mL).
C CONSTS(2) is Vlvrs in component lpc (mL).
C CONSTS(3) is Vrvrd in component lpc (mL).
C CONSTS(4) is Vrvrs in component lpc (mL).
C CONSTS(5) is Vlard in component lpc (mL).
C CONSTS(6) is Vlars in component lpc (mL).
C CONSTS(7) is Vrard in component lpc (mL).
C CONSTS(8) is Vrars in component lpc (mL).
C CONSTS(9) is Rra in component lpc (resistance).
C CONSTS(10) is Rla in component lpc (resistance).
C CONSTS(11) is PRint in component lpc (s).
C CONSTS(12) is Emaxlv in component lpc (elastance).
C CONSTS(13) is Eminlv in component lpc (elastance).
C CONSTS(14) is Emaxrv in component lpc (elastance).
C CONSTS(15) is Eminrv in component lpc (elastance).
C CONSTS(16) is Emaxra in component lpc (elastance).
C CONSTS(17) is Eminra in component lpc (elastance).
C CONSTS(18) is Emaxla in component lpc (elastance).
C CONSTS(19) is Eminla in component lpc (elastance).
C CONSTS(20) is Pbs in component lpc (mmHg).
C CONSTS(21) is Vmyo in component lpc (mL).
C CONSTS(22) is TsvK in component lpc (s).
C CONSTS(23) is TsaK in component lpc (s).
C CONSTS(24) is HR in component lpc (ratepm).
C ALGBRC(18) is Elv in component lpc (elastance).
C ALGBRC(16) is Erv in component lpc (elastance).
C ALGBRC(15) is Era in component lpc (elastance).
C ALGBRC(17) is Ela in component lpc (elastance).
C ALGBRC(59) is Vlvr in component lpc (mL).
C ALGBRC(82) is Vrvr in component lpc (mL).
C ALGBRC(56) is Vlar in component lpc (mL).
C ALGBRC(71) is Vrar in component lpc (mL).
C ALGBRC(89) is Fra in component lpc (flow).
C ALGBRC(95) is Frv in component lpc (flow).
C ALGBRC(81) is Fla in component lpc (flow).
C ALGBRC(84) is Flv in component lpc (flow).
C ALGBRC(73) is Pra in component lpc (mmHg).
C ALGBRC(75) is Prac in component lpc (mmHg).
C ALGBRC(85) is Prv in component lpc (mmHg).
C ALGBRC(87) is Prvc in component lpc (mmHg).
C ALGBRC(58) is Pla in component lpc (mmHg).
C ALGBRC(76) is Plac in component lpc (mmHg).
C ALGBRC(60) is Plv in component lpc (mmHg).
C ALGBRC(80) is Plvc in component lpc (mmHg).
C STATES(1) is Vra in component lpc (mL).
C STATES(2) is Vrv in component lpc (mL).
C STATES(3) is Vla in component lpc (mL).
C STATES(4) is Vlv in component lpc (mL).
C CONSTS(25) is COutput in component lpc (flow).
C CONSTS(104) is SV in component lpc (mL).
C ALGBRC(7) is P_QRSwave in component lpc (mV).
C CONSTS(106) is Tsv in component lpc (s).
C CONSTS(108) is Tsa in component lpc (s).
C ALGBRC(10) is trela in component lpc (s).
C ALGBRC(8) is trelv in component lpc (s).
C CONSTS(26) is Rav in component lpc (resistance).
C CONSTS(27) is Raop in component lpc (resistance).
C CONSTS(28) is Rcrb in component lpc (resistance).
C CONSTS(29) is Raod in component lpc (resistance).
C CONSTS(30) is Rtaop in component lpc (resistance).
C CONSTS(31) is Rtaod in component lpc (resistance).
C CONSTS(32) is Rsap in component lpc (resistance).
C CONSTS(33) is Rsc in component lpc (resistance).
C CONSTS(34) is Rsv in component lpc (resistance).
C CONSTS(35) is Rsao in component lpc (resistance).
C CONSTS(36) is Caop in component lpc (capacitance).
C CONSTS(37) is Caod in component lpc (capacitance).
C CONSTS(38) is Csap in component lpc (capacitance).
C CONSTS(39) is Csc in component lpc (capacitance).
C CONSTS(40) is Laop in component lpc (inductance).
C CONSTS(41) is Laod in component lpc (inductance).
C CONSTS(42) is Kc in component lpc (mmHg).
C CONSTS(43) is Do in component lpc (mL).
C CONSTS(44) is Vsa_o in component lpc (mL).
C CONSTS(45) is Vsa_max in component lpc (mL).
C CONSTS(46) is Kp1 in component lpc (mmHg).
C CONSTS(47) is Kp2 in component lpc (mmHgmlm2).
C CONSTS(48) is Kr in component lpc (resistance).
C CONSTS(49) is tau_p in component lpc (perml).
C CONSTS(50) is Kv in component lpc (mmHg).
C CONSTS(51) is Vmax_sv in component lpc (mL).
C CONSTS(52) is D2 in component lpc (mmHg).
C CONSTS(53) is K1 in component lpc (dimensionless).
C CONSTS(54) is K2 in component lpc (mmHg).
C CONSTS(55) is KR in component lpc (resistance).
C CONSTS(56) is Ro in component lpc (resistance).
C CONSTS(57) is Vo in component lpc (mL).
C CONSTS(58) is Vmax_vc in component lpc (mL).
C CONSTS(59) is Vmin_vc in component lpc (mL).
C CONSTS(60) is COtau in component lpc (s).
C CONSTS(61) is Pplc in component lpc (mmHg).
C CONSTS(62) is Px2 in component lpc (mmHg).
C CONSTS(63) is Vx8 in component lpc (mL).
C CONSTS(64) is Vx75 in component lpc (mL).
C CONSTS(65) is Vx1 in component lpc (mL).
C CONSTS(66) is Px1 in component lpc (mmHg).
C CONSTS(67) is F_vaso in component lpc (dimensionless).
C ALGBRC(53) is Rvc in component lpc (resistance).
C ALGBRC(45) is Rsa in component lpc (resistance).
C ALGBRC(1) is Paop in component lpc (mmHg).
C STATES(5) is Paopc in component lpc (mmHg).
C ALGBRC(51) is Paod in component lpc (mmHg).
C ALGBRC(50) is Paodc in component lpc (mmHg).
C ALGBRC(37) is Psa_a in component lpc (mmHg).
C ALGBRC(38) is Psa_p in component lpc (mmHg).
C ALGBRC(39) is Psa in component lpc (mmHg).
C ALGBRC(40) is Psap in component lpc (mmHg).
C ALGBRC(42) is Psc in component lpc (mmHg).
C ALGBRC(43) is Psv in component lpc (mmHg).
C ALGBRC(46) is Pvc in component lpc (mmHg).
C ALGBRC(48) is Pvcc in component lpc (mmHg).
C STATES(6) is MAP in component lpc (mmHg).
C STATES(7) is Faop in component lpc (flowLm).
C STATES(8) is Faod in component lpc (flow).
C ALGBRC(41) is Fsap in component lpc (flow).
C ALGBRC(47) is Fsa in component lpc (flow).
C ALGBRC(44) is Fsc in component lpc (flow).
C ALGBRC(49) is Fsv in component lpc (flow).
C ALGBRC(79) is Fvc in component lpc (flow).
C ALGBRC(52) is Fcrb in component lpc (flow).
C STATES(9) is Vaop in component lpc (mL).
C STATES(10) is Vaod in component lpc (mL).
C STATES(11) is Vsa in component lpc (mL).
C STATES(12) is Vsap in component lpc (mL).
C STATES(13) is Vsc in component lpc (mL).
C STATES(14) is Vsv in component lpc (mL).
C STATES(15) is Vvc in component lpc (mL).
C ALGBRC(55) is Vtot in component lpc (mL).
C ALGBRC(2) is SysArtVol in component lpc (mL).
C ALGBRC(3) is SysVenVol in component lpc (mL).
C ALGBRC(4) is PulArtVol in component lpc (mL).
C ALGBRC(5) is PulVenVol in component lpc (mL).
C ALGBRC(57) is VBcirc in component lpc (mL).
C CONSTS(105) is Ppl in component lpc (mmHg).
C CONSTS(68) is Rpuv in component lpc (resistance).
C CONSTS(69) is Rtpap in component lpc (resistance).
C CONSTS(70) is Rtpad in component lpc (resistance).
C CONSTS(71) is Rpap in component lpc (resistance).
C CONSTS(72) is Rpad in component lpc (resistance).
C CONSTS(73) is Rps in component lpc (resistance).
C CONSTS(74) is Rpa in component lpc (resistance).
C CONSTS(75) is Rpc in component lpc (resistance).
C CONSTS(76) is Rpv in component lpc (resistance).
C CONSTS(77) is Ctpap in component lpc (capacitance).
C CONSTS(78) is Ctpad in component lpc (capacitance).
C CONSTS(79) is Cpa in component lpc (capacitance).
C CONSTS(80) is Cpc in component lpc (capacitance).
C CONSTS(81) is Cpv in component lpc (capacitance).
C CONSTS(82) is Lpa in component lpc (inductance).
C CONSTS(83) is Lpad in component lpc (inductance).
C ALGBRC(93) is Ppapc in component lpc (mmHg).
C ALGBRC(91) is Ppapc1 in component lpc (mmHg).
C ALGBRC(29) is Ppapc2 in component lpc (mmHg).
C ALGBRC(96) is Ppap in component lpc (mmHg).
C ALGBRC(62) is Ppad in component lpc (mmHg).
C ALGBRC(61) is Ppadc in component lpc (mmHg).
C ALGBRC(27) is Ppa in component lpc (mmHg).
C ALGBRC(28) is Ppac in component lpc (mmHg).
C ALGBRC(30) is Ppc in component lpc (mmHg).
C ALGBRC(31) is Ppcc in component lpc (mmHg).
C ALGBRC(33) is Ppv in component lpc (mmHg).
C ALGBRC(34) is Ppvc in component lpc (mmHg).
C STATES(16) is Vpap in component lpc (mL).
C STATES(17) is Vpad in component lpc (mL).
C STATES(18) is Vpa in component lpc (mL).
C STATES(19) is Vpc in component lpc (mL).
C STATES(20) is Vpv in component lpc (mL).
C STATES(21) is Fpap in component lpc (flowLm).
C STATES(22) is Fpad in component lpc (flow).
C ALGBRC(36) is Fps in component lpc (flow).
C ALGBRC(32) is Fpa in component lpc (flow).
C ALGBRC(35) is Fpc in component lpc (flow).
C ALGBRC(78) is Fpv in component lpc (flow).
C CONSTS(84) is K_pcd in component lpc (mmHg).
C CONSTS(85) is phi_pcd in component lpc (mL).
C CONSTS(86) is Vpcd_o in component lpc (mL).
C CONSTS(87) is perifl in component lpc (mL).
C ALGBRC(64) is Ppcd in component lpc (mmHg).
C ALGBRC(65) is Ppcdc in component lpc (mmHg).
C ALGBRC(63) is Vpcd in component lpc (mL).
C CONSTS(88) is Rcorao in component lpc (resistance).
C CONSTS(89) is Rcorea in component lpc (resistance).
C CONSTS(90) is Rcorla in component lpc (resistance).
C CONSTS(91) is Rcorsa in component lpc (resistance).
C CONSTS(92) is Rcorcap in component lpc (resistance).
C CONSTS(93) is Rcorsv in component lpc (resistance).
C CONSTS(94) is Rcorlv in component lpc (resistance).
C CONSTS(95) is Rcorev in component lpc (resistance).
C CONSTS(96) is Ccorao in component lpc (capacitance).
C CONSTS(97) is Ccorea in component lpc (capacitance).
C CONSTS(98) is Ccorla in component lpc (capacitance).
C CONSTS(99) is Ccorsa in component lpc (capacitance).
C CONSTS(100) is Ccorcap in component lpc (capacitance).
C CONSTS(101) is Ccorsv in component lpc (capacitance).
C CONSTS(102) is Ccorlv in component lpc (capacitance).
C CONSTS(103) is Ccorev in component lpc (capacitance).
C ALGBRC(83) is Pcorisfc in component lpc (mmHg).
C ALGBRC(19) is Pcoraoc in component lpc (mmHg).
C ALGBRC(66) is Pcoreac in component lpc (mmHg).
C ALGBRC(90) is Pcorlac in component lpc (mmHg).
C ALGBRC(94) is Pcorsac in component lpc (mmHg).
C ALGBRC(86) is Pcorcapc in component lpc (mmHg).
C ALGBRC(72) is Pcorsvc in component lpc (mmHg).
C ALGBRC(69) is Pcorlvc in component lpc (mmHg).
C ALGBRC(68) is Pcorevc in component lpc (mmHg).
C ALGBRC(9) is Pcorao in component lpc (mmHg).
C ALGBRC(21) is Pcorea in component lpc (mmHg).
C ALGBRC(23) is Pcorla in component lpc (mmHg).
C ALGBRC(25) is Pcorsa in component lpc (mmHg).
C ALGBRC(20) is Pcorcap in component lpc (mmHg).
C ALGBRC(26) is Pcorsv in component lpc (mmHg).
C ALGBRC(24) is Pcorlv in component lpc (mmHg).
C ALGBRC(22) is Pcorev in component lpc (mmHg).
C STATES(23) is Vcorao in component lpc (mL).
C STATES(24) is Vcorea in component lpc (mL).
C STATES(25) is Vcorla in component lpc (mL).
C STATES(26) is Vcorsa in component lpc (mL).
C STATES(27) is Vcorcap in component lpc (mL).
C STATES(28) is Vcorsv in component lpc (mL).
C STATES(29) is Vcorlv in component lpc (mL).
C STATES(30) is Vcorev in component lpc (mL).
C ALGBRC(54) is Vcorcirc in component lpc (mL).
C ALGBRC(67) is Fcorao in component lpc (flow).
C ALGBRC(92) is Fcorea in component lpc (flow).
C ALGBRC(97) is Fcorla in component lpc (flow).
C ALGBRC(98) is Fcorsa in component lpc (flow).
C ALGBRC(88) is Fcorcap in component lpc (flow).
C ALGBRC(74) is Fcorsv in component lpc (flow).
C ALGBRC(70) is Fcorlv in component lpc (flow).
C ALGBRC(77) is Fcorev in component lpc (flow).
C ALGBRC(6) is beattime in component lpc (s).
C ALGBRC(11) is yla in component lpc (dimensionless).
C ALGBRC(12) is yra in component lpc (dimensionless).
C ALGBRC(13) is ylv in component lpc (dimensionless).
C ALGBRC(14) is yrv in component lpc (dimensionless).
C CONSTS(107) is hrf in component lpc (Hz).
C RATES(8) is d/dt Faod in component lpc (flow).
C RATES(7) is d/dt Faop in component lpc (flowLm).
C RATES(22) is d/dt Fpad in component lpc (flow).
C RATES(21) is d/dt Fpap in component lpc (flowLm).
C RATES(6) is d/dt MAP in component lpc (mmHg).
C RATES(5) is d/dt Paopc in component lpc (mmHg).
C RATES(9) is d/dt Vaop in component lpc (mL).
C RATES(17) is d/dt Vpad in component lpc (mL).
C RATES(10) is d/dt Vaod in component lpc (mL).
C RATES(23) is d/dt Vcorao in component lpc (mL).
C RATES(27) is d/dt Vcorcap in component lpc (mL).
C RATES(24) is d/dt Vcorea in component lpc (mL).
C RATES(30) is d/dt Vcorev in component lpc (mL).
C RATES(25) is d/dt Vcorla in component lpc (mL).
C RATES(29) is d/dt Vcorlv in component lpc (mL).
C RATES(26) is d/dt Vcorsa in component lpc (mL).
C RATES(28) is d/dt Vcorsv in component lpc (mL).
C RATES(3) is d/dt Vla in component lpc (mL).
C RATES(4) is d/dt Vlv in component lpc (mL).
C RATES(18) is d/dt Vpa in component lpc (mL).
C RATES(16) is d/dt Vpap in component lpc (mL).
C RATES(19) is d/dt Vpc in component lpc (mL).
C RATES(20) is d/dt Vpv in component lpc (mL).
C RATES(1) is d/dt Vra in component lpc (mL).
C RATES(2) is d/dt Vrv in component lpc (mL).
C RATES(11) is d/dt Vsa in component lpc (mL).
C RATES(12) is d/dt Vsap in component lpc (mL).
C RATES(13) is d/dt Vsc in component lpc (mL).
C RATES(14) is d/dt Vsv in component lpc (mL).
C RATES(15) is d/dt Vvc in component lpc (mL).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 72
      CONSTS(2) = 23
      CONSTS(3) = 103
      CONSTS(4) = 53
      CONSTS(5) = 10
      CONSTS(6) = 8
      CONSTS(7) = 10
      CONSTS(8) = 8
      CONSTS(9) = 0.001
      CONSTS(10) = 0.001
      CONSTS(11) = 0.12
      CONSTS(12) = 5.6
      CONSTS(13) = 0.186874659
      CONSTS(14) = 0.67
      CONSTS(15) = 0.1041640922
      CONSTS(16) = 0.1091675077
      CONSTS(17) = 0.0992431888
      CONSTS(18) = 0.1446191772
      CONSTS(19) = 0.1314719793
      CONSTS(20) = 0
      CONSTS(21) = 238
      CONSTS(22) = 0.35
      CONSTS(23) = 0.2
      CONSTS(24) = 77
      STATES(1) = 78.2537
      STATES(2) = 167.4806
      STATES(3) = 85.9126
      STATES(4) = 125.360568
      CONSTS(25) = 108.56912706
      CONSTS(26) = 1e-4
      CONSTS(27) = 1e-4
      CONSTS(28) = 6.8284472205
      CONSTS(29) = 0.025
      CONSTS(30) = 0.2
      CONSTS(31) = 0.3
      CONSTS(32) = 0.025
      CONSTS(33) = 0.1545054945
      CONSTS(34) = 0.1381298227
      CONSTS(35) = 0.5508058134
      CONSTS(36) = 0.3445734208
      CONSTS(37) = 1.4544677036
      CONSTS(38) = 1.4843409851
      CONSTS(39) = 7.9822364317
      CONSTS(40) = 3.5e-4
      CONSTS(41) = 3.5e-4
      CONSTS(42) = 497.7852450367
      CONSTS(43) = 50
      CONSTS(44) = 485.7624931891
      CONSTS(45) = 577.7106000108
      CONSTS(46) = 0.03
      CONSTS(47) = 0.05
      CONSTS(48) = 0.01
      CONSTS(49) = 0.1
      CONSTS(50) = 21.83
      CONSTS(51) = 3379.545
      CONSTS(52) = -5
      CONSTS(53) = 0.0968305478
      CONSTS(54) = 0.4
      CONSTS(55) = 0.001
      CONSTS(56) = 0.025
      CONSTS(57) = 129.6486
      CONSTS(58) = 350.5314
      CONSTS(59) = 50.010747
      CONSTS(60) = 15
      CONSTS(61) = -5.6
      CONSTS(62) = 2
      CONSTS(63) = 8
      CONSTS(64) = 75
      CONSTS(65) = 1
      CONSTS(66) = 1
      CONSTS(67) = 0.5
      STATES(5) = 87.93968
      STATES(6) = 90.6179
      STATES(7) = 0.698577
      STATES(8) = 23.5957
      STATES(9) = 31.1705
      STATES(10) = 138.4476
      STATES(11) = 519.7915
      STATES(12) = 129.6439
      STATES(13) = 256.8555
      STATES(14) = 2961.6507
      STATES(15) = 232.46638962
      CONSTS(68) = 1e-4
      CONSTS(69) = 0.05
      CONSTS(70) = 0.05
      CONSTS(71) = 1e-4
      CONSTS(72) = 0.03
      CONSTS(73) = 4.2958026137
      CONSTS(74) = 0.0565149137
      CONSTS(75) = 0.0309026688
      CONSTS(76) = 1e-4
      CONSTS(77) = 1.5365929068
      CONSTS(78) = 2.6893667388
      CONSTS(79) = 3.1321449506
      CONSTS(80) = 7.7147
      CONSTS(81) = 27.87028922
      CONSTS(82) = 1.801907e-4
      CONSTS(83) = 1.932239e-4
      STATES(16) = 33.1398
      STATES(17) = 60.11203897
      STATES(18) = 58.926
      STATES(19) = 107.57022
      STATES(20) = 293.0398
      STATES(21) = 1.2282
      STATES(22) = 57.1876
      CONSTS(84) = 1
      CONSTS(85) = 40
      CONSTS(86) = 785
      CONSTS(87) = 15
      CONSTS(88) = 2.642367
      CONSTS(89) = 2.642367
      CONSTS(90) = 5.073345
      CONSTS(91) = 5.073345
      CONSTS(92) = 4.227788
      CONSTS(93) = 0.4932479
      CONSTS(94) = 0.4932479
      CONSTS(95) = 0.4932479
      CONSTS(96) = 0.13
      CONSTS(97) = 0.05507
      CONSTS(98) = 0.09129
      CONSTS(99) = 0.15602
      CONSTS(100) = 1.8
      CONSTS(101) = 0.58155
      CONSTS(102) = 0.68372
      CONSTS(103) = 0.832299
      STATES(23) = 2.76087
      STATES(24) = 4.41135
      STATES(25) = 4.992799
      STATES(26) = 4.22047
      STATES(27) = 8.55228
      STATES(28) = 7.8362
      STATES(29) = 8.213955
      STATES(30) = 8.76758
      CONSTS(104) =  (CONSTS(25)/CONSTS(24))*60.0000
      CONSTS(105) = CONSTS(61) - CONSTS(20)
      CONSTS(106) =  CONSTS(22)* (( 1.00000*1.00000)/CONSTS(24)) ** (1.0 / 2)*7.74597
      CONSTS(107) = 60.0000/CONSTS(24)
      CONSTS(108) =  CONSTS(23)* (1.00000/CONSTS(24)) ** (1.0 / 2)*7.74597
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(27) = STATES(18)/CONSTS(79) - ( CONSTS(62)*1.00000)/(EXP(STATES(18)/CONSTS(63)) - 1.00000)
      ALGBRC(28) = ALGBRC(27)+CONSTS(61)
      ALGBRC(30) = STATES(19)/CONSTS(80) - ( CONSTS(62)*1.00000)/(EXP(STATES(19)/CONSTS(63)) - 1.00000)
      ALGBRC(31) = ALGBRC(30)+CONSTS(61)
      ALGBRC(32) = (ALGBRC(28) - ALGBRC(31))/CONSTS(74)
      ALGBRC(33) = STATES(20)/CONSTS(81) - ( CONSTS(62)*1.00000)/(EXP(STATES(20)/CONSTS(63)) - 1.00000)
      ALGBRC(34) = ALGBRC(33)+CONSTS(61)
      ALGBRC(35) = (ALGBRC(31) - ALGBRC(34))/CONSTS(75)
      RATES(19) = ALGBRC(32) - ALGBRC(35)
      ALGBRC(36) = (ALGBRC(28) - ALGBRC(34))/CONSTS(73)
      RATES(18) = (STATES(22) - ALGBRC(36)) - ALGBRC(32)
      ALGBRC(37) =  CONSTS(42)*arbitrary_log((STATES(11) - CONSTS(44))/CONSTS(43)+1.00000, 10)
      ALGBRC(38) =  CONSTS(46)*EXP( CONSTS(49)*(STATES(11) - CONSTS(44)))+ CONSTS(47)*STATES(11) - CONSTS(44) ** 2.00000
      ALGBRC(39) =  CONSTS(67)*ALGBRC(37)+ (1.00000 - CONSTS(67))*ALGBRC(38)
      RATES(6) = (ALGBRC(39) - STATES(6))/CONSTS(60)
      ALGBRC(40) = STATES(12)/CONSTS(38) - ( CONSTS(62)*1.00000)/(EXP(STATES(12)/CONSTS(63)) - 1.00000)
      ALGBRC(41) = (ALGBRC(40) - ALGBRC(39))/CONSTS(32)
      RATES(12) = STATES(8) - ALGBRC(41)
      ALGBRC(45) = CONSTS(35)+ CONSTS(48)*EXP( 4.00000*CONSTS(67))+ CONSTS(48)*CONSTS(45)/STATES(11) ** 2.00000
      ALGBRC(42) = STATES(13)/CONSTS(39) - ( CONSTS(62)*1.00000)/(EXP(STATES(13)/CONSTS(63)) - 1.00000)
      ALGBRC(47) = (ALGBRC(39) - ALGBRC(42))/ALGBRC(45)
      RATES(11) = ALGBRC(41) - ALGBRC(47)
      ALGBRC(43) =  CONSTS(50)*arbitrary_log(CONSTS(51)/STATES(14) - 0.990000, 10)*- 1.00000
      ALGBRC(44) = (ALGBRC(42) - ALGBRC(43))/CONSTS(33)
      RATES(13) = ALGBRC(47) - ALGBRC(44)
      ALGBRC(46) = TERNRY(STATES(15).GT.CONSTS(57), (CONSTS(52)+ CONSTS(54)*EXP(CONSTS(57)/CONSTS(59))+( CONSTS(53)*(STATES(15) - CONSTS(57)))/1.00000) - CONSTS(62)/(EXP(STATES(15)/CONSTS(63)) - 1.00000), (CONSTS(52)+ CONSTS(54)*EXP(STATES(15)/CONSTS(59))) - CONSTS(62)/(EXP(STATES(15)/CONSTS(63)) - 1.00000))
      ALGBRC(48) = ALGBRC(46)+CONSTS(61)
      ALGBRC(49) = (ALGBRC(43) - ALGBRC(48))/CONSTS(34)
      RATES(14) = ALGBRC(44) - ALGBRC(49)
      ALGBRC(50) =  ((((( CONSTS(31)*CONSTS(28)*STATES(7) -  CONSTS(31)*CONSTS(28)*STATES(8)*0.0600000)+ (( STATES(10)*CONSTS(28))/CONSTS(37))*0.0600000) -  (( CONSTS(28)*CONSTS(62))/(EXP(STATES(10)/CONSTS(63)) - 1.00000))*0.0600000)+ ALGBRC(48)*CONSTS(31)*0.0600000)/(CONSTS(28)+CONSTS(31)))*16.6667
      RATES(8) = ((ALGBRC(50) -  STATES(8)*CONSTS(29)) - ALGBRC(40))/CONSTS(41)
      RATES(7) =  (((STATES(5) -  STATES(7)*CONSTS(27)*16.6667) - ALGBRC(50))/CONSTS(40))*0.0600000
      ALGBRC(52) = (ALGBRC(50) - ALGBRC(48))/CONSTS(28)
      RATES(10) =  ((STATES(7) -  STATES(8)*0.0600000) -  ALGBRC(52)*0.0600000)*16.6667
      RATES(9) = ((STATES(5) - STATES(9)/CONSTS(36))+( CONSTS(62)*1.00000)/(EXP(STATES(9)/CONSTS(63)) - 1.00000))/CONSTS(30)
      RATES(17) =  (STATES(21) -  STATES(22)*0.0600000)*16.6667
      ALGBRC(61) = ( RATES(17)*CONSTS(70)+CONSTS(61)+STATES(17)/CONSTS(78)) - ( CONSTS(62)*1.00000)/(EXP(STATES(17)/CONSTS(63)) - 1.00000)
      RATES(22) = ((ALGBRC(61) - ALGBRC(28)) -  STATES(22)*CONSTS(72))/CONSTS(83)
      ALGBRC(54) = STATES(23)+STATES(24)+STATES(25)+STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(30)
      ALGBRC(63) = STATES(2)+STATES(1)+STATES(4)+STATES(3)+CONSTS(87)+CONSTS(21)+ALGBRC(54)
      ALGBRC(64) =  CONSTS(84)*EXP((ALGBRC(63) - CONSTS(86))/CONSTS(85)) - ( CONSTS(62)*1.00000)/(EXP(ALGBRC(63)/CONSTS(64)) - 1.00000)
      ALGBRC(65) = ALGBRC(64)+CONSTS(61)
      ALGBRC(26) = STATES(28)/CONSTS(101) - ( CONSTS(66)*1.00000)/(EXP(STATES(28)/CONSTS(65)) - 1.00000)
      ALGBRC(72) = ALGBRC(26)+ALGBRC(65)
      ALGBRC(24) = STATES(29)/CONSTS(102) - ( CONSTS(66)*1.00000)/(EXP(STATES(29)/CONSTS(65)) - 1.00000)
      ALGBRC(69) = ALGBRC(24)+ALGBRC(65)
      ALGBRC(74) = (ALGBRC(72) - ALGBRC(69))/CONSTS(93)
      ALGBRC(22) = STATES(30)/CONSTS(103) - ( CONSTS(62)*1.00000)/(EXP(STATES(30)/CONSTS(63)) - 1.00000)
      ALGBRC(68) = ALGBRC(22)+ALGBRC(65)
      ALGBRC(70) = (ALGBRC(69) - ALGBRC(68))/CONSTS(94)
      RATES(29) = ALGBRC(74) - ALGBRC(70)
      ALGBRC(6) = VOI -  INT(VOI/CONSTS(107))*CONSTS(107)
      ALGBRC(10) = ALGBRC(6)
      ALGBRC(11) = TERNRY(ALGBRC(10).GE.0.00000.AND.ALGBRC(10).LE.CONSTS(108), (1.00000 - cos((  3.14159265358979*ALGBRC(10))/CONSTS(108)))/2.00000, TERNRY(ALGBRC(10).LT. 1.50000*CONSTS(108).AND.ALGBRC(10).GE.CONSTS(108), (1.00000+cos(( 2.00000* 3.14159265358979*(ALGBRC(10) - CONSTS(108)))/CONSTS(108)))/2.00000, 0.00000)
      ALGBRC(12) = ALGBRC(11)
      ALGBRC(15) =  (CONSTS(16) - CONSTS(17))*ALGBRC(12)+CONSTS(17)
      ALGBRC(71) =  (1.00000 - ALGBRC(12))*(CONSTS(7) - CONSTS(8))+CONSTS(8)
      ALGBRC(73) =  (STATES(1) - ALGBRC(71))*ALGBRC(15) - ( CONSTS(62)*1.00000)/(EXP(STATES(1)/CONSTS(63)) - 1.00000)
      ALGBRC(75) = ALGBRC(73)+ALGBRC(65)
      ALGBRC(77) = (ALGBRC(68) - ALGBRC(75))/CONSTS(95)
      RATES(30) = ALGBRC(70) - ALGBRC(77)
      ALGBRC(17) =  (CONSTS(18) - CONSTS(19))*ALGBRC(11)+CONSTS(19)
      ALGBRC(56) =  (1.00000 - ALGBRC(11))*(CONSTS(5) - CONSTS(6))+CONSTS(6)
      ALGBRC(58) =  (STATES(3) - ALGBRC(56))*ALGBRC(17) - ( CONSTS(62)*1.00000)/(EXP(STATES(3)/CONSTS(63)) - 1.00000)
      ALGBRC(76) = ALGBRC(58)+ALGBRC(65)
      ALGBRC(78) = (ALGBRC(34) - ALGBRC(76))/CONSTS(76)
      RATES(20) = (ALGBRC(35)+ALGBRC(36)) - ALGBRC(78)
      ALGBRC(53) =  CONSTS(55)*CONSTS(58)/STATES(15) ** 2.00000+CONSTS(56)
      ALGBRC(79) = (ALGBRC(48) - ALGBRC(75))/ALGBRC(53)
      RATES(15) = (ALGBRC(49)+ALGBRC(52)) - ALGBRC(79)
      ALGBRC(8) = ALGBRC(6) - CONSTS(11)
      ALGBRC(13) = TERNRY(ALGBRC(8).GE.0.00000.AND.ALGBRC(8).LE.CONSTS(106), (1.00000 - cos((  3.14159265358979*ALGBRC(8))/CONSTS(106)))/2.00000, TERNRY(ALGBRC(8).LT. 1.50000*CONSTS(106).AND.ALGBRC(8).GE.CONSTS(106), (1.00000+cos(( 2.00000* 3.14159265358979*(ALGBRC(8) - CONSTS(106)))/CONSTS(106)))/2.00000, 0.00000)
      ALGBRC(18) =  (CONSTS(12) - CONSTS(13))*ALGBRC(13)+CONSTS(13)
      ALGBRC(59) =  (1.00000 - ALGBRC(13))*(CONSTS(1) - CONSTS(2))+CONSTS(2)
      ALGBRC(60) =  (STATES(4) - ALGBRC(59))*ALGBRC(18) - ( CONSTS(62)*1.00000)/(EXP(STATES(4)/CONSTS(63)) - 1.00000)
      ALGBRC(80) = ALGBRC(60)+ALGBRC(65)
      ALGBRC(81) = TERNRY(ALGBRC(76).GT.ALGBRC(80), (ALGBRC(76) - ALGBRC(80))/CONSTS(10), 0.00000)
      RATES(3) = ALGBRC(78) - ALGBRC(81)
      ALGBRC(84) = TERNRY(ALGBRC(80).GT.STATES(5), (ALGBRC(80) - STATES(5))/CONSTS(26), 0.00000)
      ALGBRC(19) = STATES(5)
      ALGBRC(21) = STATES(24)/CONSTS(97) - ( CONSTS(66)*1.00000)/(EXP(STATES(24)/CONSTS(65)) - 1.00000)
      ALGBRC(66) = ALGBRC(21)+ALGBRC(65)
      ALGBRC(67) = (ALGBRC(19) - ALGBRC(66))/CONSTS(88)
      RATES(5) =  (((ALGBRC(84) - RATES(9)) -  STATES(7)*16.6667) - ALGBRC(67))*(1.00000/CONSTS(96)+( (CONSTS(62)/1.00000)*EXP(STATES(23)/CONSTS(65)))/EXP(STATES(23)/CONSTS(65)) - 1.00000 ** 2.00000)
      RATES(23) = ((ALGBRC(84) - RATES(9)) -  STATES(7)*16.6667) - ALGBRC(67)
      RATES(4) = ALGBRC(81) - ALGBRC(84)
      ALGBRC(83) = ABS((ALGBRC(80) - ALGBRC(65))/2.00000)
      ALGBRC(20) = STATES(27)/CONSTS(100) - ( CONSTS(66)*1.00000)/(EXP(STATES(27)/CONSTS(65)) - 1.00000)
      ALGBRC(86) = ALGBRC(20)+ALGBRC(83)
      ALGBRC(88) = (ALGBRC(86) - ALGBRC(72))/CONSTS(92)
      RATES(28) = ALGBRC(88) - ALGBRC(74)
      ALGBRC(14) = ALGBRC(13)
      ALGBRC(16) =  (CONSTS(14) - CONSTS(15))*ALGBRC(14)+CONSTS(15)
      ALGBRC(82) =  (1.00000 - ALGBRC(14))*(CONSTS(3) - CONSTS(4))+CONSTS(4)
      ALGBRC(85) =  (STATES(2) - ALGBRC(82))*ALGBRC(16) - ( CONSTS(62)*1.00000)/(EXP(STATES(2)/CONSTS(63)) - 1.00000)
      ALGBRC(87) = ALGBRC(85)+ALGBRC(65)
      ALGBRC(89) = TERNRY(ALGBRC(75).GT.ALGBRC(87), (ALGBRC(75) - ALGBRC(87))/CONSTS(9), 0.00000)
      RATES(1) = (ALGBRC(79) - ALGBRC(89))+ALGBRC(77)
      ALGBRC(23) = STATES(25)/CONSTS(98) - ( CONSTS(66)*1.00000)/(EXP(STATES(25)/CONSTS(65)) - 1.00000)
      ALGBRC(90) = ALGBRC(23)+ALGBRC(83)
      ALGBRC(92) = (ALGBRC(66) - ALGBRC(90))/CONSTS(89)
      RATES(24) = ALGBRC(67) - ALGBRC(92)
      ALGBRC(91) = ((((( CONSTS(68)*STATES(16))/CONSTS(77) - ( CONSTS(68)*CONSTS(62)*1.00000)/(EXP(STATES(16)/CONSTS(63)) - 1.00000))+ ALGBRC(87)*CONSTS(69)) -  CONSTS(68)*CONSTS(69)*STATES(21)*16.6667)+ CONSTS(61)*CONSTS(68))/(CONSTS(69)+CONSTS(68))
      ALGBRC(29) = ((STATES(16)/CONSTS(77)+CONSTS(61)) -  CONSTS(69)*STATES(21)*16.6667) - ( CONSTS(62)*1.00000)/(EXP(STATES(16)/CONSTS(63)) - 1.00000)
      ALGBRC(93) = TERNRY(ALGBRC(87).GT.ALGBRC(91), ALGBRC(91), ALGBRC(29))
      RATES(21) =  (((ALGBRC(93) - ALGBRC(61)) -  STATES(21)*CONSTS(71)*16.6667)/CONSTS(82))*0.0600000
      ALGBRC(25) = STATES(26)/CONSTS(99) - ( CONSTS(66)*1.00000)/(EXP(STATES(26)/CONSTS(65)) - 1.00000)
      ALGBRC(94) = ALGBRC(25)+ALGBRC(83)
      ALGBRC(97) = (ALGBRC(90) - ALGBRC(94))/CONSTS(90)
      RATES(25) = ALGBRC(92) - ALGBRC(97)
      ALGBRC(95) = TERNRY(ALGBRC(87).GT.ALGBRC(93), (ALGBRC(87) - ALGBRC(93))/CONSTS(68), 0.00000)
      RATES(16) = ALGBRC(95) -  STATES(21)*16.6667
      RATES(2) = ALGBRC(89) - ALGBRC(95)
      ALGBRC(98) = (ALGBRC(94) - ALGBRC(86))/CONSTS(91)
      RATES(27) = ALGBRC(98) - ALGBRC(88)
      RATES(26) = ALGBRC(97) - ALGBRC(98)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(27) = STATES(18)/CONSTS(79) - ( CONSTS(62)*1.00000)/(EXP(STATES(18)/CONSTS(63)) - 1.00000)
      ALGBRC(28) = ALGBRC(27)+CONSTS(61)
      ALGBRC(30) = STATES(19)/CONSTS(80) - ( CONSTS(62)*1.00000)/(EXP(STATES(19)/CONSTS(63)) - 1.00000)
      ALGBRC(31) = ALGBRC(30)+CONSTS(61)
      ALGBRC(32) = (ALGBRC(28) - ALGBRC(31))/CONSTS(74)
      ALGBRC(33) = STATES(20)/CONSTS(81) - ( CONSTS(62)*1.00000)/(EXP(STATES(20)/CONSTS(63)) - 1.00000)
      ALGBRC(34) = ALGBRC(33)+CONSTS(61)
      ALGBRC(35) = (ALGBRC(31) - ALGBRC(34))/CONSTS(75)
      ALGBRC(36) = (ALGBRC(28) - ALGBRC(34))/CONSTS(73)
      ALGBRC(37) =  CONSTS(42)*arbitrary_log((STATES(11) - CONSTS(44))/CONSTS(43)+1.00000, 10)
      ALGBRC(38) =  CONSTS(46)*EXP( CONSTS(49)*(STATES(11) - CONSTS(44)))+ CONSTS(47)*STATES(11) - CONSTS(44) ** 2.00000
      ALGBRC(39) =  CONSTS(67)*ALGBRC(37)+ (1.00000 - CONSTS(67))*ALGBRC(38)
      ALGBRC(40) = STATES(12)/CONSTS(38) - ( CONSTS(62)*1.00000)/(EXP(STATES(12)/CONSTS(63)) - 1.00000)
      ALGBRC(41) = (ALGBRC(40) - ALGBRC(39))/CONSTS(32)
      ALGBRC(45) = CONSTS(35)+ CONSTS(48)*EXP( 4.00000*CONSTS(67))+ CONSTS(48)*CONSTS(45)/STATES(11) ** 2.00000
      ALGBRC(42) = STATES(13)/CONSTS(39) - ( CONSTS(62)*1.00000)/(EXP(STATES(13)/CONSTS(63)) - 1.00000)
      ALGBRC(47) = (ALGBRC(39) - ALGBRC(42))/ALGBRC(45)
      ALGBRC(43) =  CONSTS(50)*arbitrary_log(CONSTS(51)/STATES(14) - 0.990000, 10)*- 1.00000
      ALGBRC(44) = (ALGBRC(42) - ALGBRC(43))/CONSTS(33)
      ALGBRC(46) = TERNRY(STATES(15).GT.CONSTS(57), (CONSTS(52)+ CONSTS(54)*EXP(CONSTS(57)/CONSTS(59))+( CONSTS(53)*(STATES(15) - CONSTS(57)))/1.00000) - CONSTS(62)/(EXP(STATES(15)/CONSTS(63)) - 1.00000), (CONSTS(52)+ CONSTS(54)*EXP(STATES(15)/CONSTS(59))) - CONSTS(62)/(EXP(STATES(15)/CONSTS(63)) - 1.00000))
      ALGBRC(48) = ALGBRC(46)+CONSTS(61)
      ALGBRC(49) = (ALGBRC(43) - ALGBRC(48))/CONSTS(34)
      ALGBRC(50) =  ((((( CONSTS(31)*CONSTS(28)*STATES(7) -  CONSTS(31)*CONSTS(28)*STATES(8)*0.0600000)+ (( STATES(10)*CONSTS(28))/CONSTS(37))*0.0600000) -  (( CONSTS(28)*CONSTS(62))/(EXP(STATES(10)/CONSTS(63)) - 1.00000))*0.0600000)+ ALGBRC(48)*CONSTS(31)*0.0600000)/(CONSTS(28)+CONSTS(31)))*16.6667
      ALGBRC(52) = (ALGBRC(50) - ALGBRC(48))/CONSTS(28)
      ALGBRC(61) = ( RATES(17)*CONSTS(70)+CONSTS(61)+STATES(17)/CONSTS(78)) - ( CONSTS(62)*1.00000)/(EXP(STATES(17)/CONSTS(63)) - 1.00000)
      ALGBRC(54) = STATES(23)+STATES(24)+STATES(25)+STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(30)
      ALGBRC(63) = STATES(2)+STATES(1)+STATES(4)+STATES(3)+CONSTS(87)+CONSTS(21)+ALGBRC(54)
      ALGBRC(64) =  CONSTS(84)*EXP((ALGBRC(63) - CONSTS(86))/CONSTS(85)) - ( CONSTS(62)*1.00000)/(EXP(ALGBRC(63)/CONSTS(64)) - 1.00000)
      ALGBRC(65) = ALGBRC(64)+CONSTS(61)
      ALGBRC(26) = STATES(28)/CONSTS(101) - ( CONSTS(66)*1.00000)/(EXP(STATES(28)/CONSTS(65)) - 1.00000)
      ALGBRC(72) = ALGBRC(26)+ALGBRC(65)
      ALGBRC(24) = STATES(29)/CONSTS(102) - ( CONSTS(66)*1.00000)/(EXP(STATES(29)/CONSTS(65)) - 1.00000)
      ALGBRC(69) = ALGBRC(24)+ALGBRC(65)
      ALGBRC(74) = (ALGBRC(72) - ALGBRC(69))/CONSTS(93)
      ALGBRC(22) = STATES(30)/CONSTS(103) - ( CONSTS(62)*1.00000)/(EXP(STATES(30)/CONSTS(63)) - 1.00000)
      ALGBRC(68) = ALGBRC(22)+ALGBRC(65)
      ALGBRC(70) = (ALGBRC(69) - ALGBRC(68))/CONSTS(94)
      ALGBRC(6) = VOI -  INT(VOI/CONSTS(107))*CONSTS(107)
      ALGBRC(10) = ALGBRC(6)
      ALGBRC(11) = TERNRY(ALGBRC(10).GE.0.00000.AND.ALGBRC(10).LE.CONSTS(108), (1.00000 - cos((  3.14159265358979*ALGBRC(10))/CONSTS(108)))/2.00000, TERNRY(ALGBRC(10).LT. 1.50000*CONSTS(108).AND.ALGBRC(10).GE.CONSTS(108), (1.00000+cos(( 2.00000* 3.14159265358979*(ALGBRC(10) - CONSTS(108)))/CONSTS(108)))/2.00000, 0.00000)
      ALGBRC(12) = ALGBRC(11)
      ALGBRC(15) =  (CONSTS(16) - CONSTS(17))*ALGBRC(12)+CONSTS(17)
      ALGBRC(71) =  (1.00000 - ALGBRC(12))*(CONSTS(7) - CONSTS(8))+CONSTS(8)
      ALGBRC(73) =  (STATES(1) - ALGBRC(71))*ALGBRC(15) - ( CONSTS(62)*1.00000)/(EXP(STATES(1)/CONSTS(63)) - 1.00000)
      ALGBRC(75) = ALGBRC(73)+ALGBRC(65)
      ALGBRC(77) = (ALGBRC(68) - ALGBRC(75))/CONSTS(95)
      ALGBRC(17) =  (CONSTS(18) - CONSTS(19))*ALGBRC(11)+CONSTS(19)
      ALGBRC(56) =  (1.00000 - ALGBRC(11))*(CONSTS(5) - CONSTS(6))+CONSTS(6)
      ALGBRC(58) =  (STATES(3) - ALGBRC(56))*ALGBRC(17) - ( CONSTS(62)*1.00000)/(EXP(STATES(3)/CONSTS(63)) - 1.00000)
      ALGBRC(76) = ALGBRC(58)+ALGBRC(65)
      ALGBRC(78) = (ALGBRC(34) - ALGBRC(76))/CONSTS(76)
      ALGBRC(53) =  CONSTS(55)*CONSTS(58)/STATES(15) ** 2.00000+CONSTS(56)
      ALGBRC(79) = (ALGBRC(48) - ALGBRC(75))/ALGBRC(53)
      ALGBRC(8) = ALGBRC(6) - CONSTS(11)
      ALGBRC(13) = TERNRY(ALGBRC(8).GE.0.00000.AND.ALGBRC(8).LE.CONSTS(106), (1.00000 - cos((  3.14159265358979*ALGBRC(8))/CONSTS(106)))/2.00000, TERNRY(ALGBRC(8).LT. 1.50000*CONSTS(106).AND.ALGBRC(8).GE.CONSTS(106), (1.00000+cos(( 2.00000* 3.14159265358979*(ALGBRC(8) - CONSTS(106)))/CONSTS(106)))/2.00000, 0.00000)
      ALGBRC(18) =  (CONSTS(12) - CONSTS(13))*ALGBRC(13)+CONSTS(13)
      ALGBRC(59) =  (1.00000 - ALGBRC(13))*(CONSTS(1) - CONSTS(2))+CONSTS(2)
      ALGBRC(60) =  (STATES(4) - ALGBRC(59))*ALGBRC(18) - ( CONSTS(62)*1.00000)/(EXP(STATES(4)/CONSTS(63)) - 1.00000)
      ALGBRC(80) = ALGBRC(60)+ALGBRC(65)
      ALGBRC(81) = TERNRY(ALGBRC(76).GT.ALGBRC(80), (ALGBRC(76) - ALGBRC(80))/CONSTS(10), 0.00000)
      ALGBRC(84) = TERNRY(ALGBRC(80).GT.STATES(5), (ALGBRC(80) - STATES(5))/CONSTS(26), 0.00000)
      ALGBRC(19) = STATES(5)
      ALGBRC(21) = STATES(24)/CONSTS(97) - ( CONSTS(66)*1.00000)/(EXP(STATES(24)/CONSTS(65)) - 1.00000)
      ALGBRC(66) = ALGBRC(21)+ALGBRC(65)
      ALGBRC(67) = (ALGBRC(19) - ALGBRC(66))/CONSTS(88)
      ALGBRC(83) = ABS((ALGBRC(80) - ALGBRC(65))/2.00000)
      ALGBRC(20) = STATES(27)/CONSTS(100) - ( CONSTS(66)*1.00000)/(EXP(STATES(27)/CONSTS(65)) - 1.00000)
      ALGBRC(86) = ALGBRC(20)+ALGBRC(83)
      ALGBRC(88) = (ALGBRC(86) - ALGBRC(72))/CONSTS(92)
      ALGBRC(14) = ALGBRC(13)
      ALGBRC(16) =  (CONSTS(14) - CONSTS(15))*ALGBRC(14)+CONSTS(15)
      ALGBRC(82) =  (1.00000 - ALGBRC(14))*(CONSTS(3) - CONSTS(4))+CONSTS(4)
      ALGBRC(85) =  (STATES(2) - ALGBRC(82))*ALGBRC(16) - ( CONSTS(62)*1.00000)/(EXP(STATES(2)/CONSTS(63)) - 1.00000)
      ALGBRC(87) = ALGBRC(85)+ALGBRC(65)
      ALGBRC(89) = TERNRY(ALGBRC(75).GT.ALGBRC(87), (ALGBRC(75) - ALGBRC(87))/CONSTS(9), 0.00000)
      ALGBRC(23) = STATES(25)/CONSTS(98) - ( CONSTS(66)*1.00000)/(EXP(STATES(25)/CONSTS(65)) - 1.00000)
      ALGBRC(90) = ALGBRC(23)+ALGBRC(83)
      ALGBRC(92) = (ALGBRC(66) - ALGBRC(90))/CONSTS(89)
      ALGBRC(91) = ((((( CONSTS(68)*STATES(16))/CONSTS(77) - ( CONSTS(68)*CONSTS(62)*1.00000)/(EXP(STATES(16)/CONSTS(63)) - 1.00000))+ ALGBRC(87)*CONSTS(69)) -  CONSTS(68)*CONSTS(69)*STATES(21)*16.6667)+ CONSTS(61)*CONSTS(68))/(CONSTS(69)+CONSTS(68))
      ALGBRC(29) = ((STATES(16)/CONSTS(77)+CONSTS(61)) -  CONSTS(69)*STATES(21)*16.6667) - ( CONSTS(62)*1.00000)/(EXP(STATES(16)/CONSTS(63)) - 1.00000)
      ALGBRC(93) = TERNRY(ALGBRC(87).GT.ALGBRC(91), ALGBRC(91), ALGBRC(29))
      ALGBRC(25) = STATES(26)/CONSTS(99) - ( CONSTS(66)*1.00000)/(EXP(STATES(26)/CONSTS(65)) - 1.00000)
      ALGBRC(94) = ALGBRC(25)+ALGBRC(83)
      ALGBRC(97) = (ALGBRC(90) - ALGBRC(94))/CONSTS(90)
      ALGBRC(95) = TERNRY(ALGBRC(87).GT.ALGBRC(93), (ALGBRC(87) - ALGBRC(93))/CONSTS(68), 0.00000)
      ALGBRC(98) = (ALGBRC(94) - ALGBRC(86))/CONSTS(91)
      ALGBRC(1) = STATES(5) - CONSTS(61)
      ALGBRC(2) = STATES(9)+STATES(10)+STATES(11)+STATES(12)
      ALGBRC(3) = STATES(14)+STATES(15)
      ALGBRC(4) = STATES(16)+STATES(17)+STATES(18)
      ALGBRC(5) = STATES(20)
      ALGBRC(7) = TERNRY(VOI.GT.ALGBRC(6).AND.VOI.LT.ALGBRC(6)+0.0200000, 5.00000, TERNRY(VOI.GT.ALGBRC(6)+CONSTS(11).AND.VOI.LT.ALGBRC(6)+CONSTS(11)+0.0200000, 10.0000, 0.00000)
      ALGBRC(9) = ALGBRC(1)
      ALGBRC(51) = ALGBRC(50) - CONSTS(20)
      ALGBRC(55) = STATES(1)+STATES(2)+STATES(16)+STATES(17)+STATES(18)+STATES(19)+STATES(20)+STATES(3)+STATES(4)+STATES(9)+STATES(10)+STATES(11)+STATES(12)+STATES(13)+STATES(14)+STATES(15)+CONSTS(87)+ALGBRC(54)
      ALGBRC(57) = ALGBRC(55) - CONSTS(87)
      ALGBRC(62) = ALGBRC(61) - CONSTS(61)
      ALGBRC(96) = ALGBRC(93) - CONSTS(61)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END