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 75 entries in the algebraic variable array.
C There are a total of 27 entries in each of the rate and state variable arrays.
C There are a total of 66 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is V in component membrane (millivolt).
C CONSTS(1) is R in component membrane (attojoule_per_millimole_kelvin).
C CONSTS(2) is T in component membrane (kelvin).
C CONSTS(3) is F in component membrane (femtocoulomb_per_millimole).
C CONSTS(4) is Cm in component membrane (picoF).
C ALGBRC(40) is i_CaT in component T_type_calcium_channel_current (picoA).
C ALGBRC(35) is i_CaL in component L_type_calcium_channel_current (picoA).
C ALGBRC(41) is i_Kr in component rapidly_activating_delayed_rectifier_potassium_current (picoA).
C ALGBRC(43) is i_Ks in component slowly_activating_delayed_rectifier_potassium_current (picoA).
C ALGBRC(44) is i_to in component AP_sensitive_currents (picoA).
C ALGBRC(45) is i_sus in component AP_sensitive_currents (picoA).
C ALGBRC(48) is i_h in component hyperpolarisation_activated_current (picoA).
C ALGBRC(49) is i_st in component sustained_inward_current (picoA).
C ALGBRC(50) is i_b_Na in component sodium_dependent_background_current (picoA).
C ALGBRC(66) is i_NaCa in component sodium_calcium_exchange_current (picoA).
C ALGBRC(52) is i_NaK in component sodium_potassium_pump_current (picoA).
C ALGBRC(51) is i_K_ACh in component background_muscarinic_potassium_channel_current (picoA).
C ALGBRC(1) is E_Na in component reversal_potentials (millivolt).
C ALGBRC(30) is E_K in component reversal_potentials (millivolt).
C STATES(2) is Nai in component intracellular_ion_concentrations (millimolar).
C CONSTS(5) is Nao in component intracellular_ion_concentrations (millimolar).
C STATES(3) is Ki in component intracellular_ion_concentrations (millimolar).
C CONSTS(6) is Ko in component intracellular_ion_concentrations (millimolar).
C CONSTS(7) is g_CaL in component L_type_calcium_channel_current (nanoS).
C CONSTS(8) is E_CaL in component L_type_calcium_channel_current (millivolt).
C STATES(4) is Ca_sub in component intracellular_ion_concentrations (millimolar).
C STATES(5) is d in component L_type_calcium_channel_current_d_gate (dimensionless).
C STATES(6) is f in component L_type_calcium_channel_current_f_gate (dimensionless).
C STATES(7) is fCa in component L_type_calcium_channel_current_fCa_gate (dimensionless).
C ALGBRC(2) is d_infinity in component L_type_calcium_channel_current_d_gate (dimensionless).
C ALGBRC(36) is tau_d in component L_type_calcium_channel_current_d_gate (millisecond).
C ALGBRC(16) is alpha_d in component L_type_calcium_channel_current_d_gate (per_millisecond).
C ALGBRC(31) is beta_d in component L_type_calcium_channel_current_d_gate (per_millisecond).
C ALGBRC(3) is f_infinity in component L_type_calcium_channel_current_f_gate (dimensionless).
C ALGBRC(17) is tau_f in component L_type_calcium_channel_current_f_gate (millisecond).
C CONSTS(63) is alpha_fCa in component L_type_calcium_channel_current_fCa_gate (per_millisecond).
C CONSTS(9) is beta_fCa in component L_type_calcium_channel_current_fCa_gate (per_millimolar_millisecond).
C ALGBRC(4) is fCa_infinity in component L_type_calcium_channel_current_fCa_gate (dimensionless).
C ALGBRC(18) is tau_fCa in component L_type_calcium_channel_current_fCa_gate (millisecond).
C CONSTS(10) is Km_fCa in component L_type_calcium_channel_current_fCa_gate (millimolar).
C CONSTS(11) is g_CaT in component T_type_calcium_channel_current (nanoS).
C CONSTS(12) is E_CaT in component T_type_calcium_channel_current (millivolt).
C STATES(8) is d in component T_type_calcium_channel_current_d_gate (dimensionless).
C STATES(9) is f in component T_type_calcium_channel_current_f_gate (dimensionless).
C ALGBRC(5) is d_infinity in component T_type_calcium_channel_current_d_gate (dimensionless).
C ALGBRC(19) is tau_d in component T_type_calcium_channel_current_d_gate (millisecond).
C ALGBRC(6) is f_infinity in component T_type_calcium_channel_current_f_gate (dimensionless).
C ALGBRC(20) is tau_f in component T_type_calcium_channel_current_f_gate (millisecond).
C CONSTS(64) is g_Kr in component rapidly_activating_delayed_rectifier_potassium_current (nanoS).
C STATES(10) is paS in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
C STATES(11) is paF in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
C STATES(12) is piy in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (dimensionless).
C ALGBRC(7) is pa_infinity in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
C ALGBRC(21) is tau_paS in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (millisecond).
C ALGBRC(22) is tau_paF in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (millisecond).
C ALGBRC(8) is pi_infinity in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (dimensionless).
C ALGBRC(23) is tau_pi in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (millisecond).
C CONSTS(13) is g_Ks in component slowly_activating_delayed_rectifier_potassium_current (nanoS).
C ALGBRC(42) is E_Ks in component slowly_activating_delayed_rectifier_potassium_current (millivolt).
C STATES(13) is n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (dimensionless).
C ALGBRC(32) is n_infinity in component slowly_activating_delayed_rectifier_potassium_current_n_gate (dimensionless).
C ALGBRC(37) is tau_n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (millisecond).
C ALGBRC(9) is alpha_n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (per_millisecond).
C ALGBRC(24) is beta_n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (per_millisecond).
C CONSTS(14) is g_to in component AP_sensitive_currents (nanoS).
C CONSTS(15) is g_sus in component AP_sensitive_currents (nanoS).
C STATES(14) is q in component AP_sensitive_currents_q_gate (dimensionless).
C STATES(15) is r in component AP_sensitive_currents_r_gate (dimensionless).
C ALGBRC(10) is q_infinity in component AP_sensitive_currents_q_gate (dimensionless).
C ALGBRC(25) is tau_q in component AP_sensitive_currents_q_gate (millisecond).
C ALGBRC(11) is r_infinity in component AP_sensitive_currents_r_gate (dimensionless).
C ALGBRC(26) is tau_r in component AP_sensitive_currents_r_gate (millisecond).
C ALGBRC(46) is i_h_Na in component hyperpolarisation_activated_current (picoA).
C ALGBRC(47) is i_h_K in component hyperpolarisation_activated_current (picoA).
C CONSTS(16) is g_h_Na in component hyperpolarisation_activated_current (nanoS).
C CONSTS(17) is g_h_K in component hyperpolarisation_activated_current (nanoS).
C STATES(16) is y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C ALGBRC(12) is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless).
C ALGBRC(27) is tau_y in component hyperpolarisation_activated_current_y_gate (millisecond).
C CONSTS(18) is g_st in component sustained_inward_current (nanoS).
C CONSTS(19) is E_st in component sustained_inward_current (millivolt).
C STATES(17) is qa in component sustained_inward_current_qa_gate (dimensionless).
C STATES(18) is qi in component sustained_inward_current_qi_gate (dimensionless).
C ALGBRC(13) is qa_infinity in component sustained_inward_current_qa_gate (dimensionless).
C ALGBRC(38) is tau_qa in component sustained_inward_current_qa_gate (millisecond).
C ALGBRC(28) is alpha_qa in component sustained_inward_current_qa_gate (per_millisecond).
C ALGBRC(33) is beta_qa in component sustained_inward_current_qa_gate (per_millisecond).
C ALGBRC(34) is qi_infinity in component sustained_inward_current_qi_gate (dimensionless).
C ALGBRC(39) is tau_qi in component sustained_inward_current_qi_gate (millisecond).
C ALGBRC(14) is alpha_qi in component sustained_inward_current_qi_gate (per_millisecond).
C ALGBRC(29) is beta_qi in component sustained_inward_current_qi_gate (per_millisecond).
C CONSTS(20) is g_b_Na in component sodium_dependent_background_current (nanoS).
C CONSTS(65) is g_K_ACh in component background_muscarinic_potassium_channel_current (picoA).
C CONSTS(21) is Km_Kp in component sodium_potassium_pump_current (millimolar).
C CONSTS(22) is Km_Nap in component sodium_potassium_pump_current (millimolar).
C CONSTS(23) is i_NaK_max in component sodium_potassium_pump_current (picoA).
C CONSTS(24) is kNaCa in component sodium_calcium_exchange_current (picoA).
C ALGBRC(63) is x1 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(59) is x2 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(64) is x3 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(65) is x4 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(54) is k41 in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(66) is k34 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(62) is k23 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(61) is k21 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(58) is k32 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(53) is k43 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(56) is k12 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(57) is k14 in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(25) is Qci in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(26) is Qn in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(27) is Qco in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(28) is K3ni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(29) is Kci in component sodium_calcium_exchange_current (millimolar).
C CONSTS(30) is K1ni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(31) is K2ni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(32) is Kcni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(33) is K3no in component sodium_calcium_exchange_current (millimolar).
C CONSTS(34) is K1no in component sodium_calcium_exchange_current (millimolar).
C CONSTS(35) is K2no in component sodium_calcium_exchange_current (millimolar).
C CONSTS(36) is Kco in component sodium_calcium_exchange_current (millimolar).
C ALGBRC(60) is do in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(55) is di in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(37) is Cao in component intracellular_ion_concentrations (millimolar).
C ALGBRC(67) is j_Ca_dif in component intracellular_calcium_dynamics (millimolar_per_millisecond).
C ALGBRC(68) is j_rel in component intracellular_calcium_dynamics (millimolar_per_millisecond).
C ALGBRC(69) is j_up in component intracellular_calcium_dynamics (millimolar_per_millisecond).
C ALGBRC(71) is j_tr in component intracellular_calcium_dynamics (millimolar_per_millisecond).
C CONSTS(38) is tau_dif_Ca in component intracellular_calcium_dynamics (millisecond).
C CONSTS(39) is tau_tr in component intracellular_calcium_dynamics (millisecond).
C CONSTS(40) is K_rel in component intracellular_calcium_dynamics (millimolar).
C CONSTS(41) is P_up in component intracellular_calcium_dynamics (millimolar_per_millisecond).
C CONSTS(42) is P_rel in component intracellular_calcium_dynamics (per_millisecond).
C CONSTS(43) is K_up in component intracellular_calcium_dynamics (millimolar).
C STATES(19) is Ca_up in component intracellular_ion_concentrations (millimolar).
C STATES(20) is Cai in component intracellular_ion_concentrations (millimolar).
C STATES(21) is Ca_rel in component intracellular_ion_concentrations (millimolar).
C CONSTS(44) is Mgi in component intracellular_ion_concentrations (millimolar).
C CONSTS(45) is V_i in component intracellular_ion_concentrations (litre).
C CONSTS(46) is V_rel in component intracellular_ion_concentrations (litre).
C CONSTS(47) is V_up in component intracellular_ion_concentrations (litre).
C CONSTS(48) is V_sub in component intracellular_ion_concentrations (litre).
C CONSTS(49) is TMC_tot in component calcium_buffering (dimensionless).
C CONSTS(50) is CM_tot in component calcium_buffering (dimensionless).
C CONSTS(51) is TC_tot in component calcium_buffering (dimensionless).
C CONSTS(52) is CQ_tot in component calcium_buffering (dimensionless).
C ALGBRC(73) is delta_fTMC in component calcium_buffering (millimolar_per_millisecond).
C ALGBRC(75) is delta_fCMi in component calcium_buffering (millimolar_per_millisecond).
C ALGBRC(70) is delta_fCMs in component calcium_buffering (millimolar_per_millisecond).
C ALGBRC(72) is delta_fTC in component calcium_buffering (millimolar_per_millisecond).
C ALGBRC(74) is delta_fCQ in component calcium_buffering (millimolar_per_millisecond).
C ALGBRC(15) is delta_fTMM in component calcium_buffering (millimolar_per_millisecond).
C STATES(22) is fTMM in component calcium_buffering (millimolar).
C STATES(23) is fCMi in component calcium_buffering (millimolar).
C STATES(24) is fCMs in component calcium_buffering (millimolar).
C STATES(25) is fTC in component calcium_buffering (millimolar).
C STATES(26) is fTMC in component calcium_buffering (millimolar).
C STATES(27) is fCQ in component calcium_buffering (millimolar).
C CONSTS(53) is kf_TC in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(54) is kf_TMM in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(55) is kf_TMC in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(56) is kf_CM in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(57) is kf_CQ in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(58) is kb_TC in component calcium_buffering (per_millisecond).
C CONSTS(59) is kb_TMC in component calcium_buffering (per_millisecond).
C CONSTS(60) is kb_TMM in component calcium_buffering (per_millisecond).
C CONSTS(61) is kb_CM in component calcium_buffering (per_millisecond).
C CONSTS(62) is kb_CQ in component calcium_buffering (per_millisecond).
C RATES(1) is d/dt V in component membrane (millivolt).
C RATES(5) is d/dt d in component L_type_calcium_channel_current_d_gate (dimensionless).
C RATES(6) is d/dt f in component L_type_calcium_channel_current_f_gate (dimensionless).
C RATES(7) is d/dt fCa in component L_type_calcium_channel_current_fCa_gate (dimensionless).
C RATES(8) is d/dt d in component T_type_calcium_channel_current_d_gate (dimensionless).
C RATES(9) is d/dt f in component T_type_calcium_channel_current_f_gate (dimensionless).
C RATES(10) is d/dt paS in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
C RATES(11) is d/dt paF in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
C RATES(12) is d/dt piy in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (dimensionless).
C RATES(13) is d/dt n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (dimensionless).
C RATES(14) is d/dt q in component AP_sensitive_currents_q_gate (dimensionless).
C RATES(15) is d/dt r in component AP_sensitive_currents_r_gate (dimensionless).
C RATES(16) is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C RATES(17) is d/dt qa in component sustained_inward_current_qa_gate (dimensionless).
C RATES(18) is d/dt qi in component sustained_inward_current_qi_gate (dimensionless).
C RATES(2) is d/dt Nai in component intracellular_ion_concentrations (millimolar).
C RATES(3) is d/dt Ki in component intracellular_ion_concentrations (millimolar).
C RATES(20) is d/dt Cai in component intracellular_ion_concentrations (millimolar).
C RATES(4) is d/dt Ca_sub in component intracellular_ion_concentrations (millimolar).
C RATES(19) is d/dt Ca_up in component intracellular_ion_concentrations (millimolar).
C RATES(21) is d/dt Ca_rel in component intracellular_ion_concentrations (millimolar).
C RATES(25) is d/dt fTC in component calcium_buffering (millimolar).
C RATES(26) is d/dt fTMC in component calcium_buffering (millimolar).
C RATES(22) is d/dt fTMM in component calcium_buffering (millimolar).
C RATES(23) is d/dt fCMi in component calcium_buffering (millimolar).
C RATES(24) is d/dt fCMs in component calcium_buffering (millimolar).
C RATES(27) is d/dt fCQ in component calcium_buffering (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -58.600291137693
      CONSTS(1) = 8314400000000000
      CONSTS(2) = 310.15
      CONSTS(3) = 96485000000000000
      CONSTS(4) = 32
      STATES(2) = 9.438646305915
      CONSTS(5) = 140
      STATES(3) = 139.984146485614
      CONSTS(6) = 5.4
      CONSTS(7) = 0.58
      CONSTS(8) = 45
      STATES(4) = 0.00019074741
      STATES(5) = 0.000602055134
      STATES(6) = 0.626999773853
      STATES(7) = 0.589580408056
      CONSTS(9) = 60
      CONSTS(10) = 0.00035
      CONSTS(11) = 0.458
      CONSTS(12) = 45
      STATES(8) = 0.004571884917
      STATES(9) = 0.249637570396
      STATES(10) = 0.629323128348
      STATES(11) = 0.3493633709533
      STATES(12) = 0.852396631172
      CONSTS(13) = 0.0259
      STATES(13) = 0.054409723782
      CONSTS(14) = 0.18
      CONSTS(15) = 0.02
      STATES(14) = 0.531446952485
      STATES(15) = 0.005550489445
      CONSTS(16) = 0.1437375
      CONSTS(17) = 0.2312625
      STATES(16) = 0.067156687129
      CONSTS(18) = 0.015
      CONSTS(19) = 37.4
      STATES(17) = 0.426018100136
      STATES(18) = 0.333330378068
      CONSTS(20) = 0.0054
      CONSTS(21) = 1.4
      CONSTS(22) = 14
      CONSTS(23) = 3.6
      CONSTS(24) = 125
      CONSTS(25) = 0.1369
      CONSTS(26) = 0.4315
      CONSTS(27) = 0
      CONSTS(28) = 26.44
      CONSTS(29) = 0.0207
      CONSTS(30) = 395.3
      CONSTS(31) = 2.289
      CONSTS(32) = 26.44
      CONSTS(33) = 4.663
      CONSTS(34) = 1628
      CONSTS(35) = 561.4
      CONSTS(36) = 3.663
      CONSTS(37) = 2
      CONSTS(38) = 0.04
      CONSTS(39) = 60
      CONSTS(40) = 0.0012
      CONSTS(41) = 0.005
      CONSTS(42) = 0.5
      CONSTS(43) = 0.0006
      STATES(19) = 1.462338380106
      STATES(20) = 0.000312494921
      STATES(21) = 0.296742023718
      CONSTS(44) = 2.5
      CONSTS(45) = 0.0000000000015835
      CONSTS(46) = 0.0000000000000042223
      CONSTS(47) = 0.000000000000040816
      CONSTS(48) = 0.000000000000035098
      CONSTS(49) = 0.062
      CONSTS(50) = 0.045
      CONSTS(51) = 0.031
      CONSTS(52) = 10
      STATES(22) = 0.350600895635
      STATES(23) = 0.116947220413
      STATES(24) = 0.074631965653
      STATES(25) = 0.059206293446
      STATES(26) = 0.602955114871
      STATES(27) = 0.260317260703
      CONSTS(53) = 88.8
      CONSTS(54) = 2.277
      CONSTS(55) = 227.7
      CONSTS(56) = 227.7
      CONSTS(57) = 0.534
      CONSTS(58) = 0.446
      CONSTS(59) = 0.00751
      CONSTS(60) = 0.751
      CONSTS(61) = 0.542
      CONSTS(62) = 0.445
      CONSTS(63) =  CONSTS(10)*CONSTS(9)
      CONSTS(64) =  0.0250000*CONSTS(6)/1.00000 ** 0.590000
      CONSTS(65) =  0.00110000*CONSTS(6)/1.00000 ** 0.410000
      CONSTS(66) = CONSTS(5)/(CONSTS(33)+CONSTS(5))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(15) =  CONSTS(54)*CONSTS(44)*(1.00000 - (STATES(26)+STATES(22))) -  CONSTS(60)*STATES(22)
      RATES(22) = ALGBRC(15)
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/5.00000))
      ALGBRC(17) = 44.3000+ 257.100*EXP(- (STATES(1)+32.5000)/13.9000 ** 2.00000)
      RATES(6) = (ALGBRC(3) - STATES(6))/ALGBRC(17)
      ALGBRC(4) = CONSTS(10)/(CONSTS(10)+STATES(4))
      ALGBRC(18) = ALGBRC(4)/CONSTS(63)
      RATES(7) = (ALGBRC(4) - STATES(7))/ALGBRC(18)
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+26.3000)/6.00000))
      ALGBRC(19) = 1.00000/( 1.06800*EXP((STATES(1)+26.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+26.3000)/30.0000))
      RATES(8) = (ALGBRC(5) - STATES(8))/ALGBRC(19)
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+61.7000)/5.60000))
      ALGBRC(20) = 1.00000/( 0.0153000*EXP(- (STATES(1)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(1)+61.7000)/15.3800))
      RATES(9) = (ALGBRC(6) - STATES(9))/ALGBRC(20)
      ALGBRC(7) = 1.00000/(1.00000+EXP(- (STATES(1)+23.2000)/10.6000))
      ALGBRC(21) = 0.846554/( 0.00420000*EXP(STATES(1)/17.0000)+ 0.000150000*EXP(- STATES(1)/21.6000))
      RATES(10) = (ALGBRC(7) - STATES(10))/ALGBRC(21)
      ALGBRC(22) = 0.846554/( 0.0372000*EXP(STATES(1)/15.9000)+ 0.000960000*EXP(- STATES(1)/22.5000))
      RATES(11) = (ALGBRC(7) - STATES(11))/ALGBRC(22)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+28.6000)/17.1000))
      ALGBRC(23) = 1.00000/( 0.100000*EXP(- STATES(1)/54.6450)+ 0.656000*EXP(STATES(1)/106.157))
      RATES(12) = (ALGBRC(8) - STATES(12))/ALGBRC(23)
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+49.0000)/13.0000))
      ALGBRC(25) =  0.600000*(65.1700/( 0.570000*EXP( - 0.0800000*(STATES(1)+44.0000))+ 0.0650000*EXP( 0.100000*(STATES(1)+45.9300)))+10.1000)
      RATES(14) = (ALGBRC(10) - STATES(14))/ALGBRC(25)
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.3000)/15.0000))
      ALGBRC(26) =  0.660000*1.40000*(15.5900/( 1.03700*EXP( 0.0900000*(STATES(1)+30.6100))+ 0.369000*EXP( - 0.120000*(STATES(1)+23.8400)))+2.98000)
      RATES(15) = (ALGBRC(11) - STATES(15))/ALGBRC(26)
      ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(1)+64.0000)/13.5000))
      ALGBRC(27) = 0.716653/(EXP(- (STATES(1)+386.900)/45.3020)+EXP((STATES(1) - 73.0800)/19.2310))
      RATES(16) = (ALGBRC(12) - STATES(16))/ALGBRC(27)
      ALGBRC(2) = 1.00000/(1.00000+EXP(- (STATES(1)+14.1000)/6.00000))
      ALGBRC(16) = ( - 0.0283900*(STATES(1)+35.0000))/(EXP(- (STATES(1)+35.0000)/2.50000) - 1.00000) - ( 0.0849000*STATES(1))/(EXP(- STATES(1)/4.80800) - 1.00000)
      ALGBRC(31) = ( 0.0114300*(STATES(1) - 5.00000))/(EXP((STATES(1) - 5.00000)/2.50000) - 1.00000)
      ALGBRC(36) = 1.00000/(ALGBRC(16)+ALGBRC(31))
      RATES(5) = (ALGBRC(2) - STATES(5))/ALGBRC(36)
      ALGBRC(9) = 0.0140000/(1.00000+EXP(- (STATES(1) - 40.0000)/9.00000))
      ALGBRC(24) =  0.00100000*EXP(- STATES(1)/45.0000)
      ALGBRC(32) = ALGBRC(9)/(ALGBRC(9)+ALGBRC(24))
      ALGBRC(37) = 1.00000/(ALGBRC(9)+ALGBRC(24))
      RATES(13) = (ALGBRC(32) - STATES(13))/ALGBRC(37)
      ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+57.0000)/5.00000))
      ALGBRC(28) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000))
      ALGBRC(33) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000))
      ALGBRC(38) = 1.00000/(ALGBRC(28)+ALGBRC(33))
      RATES(17) = (ALGBRC(13) - STATES(17))/ALGBRC(38)
      ALGBRC(14) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000))
      ALGBRC(29) = 1.00000/( 95.0000*EXP(- STATES(1)/10.0000)+ 50.0000*EXP(- STATES(1)/700.000))+0.000229000/(1.00000+EXP(- STATES(1)/5.00000))
      ALGBRC(34) = ALGBRC(14)/(ALGBRC(14)+ALGBRC(29))
      ALGBRC(39) = 6.65000/(ALGBRC(14)+ALGBRC(29))
      RATES(18) = (ALGBRC(34) - STATES(18))/ALGBRC(39)
      ALGBRC(30) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(6)/STATES(3))
      ALGBRC(41) =  CONSTS(64)*(STATES(1) - ALGBRC(30))*( 0.600000*STATES(11)+ 0.400000*STATES(10))*STATES(12)
      ALGBRC(42) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log((CONSTS(6)+ 0.120000*CONSTS(5))/(STATES(3)+ 0.120000*STATES(2)))
      ALGBRC(43) =  CONSTS(13)*(STATES(1) - ALGBRC(42))*STATES(13) ** 2.00000
      ALGBRC(44) =  CONSTS(14)*(STATES(1) - ALGBRC(30))*STATES(14)*STATES(15)
      ALGBRC(45) =  CONSTS(15)*(STATES(1) - ALGBRC(30))*STATES(15)
      ALGBRC(1) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(5)/STATES(2))
      ALGBRC(52) =  CONSTS(23)*1.00000+CONSTS(21)/CONSTS(6) ** 1.20000 ** - 1.00000*1.00000+CONSTS(22)/STATES(2) ** 1.30000 ** - 1.00000*1.00000+EXP(- ((STATES(1) - ALGBRC(1))+120.000)/30.0000) ** - 1.00000
      ALGBRC(51) = ( CONSTS(65)*(STATES(3) -  CONSTS(6)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))/1.00000
      ALGBRC(47) =  CONSTS(17)*(STATES(1) - ALGBRC(30))*STATES(16) ** 2.00000
      RATES(3) = ( - (ALGBRC(41)+ALGBRC(43)+ALGBRC(44)+ALGBRC(45)+ALGBRC(47)+ALGBRC(51)+ - 2.00000*ALGBRC(52))*CONSTS(4))/( 1.00000*CONSTS(3)*(CONSTS(45)+CONSTS(48)))
      ALGBRC(40) =  CONSTS(11)*(STATES(1) - CONSTS(12))*STATES(8)*STATES(9)
      ALGBRC(35) =  CONSTS(7)*(STATES(1) - CONSTS(8))*STATES(5)*STATES(6)*STATES(7)
      ALGBRC(46) =  CONSTS(16)*(STATES(1) - ALGBRC(1))*STATES(16) ** 2.00000
      ALGBRC(48) = ALGBRC(46)+ALGBRC(47)
      ALGBRC(49) =  CONSTS(18)*(STATES(1) - CONSTS(19))*STATES(17)*STATES(18)
      ALGBRC(50) =  CONSTS(20)*(STATES(1) - ALGBRC(1))
      ALGBRC(54) = EXP(( - CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(60) = 1.00000+ (CONSTS(37)/CONSTS(36))*(1.00000+EXP(( CONSTS(27)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))+ (CONSTS(5)/CONSTS(34))*(1.00000+CONSTS(5)/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(33))
      ALGBRC(62) = ( (( (CONSTS(5)/CONSTS(34))*CONSTS(5))/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(33))*EXP(( - CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(60)
      ALGBRC(61) = ( (CONSTS(37)/CONSTS(36))*EXP(( CONSTS(27)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(60)
      ALGBRC(58) = EXP(( CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(53) = STATES(2)/(CONSTS(28)+STATES(2))
      ALGBRC(63) =  ALGBRC(54)*CONSTS(66)*(ALGBRC(62)+ALGBRC(61))+ ALGBRC(61)*ALGBRC(58)*(ALGBRC(53)+ALGBRC(54))
      ALGBRC(55) = 1.00000+ (STATES(4)/CONSTS(29))*(1.00000+EXP(( - CONSTS(25)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+STATES(2)/CONSTS(32))+ (STATES(2)/CONSTS(30))*(1.00000+ (STATES(2)/CONSTS(31))*(1.00000+STATES(2)/CONSTS(28)))
      ALGBRC(56) = ( (STATES(4)/CONSTS(29))*EXP(( - CONSTS(25)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(55)
      ALGBRC(57) = ( (( (STATES(2)/CONSTS(30))*STATES(2))/CONSTS(31))*(1.00000+STATES(2)/CONSTS(28))*EXP(( CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(55)
      ALGBRC(59) =  ALGBRC(58)*ALGBRC(53)*(ALGBRC(57)+ALGBRC(56))+ ALGBRC(54)*ALGBRC(56)*(CONSTS(66)+ALGBRC(58))
      ALGBRC(64) =  ALGBRC(57)*ALGBRC(53)*(ALGBRC(62)+ALGBRC(61))+ ALGBRC(56)*ALGBRC(62)*(ALGBRC(53)+ALGBRC(54))
      ALGBRC(65) =  ALGBRC(62)*CONSTS(66)*(ALGBRC(57)+ALGBRC(56))+ ALGBRC(57)*ALGBRC(61)*(CONSTS(66)+ALGBRC(58))
      ALGBRC(66) = ( CONSTS(24)*( ALGBRC(59)*ALGBRC(61) -  ALGBRC(63)*ALGBRC(56)))/(ALGBRC(63)+ALGBRC(59)+ALGBRC(64)+ALGBRC(65))
      RATES(1) = - (ALGBRC(35)+ALGBRC(40)+ALGBRC(41)+ALGBRC(43)+ALGBRC(44)+ALGBRC(45)+ALGBRC(48)+ALGBRC(49)+ALGBRC(50)+ALGBRC(51)+ALGBRC(52)+ALGBRC(66))/1.00000
      RATES(2) = ( - (ALGBRC(46)+ALGBRC(49)+ALGBRC(50)+ 3.00000*ALGBRC(52)+ 3.00000*ALGBRC(66))*CONSTS(4))/( 1.00000*CONSTS(3)*(CONSTS(45)+CONSTS(48)))
      ALGBRC(67) = (STATES(4) - STATES(20))/CONSTS(38)
      ALGBRC(68) = ( CONSTS(42)*(STATES(21) - STATES(4)))/(1.00000+CONSTS(40)/STATES(4) ** 2.00000)
      ALGBRC(70) =  CONSTS(56)*STATES(4)*(1.00000 - STATES(24)) -  CONSTS(61)*STATES(24)
      RATES(4) = (( - ((ALGBRC(35)+ALGBRC(40)) -  2.00000*ALGBRC(66))*CONSTS(4))/( 1.00000*2.00000*CONSTS(3))+ ALGBRC(68)*CONSTS(46))/CONSTS(48) - (ALGBRC(67)+ CONSTS(50)*ALGBRC(70))
      RATES(24) = ALGBRC(70)
      ALGBRC(69) = CONSTS(41)/(1.00000+CONSTS(43)/STATES(20))
      ALGBRC(71) = (STATES(19) - STATES(21))/CONSTS(39)
      RATES(19) = ALGBRC(69) - ( ALGBRC(71)*CONSTS(46))/CONSTS(47)
      ALGBRC(72) =  CONSTS(53)*STATES(20)*(1.00000 - STATES(25)) -  CONSTS(58)*STATES(25)
      RATES(25) = ALGBRC(72)
      ALGBRC(74) =  CONSTS(57)*STATES(21)*(1.00000 - STATES(27)) -  CONSTS(62)*STATES(27)
      RATES(21) = ALGBRC(71) - (ALGBRC(68)+ CONSTS(52)*ALGBRC(74))
      ALGBRC(73) =  CONSTS(55)*STATES(20)*(1.00000 - (STATES(26)+STATES(22))) -  CONSTS(59)*STATES(26)
      RATES(26) = ALGBRC(73)
      RATES(27) = ALGBRC(74)
      ALGBRC(75) =  CONSTS(56)*STATES(20)*(1.00000 - STATES(23)) -  CONSTS(61)*STATES(23)
      RATES(20) = ( ALGBRC(67)*CONSTS(48) -  ALGBRC(69)*CONSTS(47))/CONSTS(45) - ( CONSTS(50)*ALGBRC(75)+ CONSTS(51)*ALGBRC(72)+ CONSTS(49)*ALGBRC(73))
      RATES(23) = ALGBRC(75)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(15) =  CONSTS(54)*CONSTS(44)*(1.00000 - (STATES(26)+STATES(22))) -  CONSTS(60)*STATES(22)
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/5.00000))
      ALGBRC(17) = 44.3000+ 257.100*EXP(- (STATES(1)+32.5000)/13.9000 ** 2.00000)
      ALGBRC(4) = CONSTS(10)/(CONSTS(10)+STATES(4))
      ALGBRC(18) = ALGBRC(4)/CONSTS(63)
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+26.3000)/6.00000))
      ALGBRC(19) = 1.00000/( 1.06800*EXP((STATES(1)+26.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+26.3000)/30.0000))
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+61.7000)/5.60000))
      ALGBRC(20) = 1.00000/( 0.0153000*EXP(- (STATES(1)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(1)+61.7000)/15.3800))
      ALGBRC(7) = 1.00000/(1.00000+EXP(- (STATES(1)+23.2000)/10.6000))
      ALGBRC(21) = 0.846554/( 0.00420000*EXP(STATES(1)/17.0000)+ 0.000150000*EXP(- STATES(1)/21.6000))
      ALGBRC(22) = 0.846554/( 0.0372000*EXP(STATES(1)/15.9000)+ 0.000960000*EXP(- STATES(1)/22.5000))
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+28.6000)/17.1000))
      ALGBRC(23) = 1.00000/( 0.100000*EXP(- STATES(1)/54.6450)+ 0.656000*EXP(STATES(1)/106.157))
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+49.0000)/13.0000))
      ALGBRC(25) =  0.600000*(65.1700/( 0.570000*EXP( - 0.0800000*(STATES(1)+44.0000))+ 0.0650000*EXP( 0.100000*(STATES(1)+45.9300)))+10.1000)
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.3000)/15.0000))
      ALGBRC(26) =  0.660000*1.40000*(15.5900/( 1.03700*EXP( 0.0900000*(STATES(1)+30.6100))+ 0.369000*EXP( - 0.120000*(STATES(1)+23.8400)))+2.98000)
      ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(1)+64.0000)/13.5000))
      ALGBRC(27) = 0.716653/(EXP(- (STATES(1)+386.900)/45.3020)+EXP((STATES(1) - 73.0800)/19.2310))
      ALGBRC(2) = 1.00000/(1.00000+EXP(- (STATES(1)+14.1000)/6.00000))
      ALGBRC(16) = ( - 0.0283900*(STATES(1)+35.0000))/(EXP(- (STATES(1)+35.0000)/2.50000) - 1.00000) - ( 0.0849000*STATES(1))/(EXP(- STATES(1)/4.80800) - 1.00000)
      ALGBRC(31) = ( 0.0114300*(STATES(1) - 5.00000))/(EXP((STATES(1) - 5.00000)/2.50000) - 1.00000)
      ALGBRC(36) = 1.00000/(ALGBRC(16)+ALGBRC(31))
      ALGBRC(9) = 0.0140000/(1.00000+EXP(- (STATES(1) - 40.0000)/9.00000))
      ALGBRC(24) =  0.00100000*EXP(- STATES(1)/45.0000)
      ALGBRC(32) = ALGBRC(9)/(ALGBRC(9)+ALGBRC(24))
      ALGBRC(37) = 1.00000/(ALGBRC(9)+ALGBRC(24))
      ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+57.0000)/5.00000))
      ALGBRC(28) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000))
      ALGBRC(33) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000))
      ALGBRC(38) = 1.00000/(ALGBRC(28)+ALGBRC(33))
      ALGBRC(14) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000))
      ALGBRC(29) = 1.00000/( 95.0000*EXP(- STATES(1)/10.0000)+ 50.0000*EXP(- STATES(1)/700.000))+0.000229000/(1.00000+EXP(- STATES(1)/5.00000))
      ALGBRC(34) = ALGBRC(14)/(ALGBRC(14)+ALGBRC(29))
      ALGBRC(39) = 6.65000/(ALGBRC(14)+ALGBRC(29))
      ALGBRC(30) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(6)/STATES(3))
      ALGBRC(41) =  CONSTS(64)*(STATES(1) - ALGBRC(30))*( 0.600000*STATES(11)+ 0.400000*STATES(10))*STATES(12)
      ALGBRC(42) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log((CONSTS(6)+ 0.120000*CONSTS(5))/(STATES(3)+ 0.120000*STATES(2)))
      ALGBRC(43) =  CONSTS(13)*(STATES(1) - ALGBRC(42))*STATES(13) ** 2.00000
      ALGBRC(44) =  CONSTS(14)*(STATES(1) - ALGBRC(30))*STATES(14)*STATES(15)
      ALGBRC(45) =  CONSTS(15)*(STATES(1) - ALGBRC(30))*STATES(15)
      ALGBRC(1) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(5)/STATES(2))
      ALGBRC(52) =  CONSTS(23)*1.00000+CONSTS(21)/CONSTS(6) ** 1.20000 ** - 1.00000*1.00000+CONSTS(22)/STATES(2) ** 1.30000 ** - 1.00000*1.00000+EXP(- ((STATES(1) - ALGBRC(1))+120.000)/30.0000) ** - 1.00000
      ALGBRC(51) = ( CONSTS(65)*(STATES(3) -  CONSTS(6)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))/1.00000
      ALGBRC(47) =  CONSTS(17)*(STATES(1) - ALGBRC(30))*STATES(16) ** 2.00000
      ALGBRC(40) =  CONSTS(11)*(STATES(1) - CONSTS(12))*STATES(8)*STATES(9)
      ALGBRC(35) =  CONSTS(7)*(STATES(1) - CONSTS(8))*STATES(5)*STATES(6)*STATES(7)
      ALGBRC(46) =  CONSTS(16)*(STATES(1) - ALGBRC(1))*STATES(16) ** 2.00000
      ALGBRC(48) = ALGBRC(46)+ALGBRC(47)
      ALGBRC(49) =  CONSTS(18)*(STATES(1) - CONSTS(19))*STATES(17)*STATES(18)
      ALGBRC(50) =  CONSTS(20)*(STATES(1) - ALGBRC(1))
      ALGBRC(54) = EXP(( - CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(60) = 1.00000+ (CONSTS(37)/CONSTS(36))*(1.00000+EXP(( CONSTS(27)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))+ (CONSTS(5)/CONSTS(34))*(1.00000+CONSTS(5)/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(33))
      ALGBRC(62) = ( (( (CONSTS(5)/CONSTS(34))*CONSTS(5))/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(33))*EXP(( - CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(60)
      ALGBRC(61) = ( (CONSTS(37)/CONSTS(36))*EXP(( CONSTS(27)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(60)
      ALGBRC(58) = EXP(( CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(53) = STATES(2)/(CONSTS(28)+STATES(2))
      ALGBRC(63) =  ALGBRC(54)*CONSTS(66)*(ALGBRC(62)+ALGBRC(61))+ ALGBRC(61)*ALGBRC(58)*(ALGBRC(53)+ALGBRC(54))
      ALGBRC(55) = 1.00000+ (STATES(4)/CONSTS(29))*(1.00000+EXP(( - CONSTS(25)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+STATES(2)/CONSTS(32))+ (STATES(2)/CONSTS(30))*(1.00000+ (STATES(2)/CONSTS(31))*(1.00000+STATES(2)/CONSTS(28)))
      ALGBRC(56) = ( (STATES(4)/CONSTS(29))*EXP(( - CONSTS(25)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(55)
      ALGBRC(57) = ( (( (STATES(2)/CONSTS(30))*STATES(2))/CONSTS(31))*(1.00000+STATES(2)/CONSTS(28))*EXP(( CONSTS(26)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(55)
      ALGBRC(59) =  ALGBRC(58)*ALGBRC(53)*(ALGBRC(57)+ALGBRC(56))+ ALGBRC(54)*ALGBRC(56)*(CONSTS(66)+ALGBRC(58))
      ALGBRC(64) =  ALGBRC(57)*ALGBRC(53)*(ALGBRC(62)+ALGBRC(61))+ ALGBRC(56)*ALGBRC(62)*(ALGBRC(53)+ALGBRC(54))
      ALGBRC(65) =  ALGBRC(62)*CONSTS(66)*(ALGBRC(57)+ALGBRC(56))+ ALGBRC(57)*ALGBRC(61)*(CONSTS(66)+ALGBRC(58))
      ALGBRC(66) = ( CONSTS(24)*( ALGBRC(59)*ALGBRC(61) -  ALGBRC(63)*ALGBRC(56)))/(ALGBRC(63)+ALGBRC(59)+ALGBRC(64)+ALGBRC(65))
      ALGBRC(67) = (STATES(4) - STATES(20))/CONSTS(38)
      ALGBRC(68) = ( CONSTS(42)*(STATES(21) - STATES(4)))/(1.00000+CONSTS(40)/STATES(4) ** 2.00000)
      ALGBRC(70) =  CONSTS(56)*STATES(4)*(1.00000 - STATES(24)) -  CONSTS(61)*STATES(24)
      ALGBRC(69) = CONSTS(41)/(1.00000+CONSTS(43)/STATES(20))
      ALGBRC(71) = (STATES(19) - STATES(21))/CONSTS(39)
      ALGBRC(72) =  CONSTS(53)*STATES(20)*(1.00000 - STATES(25)) -  CONSTS(58)*STATES(25)
      ALGBRC(74) =  CONSTS(57)*STATES(21)*(1.00000 - STATES(27)) -  CONSTS(62)*STATES(27)
      ALGBRC(73) =  CONSTS(55)*STATES(20)*(1.00000 - (STATES(26)+STATES(22))) -  CONSTS(59)*STATES(26)
      ALGBRC(75) =  CONSTS(56)*STATES(20)*(1.00000 - STATES(23)) -  CONSTS(61)*STATES(23)
      RETURN
      END