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