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 122 entries in the algebraic variable array.
C There are a total of 2 entries in each of the rate and state variable arrays.
C There are a total of 34 entries in the constant variable array.
C
C
C CONSTS(1) is k1_m in component gas_exchange (dimensionless).
C CONSTS(2) is k1_81m in component gas_exchange (dimensionless).
C CONSTS(3) is k2_m in component gas_exchange (dimensionless).
C CONSTS(4) is k3_m in component gas_exchange (dimensionless).
C CONSTS(5) is K3_81m in component gas_exchange (dimensionless).
C CONSTS(6) is k4_m in component gas_exchange (dimensionless).
C CONSTS(7) is pH_m in component gas_exchange (dimensionless).
C CONSTS(8) is k1_f in component gas_exchange (dimensionless).
C CONSTS(9) is k2_f in component gas_exchange (dimensionless).
C CONSTS(10) is k3_f in component gas_exchange (dimensionless).
C CONSTS(11) is k4_f in component gas_exchange (dimensionless).
C CONSTS(12) is pH_f in component gas_exchange (dimensionless).
C STATES(1) is Pin_m in component gas_exchange (mmHg).
C STATES(2) is Pin_f in component gas_exchange (mmHg).
C ALGBRC(1) is x_m in component gas_exchange (dimensionless).
C ALGBRC(2) is x_m1 in component gas_exchange (dimensionless).
C ALGBRC(4) is x_f in component gas_exchange (dimensionless).
C ALGBRC(5) is x_f1 in component gas_exchange (dimensionless).
C ALGBRC(3) is S_m in component gas_exchange (dimensionless).
C ALGBRC(6) is S_f in component gas_exchange (dimensionless).
C VOI is t in component gas_exchange (second).
C ALGBRC(7) is O2_m in component gas_exchange (ml_per_ml).
C ALGBRC(8) is O2_f in component gas_exchange (ml_per_ml).
C CONSTS(13) is Hb_m in component gas_exchange (dimensionless).
C CONSTS(14) is Hb_f in component gas_exchange (dimensionless).
C CONSTS(15) is a in component gas_exchange (dimensionless).
C CONSTS(16) is ks_m in component gas_exchange (dimensionless).
C CONSTS(17) is ks_f in component gas_exchange (dimensionless).
C CONSTS(18) is kh in component gas_exchange (dimensionless).
C CONSTS(19) is B in component gas_exchange (dimensionless).
C CONSTS(20) is V_m in component gas_exchange (ml).
C CONSTS(21) is V_f in component gas_exchange (ml).
C CONSTS(22) is Q_m in component gas_exchange (ml_per_min).
C CONSTS(23) is Q_f in component gas_exchange (ml_per_min).
C CONSTS(24) is Dm in component gas_exchange (ml_per_ml_min_mmHg).
C CONSTS(25) is c1 in component gas_exchange (dimensionless).
C CONSTS(26) is c2 in component gas_exchange (dimensionless).
C CONSTS(27) is c3 in component gas_exchange (dimensionless).
C CONSTS(28) is c4 in component gas_exchange (dimensionless).
C ALGBRC(11) is theta_m in component gas_exchange (dimensionless).
C ALGBRC(14) is theta_f in component gas_exchange (dimensionless).
C ALGBRC(9) is x1_m in component gas_exchange (dimensionless).
C ALGBRC(12) is x1_f in component gas_exchange (dimensionless).
C ALGBRC(15) is R_m in component gas_exchange (dimensionless).
C ALGBRC(16) is R_f in component gas_exchange (dimensionless).
C CONSTS(29) is R_p in component gas_exchange (dimensionless).
C ALGBRC(17) is R_t in component gas_exchange (dimensionless).
C ALGBRC(18) is Dp in component gas_exchange (dimensionless).
C ALGBRC(19) is dodt_m in component gas_exchange (dimensionless).
C ALGBRC(20) is dodt_f in component gas_exchange (dimensionless).
C ALGBRC(21) is O2_mnew in component gas_exchange (ml_per_ml).
C ALGBRC(22) is O2_fnew in component gas_exchange (ml_per_ml).
C ALGBRC(23) is f_m in component gas_exchange (dimensionless).
C ALGBRC(24) is f_f in component gas_exchange (dimensionless).
C ALGBRC(25) is upper_m in component gas_exchange (dimensionless).
C ALGBRC(26) is lower_m in component gas_exchange (dimensionless).
C ALGBRC(28) is upper_f in component gas_exchange (dimensionless).
C ALGBRC(29) is lower_f in component gas_exchange (dimensionless).
C ALGBRC(30) is dfdP_f in component gas_exchange (dimensionless).
C ALGBRC(27) is dfdP_m in component gas_exchange (dimensionless).
C CONSTS(30) is tt in component gas_exchange (second).
C CONSTS(31) is dt in component gas_exchange (second).
C ALGBRC(31) is Pend1_m in component gas_exchange (mmHg).
C ALGBRC(32) is Pend1_f in component gas_exchange (mmHg).
C ALGBRC(10) is theta_96_m in component gas_exchange (dimensionless).
C ALGBRC(13) is theta_96_f in component gas_exchange (dimensionless).
C ALGBRC(33) is x_m_sec in component gas_exchange (dimensionless).
C ALGBRC(34) is x_m1_sec in component gas_exchange (dimensionless).
C ALGBRC(35) is S_m_sec in component gas_exchange (dimensionless).
C ALGBRC(36) is x_f_sec in component gas_exchange (dimensionless).
C ALGBRC(37) is x_f1_sec in component gas_exchange (dimensionless).
C ALGBRC(38) is S_f_sec in component gas_exchange (dimensionless).
C ALGBRC(40) is theta_96_m_sec in component gas_exchange (dimensionless).
C ALGBRC(43) is theta_96_f_sec in component gas_exchange (dimensionless).
C ALGBRC(39) is x1_m_sec in component gas_exchange (dimensionless).
C ALGBRC(41) is theta_m_sec in component gas_exchange (dimensionless).
C ALGBRC(42) is x1_f_sec in component gas_exchange (dimensionless).
C ALGBRC(44) is theta_f_sec in component gas_exchange (dimensionless).
C ALGBRC(45) is R_m_sec in component gas_exchange (dimensionless).
C ALGBRC(46) is R_f_sec in component gas_exchange (dimensionless).
C CONSTS(32) is R_p_sec in component gas_exchange (dimensionless).
C ALGBRC(47) is R_t_sec in component gas_exchange (dimensionless).
C ALGBRC(48) is Dp_sec in component gas_exchange (dimensionless).
C ALGBRC(49) is dodt_m_sec in component gas_exchange (dimensionless).
C ALGBRC(50) is dodt_f_sec in component gas_exchange (dimensionless).
C ALGBRC(51) is O2_mnew_sec in component gas_exchange (ml_per_ml).
C ALGBRC(52) is O2_fnew_sec in component gas_exchange (ml_per_ml).
C ALGBRC(53) is f_m_sec in component gas_exchange (dimensionless).
C ALGBRC(54) is f_f_sec in component gas_exchange (dimensionless).
C ALGBRC(55) is upper_m_sec in component gas_exchange (dimensionless).
C ALGBRC(56) is lower_m_sec in component gas_exchange (dimensionless).
C ALGBRC(57) is dfdP_m_sec in component gas_exchange (dimensionless).
C ALGBRC(58) is upper_f_sec in component gas_exchange (dimensionless).
C ALGBRC(59) is lower_f_sec in component gas_exchange (dimensionless).
C ALGBRC(60) is dfdP_f_sec in component gas_exchange (dimensionless).
C ALGBRC(61) is Pend2_m in component gas_exchange (mmHg).
C ALGBRC(62) is Pend2_f in component gas_exchange (mmHg).
C ALGBRC(63) is x_m_trd in component gas_exchange (dimensionless).
C ALGBRC(64) is x_m1_trd in component gas_exchange (dimensionless).
C ALGBRC(65) is S_m_trd in component gas_exchange (dimensionless).
C ALGBRC(66) is x_f_trd in component gas_exchange (dimensionless).
C ALGBRC(67) is x_f1_trd in component gas_exchange (dimensionless).
C ALGBRC(68) is S_f_trd in component gas_exchange (dimensionless).
C ALGBRC(70) is theta_96_m_trd in component gas_exchange (dimensionless).
C ALGBRC(73) is theta_96_f_trd in component gas_exchange (dimensionless).
C ALGBRC(69) is x1_m_trd in component gas_exchange (dimensionless).
C ALGBRC(71) is theta_m_trd in component gas_exchange (dimensionless).
C ALGBRC(72) is x1_f_trd in component gas_exchange (dimensionless).
C ALGBRC(74) is theta_f_trd in component gas_exchange (dimensionless).
C ALGBRC(75) is R_m_trd in component gas_exchange (dimensionless).
C ALGBRC(76) is R_f_trd in component gas_exchange (dimensionless).
C CONSTS(33) is R_p_trd in component gas_exchange (dimensionless).
C ALGBRC(77) is R_t_trd in component gas_exchange (dimensionless).
C ALGBRC(78) is Dp_trd in component gas_exchange (dimensionless).
C ALGBRC(79) is dodt_m_trd in component gas_exchange (dimensionless).
C ALGBRC(80) is dodt_f_trd in component gas_exchange (dimensionless).
C ALGBRC(81) is O2_mnew_trd in component gas_exchange (ml_per_ml).
C ALGBRC(82) is O2_fnew_trd in component gas_exchange (ml_per_ml).
C ALGBRC(83) is f_m_trd in component gas_exchange (dimensionless).
C ALGBRC(84) is f_f_trd in component gas_exchange (dimensionless).
C ALGBRC(85) is upper_m_trd in component gas_exchange (dimensionless).
C ALGBRC(86) is lower_m_trd in component gas_exchange (dimensionless).
C ALGBRC(87) is dfdP_m_trd in component gas_exchange (dimensionless).
C ALGBRC(88) is upper_f_trd in component gas_exchange (dimensionless).
C ALGBRC(89) is lower_f_trd in component gas_exchange (dimensionless).
C ALGBRC(90) is dfdP_f_trd in component gas_exchange (dimensionless).
C ALGBRC(91) is Pend3_m in component gas_exchange (mmHg).
C ALGBRC(92) is Pend3_f in component gas_exchange (mmHg).
C ALGBRC(93) is x_m_for in component gas_exchange (dimensionless).
C ALGBRC(94) is x_m1_for in component gas_exchange (dimensionless).
C ALGBRC(95) is S_m_for in component gas_exchange (dimensionless).
C ALGBRC(96) is x_f_for in component gas_exchange (dimensionless).
C ALGBRC(97) is x_f1_for in component gas_exchange (dimensionless).
C ALGBRC(98) is S_f_for in component gas_exchange (dimensionless).
C ALGBRC(100) is theta_96_m_for in component gas_exchange (dimensionless).
C ALGBRC(103) is theta_96_f_for in component gas_exchange (dimensionless).
C ALGBRC(99) is x1_m_for in component gas_exchange (dimensionless).
C ALGBRC(101) is theta_m_for in component gas_exchange (dimensionless).
C ALGBRC(102) is x1_f_for in component gas_exchange (dimensionless).
C ALGBRC(104) is theta_f_for in component gas_exchange (dimensionless).
C ALGBRC(105) is R_m_for in component gas_exchange (dimensionless).
C ALGBRC(106) is R_f_for in component gas_exchange (dimensionless).
C CONSTS(34) is R_p_for in component gas_exchange (dimensionless).
C ALGBRC(107) is R_t_for in component gas_exchange (dimensionless).
C ALGBRC(108) is Dp_for in component gas_exchange (dimensionless).
C ALGBRC(109) is dodt_m_for in component gas_exchange (dimensionless).
C ALGBRC(111) is dodt_f_for in component gas_exchange (dimensionless).
C ALGBRC(110) is O2_mnew_final in component gas_exchange (ml_per_ml).
C ALGBRC(113) is O2_fnew_final in component gas_exchange (ml_per_ml).
C ALGBRC(112) is f_m_final in component gas_exchange (dimensionless).
C ALGBRC(115) is f_f_final in component gas_exchange (dimensionless).
C ALGBRC(114) is upper_m_final in component gas_exchange (dimensionless).
C ALGBRC(116) is lower_m_final in component gas_exchange (dimensionless).
C ALGBRC(118) is dfdP_m_final in component gas_exchange (dimensionless).
C ALGBRC(117) is upper_f_final in component gas_exchange (dimensionless).
C ALGBRC(119) is lower_f_final in component gas_exchange (dimensionless).
C ALGBRC(120) is dfdP_f_final in component gas_exchange (dimensionless).
C ALGBRC(121) is Pfinal_m in component gas_exchange (mmHg).
C ALGBRC(122) is Pfinal_f in component gas_exchange (mmHg).
C RATES(1) is d/dt Pin_m in component gas_exchange (mmHg).
C RATES(2) is d/dt Pin_f in component gas_exchange (mmHg).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 1.4452
CONSTS(2) = 1.4452
CONSTS(3) = 0.456
CONSTS(4) = 0.3711
CONSTS(5) = 0.3711
CONSTS(6) = 7.8322e+03
CONSTS(7) = 7.4
CONSTS(8) = 1.302
CONSTS(9) = 0.464
CONSTS(10) = 0.395
CONSTS(11) = 2.2643e+03
CONSTS(12) = 7.35
STATES(1) = 33
STATES(2) = 15
CONSTS(13) = 12.5
CONSTS(14) = 15.5
CONSTS(15) = 3.0263e-05
CONSTS(16) = 164
CONSTS(17) = 164
CONSTS(18) = 164
CONSTS(19) = 1
CONSTS(20) = 4.89e-7
CONSTS(21) = 4.89e-7
CONSTS(22) = 1.2014e-04
CONSTS(23) = 1.2014e-04
CONSTS(24) = 6.21
CONSTS(25) = 3.287
CONSTS(26) = 0.1117
CONSTS(27) = 7.05e-3
CONSTS(28) = 0.8142
CONSTS(29) = 1.00000/CONSTS(24)
CONSTS(30) = CONSTS(21)/(CONSTS(23)/60.0000)
CONSTS(31) = CONSTS(30)/100.000
CONSTS(32) = 1.00000/CONSTS(24)
CONSTS(33) = 1.00000/CONSTS(24)
CONSTS(34) = 1.00000/CONSTS(24)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(7) = (( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*STATES(1)
ALGBRC(1) = ((arbitrary_log(STATES(1), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(2) = 10.0000 ** ALGBRC(1)
ALGBRC(3) = ( 100.000*ALGBRC(2))/(1.00000+ALGBRC(2))
ALGBRC(9) = TERNRY(ALGBRC(3).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(3))))
ALGBRC(10) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(9))+ CONSTS(27)*ALGBRC(3)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(11) = TERNRY(ALGBRC(3).GT.96.0000, ALGBRC(10) - ( ALGBRC(10)*(96.0000 - ALGBRC(3)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(9))+ CONSTS(27)*ALGBRC(3)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(15) = 1.00000/( ALGBRC(11)*CONSTS(20))
ALGBRC(4) = ((arbitrary_log(STATES(2), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(5) = 10.0000 ** ALGBRC(4)
ALGBRC(6) = ( 100.000*ALGBRC(5))/(1.00000+ALGBRC(5))
ALGBRC(12) = TERNRY(ALGBRC(6).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(6))))
ALGBRC(13) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(12))+ CONSTS(27)*ALGBRC(6)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(14) = TERNRY(ALGBRC(6).GT.96.0000, ALGBRC(13) - ( ALGBRC(13)*(96.0000 - ALGBRC(6)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(12))+ CONSTS(27)*ALGBRC(6)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(16) = 1.00000/( ALGBRC(14)*CONSTS(21))
ALGBRC(17) = ALGBRC(15)+ALGBRC(16)+CONSTS(29)
ALGBRC(18) = (1.00000/ALGBRC(17))/60.0000
ALGBRC(19) = ( CONSTS(19)*- 1.00000*ALGBRC(18)*(STATES(1) - STATES(2)))/CONSTS(20)
ALGBRC(21) = ALGBRC(7)+( ALGBRC(19)*CONSTS(31))/2.00000
ALGBRC(23) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*STATES(1)) - ALGBRC(21)
ALGBRC(25) = 1.34000*CONSTS(13)*CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(26) = 100.000*CONSTS(4)* CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(27) = ALGBRC(25)/ALGBRC(26)+CONSTS(15)
ALGBRC(31) = STATES(1) - ALGBRC(23)/ALGBRC(27)
ALGBRC(8) = (( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*STATES(2)
ALGBRC(20) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(19)
ALGBRC(22) = ALGBRC(8)+( ALGBRC(20)*CONSTS(31))/2.00000
ALGBRC(24) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*STATES(2)) - ALGBRC(22)
ALGBRC(28) = 1.34000*CONSTS(14)*CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(29) = 100.000*CONSTS(10)* CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(30) = ALGBRC(28)/ALGBRC(29)+CONSTS(15)
ALGBRC(32) = STATES(2) - ALGBRC(24)/ALGBRC(30)
ALGBRC(33) = ((arbitrary_log(ALGBRC(31), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(34) = 10.0000 ** ALGBRC(33)
ALGBRC(35) = ( 100.000*ALGBRC(34))/(1.00000+ALGBRC(34))
ALGBRC(39) = TERNRY(ALGBRC(35).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(35))))
ALGBRC(40) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(39))+ CONSTS(27)*ALGBRC(35)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(41) = TERNRY(ALGBRC(35).GT.96.0000, ALGBRC(40) - ( ALGBRC(40)*(96.0000 - ALGBRC(35)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(39))+ CONSTS(27)*ALGBRC(35)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(45) = 1.00000/( ALGBRC(41)*CONSTS(20))
ALGBRC(36) = ((arbitrary_log(ALGBRC(32), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(37) = 10.0000 ** ALGBRC(36)
ALGBRC(38) = ( 100.000*ALGBRC(37))/(1.00000+ALGBRC(37))
ALGBRC(42) = TERNRY(ALGBRC(38).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(38))))
ALGBRC(43) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(42))+ CONSTS(27)*ALGBRC(38)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(44) = TERNRY(ALGBRC(38).GT.96.0000, ALGBRC(43) - ( ALGBRC(43)*(96.0000 - ALGBRC(38)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(42))+ CONSTS(27)*ALGBRC(38)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(46) = 1.00000/( ALGBRC(44)*CONSTS(21))
ALGBRC(47) = ALGBRC(45)+ALGBRC(46)+CONSTS(32)
ALGBRC(48) = (1.00000/ALGBRC(47))/60.0000
ALGBRC(49) = ( CONSTS(19)*- 1.00000*ALGBRC(48)*(ALGBRC(31) - ALGBRC(32)))/CONSTS(20)
ALGBRC(51) = ALGBRC(7)+( ALGBRC(49)*CONSTS(31))/2.00000
ALGBRC(53) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*ALGBRC(31) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*ALGBRC(31)) - ALGBRC(51)
ALGBRC(55) = 1.34000*CONSTS(13)*CONSTS(6)*ALGBRC(31) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(56) = 100.000*CONSTS(4)* CONSTS(6)*ALGBRC(31) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(57) = ALGBRC(55)/ALGBRC(56)+CONSTS(15)
ALGBRC(61) = ALGBRC(31) - ALGBRC(53)/ALGBRC(57)
ALGBRC(50) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(49)
ALGBRC(52) = ALGBRC(8)+( ALGBRC(50)*CONSTS(31))/2.00000
ALGBRC(54) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*ALGBRC(32) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*ALGBRC(32)) - ALGBRC(52)
ALGBRC(58) = 1.34000*CONSTS(14)*CONSTS(11)*ALGBRC(32) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(59) = 100.000*CONSTS(10)* CONSTS(11)*ALGBRC(32) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(60) = ALGBRC(58)/ALGBRC(59)+CONSTS(15)
ALGBRC(62) = ALGBRC(32) - ALGBRC(54)/ALGBRC(60)
ALGBRC(63) = ((arbitrary_log(ALGBRC(61), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(64) = 10.0000 ** ALGBRC(63)
ALGBRC(65) = ( 100.000*ALGBRC(64))/(1.00000+ALGBRC(64))
ALGBRC(69) = TERNRY(ALGBRC(65).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(65))))
ALGBRC(70) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(69))+ CONSTS(27)*ALGBRC(65)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(71) = TERNRY(ALGBRC(65).GT.96.0000, ALGBRC(70) - ( ALGBRC(70)*(96.0000 - ALGBRC(65)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(69))+ CONSTS(27)*ALGBRC(65)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(75) = 1.00000/( ALGBRC(71)*CONSTS(20))
ALGBRC(66) = ((arbitrary_log(ALGBRC(62), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(67) = 10.0000 ** ALGBRC(66)
ALGBRC(68) = ( 100.000*ALGBRC(67))/(1.00000+ALGBRC(67))
ALGBRC(72) = TERNRY(ALGBRC(68).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(68))))
ALGBRC(73) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(72))+ CONSTS(27)*ALGBRC(68)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(74) = TERNRY(ALGBRC(68).GT.96.0000, ALGBRC(73) - ( ALGBRC(73)*(96.0000 - ALGBRC(68)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(72))+ CONSTS(27)*ALGBRC(68)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(76) = 1.00000/( ALGBRC(74)*CONSTS(21))
ALGBRC(77) = ALGBRC(75)+ALGBRC(76)+CONSTS(33)
ALGBRC(78) = (1.00000/ALGBRC(77))/60.0000
ALGBRC(79) = ( CONSTS(19)*- 1.00000*ALGBRC(78)*(ALGBRC(61) - ALGBRC(62)))/CONSTS(20)
ALGBRC(81) = ALGBRC(7)+ ALGBRC(79)*CONSTS(31)
ALGBRC(83) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*ALGBRC(61) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*ALGBRC(61)) - ALGBRC(81)
ALGBRC(85) = 1.34000*CONSTS(13)*CONSTS(6)*ALGBRC(61) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(86) = 100.000*CONSTS(4)* CONSTS(6)*ALGBRC(61) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(87) = ALGBRC(85)/ALGBRC(86)+CONSTS(15)
ALGBRC(91) = ALGBRC(61) - ALGBRC(83)/ALGBRC(87)
ALGBRC(80) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(79)
ALGBRC(82) = ALGBRC(8)+ ALGBRC(80)*CONSTS(31)
ALGBRC(84) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*ALGBRC(62) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*ALGBRC(62)) - ALGBRC(82)
ALGBRC(88) = 1.34000*CONSTS(14)*CONSTS(11)*ALGBRC(62) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(89) = 100.000*CONSTS(10)* CONSTS(11)*ALGBRC(62) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(90) = ALGBRC(88)/ALGBRC(89)+CONSTS(15)
ALGBRC(92) = ALGBRC(62) - ALGBRC(84)/ALGBRC(90)
ALGBRC(93) = ((arbitrary_log(ALGBRC(91), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(94) = 10.0000 ** ALGBRC(93)
ALGBRC(95) = ( 100.000*ALGBRC(94))/(1.00000+ALGBRC(94))
ALGBRC(99) = TERNRY(ALGBRC(95).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(95))))
ALGBRC(100) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(99))+ CONSTS(27)*ALGBRC(95)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(101) = TERNRY(ALGBRC(95).GT.96.0000, ALGBRC(100) - ( ALGBRC(100)*(96.0000 - ALGBRC(95)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(99))+ CONSTS(27)*ALGBRC(95)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(105) = 1.00000/( ALGBRC(101)*CONSTS(20))
ALGBRC(96) = ((arbitrary_log(ALGBRC(92), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(97) = 10.0000 ** ALGBRC(96)
ALGBRC(98) = ( 100.000*ALGBRC(97))/(1.00000+ALGBRC(97))
ALGBRC(102) = TERNRY(ALGBRC(98).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(98))))
ALGBRC(103) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(102))+ CONSTS(27)*ALGBRC(98)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(104) = TERNRY(ALGBRC(98).GT.96.0000, ALGBRC(103) - ( ALGBRC(103)*(96.0000 - ALGBRC(98)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(102))+ CONSTS(27)*ALGBRC(98)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(106) = 1.00000/( ALGBRC(104)*CONSTS(21))
ALGBRC(107) = ALGBRC(105)+ALGBRC(106)+CONSTS(34)
ALGBRC(108) = (1.00000/ALGBRC(107))/60.0000
ALGBRC(109) = ( CONSTS(19)*- 1.00000*ALGBRC(108)*(ALGBRC(91) - ALGBRC(92)))/CONSTS(20)
ALGBRC(110) = ALGBRC(7)+ (1.00000/6.00000)*(ALGBRC(19)+ 2.00000*ALGBRC(49)+ 2.00000*ALGBRC(79)+ALGBRC(109))*CONSTS(31)
ALGBRC(112) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*ALGBRC(91) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*ALGBRC(91)) - ALGBRC(110)
ALGBRC(114) = 1.34000*CONSTS(13)*CONSTS(6)*ALGBRC(91) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(116) = 100.000*CONSTS(4)* CONSTS(6)*ALGBRC(91) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(118) = ALGBRC(114)/ALGBRC(116)+CONSTS(15)
ALGBRC(121) = ALGBRC(91) - ALGBRC(112)/ALGBRC(118)
RATES(1) = (ALGBRC(121) - STATES(1))/CONSTS(31)
ALGBRC(111) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(109)
ALGBRC(113) = ALGBRC(8)+ (1.00000/6.00000)*(ALGBRC(20)+ 2.00000*ALGBRC(50)+ 2.00000*ALGBRC(80)+ALGBRC(111))*CONSTS(31)
ALGBRC(115) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*ALGBRC(92) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*ALGBRC(92)) - ALGBRC(113)
ALGBRC(117) = 1.34000*CONSTS(14)*CONSTS(11)*ALGBRC(92) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(119) = 100.000*CONSTS(10)* CONSTS(11)*ALGBRC(92) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(120) = ALGBRC(117)/ALGBRC(119)+CONSTS(15)
ALGBRC(122) = ALGBRC(92) - ALGBRC(115)/ALGBRC(120)
RATES(2) = (ALGBRC(122) - STATES(2))/CONSTS(31)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(7) = (( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*STATES(1)
ALGBRC(1) = ((arbitrary_log(STATES(1), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(2) = 10.0000 ** ALGBRC(1)
ALGBRC(3) = ( 100.000*ALGBRC(2))/(1.00000+ALGBRC(2))
ALGBRC(9) = TERNRY(ALGBRC(3).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(3))))
ALGBRC(10) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(9))+ CONSTS(27)*ALGBRC(3)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(11) = TERNRY(ALGBRC(3).GT.96.0000, ALGBRC(10) - ( ALGBRC(10)*(96.0000 - ALGBRC(3)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(9))+ CONSTS(27)*ALGBRC(3)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(15) = 1.00000/( ALGBRC(11)*CONSTS(20))
ALGBRC(4) = ((arbitrary_log(STATES(2), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(5) = 10.0000 ** ALGBRC(4)
ALGBRC(6) = ( 100.000*ALGBRC(5))/(1.00000+ALGBRC(5))
ALGBRC(12) = TERNRY(ALGBRC(6).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(6))))
ALGBRC(13) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(12))+ CONSTS(27)*ALGBRC(6)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(14) = TERNRY(ALGBRC(6).GT.96.0000, ALGBRC(13) - ( ALGBRC(13)*(96.0000 - ALGBRC(6)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(12))+ CONSTS(27)*ALGBRC(6)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(16) = 1.00000/( ALGBRC(14)*CONSTS(21))
ALGBRC(17) = ALGBRC(15)+ALGBRC(16)+CONSTS(29)
ALGBRC(18) = (1.00000/ALGBRC(17))/60.0000
ALGBRC(19) = ( CONSTS(19)*- 1.00000*ALGBRC(18)*(STATES(1) - STATES(2)))/CONSTS(20)
ALGBRC(21) = ALGBRC(7)+( ALGBRC(19)*CONSTS(31))/2.00000
ALGBRC(23) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*STATES(1)) - ALGBRC(21)
ALGBRC(25) = 1.34000*CONSTS(13)*CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(26) = 100.000*CONSTS(4)* CONSTS(6)*STATES(1) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(27) = ALGBRC(25)/ALGBRC(26)+CONSTS(15)
ALGBRC(31) = STATES(1) - ALGBRC(23)/ALGBRC(27)
ALGBRC(8) = (( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*STATES(2)
ALGBRC(20) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(19)
ALGBRC(22) = ALGBRC(8)+( ALGBRC(20)*CONSTS(31))/2.00000
ALGBRC(24) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*STATES(2)) - ALGBRC(22)
ALGBRC(28) = 1.34000*CONSTS(14)*CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(29) = 100.000*CONSTS(10)* CONSTS(11)*STATES(2) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(30) = ALGBRC(28)/ALGBRC(29)+CONSTS(15)
ALGBRC(32) = STATES(2) - ALGBRC(24)/ALGBRC(30)
ALGBRC(33) = ((arbitrary_log(ALGBRC(31), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(34) = 10.0000 ** ALGBRC(33)
ALGBRC(35) = ( 100.000*ALGBRC(34))/(1.00000+ALGBRC(34))
ALGBRC(39) = TERNRY(ALGBRC(35).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(35))))
ALGBRC(40) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(39))+ CONSTS(27)*ALGBRC(35)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(41) = TERNRY(ALGBRC(35).GT.96.0000, ALGBRC(40) - ( ALGBRC(40)*(96.0000 - ALGBRC(35)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(39))+ CONSTS(27)*ALGBRC(35)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(45) = 1.00000/( ALGBRC(41)*CONSTS(20))
ALGBRC(36) = ((arbitrary_log(ALGBRC(32), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(37) = 10.0000 ** ALGBRC(36)
ALGBRC(38) = ( 100.000*ALGBRC(37))/(1.00000+ALGBRC(37))
ALGBRC(42) = TERNRY(ALGBRC(38).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(38))))
ALGBRC(43) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(42))+ CONSTS(27)*ALGBRC(38)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(44) = TERNRY(ALGBRC(38).GT.96.0000, ALGBRC(43) - ( ALGBRC(43)*(96.0000 - ALGBRC(38)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(42))+ CONSTS(27)*ALGBRC(38)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(46) = 1.00000/( ALGBRC(44)*CONSTS(21))
ALGBRC(47) = ALGBRC(45)+ALGBRC(46)+CONSTS(32)
ALGBRC(48) = (1.00000/ALGBRC(47))/60.0000
ALGBRC(49) = ( CONSTS(19)*- 1.00000*ALGBRC(48)*(ALGBRC(31) - ALGBRC(32)))/CONSTS(20)
ALGBRC(51) = ALGBRC(7)+( ALGBRC(49)*CONSTS(31))/2.00000
ALGBRC(53) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*ALGBRC(31) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*ALGBRC(31)) - ALGBRC(51)
ALGBRC(55) = 1.34000*CONSTS(13)*CONSTS(6)*ALGBRC(31) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(56) = 100.000*CONSTS(4)* CONSTS(6)*ALGBRC(31) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(57) = ALGBRC(55)/ALGBRC(56)+CONSTS(15)
ALGBRC(61) = ALGBRC(31) - ALGBRC(53)/ALGBRC(57)
ALGBRC(50) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(49)
ALGBRC(52) = ALGBRC(8)+( ALGBRC(50)*CONSTS(31))/2.00000
ALGBRC(54) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*ALGBRC(32) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*ALGBRC(32)) - ALGBRC(52)
ALGBRC(58) = 1.34000*CONSTS(14)*CONSTS(11)*ALGBRC(32) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(59) = 100.000*CONSTS(10)* CONSTS(11)*ALGBRC(32) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(60) = ALGBRC(58)/ALGBRC(59)+CONSTS(15)
ALGBRC(62) = ALGBRC(32) - ALGBRC(54)/ALGBRC(60)
ALGBRC(63) = ((arbitrary_log(ALGBRC(61), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(64) = 10.0000 ** ALGBRC(63)
ALGBRC(65) = ( 100.000*ALGBRC(64))/(1.00000+ALGBRC(64))
ALGBRC(69) = TERNRY(ALGBRC(65).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(65))))
ALGBRC(70) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(69))+ CONSTS(27)*ALGBRC(65)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(71) = TERNRY(ALGBRC(65).GT.96.0000, ALGBRC(70) - ( ALGBRC(70)*(96.0000 - ALGBRC(65)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(69))+ CONSTS(27)*ALGBRC(65)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(75) = 1.00000/( ALGBRC(71)*CONSTS(20))
ALGBRC(66) = ((arbitrary_log(ALGBRC(62), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(67) = 10.0000 ** ALGBRC(66)
ALGBRC(68) = ( 100.000*ALGBRC(67))/(1.00000+ALGBRC(67))
ALGBRC(72) = TERNRY(ALGBRC(68).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(68))))
ALGBRC(73) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(72))+ CONSTS(27)*ALGBRC(68)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(74) = TERNRY(ALGBRC(68).GT.96.0000, ALGBRC(73) - ( ALGBRC(73)*(96.0000 - ALGBRC(68)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(72))+ CONSTS(27)*ALGBRC(68)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(76) = 1.00000/( ALGBRC(74)*CONSTS(21))
ALGBRC(77) = ALGBRC(75)+ALGBRC(76)+CONSTS(33)
ALGBRC(78) = (1.00000/ALGBRC(77))/60.0000
ALGBRC(79) = ( CONSTS(19)*- 1.00000*ALGBRC(78)*(ALGBRC(61) - ALGBRC(62)))/CONSTS(20)
ALGBRC(81) = ALGBRC(7)+ ALGBRC(79)*CONSTS(31)
ALGBRC(83) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*ALGBRC(61) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*ALGBRC(61)) - ALGBRC(81)
ALGBRC(85) = 1.34000*CONSTS(13)*CONSTS(6)*ALGBRC(61) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(86) = 100.000*CONSTS(4)* CONSTS(6)*ALGBRC(61) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(87) = ALGBRC(85)/ALGBRC(86)+CONSTS(15)
ALGBRC(91) = ALGBRC(61) - ALGBRC(83)/ALGBRC(87)
ALGBRC(80) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(79)
ALGBRC(82) = ALGBRC(8)+ ALGBRC(80)*CONSTS(31)
ALGBRC(84) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*ALGBRC(62) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*ALGBRC(62)) - ALGBRC(82)
ALGBRC(88) = 1.34000*CONSTS(14)*CONSTS(11)*ALGBRC(62) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(89) = 100.000*CONSTS(10)* CONSTS(11)*ALGBRC(62) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(90) = ALGBRC(88)/ALGBRC(89)+CONSTS(15)
ALGBRC(92) = ALGBRC(62) - ALGBRC(84)/ALGBRC(90)
ALGBRC(93) = ((arbitrary_log(ALGBRC(91), 10) - CONSTS(1))+ CONSTS(3)*(CONSTS(7) - 7.40000))/CONSTS(4)
ALGBRC(94) = 10.0000 ** ALGBRC(93)
ALGBRC(95) = ( 100.000*ALGBRC(94))/(1.00000+ALGBRC(94))
ALGBRC(99) = TERNRY(ALGBRC(95).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(95))))
ALGBRC(100) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(99))+ CONSTS(27)*ALGBRC(95)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18)
ALGBRC(101) = TERNRY(ALGBRC(95).GT.96.0000, ALGBRC(100) - ( ALGBRC(100)*(96.0000 - ALGBRC(95)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(99))+ CONSTS(27)*ALGBRC(95)) - CONSTS(28))*1.34000)/100.000)*CONSTS(13))/0.200000)*CONSTS(16))/CONSTS(18))
ALGBRC(105) = 1.00000/( ALGBRC(101)*CONSTS(20))
ALGBRC(96) = ((arbitrary_log(ALGBRC(92), 10) - CONSTS(8))+ CONSTS(9)*(CONSTS(12) - 7.40000))/CONSTS(10)
ALGBRC(97) = 10.0000 ** ALGBRC(96)
ALGBRC(98) = ( 100.000*ALGBRC(97))/(1.00000+ALGBRC(97))
ALGBRC(102) = TERNRY(ALGBRC(98).GT.96.0000, EXP( - 1.00000*CONSTS(26)*(100.000 - 96.0000)), EXP( - 1.00000*CONSTS(26)*(100.000 - ALGBRC(98))))
ALGBRC(103) = ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(102))+ CONSTS(27)*ALGBRC(98)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18)
ALGBRC(104) = TERNRY(ALGBRC(98).GT.96.0000, ALGBRC(103) - ( ALGBRC(103)*(96.0000 - ALGBRC(98)))/- 4.00000, ( (( (( (( CONSTS(25)*(1.00000 - ALGBRC(102))+ CONSTS(27)*ALGBRC(98)) - CONSTS(28))*1.34000)/100.000)*CONSTS(14))/0.200000)*CONSTS(17))/CONSTS(18))
ALGBRC(106) = 1.00000/( ALGBRC(104)*CONSTS(21))
ALGBRC(107) = ALGBRC(105)+ALGBRC(106)+CONSTS(34)
ALGBRC(108) = (1.00000/ALGBRC(107))/60.0000
ALGBRC(109) = ( CONSTS(19)*- 1.00000*ALGBRC(108)*(ALGBRC(91) - ALGBRC(92)))/CONSTS(20)
ALGBRC(110) = ALGBRC(7)+ (1.00000/6.00000)*(ALGBRC(19)+ 2.00000*ALGBRC(49)+ 2.00000*ALGBRC(79)+ALGBRC(109))*CONSTS(31)
ALGBRC(112) = ((( 1.34000*CONSTS(13))/100.000)/( CONSTS(6)*ALGBRC(91) ** - 1.00000/CONSTS(4)+1.00000)+ CONSTS(15)*ALGBRC(91)) - ALGBRC(110)
ALGBRC(114) = 1.34000*CONSTS(13)*CONSTS(6)*ALGBRC(91) ** - 1.00000/CONSTS(4) - 1.00000
ALGBRC(116) = 100.000*CONSTS(4)* CONSTS(6)*ALGBRC(91) ** - 1.00000/CONSTS(4)+1.00000 ** 2.00000
ALGBRC(118) = ALGBRC(114)/ALGBRC(116)+CONSTS(15)
ALGBRC(121) = ALGBRC(91) - ALGBRC(112)/ALGBRC(118)
ALGBRC(111) = (( (- 1.00000/CONSTS(19))*CONSTS(20))/CONSTS(21))*ALGBRC(109)
ALGBRC(113) = ALGBRC(8)+ (1.00000/6.00000)*(ALGBRC(20)+ 2.00000*ALGBRC(50)+ 2.00000*ALGBRC(80)+ALGBRC(111))*CONSTS(31)
ALGBRC(115) = ((( 1.34000*CONSTS(14))/100.000)/( CONSTS(11)*ALGBRC(92) ** - 1.00000/CONSTS(10)+1.00000)+ CONSTS(15)*ALGBRC(92)) - ALGBRC(113)
ALGBRC(117) = 1.34000*CONSTS(14)*CONSTS(11)*ALGBRC(92) ** - 1.00000/CONSTS(10) - 1.00000
ALGBRC(119) = 100.000*CONSTS(10)* CONSTS(11)*ALGBRC(92) ** - 1.00000/CONSTS(10)+1.00000 ** 2.00000
ALGBRC(120) = ALGBRC(117)/ALGBRC(119)+CONSTS(15)
ALGBRC(122) = ALGBRC(92) - ALGBRC(115)/ALGBRC(120)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END