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 41 entries in the algebraic variable array.
C There are a total of 16 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 STATES(1) is V in component membrane (millivolt).
C CONSTS(1) is R in component membrane (joule_per_kilomole_kelvin).
C CONSTS(2) is T in component membrane (kelvin).
C CONSTS(3) is F in component membrane (coulomb_per_mole).
C CONSTS(50) is RTONF in component membrane (millivolt).
C CONSTS(4) is C_m in component membrane (microF).
C ALGBRC(27) is i_b_K in component potassium_background_current (nanoA).
C ALGBRC(28) is i_K1 in component time_independent_potassium_current (nanoA).
C ALGBRC(19) is i_to in component transient_outward_current (nanoA).
C ALGBRC(23) is i_b_Na in component sodium_background_current (nanoA).
C ALGBRC(25) is i_b_Ca in component calcium_background_current (nanoA).
C ALGBRC(21) is i_NaK in component sodium_potassium_pump (nanoA).
C ALGBRC(26) is i_NaCa in component Na_Ca_exchanger (nanoA).
C ALGBRC(15) is i_Na in component fast_sodium_current (nanoA).
C ALGBRC(39) is i_Ca_L in component L_type_calcium_current (nanoA).
C ALGBRC(5) is i_Stim in component membrane (nanoA).
C CONSTS(5) is stim_start in component membrane (second).
C CONSTS(6) is stim_end in component membrane (second).
C CONSTS(7) is stim_period in component membrane (second).
C CONSTS(8) is stim_duration in component membrane (second).
C CONSTS(9) is stim_amplitude in component membrane (nanoA).
C CONSTS(10) is g_Na in component fast_sodium_current (microS).
C ALGBRC(11) is E_mh in component fast_sodium_current (millivolt).
C CONSTS(11) is Na_o in component extracellular_sodium_concentration (millimolar).
C STATES(2) is Na_i in component intracellular_sodium_concentration (millimolar).
C CONSTS(12) is K_c in component extracellular_potassium_concentration (millimolar).
C STATES(3) is K_i in component intracellular_potassium_concentration (millimolar).
C STATES(4) is m in component fast_sodium_current_m_gate (dimensionless).
C STATES(5) is h in component fast_sodium_current_h_gate (dimensionless).
C ALGBRC(7) is alpha_m in component fast_sodium_current_m_gate (per_second).
C ALGBRC(13) is beta_m in component fast_sodium_current_m_gate (per_second).
C CONSTS(13) is delta_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(1) is E0_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(2) is alpha_h in component fast_sodium_current_h_gate (per_second).
C ALGBRC(8) is beta_h in component fast_sodium_current_h_gate (per_second).
C CONSTS(14) is g_to in component transient_outward_current (microS).
C ALGBRC(17) is E_K in component transient_outward_current (millivolt).
C CONSTS(15) is g_to_s in component transient_outward_current (dimensionless).
C STATES(6) is r in component transient_outward_current_r_gate (dimensionless).
C STATES(7) is s in component transient_outward_current_s_gate (dimensionless).
C ALGBRC(3) is alpha_s in component transient_outward_current_s_gate (per_second).
C ALGBRC(9) is beta_s in component transient_outward_current_s_gate (per_second).
C CONSTS(16) is i_NaK_max in component sodium_potassium_pump (nanoA).
C CONSTS(17) is K_mK in component sodium_potassium_pump (millimolar).
C CONSTS(18) is K_mNa in component sodium_potassium_pump (millimolar).
C ALGBRC(22) is E_Na in component sodium_background_current (millivolt).
C CONSTS(19) is g_b_Na in component sodium_background_current (microS).
C ALGBRC(24) is E_Ca in component calcium_background_current (millivolt).
C CONSTS(20) is g_b_Ca in component calcium_background_current (microS).
C CONSTS(21) is Ca_o in component extracellular_calcium_concentration (millimolar).
C STATES(8) is Ca_i in component intracellular_calcium_concentration (millimolar).
C CONSTS(22) is k_NaCa in component Na_Ca_exchanger (nanoA).
C CONSTS(23) is n_NaCa in component Na_Ca_exchanger (dimensionless).
C CONSTS(24) is d_NaCa in component Na_Ca_exchanger (dimensionless).
C CONSTS(25) is gamma in component Na_Ca_exchanger (dimensionless).
C CONSTS(26) is g_b_K in component potassium_background_current (microS).
C CONSTS(27) is g_K1 in component time_independent_potassium_current (microS).
C CONSTS(28) is K_m_K1 in component time_independent_potassium_current (millimolar).
C ALGBRC(34) is i_Ca_L_Ca in component L_type_calcium_current (nanoA).
C ALGBRC(35) is i_Ca_L_K in component L_type_calcium_current (nanoA).
C ALGBRC(37) is i_Ca_L_Na in component L_type_calcium_current (nanoA).
C CONSTS(29) is P_Ca_L in component L_type_calcium_current (nanoA_per_millimolar).
C STATES(9) is d in component L_type_calcium_current_d_gate (dimensionless).
C STATES(10) is f_Ca in component L_type_calcium_current_f_Ca_gate (dimensionless).
C ALGBRC(33) is CaChon in component L_type_calcium_current_f_Ca_gate (dimensionless).
C ALGBRC(10) is alpha_d in component L_type_calcium_current_d_gate (per_second).
C ALGBRC(14) is beta_d in component L_type_calcium_current_d_gate (per_second).
C ALGBRC(4) is E0_d in component L_type_calcium_current_d_gate (millivolt).
C ALGBRC(30) is alpha_f_Ca in component L_type_calcium_current_f_Ca_gate (per_second).
C ALGBRC(31) is beta_f_Ca in component L_type_calcium_current_f_Ca_gate (per_second).
C ALGBRC(32) is CaChoff in component L_type_calcium_current_f_Ca_gate (dimensionless).
C ALGBRC(29) is E0_f in component L_type_calcium_current_f_Ca_gate (millivolt).
C ALGBRC(38) is i_up in component sarcoplasmic_reticulum_calcium_pump (millimolar_per_second).
C CONSTS(52) is K_1 in component sarcoplasmic_reticulum_calcium_pump (dimensionless).
C ALGBRC(36) is K_2 in component sarcoplasmic_reticulum_calcium_pump (millimolar).
C CONSTS(30) is K_cyca in component sarcoplasmic_reticulum_calcium_pump (millimolar).
C CONSTS(31) is K_xcs in component sarcoplasmic_reticulum_calcium_pump (dimensionless).
C CONSTS(32) is K_srca in component sarcoplasmic_reticulum_calcium_pump (millimolar).
C CONSTS(33) is alpha_up in component sarcoplasmic_reticulum_calcium_pump (millimolar_per_second).
C CONSTS(34) is beta_up in component sarcoplasmic_reticulum_calcium_pump (millimolar_per_second).
C STATES(11) is Ca_up in component intracellular_calcium_concentration (millimolar).
C ALGBRC(40) is i_rel in component calcium_release (millimolar_per_second).
C ALGBRC(12) is VoltDep in component calcium_release (dimensionless).
C ALGBRC(16) is RegBindSite in component calcium_release (dimensionless).
C ALGBRC(18) is ActRate in component calcium_release (per_second).
C ALGBRC(20) is InactRate in component calcium_release (per_second).
C CONSTS(35) is K_leak_rate in component calcium_release (per_second).
C CONSTS(36) is K_m_rel in component calcium_release (per_second).
C ALGBRC(6) is PrecFrac in component calcium_release (dimensionless).
C STATES(12) is ActFrac in component calcium_release (dimensionless).
C STATES(13) is ProdFrac in component calcium_release (dimensionless).
C CONSTS(37) is ProdFracRate in component calcium_release (per_second).
C STATES(14) is Ca_rel in component intracellular_calcium_concentration (millimolar).
C ALGBRC(41) is i_trans in component calcium_translocation (millimolar_per_second).
C CONSTS(38) is alpha_tr in component calcium_translocation (per_second).
C CONSTS(54) is V_i in component intracellular_calcium_concentration (micrometre3).
C STATES(15) is Ca_Calmod in component intracellular_calcium_concentration (millimolar).
C STATES(16) is Ca_Trop in component intracellular_calcium_concentration (millimolar).
C CONSTS(39) is Calmod in component intracellular_calcium_concentration (millimolar).
C CONSTS(40) is Trop in component intracellular_calcium_concentration (millimolar).
C CONSTS(41) is alpha_Calmod in component intracellular_calcium_concentration (per_millimolar_second).
C CONSTS(42) is beta_Calmod in component intracellular_calcium_concentration (per_second).
C CONSTS(43) is alpha_Trop in component intracellular_calcium_concentration (per_millimolar_second).
C CONSTS(44) is beta_Trop in component intracellular_calcium_concentration (per_second).
C CONSTS(45) is radius in component intracellular_calcium_concentration (micrometre).
C CONSTS(46) is length in component intracellular_calcium_concentration (micrometre).
C CONSTS(51) is V_Cell in component intracellular_calcium_concentration (micrometre3).
C CONSTS(53) is V_i_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(47) is V_rel_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(48) is V_e_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(49) is V_up_ratio in component intracellular_calcium_concentration (dimensionless).
C RATES(1) is d/dt V in component membrane (millivolt).
C RATES(4) is d/dt m in component fast_sodium_current_m_gate (dimensionless).
C RATES(5) is d/dt h in component fast_sodium_current_h_gate (dimensionless).
C RATES(6) is d/dt r in component transient_outward_current_r_gate (dimensionless).
C RATES(7) is d/dt s in component transient_outward_current_s_gate (dimensionless).
C RATES(9) is d/dt d in component L_type_calcium_current_d_gate (dimensionless).
C RATES(10) is d/dt f_Ca in component L_type_calcium_current_f_Ca_gate (dimensionless).
C RATES(12) is d/dt ActFrac in component calcium_release (dimensionless).
C RATES(13) is d/dt ProdFrac in component calcium_release (dimensionless).
C RATES(2) is d/dt Na_i in component intracellular_sodium_concentration (millimolar).
C RATES(3) is d/dt K_i in component intracellular_potassium_concentration (millimolar).
C RATES(8) is d/dt Ca_i in component intracellular_calcium_concentration (millimolar).
C RATES(15) is d/dt Ca_Calmod in component intracellular_calcium_concentration (millimolar).
C RATES(16) is d/dt Ca_Trop in component intracellular_calcium_concentration (millimolar).
C RATES(11) is d/dt Ca_up in component intracellular_calcium_concentration (millimolar).
C RATES(14) is d/dt Ca_rel in component intracellular_calcium_concentration (millimolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = -91.6
CONSTS(1) = 8314.472
CONSTS(2) = 310
CONSTS(3) = 96485.3415
CONSTS(4) = 4e-5
CONSTS(5) = 0.1
CONSTS(6) = 100000
CONSTS(7) = 1
CONSTS(8) = 0.002
CONSTS(9) = -1.3
CONSTS(10) = 0.5
CONSTS(11) = 140
STATES(2) = 6.48
CONSTS(12) = 4
STATES(3) = 140
STATES(4) = 0.076
STATES(5) = 0.015
CONSTS(13) = 1e-5
CONSTS(14) = 0.01
CONSTS(15) = 0
STATES(6) = 0
STATES(7) = 1
CONSTS(16) = 0.14
CONSTS(17) = 1
CONSTS(18) = 40
CONSTS(19) = 0.00012
CONSTS(20) = 5e-5
CONSTS(21) = 2
STATES(8) = 1e-5
CONSTS(22) = 0.0001
CONSTS(23) = 3
CONSTS(24) = 0.0001
CONSTS(25) = 0.5
CONSTS(26) = 0.0017
CONSTS(27) = 0.017
CONSTS(28) = 10
CONSTS(29) = 0.05
STATES(9) = 0.0011
STATES(10) = 0.785
CONSTS(30) = 0.0003
CONSTS(31) = 0.4
CONSTS(32) = 0.5
CONSTS(33) = 3
CONSTS(34) = 0.23
STATES(11) = 0.3
CONSTS(35) = 0
CONSTS(36) = 250
STATES(12) = 0
STATES(13) = 0
CONSTS(37) = 1
STATES(14) = 0.3
CONSTS(38) = 50
STATES(15) = 0.0005
STATES(16) = 0.0015
CONSTS(39) = 0.02
CONSTS(40) = 0.15
CONSTS(41) = 100000
CONSTS(42) = 50
CONSTS(43) = 100000
CONSTS(44) = 200
CONSTS(45) = 0.01
CONSTS(46) = 0.08
CONSTS(47) = 0.1
CONSTS(48) = 0.4
CONSTS(49) = 0.01
CONSTS(50) = ( CONSTS(1)*CONSTS(2))/CONSTS(3)
CONSTS(51) = 3.14159*CONSTS(45) ** 2.00000*CONSTS(46)
CONSTS(52) = ( CONSTS(30)*CONSTS(31))/CONSTS(32)
CONSTS(53) = ((1.00000 - CONSTS(48)) - CONSTS(49)) - CONSTS(47)
CONSTS(54) = CONSTS(51)*CONSTS(53)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(6) = 333.000*(1.00000/(1.00000+EXP(- (STATES(1)+4.00000)/5.00000)) - STATES(6))
ALGBRC(2) = 20.0000*EXP( - 0.125000*(STATES(1)+75.0000))
ALGBRC(8) = 2000.00/(1.00000+ 320.000*EXP( - 0.100000*(STATES(1)+75.0000)))
RATES(5) = ALGBRC(2)*(1.00000 - STATES(5)) - ALGBRC(8)*STATES(5)
ALGBRC(3) = 0.0330000*EXP(- STATES(1)/17.0000)
ALGBRC(9) = 33.0000/(1.00000+EXP( - 0.125000*(STATES(1)+10.0000)))
RATES(7) = ALGBRC(3)*(1.00000 - STATES(7)) - ALGBRC(9)*STATES(7)
ALGBRC(1) = STATES(1)+41.0000
ALGBRC(7) = TERNRY(ABS(ALGBRC(1)).LT.CONSTS(13), 2000.00, ( 200.000*ALGBRC(1))/(1.00000 - EXP( - 0.100000*ALGBRC(1))))
ALGBRC(13) = 8000.00*EXP( - 0.0560000*(STATES(1)+66.0000))
RATES(4) = ALGBRC(7)*(1.00000 - STATES(4)) - ALGBRC(13)*STATES(4)
ALGBRC(4) = STATES(1)+19.0000
ALGBRC(10) = TERNRY(ABS(ALGBRC(4)).LT.0.000100000, 120.000, ( 30.0000*ALGBRC(4))/(1.00000 - EXP(- ALGBRC(4)/4.00000)))
ALGBRC(14) = TERNRY(ABS(ALGBRC(4)).LT.0.000100000, 120.000, ( 12.0000*ALGBRC(4))/(EXP(ALGBRC(4)/10.0000) - 1.00000))
RATES(9) = ALGBRC(10)*(1.00000 - STATES(9)) - ALGBRC(14)*STATES(9)
ALGBRC(12) = EXP( 0.0800000*(STATES(1) - 40.0000))
ALGBRC(16) = STATES(8)/(STATES(8)+0.000500000) ** 2.00000
ALGBRC(18) = 600.000*ALGBRC(12)+ 500.000*ALGBRC(16)
ALGBRC(20) = 60.0000+ 500.000*ALGBRC(16)
ALGBRC(6) = (1.00000 - STATES(12)) - STATES(13)
RATES(12) = ALGBRC(6)*ALGBRC(18) - STATES(12)*ALGBRC(20)
RATES(13) = STATES(12)*ALGBRC(20) - CONSTS(37)*STATES(13)
ALGBRC(29) = STATES(1)+34.0000
ALGBRC(30) = TERNRY(ABS(ALGBRC(29)).LT.0.000100000, 25.0000, ( 6.25000*ALGBRC(29))/(EXP(ALGBRC(29)/4.00000) - 1.00000))
ALGBRC(31) = 12.0000/(1.00000+EXP(- ALGBRC(29)/4.00000))
ALGBRC(32) = STATES(8)/(0.00100000+STATES(8))
RATES(10) = ( 120.000*(1.00000 - STATES(10))*ALGBRC(32)+ (1.00000 - STATES(10))*(1.00000 - ALGBRC(32)))*ALGBRC(31) - ALGBRC(30)*STATES(10)
ALGBRC(17) = CONSTS(50)*log(CONSTS(12)/STATES(3))
ALGBRC(27) = CONSTS(26)*(STATES(1) - ALGBRC(17))
ALGBRC(28) = ( (( CONSTS(27)*CONSTS(12))/(CONSTS(12)+CONSTS(28)))*(STATES(1) - ALGBRC(17)))/(1.00000+EXP(( ((STATES(1) - ALGBRC(17)) - 10.0000)*2.00000)/CONSTS(50)))
ALGBRC(19) = CONSTS(14)*(CONSTS(15)+ STATES(7)*(1.00000 - CONSTS(15)))*STATES(6)*(STATES(1) - ALGBRC(17))
ALGBRC(21) = ( (( CONSTS(16)*CONSTS(12))/(CONSTS(17)+CONSTS(12)))*STATES(2))/(CONSTS(18)+STATES(2))
ALGBRC(33) = (1.00000 - STATES(10))*(1.00000 - ALGBRC(32))
ALGBRC(35) = ((( 0.00200000*CONSTS(29)*STATES(9)*ALGBRC(33)*(STATES(1) - 50.0000))/CONSTS(50))/(1.00000 - EXP(- (STATES(1) - 50.0000)/CONSTS(50))))*( STATES(3)*EXP(50.0000/CONSTS(50)) - CONSTS(12)*EXP(- (STATES(1) - 50.0000)/CONSTS(50)))
RATES(3) = (- 1.00000/( 1.00000*CONSTS(54)*CONSTS(3)))*((ALGBRC(28)+ALGBRC(35)+ALGBRC(19)+ALGBRC(27)) - 2.00000*ALGBRC(21))
ALGBRC(22) = CONSTS(50)*log(CONSTS(11)/STATES(2))
ALGBRC(23) = CONSTS(19)*(STATES(1) - ALGBRC(22))
ALGBRC(26) = ( CONSTS(22)*( EXP(( CONSTS(25)*(CONSTS(23) - 2.00000)*STATES(1))/CONSTS(50))*STATES(2) ** CONSTS(23)*CONSTS(21) - EXP(( (CONSTS(25) - 1.00000)*(CONSTS(23) - 2.00000)*STATES(1))/CONSTS(50))*CONSTS(11) ** CONSTS(23)*STATES(8)))/( (1.00000+ CONSTS(24)*( STATES(8)*CONSTS(11) ** CONSTS(23)+ CONSTS(21)*STATES(2) ** CONSTS(23)))*(1.00000+STATES(8)/0.00690000))
ALGBRC(11) = CONSTS(50)*log((CONSTS(11)+ 0.120000*CONSTS(12))/(STATES(2)+ 0.120000*STATES(3)))
ALGBRC(15) = CONSTS(10)*STATES(4) ** 3.00000*STATES(5)*(STATES(1) - ALGBRC(11))
ALGBRC(37) = ((( 0.0100000*CONSTS(29)*STATES(9)*ALGBRC(33)*(STATES(1) - 50.0000))/CONSTS(50))/(1.00000 - EXP(- (STATES(1) - 50.0000)/CONSTS(50))))*( STATES(2)*EXP(50.0000/CONSTS(50)) - CONSTS(11)*EXP(- (STATES(1) - 50.0000)/CONSTS(50)))
RATES(2) = (- 1.00000/( 1.00000*CONSTS(54)*CONSTS(3)))*(ALGBRC(15)+ALGBRC(23)+ 3.00000*ALGBRC(21)+ 3.00000*ALGBRC(26)+ALGBRC(37))
ALGBRC(24) = 0.500000*CONSTS(50)*log(CONSTS(21)/STATES(8))
ALGBRC(25) = CONSTS(20)*(STATES(1) - ALGBRC(24))
ALGBRC(34) = ((( 4.00000*CONSTS(29)*STATES(9)*ALGBRC(33)*(STATES(1) - 50.0000))/CONSTS(50))/(1.00000 - EXP(( - (STATES(1) - 50.0000)*2.00000)/CONSTS(50))))*( STATES(8)*EXP(100.000/CONSTS(50)) - CONSTS(21)*EXP(( - (STATES(1) - 50.0000)*2.00000)/CONSTS(50)))
ALGBRC(39) = ALGBRC(34)+ALGBRC(35)+ALGBRC(37)
ALGBRC(5) = TERNRY(VOI.GE.CONSTS(5).AND.VOI.LE.CONSTS(6).AND.(VOI - CONSTS(5)) - INT((VOI - CONSTS(5))/CONSTS(7))*CONSTS(7).LE.CONSTS(8), CONSTS(9), 0.00000)
RATES(1) = (- 1.00000/CONSTS(4))*(ALGBRC(5)+ALGBRC(27)+ALGBRC(28)+ALGBRC(19)+ALGBRC(23)+ALGBRC(25)+ALGBRC(21)+ALGBRC(26)+ALGBRC(15)+ALGBRC(39))
RATES(15) = CONSTS(41)*STATES(8)*(CONSTS(39) - STATES(15)) - CONSTS(42)*STATES(15)
RATES(16) = CONSTS(43)*STATES(8)*(CONSTS(40) - STATES(16)) - CONSTS(44)*STATES(16)
ALGBRC(36) = STATES(8)+ STATES(11)*CONSTS(52)+ CONSTS(30)*CONSTS(31)+CONSTS(30)
ALGBRC(38) = (STATES(8)/ALGBRC(36))*CONSTS(33) - (( STATES(11)*CONSTS(52))/ALGBRC(36))*CONSTS(34)
ALGBRC(41) = (STATES(11) - STATES(14))*CONSTS(38)
RATES(11) = (CONSTS(53)/CONSTS(49))*ALGBRC(38) - ALGBRC(41)
ALGBRC(40) = ( STATES(12)/(STATES(12)+0.250000) ** 2.00000*CONSTS(36)+CONSTS(35))*STATES(14)
RATES(14) = (CONSTS(49)/CONSTS(47))*ALGBRC(41) - ALGBRC(40)
RATES(8) = ((( (- 1.00000/( 2.00000*1.00000*CONSTS(54)*CONSTS(3)))*((ALGBRC(34)+ALGBRC(25)) - 2.00000*ALGBRC(26))+( ALGBRC(40)*CONSTS(47))/CONSTS(53)) - RATES(15)) - RATES(16)) - ALGBRC(38)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(2) = 20.0000*EXP( - 0.125000*(STATES(1)+75.0000))
ALGBRC(8) = 2000.00/(1.00000+ 320.000*EXP( - 0.100000*(STATES(1)+75.0000)))
ALGBRC(3) = 0.0330000*EXP(- STATES(1)/17.0000)
ALGBRC(9) = 33.0000/(1.00000+EXP( - 0.125000*(STATES(1)+10.0000)))
ALGBRC(1) = STATES(1)+41.0000
ALGBRC(7) = TERNRY(ABS(ALGBRC(1)).LT.CONSTS(13), 2000.00, ( 200.000*ALGBRC(1))/(1.00000 - EXP( - 0.100000*ALGBRC(1))))
ALGBRC(13) = 8000.00*EXP( - 0.0560000*(STATES(1)+66.0000))
ALGBRC(4) = STATES(1)+19.0000
ALGBRC(10) = TERNRY(ABS(ALGBRC(4)).LT.0.000100000, 120.000, ( 30.0000*ALGBRC(4))/(1.00000 - EXP(- ALGBRC(4)/4.00000)))
ALGBRC(14) = TERNRY(ABS(ALGBRC(4)).LT.0.000100000, 120.000, ( 12.0000*ALGBRC(4))/(EXP(ALGBRC(4)/10.0000) - 1.00000))
ALGBRC(12) = EXP( 0.0800000*(STATES(1) - 40.0000))
ALGBRC(16) = STATES(8)/(STATES(8)+0.000500000) ** 2.00000
ALGBRC(18) = 600.000*ALGBRC(12)+ 500.000*ALGBRC(16)
ALGBRC(20) = 60.0000+ 500.000*ALGBRC(16)
ALGBRC(6) = (1.00000 - STATES(12)) - STATES(13)
ALGBRC(29) = STATES(1)+34.0000
ALGBRC(30) = TERNRY(ABS(ALGBRC(29)).LT.0.000100000, 25.0000, ( 6.25000*ALGBRC(29))/(EXP(ALGBRC(29)/4.00000) - 1.00000))
ALGBRC(31) = 12.0000/(1.00000+EXP(- ALGBRC(29)/4.00000))
ALGBRC(32) = STATES(8)/(0.00100000+STATES(8))
ALGBRC(17) = CONSTS(50)*log(CONSTS(12)/STATES(3))
ALGBRC(27) = CONSTS(26)*(STATES(1) - ALGBRC(17))
ALGBRC(28) = ( (( CONSTS(27)*CONSTS(12))/(CONSTS(12)+CONSTS(28)))*(STATES(1) - ALGBRC(17)))/(1.00000+EXP(( ((STATES(1) - ALGBRC(17)) - 10.0000)*2.00000)/CONSTS(50)))
ALGBRC(19) = CONSTS(14)*(CONSTS(15)+ STATES(7)*(1.00000 - CONSTS(15)))*STATES(6)*(STATES(1) - ALGBRC(17))
ALGBRC(21) = ( (( CONSTS(16)*CONSTS(12))/(CONSTS(17)+CONSTS(12)))*STATES(2))/(CONSTS(18)+STATES(2))
ALGBRC(33) = (1.00000 - STATES(10))*(1.00000 - ALGBRC(32))
ALGBRC(35) = ((( 0.00200000*CONSTS(29)*STATES(9)*ALGBRC(33)*(STATES(1) - 50.0000))/CONSTS(50))/(1.00000 - EXP(- (STATES(1) - 50.0000)/CONSTS(50))))*( STATES(3)*EXP(50.0000/CONSTS(50)) - CONSTS(12)*EXP(- (STATES(1) - 50.0000)/CONSTS(50)))
ALGBRC(22) = CONSTS(50)*log(CONSTS(11)/STATES(2))
ALGBRC(23) = CONSTS(19)*(STATES(1) - ALGBRC(22))
ALGBRC(26) = ( CONSTS(22)*( EXP(( CONSTS(25)*(CONSTS(23) - 2.00000)*STATES(1))/CONSTS(50))*STATES(2) ** CONSTS(23)*CONSTS(21) - EXP(( (CONSTS(25) - 1.00000)*(CONSTS(23) - 2.00000)*STATES(1))/CONSTS(50))*CONSTS(11) ** CONSTS(23)*STATES(8)))/( (1.00000+ CONSTS(24)*( STATES(8)*CONSTS(11) ** CONSTS(23)+ CONSTS(21)*STATES(2) ** CONSTS(23)))*(1.00000+STATES(8)/0.00690000))
ALGBRC(11) = CONSTS(50)*log((CONSTS(11)+ 0.120000*CONSTS(12))/(STATES(2)+ 0.120000*STATES(3)))
ALGBRC(15) = CONSTS(10)*STATES(4) ** 3.00000*STATES(5)*(STATES(1) - ALGBRC(11))
ALGBRC(37) = ((( 0.0100000*CONSTS(29)*STATES(9)*ALGBRC(33)*(STATES(1) - 50.0000))/CONSTS(50))/(1.00000 - EXP(- (STATES(1) - 50.0000)/CONSTS(50))))*( STATES(2)*EXP(50.0000/CONSTS(50)) - CONSTS(11)*EXP(- (STATES(1) - 50.0000)/CONSTS(50)))
ALGBRC(24) = 0.500000*CONSTS(50)*log(CONSTS(21)/STATES(8))
ALGBRC(25) = CONSTS(20)*(STATES(1) - ALGBRC(24))
ALGBRC(34) = ((( 4.00000*CONSTS(29)*STATES(9)*ALGBRC(33)*(STATES(1) - 50.0000))/CONSTS(50))/(1.00000 - EXP(( - (STATES(1) - 50.0000)*2.00000)/CONSTS(50))))*( STATES(8)*EXP(100.000/CONSTS(50)) - CONSTS(21)*EXP(( - (STATES(1) - 50.0000)*2.00000)/CONSTS(50)))
ALGBRC(39) = ALGBRC(34)+ALGBRC(35)+ALGBRC(37)
ALGBRC(5) = TERNRY(VOI.GE.CONSTS(5).AND.VOI.LE.CONSTS(6).AND.(VOI - CONSTS(5)) - INT((VOI - CONSTS(5))/CONSTS(7))*CONSTS(7).LE.CONSTS(8), CONSTS(9), 0.00000)
ALGBRC(36) = STATES(8)+ STATES(11)*CONSTS(52)+ CONSTS(30)*CONSTS(31)+CONSTS(30)
ALGBRC(38) = (STATES(8)/ALGBRC(36))*CONSTS(33) - (( STATES(11)*CONSTS(52))/ALGBRC(36))*CONSTS(34)
ALGBRC(41) = (STATES(11) - STATES(14))*CONSTS(38)
ALGBRC(40) = ( STATES(12)/(STATES(12)+0.250000) ** 2.00000*CONSTS(36)+CONSTS(35))*STATES(14)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END