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