C
C There are a total of 34 entries in the algebraic variable array.
C There are a total of 9 entries in each of the rate and state variable arrays.
C There are a total of 54 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is isotonic in component parameters (dimensionless).
C CONSTS(2) is alpha_1 in component parameters (per_micrometre).
C CONSTS(3) is beta_1 in component parameters (millinewton).
C CONSTS(4) is alpha_2 in component parameters (per_micrometre).
C CONSTS(5) is beta_2 in component parameters (millinewton).
C CONSTS(6) is alpha_3 in component parameters (per_micrometre).
C CONSTS(7) is beta_3 in component parameters (millinewton).
C CONSTS(8) is lambda in component parameters (millinewton).
C CONSTS(9) is A_half in component parameters (micromole).
C CONSTS(10) is mu in component parameters (dimensionless).
C CONSTS(11) is chi in component parameters (dimensionless).
C CONSTS(12) is chi_0 in component parameters (dimensionless).
C CONSTS(13) is m_0 in component parameters (dimensionless).
C CONSTS(14) is v_max in component parameters (micrometre_per_second).
C CONSTS(15) is a in component parameters (dimensionless).
C CONSTS(16) is d_h in component parameters (dimensionless).
C CONSTS(17) is alpha_P in component parameters (dimensionless).
C ALGBRC(10) is l in component length (micrometre).
C ALGBRC(6) is F_muscle in component force (millinewton).
C ALGBRC(12) is flag in component isotonic (dimensionless).
C CONSTS(18) is F_afterload in component isotonic (millinewton).
C ALGBRC(14) is isotonic_mode in component isotonic (dimensionless).
C CONSTS(19) is l_0 in component isotonic (micrometre).
C CONSTS(20) is S_0 in component parameters_izakov_et_al_1991 (micrometre).
C ALGBRC(1) is q_v in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(21) is q_1 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(22) is q_2 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(23) is q_3 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(24) is q_4 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(25) is v_star in component parameters_izakov_et_al_1991 (micrometre_per_second).
C CONSTS(51) is v_1 in component parameters_izakov_et_al_1991 (micrometre_per_second).
C CONSTS(26) is alpha_G in component parameters_izakov_et_al_1991 (dimensionless).
C CONSTS(27) is a_on in component parameters_izakov_et_al_1991 (per_second_micromole).
C CONSTS(28) is a_off in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(29) is k_A in component parameters_izakov_et_al_1991 (per_micromole).
C STATES(1) is v in component CE_velocity (micrometre_per_second).
C CONSTS(30) is alpha_Q in component parameters_izakov_et_al_1991 (dimensionless).
C CONSTS(31) is beta_Q in component parameters_izakov_et_al_1991 (dimensionless).
C ALGBRC(32) is F_CE in component force (millinewton).
C ALGBRC(5) is F_XSE in component force (millinewton).
C ALGBRC(2) is F_SE in component force (millinewton).
C ALGBRC(3) is F_PE in component force (millinewton).
C STATES(2) is N in component crossbridge_kinetics (dimensionless).
C ALGBRC(18) is k_P_vis in component CE_velocity (millinewton_second_per_micrometre).
C ALGBRC(21) is k_S_vis in component PE_velocity (millinewton_second_per_micrometre).
C STATES(3) is w in component PE_velocity (micrometre_per_second).
C STATES(4) is l_1 in component length (micrometre).
C STATES(5) is l_2 in component length (micrometre).
C STATES(6) is l_3 in component length (micrometre).
C ALGBRC(31) is p_v in component average_crossbridge_force (dimensionless).
C ALGBRC(28) is K_chi in component crossbridge_kinetics (per_second).
C ALGBRC(7) is M_A in component crossbridge_kinetics (dimensionless).
C ALGBRC(8) is n_1 in component crossbridge_kinetics (dimensionless).
C ALGBRC(9) is L_oz in component crossbridge_kinetics (dimensionless).
C ALGBRC(26) is k_p_v in component crossbridge_kinetics (per_second).
C ALGBRC(27) is k_m_v in component crossbridge_kinetics (per_second).
C STATES(7) is A in component calcium_handling (micromole).
C ALGBRC(25) is G_star in component average_crossbridge_force (dimensionless).
C ALGBRC(22) is P_star in component average_crossbridge_force (dimensionless).
C CONSTS(32) is v_0 in component crossbridge_kinetics (micrometre_per_second).
C CONSTS(33) is q_star in component crossbridge_kinetics (per_second).
C ALGBRC(4) is dl_1_dt in component length (micrometre_per_second).
C ALGBRC(23) is dl_2_dt in component length (micrometre_per_second).
C ALGBRC(24) is dl_3_dt in component length (micrometre_per_second).
C ALGBRC(19) is phi_chi_2 in component CE_velocity (micrometre_per_second).
C ALGBRC(34) is phi_chi in component CE_velocity (micrometre_per_second2).
C ALGBRC(33) is p_prime_v in component average_crossbridge_force (second_per_micrometre).
C CONSTS(34) is alpha_P_lengthening in component CE_velocity (per_micrometre).
C CONSTS(35) is beta_P_lengthening in component CE_velocity (millinewton_second_per_micrometre).
C CONSTS(36) is alpha_P_shortening in component CE_velocity (per_micrometre).
C CONSTS(37) is beta_P_shortening in component CE_velocity (millinewton_second_per_micrometre).
C ALGBRC(16) is alp_p in component CE_velocity (per_micrometre).
C CONSTS(38) is alpha_S_lengthening in component PE_velocity (per_micrometre).
C CONSTS(39) is beta_S_lengthening in component PE_velocity (millinewton_second_per_micrometre).
C CONSTS(40) is alpha_S_shortening in component PE_velocity (per_micrometre).
C CONSTS(41) is beta_S_shortening in component PE_velocity (millinewton_second_per_micrometre).
C ALGBRC(20) is alp_s in component PE_velocity (per_micrometre).
C CONSTS(54) is gamma in component average_crossbridge_force (dimensionless).
C CONSTS(52) is case_1 in component average_crossbridge_force (second_per_micrometre).
C ALGBRC(29) is case_2 in component average_crossbridge_force (second_per_micrometre).
C CONSTS(53) is case_3 in component average_crossbridge_force (second_per_micrometre).
C ALGBRC(30) is case_4 in component average_crossbridge_force (second_per_micrometre).
C ALGBRC(15) is dA_dt in component calcium_handling (micromole_per_second).
C ALGBRC(11) is N_A in component calcium_handling (dimensionless).
C ALGBRC(13) is pi_N_A in component calcium_handling (dimensionless).
C STATES(8) is B in component calcium_handling (micromole).
C ALGBRC(17) is dB_dt in component calcium_handling (micromole_per_second).
C STATES(9) is Ca_C in component calcium_handling (micromole).
C CONSTS(42) is A_tot in component calcium_handling (micromole).
C CONSTS(43) is B_tot in component calcium_handling (micromole).
C CONSTS(44) is b_on in component calcium_handling (per_second_micromole).
C CONSTS(45) is b_off in component calcium_handling (per_second).
C CONSTS(46) is a_c in component calcium_handling (per_second2).
C CONSTS(47) is r_Ca in component calcium_handling (per_second).
C CONSTS(48) is q_Ca in component calcium_handling (per_micromole).
C CONSTS(49) is t_d in component calcium_handling (second).
C CONSTS(50) is Ca_m in component calcium_handling (micromole).
C RATES(2) is d/dt N in component crossbridge_kinetics (dimensionless).
C RATES(4) is d/dt l_1 in component length (micrometre).
C RATES(5) is d/dt l_2 in component length (micrometre).
C RATES(6) is d/dt l_3 in component length (micrometre).
C RATES(1) is d/dt v in component CE_velocity (micrometre_per_second).
C RATES(3) is d/dt w in component PE_velocity (micrometre_per_second).
C RATES(7) is d/dt A in component calcium_handling (micromole).
C RATES(8) is d/dt B in component calcium_handling (micromole).
C RATES(9) is d/dt Ca_C in component calcium_handling (micromole).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0
      CONSTS(2) = 19
      CONSTS(3) = 0.29
      CONSTS(4) = 14.6
      CONSTS(5) = 0.000924
      CONSTS(6) = 48
      CONSTS(7) = 0.01
      CONSTS(8) = 96
      CONSTS(9) = 42
      CONSTS(10) = 3
      CONSTS(11) = 0.705
      CONSTS(12) = 3
      CONSTS(13) = 0.9
      CONSTS(14) = 5.6
      CONSTS(15) = 0.25
      CONSTS(16) = 0.5
      CONSTS(17) = 4
      CONSTS(18) = 2
      CONSTS(19) = 0.527
      CONSTS(20) = 1.14
      CONSTS(21) = 17.3
      CONSTS(22) = 259
      CONSTS(23) = 17.3
      CONSTS(24) = 15
      CONSTS(25) = 5.3035675
      CONSTS(26) = 1
      CONSTS(27) = 32.85
      CONSTS(28) = 290
      CONSTS(29) = 0.04
      STATES(1) = 0
      CONSTS(30) = 10
      CONSTS(31) = 5000
      STATES(2) = 0.0001
      STATES(3) = 0
      STATES(4) = 0.437
      STATES(5) = 0.439
      STATES(6) = 0.089
      STATES(7) = 0.7
      CONSTS(32) = 560
      CONSTS(33) = 1000
      CONSTS(34) = 16
      CONSTS(35) = 0.0015
      CONSTS(36) = 16
      CONSTS(37) = 0.0015
      CONSTS(38) = 39
      CONSTS(39) = 0.008
      CONSTS(40) = 46
      CONSTS(41) = 0.006
      STATES(8) = 0
      STATES(9) = 0
      CONSTS(42) = 70
      CONSTS(43) = 28
      CONSTS(44) = 37.143
      CONSTS(45) = 182
      CONSTS(46) = 5200
      CONSTS(47) = 650
      CONSTS(48) = 0.714
      CONSTS(49) = 0.033
      CONSTS(50) = 2.1
      CONSTS(51) = CONSTS(14)/10.0000
      CONSTS(52) = ( CONSTS(15)*(0.400000+ 0.400000*CONSTS(15)))/( CONSTS(14)* (CONSTS(15)+1.00000)*0.400000 ** 2.00000)
      CONSTS(53) = ( 0.400000*CONSTS(15)+1.00000)/( CONSTS(15)*CONSTS(14))
      CONSTS(54) = ( CONSTS(15)*CONSTS(16)*CONSTS(51)/CONSTS(14) ** 2.00000)/( 3.00000*CONSTS(15)*CONSTS(16) - ( (CONSTS(15)+1.00000)*CONSTS(51))/CONSTS(14))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = STATES(1)
      RATES(4) = ALGBRC(4)
      ALGBRC(9) = (STATES(4)+CONSTS(20))/(0.460000+CONSTS(20))
      ALGBRC(11) = ( CONSTS(42)*STATES(2))/( ALGBRC(9)*STATES(7))
      ALGBRC(13) = TERNRY(ALGBRC(11).GE.1.00000, 1.00000, 0.0200000 ** ALGBRC(11))
      ALGBRC(15) =  CONSTS(27)*(CONSTS(42) - STATES(7))*STATES(9) -  CONSTS(28)*EXP( - CONSTS(29)*STATES(7))*ALGBRC(13)*STATES(7)
      RATES(7) = ALGBRC(15)
      ALGBRC(17) =  CONSTS(44)*(CONSTS(43) - STATES(8))*STATES(9) -  CONSTS(45)*STATES(8)
      RATES(8) = ALGBRC(17)
      RATES(9) = TERNRY(VOI.LT.CONSTS(49),  4.00000*CONSTS(46)*CONSTS(50)*VOI*(1.00000 - EXP( - CONSTS(46)*VOI ** 2.00000))*EXP( - CONSTS(46)*VOI ** 2.00000), (- ALGBRC(15) - ALGBRC(17)) -  CONSTS(47)*EXP( - CONSTS(48)*STATES(9))*STATES(9))
      ALGBRC(21) = TERNRY(STATES(3).LE.STATES(1),  CONSTS(39)*EXP( CONSTS(38)*(STATES(5) - STATES(4))),  CONSTS(41)*EXP( CONSTS(40)*(STATES(5) - STATES(4))))
      ALGBRC(10) = STATES(5)+STATES(6)
      ALGBRC(5) =  CONSTS(7)*(EXP( CONSTS(6)*STATES(6)) - 1.00000)
      ALGBRC(6) = ALGBRC(5)
      ALGBRC(12) = TERNRY(ALGBRC(10).GE.CONSTS(19).AND.VOI.LT.0.150000, 0.00000, 1.00000)
      ALGBRC(14) = TERNRY(CONSTS(1).EQ.0.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.AND.ALGBRC(6).GE.CONSTS(18), 1.00000, TERNRY(CONSTS(1).EQ.1.00000.AND.ALGBRC(10).GE.CONSTS(19).AND.ALGBRC(12).EQ.1.00000, 0.00000, 0.0/0.0)
      ALGBRC(19) = TERNRY(ALGBRC(14).EQ.1.00000, ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))), ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6))))
      ALGBRC(23) = TERNRY(ALGBRC(21).EQ.0.00000, ALGBRC(19), STATES(3))
      RATES(5) = ALGBRC(23)
      ALGBRC(24) = TERNRY(ALGBRC(14).EQ.1.00000, 0.00000, TERNRY(ALGBRC(14).EQ.0.00000.AND.ALGBRC(21).EQ.0.00000, - ALGBRC(19), - STATES(3))
      RATES(6) = ALGBRC(24)
      ALGBRC(7) = STATES(7) ** CONSTS(10)/(STATES(7) ** CONSTS(10)+CONSTS(9) ** CONSTS(10))
      ALGBRC(8) =  0.600000*STATES(4)+0.500000
      ALGBRC(1) = TERNRY(STATES(1).LE.0.00000, CONSTS(21) - ( CONSTS(22)*STATES(1))/CONSTS(14), TERNRY(STATES(1).LE.CONSTS(25).AND.0.00000.LT.STATES(1), ( (CONSTS(24) - CONSTS(23))*STATES(1))/CONSTS(25)+CONSTS(23), CONSTS(24)/1.00000+( CONSTS(31)*(STATES(1) - CONSTS(25)))/CONSTS(14) ** CONSTS(30))
      ALGBRC(22) = TERNRY(STATES(1).LE.0.00000, ( CONSTS(15)*(1.00000+STATES(1)/CONSTS(14)))/(CONSTS(15) - STATES(1)/CONSTS(14)), (1.00000+CONSTS(16)) - ( CONSTS(16) ** 2.00000*CONSTS(15))/( (( CONSTS(15)*CONSTS(16))/CONSTS(54))*STATES(1)/CONSTS(14) ** 2.00000+( (CONSTS(15)+1.00000)*STATES(1))/CONSTS(14)+ CONSTS(15)*CONSTS(16)))
      ALGBRC(25) = TERNRY(- CONSTS(14).LE.STATES(1).AND.STATES(1).LE.0.00000, 1.00000+( 0.600000*STATES(1))/CONSTS(14), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(51), ALGBRC(22)/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000), ( ALGBRC(22)*EXP( - CONSTS(26)*(STATES(1) - CONSTS(51))/CONSTS(14) ** CONSTS(17)))/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000))
      ALGBRC(26) =  CONSTS(11)*CONSTS(12)*ALGBRC(1)*CONSTS(13)*ALGBRC(25)
      ALGBRC(27) = TERNRY(STATES(1).LE.CONSTS(25),  CONSTS(12)*ALGBRC(1)*(1.00000 -  CONSTS(11)*CONSTS(13)*ALGBRC(25)),  CONSTS(12)*( CONSTS(24)*(1.00000 -  CONSTS(11)*CONSTS(13)*ALGBRC(25))+( CONSTS(33)*(STATES(1) - CONSTS(25)))/(CONSTS(32) - CONSTS(25))))
      ALGBRC(28) =  ALGBRC(26)*ALGBRC(7)*ALGBRC(8)*ALGBRC(9)*(1.00000 - STATES(2)) -  ALGBRC(27)*STATES(2)
      RATES(2) = ALGBRC(28)
      ALGBRC(18) = TERNRY(STATES(1).LE.0.00000,  CONSTS(35)*EXP( CONSTS(34)*STATES(4)),  CONSTS(37)*EXP( CONSTS(36)*STATES(4)))
      ALGBRC(31) = ALGBRC(22)/ALGBRC(25)
      ALGBRC(29) = ( CONSTS(15)*1.00000*(1.00000+ 0.400000*CONSTS(15)+( 1.20000*STATES(1))/CONSTS(14)+ 0.600000*STATES(1)/CONSTS(14) ** 2.00000))/( CONSTS(14)* (CONSTS(15) - STATES(1)/CONSTS(14))*(1.00000+( 0.600000*STATES(1))/CONSTS(14)) ** 2.00000)
      ALGBRC(30) =  (1.00000/CONSTS(14))*EXP( - CONSTS(26)*STATES(1)/CONSTS(14) - CONSTS(51)/CONSTS(14) ** CONSTS(17))*(( 0.400000*CONSTS(15)+1.00000)/CONSTS(15)+ CONSTS(26)*CONSTS(17)*(1.00000+( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14)))*STATES(1)/CONSTS(14) - CONSTS(51)/CONSTS(14) ** CONSTS(17) - 1.00000)
      ALGBRC(33) = TERNRY(STATES(1).LE.- CONSTS(14), CONSTS(52), TERNRY(- CONSTS(14).LT.STATES(1).AND.STATES(1).LE.0.00000, ALGBRC(29), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(51), CONSTS(53), ALGBRC(30))
      ALGBRC(16) = TERNRY(STATES(1).LE.0.00000, CONSTS(34), CONSTS(36))
      ALGBRC(34) = TERNRY(ALGBRC(14).EQ.1.00000, - ( CONSTS(8)*ALGBRC(28)*ALGBRC(31)+ ALGBRC(16)*ALGBRC(18)*STATES(1) ** 2.00000+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(33)+ALGBRC(18)), - ( CONSTS(8)*ALGBRC(28)*ALGBRC(31)+ ALGBRC(16)*ALGBRC(18)*STATES(1) ** 2.00000+ ( CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6)))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(33)+ALGBRC(18)))
      RATES(1) = TERNRY(ALGBRC(21).EQ.0.00000, ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*(ALGBRC(19) - STATES(1)) - ( CONSTS(8)*ALGBRC(28)*ALGBRC(31)+ ALGBRC(16)*ALGBRC(18)*STATES(1) ** 2.00000))/( CONSTS(8)*STATES(2)*ALGBRC(33)+ALGBRC(18)), ALGBRC(34))
      ALGBRC(20) = TERNRY(STATES(3).LE.STATES(1), CONSTS(38), CONSTS(40))
      RATES(3) = TERNRY(ALGBRC(14).EQ.1.00000, (( ALGBRC(21)*(ALGBRC(34) -  ALGBRC(20)*STATES(3) - STATES(1) ** 2.00000) -  CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*(STATES(3) - STATES(1))) -  CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))*STATES(3))/ALGBRC(21), (ALGBRC(34) -  ALGBRC(20)*STATES(3) - STATES(1) ** 2.00000) - ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*(STATES(3) - STATES(1))+ ( CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6)))*STATES(3))/ALGBRC(21))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = STATES(1)
      ALGBRC(9) = (STATES(4)+CONSTS(20))/(0.460000+CONSTS(20))
      ALGBRC(11) = ( CONSTS(42)*STATES(2))/( ALGBRC(9)*STATES(7))
      ALGBRC(13) = TERNRY(ALGBRC(11).GE.1.00000, 1.00000, 0.0200000 ** ALGBRC(11))
      ALGBRC(15) =  CONSTS(27)*(CONSTS(42) - STATES(7))*STATES(9) -  CONSTS(28)*EXP( - CONSTS(29)*STATES(7))*ALGBRC(13)*STATES(7)
      ALGBRC(17) =  CONSTS(44)*(CONSTS(43) - STATES(8))*STATES(9) -  CONSTS(45)*STATES(8)
      ALGBRC(21) = TERNRY(STATES(3).LE.STATES(1),  CONSTS(39)*EXP( CONSTS(38)*(STATES(5) - STATES(4))),  CONSTS(41)*EXP( CONSTS(40)*(STATES(5) - STATES(4))))
      ALGBRC(10) = STATES(5)+STATES(6)
      ALGBRC(5) =  CONSTS(7)*(EXP( CONSTS(6)*STATES(6)) - 1.00000)
      ALGBRC(6) = ALGBRC(5)
      ALGBRC(12) = TERNRY(ALGBRC(10).GE.CONSTS(19).AND.VOI.LT.0.150000, 0.00000, 1.00000)
      ALGBRC(14) = TERNRY(CONSTS(1).EQ.0.00000, 0.00000, TERNRY(CONSTS(1).EQ.1.00000.AND.ALGBRC(6).GE.CONSTS(18), 1.00000, TERNRY(CONSTS(1).EQ.1.00000.AND.ALGBRC(10).GE.CONSTS(19).AND.ALGBRC(12).EQ.1.00000, 0.00000, 0.0/0.0)
      ALGBRC(19) = TERNRY(ALGBRC(14).EQ.1.00000, ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))), ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6))))
      ALGBRC(23) = TERNRY(ALGBRC(21).EQ.0.00000, ALGBRC(19), STATES(3))
      ALGBRC(24) = TERNRY(ALGBRC(14).EQ.1.00000, 0.00000, TERNRY(ALGBRC(14).EQ.0.00000.AND.ALGBRC(21).EQ.0.00000, - ALGBRC(19), - STATES(3))
      ALGBRC(7) = STATES(7) ** CONSTS(10)/(STATES(7) ** CONSTS(10)+CONSTS(9) ** CONSTS(10))
      ALGBRC(8) =  0.600000*STATES(4)+0.500000
      ALGBRC(1) = TERNRY(STATES(1).LE.0.00000, CONSTS(21) - ( CONSTS(22)*STATES(1))/CONSTS(14), TERNRY(STATES(1).LE.CONSTS(25).AND.0.00000.LT.STATES(1), ( (CONSTS(24) - CONSTS(23))*STATES(1))/CONSTS(25)+CONSTS(23), CONSTS(24)/1.00000+( CONSTS(31)*(STATES(1) - CONSTS(25)))/CONSTS(14) ** CONSTS(30))
      ALGBRC(22) = TERNRY(STATES(1).LE.0.00000, ( CONSTS(15)*(1.00000+STATES(1)/CONSTS(14)))/(CONSTS(15) - STATES(1)/CONSTS(14)), (1.00000+CONSTS(16)) - ( CONSTS(16) ** 2.00000*CONSTS(15))/( (( CONSTS(15)*CONSTS(16))/CONSTS(54))*STATES(1)/CONSTS(14) ** 2.00000+( (CONSTS(15)+1.00000)*STATES(1))/CONSTS(14)+ CONSTS(15)*CONSTS(16)))
      ALGBRC(25) = TERNRY(- CONSTS(14).LE.STATES(1).AND.STATES(1).LE.0.00000, 1.00000+( 0.600000*STATES(1))/CONSTS(14), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(51), ALGBRC(22)/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000), ( ALGBRC(22)*EXP( - CONSTS(26)*(STATES(1) - CONSTS(51))/CONSTS(14) ** CONSTS(17)))/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000))
      ALGBRC(26) =  CONSTS(11)*CONSTS(12)*ALGBRC(1)*CONSTS(13)*ALGBRC(25)
      ALGBRC(27) = TERNRY(STATES(1).LE.CONSTS(25),  CONSTS(12)*ALGBRC(1)*(1.00000 -  CONSTS(11)*CONSTS(13)*ALGBRC(25)),  CONSTS(12)*( CONSTS(24)*(1.00000 -  CONSTS(11)*CONSTS(13)*ALGBRC(25))+( CONSTS(33)*(STATES(1) - CONSTS(25)))/(CONSTS(32) - CONSTS(25))))
      ALGBRC(28) =  ALGBRC(26)*ALGBRC(7)*ALGBRC(8)*ALGBRC(9)*(1.00000 - STATES(2)) -  ALGBRC(27)*STATES(2)
      ALGBRC(18) = TERNRY(STATES(1).LE.0.00000,  CONSTS(35)*EXP( CONSTS(34)*STATES(4)),  CONSTS(37)*EXP( CONSTS(36)*STATES(4)))
      ALGBRC(31) = ALGBRC(22)/ALGBRC(25)
      ALGBRC(29) = ( CONSTS(15)*1.00000*(1.00000+ 0.400000*CONSTS(15)+( 1.20000*STATES(1))/CONSTS(14)+ 0.600000*STATES(1)/CONSTS(14) ** 2.00000))/( CONSTS(14)* (CONSTS(15) - STATES(1)/CONSTS(14))*(1.00000+( 0.600000*STATES(1))/CONSTS(14)) ** 2.00000)
      ALGBRC(30) =  (1.00000/CONSTS(14))*EXP( - CONSTS(26)*STATES(1)/CONSTS(14) - CONSTS(51)/CONSTS(14) ** CONSTS(17))*(( 0.400000*CONSTS(15)+1.00000)/CONSTS(15)+ CONSTS(26)*CONSTS(17)*(1.00000+( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14)))*STATES(1)/CONSTS(14) - CONSTS(51)/CONSTS(14) ** CONSTS(17) - 1.00000)
      ALGBRC(33) = TERNRY(STATES(1).LE.- CONSTS(14), CONSTS(52), TERNRY(- CONSTS(14).LT.STATES(1).AND.STATES(1).LE.0.00000, ALGBRC(29), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(51), CONSTS(53), ALGBRC(30))
      ALGBRC(16) = TERNRY(STATES(1).LE.0.00000, CONSTS(34), CONSTS(36))
      ALGBRC(34) = TERNRY(ALGBRC(14).EQ.1.00000, - ( CONSTS(8)*ALGBRC(28)*ALGBRC(31)+ ALGBRC(16)*ALGBRC(18)*STATES(1) ** 2.00000+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(33)+ALGBRC(18)), - ( CONSTS(8)*ALGBRC(28)*ALGBRC(31)+ ALGBRC(16)*ALGBRC(18)*STATES(1) ** 2.00000+ ( CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6)))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(33)+ALGBRC(18)))
      ALGBRC(20) = TERNRY(STATES(3).LE.STATES(1), CONSTS(38), CONSTS(40))
      ALGBRC(2) =  CONSTS(3)*(EXP( CONSTS(2)*(STATES(5) - STATES(4))) - 1.00000)
      ALGBRC(3) =  CONSTS(5)*(EXP( CONSTS(4)*STATES(5)) - 1.00000)
      ALGBRC(32) =  CONSTS(8)*ALGBRC(31)*STATES(2)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END