C
C There are a total of 76 entries in the algebraic variable array.
C There are a total of 29 entries in each of the rate and state variable arrays.
C There are a total of 86 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is Vm in component Vm (millivolt).
C CONSTS(1) is Cm in component model_parameters (picoF).
C ALGBRC(32) is i_CaT in component i_CaT (picoA).
C ALGBRC(15) is i_CaL in component i_CaL (picoA).
C ALGBRC(50) is i_f in component i_f (picoA).
C ALGBRC(51) is i_st in component i_st (picoA).
C ALGBRC(38) is i_Kr in component i_Kr (picoA).
C ALGBRC(43) is i_Ks in component i_Ks (picoA).
C ALGBRC(45) is i_to in component AP_sensitive_currents (picoA).
C ALGBRC(47) is i_sus in component AP_sensitive_currents (picoA).
C ALGBRC(53) is i_NaK in component i_NaK (picoA).
C ALGBRC(67) is i_NaCa in component i_NaCa (picoA).
C ALGBRC(54) is i_b_Ca in component i_b_Ca (picoA).
C ALGBRC(52) is i_b_Na in component i_b_Na (picoA).
C CONSTS(76) is E_Na in component electric_potentials (millivolt).
C CONSTS(77) is E_K in component electric_potentials (millivolt).
C CONSTS(78) is E_Ks in component electric_potentials (millivolt).
C CONSTS(2) is R in component model_parameters (joule_per_kilomole_kelvin).
C CONSTS(3) is T in component model_parameters (kelvin).
C CONSTS(4) is F in component model_parameters (coulomb_per_mole).
C CONSTS(5) is Nai in component model_parameters (millimolar).
C CONSTS(6) is Nao in component model_parameters (millimolar).
C CONSTS(7) is Ki in component model_parameters (millimolar).
C CONSTS(8) is Ko in component model_parameters (millimolar).
C CONSTS(9) is E_CaL in component i_CaL (millivolt).
C CONSTS(10) is g_CaL in component i_CaL (nanoS_per_picoF).
C STATES(2) is Ca_sub in component calcium_dynamics (millimolar).
C STATES(3) is dL in component i_CaL_dL_gate (dimensionless).
C STATES(4) is fL in component i_CaL_fL_gate (dimensionless).
C STATES(5) is fCa in component i_CaL_fCa_gate (dimensionless).
C ALGBRC(1) is dL_infinity in component i_CaL_dL_gate (dimensionless).
C ALGBRC(46) is tau_dL in component i_CaL_dL_gate (millisecond).
C ALGBRC(33) is alpha_dL in component i_CaL_dL_gate (per_millisecond).
C ALGBRC(44) is beta_dL in component i_CaL_dL_gate (per_millisecond).
C ALGBRC(17) is adVm in component i_CaL_dL_gate (millivolt).
C ALGBRC(39) is bdVm in component i_CaL_dL_gate (millivolt).
C ALGBRC(2) is fL_infinity in component i_CaL_fL_gate (dimensionless).
C ALGBRC(18) is tau_fL in component i_CaL_fL_gate (millisecond).
C CONSTS(11) is alpha_fCa in component i_CaL_fCa_gate (per_millisecond).
C ALGBRC(3) is fCa_infinity in component i_CaL_fCa_gate (dimensionless).
C ALGBRC(19) is tau_fCa in component i_CaL_fCa_gate (millisecond).
C CONSTS(12) is Km_fCa in component i_CaL_fCa_gate (millimolar).
C CONSTS(13) is g_CaT in component i_CaT (nanoS_per_picoF).
C CONSTS(14) is E_CaT in component i_CaT (millivolt).
C STATES(6) is dT in component i_CaT_dT_gate (dimensionless).
C STATES(7) is fT in component i_CaT_fT_gate (dimensionless).
C ALGBRC(4) is dT_infinity in component i_CaT_dT_gate (dimensionless).
C ALGBRC(20) is tau_dT in component i_CaT_dT_gate (millisecond).
C ALGBRC(5) is fT_infinity in component i_CaT_fT_gate (dimensionless).
C ALGBRC(21) is tau_fT in component i_CaT_fT_gate (millisecond).
C CONSTS(15) is g_Kr in component i_Kr (nanoS_per_picoF).
C STATES(8) is paS in component i_Kr_pa_gate (dimensionless).
C STATES(9) is paF in component i_Kr_pa_gate (dimensionless).
C STATES(10) is pi_ in component i_Kr_pi_gate (dimensionless).
C ALGBRC(6) is pa_infinity in component i_Kr_pa_gate (dimensionless).
C ALGBRC(22) is tau_paS in component i_Kr_pa_gate (millisecond).
C ALGBRC(23) is tau_paF in component i_Kr_pa_gate (millisecond).
C ALGBRC(7) is pi_infinity in component i_Kr_pi_gate (dimensionless).
C ALGBRC(24) is tau_pi in component i_Kr_pi_gate (millisecond).
C CONSTS(16) is g_Ks in component i_Ks (nanoS_per_picoF).
C STATES(11) is n in component i_Ks_n_gate (dimensionless).
C ALGBRC(34) is n_infinity in component i_Ks_n_gate (dimensionless).
C ALGBRC(40) is tau_n in component i_Ks_n_gate (millisecond).
C ALGBRC(8) is alpha_n in component i_Ks_n_gate (per_millisecond).
C ALGBRC(25) is beta_n in component i_Ks_n_gate (per_millisecond).
C CONSTS(17) is g_to in component AP_sensitive_currents (nanoS_per_picoF).
C CONSTS(18) is g_sus in component AP_sensitive_currents (nanoS_per_picoF).
C STATES(12) is q in component AP_sensitive_currents_q_gate (dimensionless).
C STATES(13) is r in component AP_sensitive_currents_r_gate (dimensionless).
C ALGBRC(9) is q_infinity in component AP_sensitive_currents_q_gate (dimensionless).
C ALGBRC(26) is tau_q in component AP_sensitive_currents_q_gate (millisecond).
C ALGBRC(10) is r_infinity in component AP_sensitive_currents_r_gate (dimensionless).
C ALGBRC(27) is tau_r in component AP_sensitive_currents_r_gate (millisecond).
C ALGBRC(48) is i_f_Na in component i_f (picoA).
C ALGBRC(49) is i_f_K in component i_f (picoA).
C CONSTS(19) is g_if in component i_f (nanoS_per_picoF).
C STATES(14) is y in component i_f_y_gate (dimensionless).
C ALGBRC(11) is y_infinity in component i_f_y_gate (dimensionless).
C ALGBRC(28) is tau_y in component i_f_y_gate (millisecond).
C CONSTS(20) is VIf_half in component i_f_y_gate (millivolt).
C CONSTS(21) is g_st in component i_st (nanoS_per_picoF).
C CONSTS(22) is E_st in component i_st (millivolt).
C STATES(15) is qa in component i_st_qa_gate (dimensionless).
C STATES(16) is qi in component i_st_qi_gate (dimensionless).
C ALGBRC(12) is qa_infinity in component i_st_qa_gate (dimensionless).
C ALGBRC(41) is tau_qa in component i_st_qa_gate (millisecond).
C ALGBRC(29) is alpha_qa in component i_st_qa_gate (per_millisecond).
C ALGBRC(35) is beta_qa in component i_st_qa_gate (per_millisecond).
C ALGBRC(36) is qi_infinity in component i_st_qi_gate (dimensionless).
C ALGBRC(42) is tau_qi in component i_st_qi_gate (millisecond).
C ALGBRC(13) is alpha_qi in component i_st_qi_gate (per_millisecond).
C ALGBRC(30) is beta_qi in component i_st_qi_gate (per_millisecond).
C CONSTS(23) is g_b_Na in component i_b_Na (nanoS_per_picoF).
C CONSTS(24) is Km_Kp in component i_NaK (millimolar).
C CONSTS(25) is Km_Nap in component i_NaK (millimolar).
C CONSTS(26) is i_NaK_max in component i_NaK (picoA_per_picoF).
C CONSTS(27) is g_b_Ca in component i_b_Ca (nanoS_per_picoF).
C CONSTS(28) is kNaCa in component i_NaCa (picoA_per_picoF).
C ALGBRC(64) is x1 in component i_NaCa (dimensionless).
C ALGBRC(60) is x2 in component i_NaCa (dimensionless).
C ALGBRC(65) is x3 in component i_NaCa (dimensionless).
C ALGBRC(66) is x4 in component i_NaCa (dimensionless).
C ALGBRC(55) is k41 in component i_NaCa (dimensionless).
C CONSTS(81) is k34 in component i_NaCa (dimensionless).
C ALGBRC(63) is k23 in component i_NaCa (dimensionless).
C ALGBRC(62) is k21 in component i_NaCa (dimensionless).
C ALGBRC(59) is k32 in component i_NaCa (dimensionless).
C CONSTS(80) is k43 in component i_NaCa (dimensionless).
C ALGBRC(57) is k12 in component i_NaCa (dimensionless).
C ALGBRC(58) is k14 in component i_NaCa (dimensionless).
C CONSTS(29) is Qci in component i_NaCa (dimensionless).
C CONSTS(30) is Qn in component i_NaCa (dimensionless).
C CONSTS(31) is Qco in component i_NaCa (dimensionless).
C CONSTS(32) is K3ni in component i_NaCa (millimolar).
C CONSTS(33) is Kci in component i_NaCa (millimolar).
C CONSTS(34) is K1ni in component i_NaCa (millimolar).
C CONSTS(35) is K2ni in component i_NaCa (millimolar).
C CONSTS(36) is Kcni in component i_NaCa (millimolar).
C CONSTS(37) is K3no in component i_NaCa (millimolar).
C CONSTS(38) is K1no in component i_NaCa (millimolar).
C CONSTS(39) is K2no in component i_NaCa (millimolar).
C CONSTS(40) is Kco in component i_NaCa (millimolar).
C CONSTS(79) is RTOnF in component i_NaCa (millivolt).
C ALGBRC(61) is do in component i_NaCa (dimensionless).
C ALGBRC(56) is di in component i_NaCa (dimensionless).
C CONSTS(41) is Cao in component model_parameters (millimolar).
C ALGBRC(68) is j_SRCarel in component j_SRCarel (millimolar_per_millisecond).
C STATES(17) is R in component j_SRCarel (dimensionless).
C STATES(18) is O in component j_SRCarel (dimensionless).
C STATES(19) is I in component j_SRCarel (dimensionless).
C STATES(20) is RI in component j_SRCarel (dimensionless).
C CONSTS(42) is ks in component j_SRCarel (per_millisecond).
C CONSTS(43) is MaxSR in component j_SRCarel (dimensionless).
C CONSTS(44) is MinSR in component j_SRCarel (dimensionless).
C CONSTS(45) is EC50_SR in component j_SRCarel (millimolar).
C CONSTS(46) is HSR in component j_SRCarel (dimensionless).
C ALGBRC(31) is koSRCa in component j_SRCarel (per_millimolar2_millisecond).
C ALGBRC(37) is kiSRCa in component j_SRCarel (per_millimolar_millisecond).
C CONSTS(47) is koCa in component j_SRCarel (per_millimolar2_millisecond).
C CONSTS(48) is kiCa in component j_SRCarel (per_millimolar_millisecond).
C ALGBRC(14) is kCaSR in component j_SRCarel (dimensionless).
C CONSTS(49) is kim in component j_SRCarel (per_millisecond).
C CONSTS(50) is kom in component j_SRCarel (per_millisecond).
C STATES(21) is Ca_jsr in component calcium_dynamics (millimolar).
C ALGBRC(69) is j_Ca_dif in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C ALGBRC(71) is j_up in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C ALGBRC(73) is j_tr in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C CONSTS(51) is tau_dif_Ca in component intracellular_calcium_fluxes (millisecond).
C CONSTS(52) is tau_tr in component intracellular_calcium_fluxes (millisecond).
C CONSTS(53) is P_up in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C CONSTS(54) is K_up in component intracellular_calcium_fluxes (millimolar).
C STATES(22) is Ca_nsr in component calcium_dynamics (millimolar).
C STATES(23) is Cai in component calcium_dynamics (millimolar).
C CONSTS(55) is TC_tot in component calcium_buffering (millimolar).
C CONSTS(56) is TMC_tot in component calcium_buffering (millimolar).
C CONSTS(57) is CM_tot in component calcium_buffering (millimolar).
C CONSTS(58) is CQ_tot in component calcium_buffering (millimolar).
C ALGBRC(72) is delta_fTC in component calcium_buffering (per_millisecond).
C ALGBRC(74) is delta_fTMC in component calcium_buffering (per_millisecond).
C ALGBRC(70) is delta_fCMs in component calcium_buffering (per_millisecond).
C ALGBRC(76) is delta_fCMi in component calcium_buffering (per_millisecond).
C ALGBRC(75) is delta_fCQ in component calcium_buffering (per_millisecond).
C ALGBRC(16) is delta_fTMM in component calcium_buffering (per_millisecond).
C STATES(24) is fTMM in component calcium_buffering (dimensionless).
C STATES(25) is fCMi in component calcium_buffering (dimensionless).
C STATES(26) is fCMs in component calcium_buffering (dimensionless).
C STATES(27) is fTC in component calcium_buffering (dimensionless).
C STATES(28) is fTMC in component calcium_buffering (dimensionless).
C STATES(29) is fCQ in component calcium_buffering (dimensionless).
C CONSTS(59) is kf_TC in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(60) is kf_TMM in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(61) is kf_TMC in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(62) is kf_CM in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(63) is kf_CQ in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(64) is kb_TC in component calcium_buffering (per_millisecond).
C CONSTS(65) is kb_TMC in component calcium_buffering (per_millisecond).
C CONSTS(66) is kb_TMM in component calcium_buffering (per_millisecond).
C CONSTS(67) is kb_CM in component calcium_buffering (per_millisecond).
C CONSTS(68) is kb_CQ in component calcium_buffering (per_millisecond).
C CONSTS(69) is Mgi in component model_parameters (millimolar).
C CONSTS(85) is V_i in component model_parameters (picolitre).
C CONSTS(84) is V_jsr in component model_parameters (picolitre).
C CONSTS(86) is V_nsr in component model_parameters (picolitre).
C CONSTS(83) is V_sub in component model_parameters (picolitre).
C CONSTS(82) is V_cell in component model_parameters (picolitre).
C CONSTS(70) is V_jsr_part in component model_parameters (dimensionless).
C CONSTS(71) is V_i_part in component model_parameters (dimensionless).
C CONSTS(72) is V_nsr_part in component model_parameters (dimensionless).
C CONSTS(73) is R_cell in component model_parameters (micrometre).
C CONSTS(74) is L_cell in component model_parameters (micrometre).
C CONSTS(75) is L_sub in component model_parameters (micrometre).
C RATES(1) is d/dt Vm in component Vm (millivolt).
C RATES(3) is d/dt dL in component i_CaL_dL_gate (dimensionless).
C RATES(4) is d/dt fL in component i_CaL_fL_gate (dimensionless).
C RATES(5) is d/dt fCa in component i_CaL_fCa_gate (dimensionless).
C RATES(6) is d/dt dT in component i_CaT_dT_gate (dimensionless).
C RATES(7) is d/dt fT in component i_CaT_fT_gate (dimensionless).
C RATES(8) is d/dt paS in component i_Kr_pa_gate (dimensionless).
C RATES(9) is d/dt paF in component i_Kr_pa_gate (dimensionless).
C RATES(10) is d/dt pi_ in component i_Kr_pi_gate (dimensionless).
C RATES(11) is d/dt n in component i_Ks_n_gate (dimensionless).
C RATES(12) is d/dt q in component AP_sensitive_currents_q_gate (dimensionless).
C RATES(13) is d/dt r in component AP_sensitive_currents_r_gate (dimensionless).
C RATES(14) is d/dt y in component i_f_y_gate (dimensionless).
C RATES(15) is d/dt qa in component i_st_qa_gate (dimensionless).
C RATES(16) is d/dt qi in component i_st_qi_gate (dimensionless).
C RATES(17) is d/dt R in component j_SRCarel (dimensionless).
C RATES(18) is d/dt O in component j_SRCarel (dimensionless).
C RATES(19) is d/dt I in component j_SRCarel (dimensionless).
C RATES(20) is d/dt RI in component j_SRCarel (dimensionless).
C RATES(27) is d/dt fTC in component calcium_buffering (dimensionless).
C RATES(28) is d/dt fTMC in component calcium_buffering (dimensionless).
C RATES(24) is d/dt fTMM in component calcium_buffering (dimensionless).
C RATES(25) is d/dt fCMi in component calcium_buffering (dimensionless).
C RATES(26) is d/dt fCMs in component calcium_buffering (dimensionless).
C RATES(29) is d/dt fCQ in component calcium_buffering (dimensionless).
C RATES(23) is d/dt Cai in component calcium_dynamics (millimolar).
C RATES(2) is d/dt Ca_sub in component calcium_dynamics (millimolar).
C RATES(22) is d/dt Ca_nsr in component calcium_dynamics (millimolar).
C RATES(21) is d/dt Ca_jsr in component calcium_dynamics (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -65
      CONSTS(1) = 32
      CONSTS(2) = 8314.4
      CONSTS(3) = 310.15
      CONSTS(4) = 96485
      CONSTS(5) = 10
      CONSTS(6) = 140
      CONSTS(7) = 140
      CONSTS(8) = 5.4
      CONSTS(9) = 45
      CONSTS(10) = 0.464
      STATES(2) = 0.000223
      STATES(3) = 0
      STATES(4) = 1
      STATES(5) = 1
      CONSTS(11) = 0.021
      CONSTS(12) = 0.00035
      CONSTS(13) = 0.1832
      CONSTS(14) = 45
      STATES(6) = 0
      STATES(7) = 1
      CONSTS(15) = 0.08113973
      STATES(8) = 0
      STATES(9) = 0
      STATES(10) = 1
      CONSTS(16) = 0.0259
      STATES(11) = 0
      CONSTS(17) = 0.252
      CONSTS(18) = 0.02
      STATES(12) = 1
      STATES(13) = 0
      CONSTS(19) = 0.15
      STATES(14) = 1
      CONSTS(20) = -64
      CONSTS(21) = 0.003
      CONSTS(22) = 37.4
      STATES(15) = 0
      STATES(16) = 1
      CONSTS(23) = 0.00486
      CONSTS(24) = 1.4
      CONSTS(25) = 14
      CONSTS(26) = 2.88
      CONSTS(27) = 0.0006
      CONSTS(28) = 187.5
      CONSTS(29) = 0.1369
      CONSTS(30) = 0.4315
      CONSTS(31) = 0
      CONSTS(32) = 26.44
      CONSTS(33) = 0.0207
      CONSTS(34) = 395.3
      CONSTS(35) = 2.289
      CONSTS(36) = 26.44
      CONSTS(37) = 4.663
      CONSTS(38) = 1628
      CONSTS(39) = 561.4
      CONSTS(40) = 3.663
      CONSTS(41) = 2
      STATES(17) = 0.7499955
      STATES(18) = 3.4e-6
      STATES(19) = 1.1e-6
      STATES(20) = 0.25
      CONSTS(42) = 250000
      CONSTS(43) = 15
      CONSTS(44) = 1
      CONSTS(45) = 0.45
      CONSTS(46) = 2.5
      CONSTS(47) = 10
      CONSTS(48) = 0.5
      CONSTS(49) = 0.005
      CONSTS(50) = 0.06
      STATES(21) = 0.029
      CONSTS(51) = 0.04
      CONSTS(52) = 40
      CONSTS(53) = 0.012
      CONSTS(54) = 0.0006
      STATES(22) = 1.35
      STATES(23) = 0.0001
      CONSTS(55) = 0.031
      CONSTS(56) = 0.062
      CONSTS(57) = 0.045
      CONSTS(58) = 10
      STATES(24) = 0.69
      STATES(25) = 0.042
      STATES(26) = 0.089
      STATES(27) = 0.02
      STATES(28) = 0.22
      STATES(29) = 0.032
      CONSTS(59) = 88.8
      CONSTS(60) = 2.277
      CONSTS(61) = 227.7
      CONSTS(62) = 227.7
      CONSTS(63) = 0.534
      CONSTS(64) = 0.446
      CONSTS(65) = 0.00751
      CONSTS(66) = 0.751
      CONSTS(67) = 0.542
      CONSTS(68) = 0.445
      CONSTS(69) = 2.5
      CONSTS(70) = 0.0012
      CONSTS(71) = 0.46
      CONSTS(72) = 0.0116
      CONSTS(73) = 4
      CONSTS(74) = 70
      CONSTS(75) = 0.02
      CONSTS(76) =  (( CONSTS(2)*CONSTS(3))/CONSTS(4))*log(CONSTS(6)/CONSTS(5))
      CONSTS(77) =  (( CONSTS(2)*CONSTS(3))/CONSTS(4))*log(CONSTS(8)/CONSTS(7))
      CONSTS(78) =  (( CONSTS(2)*CONSTS(3))/CONSTS(4))*log((CONSTS(8)+ 0.120000*CONSTS(6))/(CONSTS(7)+ 0.120000*CONSTS(5)))
      CONSTS(79) = ( CONSTS(2)*CONSTS(3))/CONSTS(4)
      CONSTS(80) = CONSTS(5)/(CONSTS(32)+CONSTS(5))
      CONSTS(81) = CONSTS(6)/(CONSTS(37)+CONSTS(6))
      CONSTS(82) =  0.00100000* 3.14159265358979*CONSTS(73) ** 2.00000*CONSTS(74)
      CONSTS(83) =  0.00100000*2.00000* 3.14159265358979*CONSTS(75)*(CONSTS(73) - CONSTS(75)/2.00000)*CONSTS(74)
      CONSTS(84) =  CONSTS(70)*CONSTS(82)
      CONSTS(85) =  CONSTS(71)*CONSTS(82) - CONSTS(83)
      CONSTS(86) =  CONSTS(72)*CONSTS(82)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(16) =  CONSTS(60)*CONSTS(69)*(1.00000 - (STATES(28)+STATES(24))) -  CONSTS(66)*STATES(24)
      RATES(24) = ALGBRC(16)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/7.30000))
      ALGBRC(18) = 44.3000+ 257.100*EXP(- (STATES(1)+32.5000)/13.9000 ** 2.00000)
      RATES(4) = (ALGBRC(2) - STATES(4))/ALGBRC(18)
      ALGBRC(3) = CONSTS(12)/(CONSTS(12)+STATES(2))
      ALGBRC(19) = ALGBRC(3)/CONSTS(11)
      RATES(5) = (ALGBRC(3) - STATES(5))/ALGBRC(19)
      ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+26.3000)/6.00000))
      ALGBRC(20) = 1.00000/( 1.06800*EXP((STATES(1)+26.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+26.3000)/30.0000))
      RATES(6) = (ALGBRC(4) - STATES(6))/ALGBRC(20)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+61.7000)/5.60000))
      ALGBRC(21) = 1.00000/( 0.0153000*EXP(- (STATES(1)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(1)+61.7000)/15.3800))
      RATES(7) = (ALGBRC(5) - STATES(7))/ALGBRC(21)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1)+23.2000)/10.6000))
      ALGBRC(22) = 0.846554/( 0.00420000*EXP(STATES(1)/17.0000)+ 0.000150000*EXP(- STATES(1)/21.6000))
      RATES(8) = (ALGBRC(6) - STATES(8))/ALGBRC(22)
      ALGBRC(23) = 0.846554/( 0.0372000*EXP(STATES(1)/15.9000)+ 0.000960000*EXP(- STATES(1)/22.5000))
      RATES(9) = (ALGBRC(6) - STATES(9))/ALGBRC(23)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+28.6000)/17.1000))
      ALGBRC(24) = 1.00000/( 0.100000*EXP(- STATES(1)/54.6450)+ 0.656000*EXP(STATES(1)/106.157))
      RATES(10) = (ALGBRC(7) - STATES(10))/ALGBRC(24)
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+49.0000)/13.0000))
      ALGBRC(26) = 6.06000+39.1020/( 0.570000*EXP( - 0.0800000*(STATES(1)+44.0000))+ 0.0650000*EXP( 0.100000*(STATES(1)+45.9300)))
      RATES(12) = (ALGBRC(9) - STATES(12))/ALGBRC(26)
      ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.3000)/15.0000))
      ALGBRC(27) = 2.75352+14.4052/( 1.03700*EXP( 0.0900000*(STATES(1)+30.6100))+ 0.369000*EXP( - 0.120000*(STATES(1)+23.8400)))
      RATES(13) = (ALGBRC(10) - STATES(13))/ALGBRC(27)
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1) - CONSTS(20))/13.5000))
      ALGBRC(28) = 0.716653/(EXP(- (STATES(1)+386.900)/45.3020)+EXP((STATES(1) - 73.0800)/19.2310))
      RATES(14) = (ALGBRC(11) - STATES(14))/ALGBRC(28)
      ALGBRC(14) = CONSTS(43) - (CONSTS(43) - CONSTS(44))/(1.00000+CONSTS(45)/STATES(21) ** CONSTS(46))
      ALGBRC(31) = CONSTS(47)/ALGBRC(14)
      ALGBRC(37) =  CONSTS(48)*ALGBRC(14)
      RATES(17) = ( CONSTS(49)*STATES(20) -  ALGBRC(37)*STATES(2)*STATES(17)) - ( ALGBRC(31)*STATES(2) ** 2.00000*STATES(17) -  CONSTS(50)*STATES(18))
      RATES(18) = ( ALGBRC(31)*STATES(2) ** 2.00000*STATES(17) -  CONSTS(50)*STATES(18)) - ( ALGBRC(37)*STATES(2)*STATES(18) -  CONSTS(49)*STATES(19))
      RATES(19) = ( ALGBRC(37)*STATES(2)*STATES(18) -  CONSTS(49)*STATES(19)) - ( CONSTS(50)*STATES(19) -  ALGBRC(31)*STATES(2) ** 2.00000*STATES(20))
      RATES(20) = ( CONSTS(50)*STATES(19) -  ALGBRC(31)*STATES(2) ** 2.00000*STATES(20)) - ( CONSTS(49)*STATES(20) -  ALGBRC(37)*STATES(2)*STATES(17))
      ALGBRC(8) = 0.0140000/(1.00000+EXP(- (STATES(1) - 40.0000)/9.00000))
      ALGBRC(25) =  0.00100000*EXP(- STATES(1)/45.0000)
      ALGBRC(34) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(25))
      ALGBRC(40) = 1.00000/(ALGBRC(8)+ALGBRC(25))
      RATES(11) = (ALGBRC(34) - STATES(11))/ALGBRC(40)
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+57.0000)/5.00000))
      ALGBRC(29) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000))
      ALGBRC(35) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000))
      ALGBRC(41) = 1.00000/(ALGBRC(29)+ALGBRC(35))
      RATES(15) = (ALGBRC(12) - STATES(15))/ALGBRC(41)
      ALGBRC(13) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000))
      ALGBRC(30) = 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(36) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(30))
      ALGBRC(42) = 6.65000/(ALGBRC(13)+ALGBRC(30))
      RATES(16) = (ALGBRC(36) - STATES(16))/ALGBRC(42)
      ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+13.5000)/6.00000))
      ALGBRC(17) = TERNRY(STATES(1).EQ.- 35.0000, - 35.0000, TERNRY(STATES(1).EQ.0.00000, 1.00000e-05, STATES(1))
      ALGBRC(33) = ( - 0.0283900*(ALGBRC(17)+35.0000))/(EXP(- (ALGBRC(17)+35.0000)/2.50000) - 1.00000) - ( 0.0849000*ALGBRC(17))/(EXP(- ALGBRC(17)/4.80000) - 1.00000)
      ALGBRC(39) = TERNRY(STATES(1).EQ.5.00000, 5.00001, STATES(1))
      ALGBRC(44) = ( 0.0114300*(ALGBRC(39) - 5.00000))/(EXP((ALGBRC(39) - 5.00000)/2.50000) - 1.00000)
      ALGBRC(46) = 1.00000/(ALGBRC(33)+ALGBRC(44))
      RATES(3) = (ALGBRC(1) - STATES(3))/ALGBRC(46)
      ALGBRC(32) =  CONSTS(1)*CONSTS(13)*(STATES(1) - CONSTS(14))*STATES(6)*STATES(7)
      ALGBRC(15) =  CONSTS(1)*CONSTS(10)*(STATES(1) - CONSTS(9))*STATES(3)*STATES(4)*STATES(5)
      ALGBRC(48) =  CONSTS(1)*0.383300*CONSTS(19)*(STATES(1) - CONSTS(76))*STATES(14) ** 2.00000
      ALGBRC(49) =  CONSTS(1)*0.616700*CONSTS(19)*(STATES(1) - CONSTS(77))*STATES(14) ** 2.00000
      ALGBRC(50) = ALGBRC(48)+ALGBRC(49)
      ALGBRC(51) =  CONSTS(1)*CONSTS(21)*(STATES(1) - CONSTS(22))*STATES(15)*STATES(16)
      ALGBRC(38) =  CONSTS(1)*CONSTS(15)*(STATES(1) - CONSTS(77))*( 0.600000*STATES(9)+ 0.400000*STATES(8))*STATES(10)
      ALGBRC(43) =  CONSTS(1)*CONSTS(16)*(STATES(1) - CONSTS(78))*STATES(11) ** 2.00000
      ALGBRC(45) =  CONSTS(1)*CONSTS(17)*(STATES(1) - CONSTS(77))*STATES(12)*STATES(13)
      ALGBRC(47) =  CONSTS(1)*CONSTS(18)*(STATES(1) - CONSTS(77))*STATES(13)
      ALGBRC(53) = ( CONSTS(1)*CONSTS(26))/( (1.00000+CONSTS(24)/CONSTS(8) ** 1.20000)*(1.00000+CONSTS(25)/CONSTS(5) ** 1.30000)*(1.00000+EXP(- ((STATES(1) - CONSTS(76))+120.000)/30.0000)))
      ALGBRC(55) = EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
      ALGBRC(61) = 1.00000+ (CONSTS(41)/CONSTS(40))*(1.00000+EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))+ (CONSTS(6)/CONSTS(38))*(1.00000+ (CONSTS(6)/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37)))
      ALGBRC(63) = ( (( (CONSTS(6)/CONSTS(38))*CONSTS(6))/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37))*EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(61)
      ALGBRC(62) = ( (CONSTS(41)/CONSTS(40))*EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))/ALGBRC(61)
      ALGBRC(59) = EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
      ALGBRC(64) =  ALGBRC(55)*CONSTS(81)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(62)*ALGBRC(59)*(CONSTS(80)+ALGBRC(55))
      ALGBRC(56) = 1.00000+ (STATES(2)/CONSTS(33))*(1.00000+EXP(( - CONSTS(29)*STATES(1))/CONSTS(79))+CONSTS(5)/CONSTS(36))+ (CONSTS(5)/CONSTS(34))*(1.00000+ (CONSTS(5)/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32)))
      ALGBRC(57) = ( (STATES(2)/CONSTS(33))*EXP(( - CONSTS(29)*STATES(1))/CONSTS(79)))/ALGBRC(56)
      ALGBRC(58) = ( (( (CONSTS(5)/CONSTS(34))*CONSTS(5))/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32))*EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(56)
      ALGBRC(60) =  ALGBRC(59)*CONSTS(80)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(55)*ALGBRC(57)*(CONSTS(81)+ALGBRC(59))
      ALGBRC(65) =  ALGBRC(58)*CONSTS(80)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(57)*ALGBRC(63)*(CONSTS(80)+ALGBRC(55))
      ALGBRC(66) =  ALGBRC(63)*CONSTS(81)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(58)*ALGBRC(62)*(CONSTS(81)+ALGBRC(59))
      ALGBRC(67) = ( CONSTS(1)*CONSTS(28)*( ALGBRC(60)*ALGBRC(62) -  ALGBRC(64)*ALGBRC(57)))/(ALGBRC(64)+ALGBRC(60)+ALGBRC(65)+ALGBRC(66))
      ALGBRC(54) =  CONSTS(1)*CONSTS(27)*(STATES(1) - CONSTS(9))
      ALGBRC(52) =  CONSTS(1)*CONSTS(23)*(STATES(1) - CONSTS(76))
      RATES(1) = - (ALGBRC(15)+ALGBRC(32)+ALGBRC(50)+ALGBRC(51)+ALGBRC(38)+ALGBRC(43)+ALGBRC(45)+ALGBRC(47)+ALGBRC(53)+ALGBRC(67)+ALGBRC(54)+ALGBRC(52))/CONSTS(1)
      ALGBRC(70) =  CONSTS(62)*STATES(2)*(1.00000 - STATES(26)) -  CONSTS(67)*STATES(26)
      RATES(26) = ALGBRC(70)
      ALGBRC(68) =  CONSTS(42)*STATES(18)*(STATES(21) - STATES(2))
      ALGBRC(69) = (STATES(2) - STATES(23))/CONSTS(51)
      RATES(2) = ( ALGBRC(68)*CONSTS(84))/CONSTS(83) - (((ALGBRC(15)+ALGBRC(32)+ALGBRC(54)) -  2.00000*ALGBRC(67))/( 2.00000*CONSTS(4)*CONSTS(83))+ALGBRC(69)+ CONSTS(57)*ALGBRC(70))
      ALGBRC(72) =  CONSTS(59)*STATES(23)*(1.00000 - STATES(27)) -  CONSTS(64)*STATES(27)
      RATES(27) = ALGBRC(72)
      ALGBRC(71) = CONSTS(53)/(1.00000+CONSTS(54)/STATES(23))
      ALGBRC(73) = (STATES(22) - STATES(21))/CONSTS(52)
      RATES(22) = ALGBRC(71) - ( ALGBRC(73)*CONSTS(84))/CONSTS(86)
      ALGBRC(74) =  CONSTS(61)*STATES(23)*(1.00000 - (STATES(28)+STATES(24))) -  CONSTS(65)*STATES(28)
      RATES(28) = ALGBRC(74)
      ALGBRC(75) =  CONSTS(63)*STATES(21)*(1.00000 - STATES(29)) -  CONSTS(68)*STATES(29)
      RATES(29) = ALGBRC(75)
      RATES(21) = ALGBRC(73) - (ALGBRC(68)+ CONSTS(58)*ALGBRC(75))
      ALGBRC(76) =  CONSTS(62)*STATES(23)*(1.00000 - STATES(25)) -  CONSTS(67)*STATES(25)
      RATES(25) = ALGBRC(76)
      RATES(23) = ( ALGBRC(69)*CONSTS(83) -  ALGBRC(71)*CONSTS(86))/CONSTS(85) - ( CONSTS(57)*ALGBRC(76)+ CONSTS(55)*ALGBRC(72)+ CONSTS(56)*ALGBRC(74))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(16) =  CONSTS(60)*CONSTS(69)*(1.00000 - (STATES(28)+STATES(24))) -  CONSTS(66)*STATES(24)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/7.30000))
      ALGBRC(18) = 44.3000+ 257.100*EXP(- (STATES(1)+32.5000)/13.9000 ** 2.00000)
      ALGBRC(3) = CONSTS(12)/(CONSTS(12)+STATES(2))
      ALGBRC(19) = ALGBRC(3)/CONSTS(11)
      ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+26.3000)/6.00000))
      ALGBRC(20) = 1.00000/( 1.06800*EXP((STATES(1)+26.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+26.3000)/30.0000))
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+61.7000)/5.60000))
      ALGBRC(21) = 1.00000/( 0.0153000*EXP(- (STATES(1)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(1)+61.7000)/15.3800))
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1)+23.2000)/10.6000))
      ALGBRC(22) = 0.846554/( 0.00420000*EXP(STATES(1)/17.0000)+ 0.000150000*EXP(- STATES(1)/21.6000))
      ALGBRC(23) = 0.846554/( 0.0372000*EXP(STATES(1)/15.9000)+ 0.000960000*EXP(- STATES(1)/22.5000))
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+28.6000)/17.1000))
      ALGBRC(24) = 1.00000/( 0.100000*EXP(- STATES(1)/54.6450)+ 0.656000*EXP(STATES(1)/106.157))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+49.0000)/13.0000))
      ALGBRC(26) = 6.06000+39.1020/( 0.570000*EXP( - 0.0800000*(STATES(1)+44.0000))+ 0.0650000*EXP( 0.100000*(STATES(1)+45.9300)))
      ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.3000)/15.0000))
      ALGBRC(27) = 2.75352+14.4052/( 1.03700*EXP( 0.0900000*(STATES(1)+30.6100))+ 0.369000*EXP( - 0.120000*(STATES(1)+23.8400)))
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1) - CONSTS(20))/13.5000))
      ALGBRC(28) = 0.716653/(EXP(- (STATES(1)+386.900)/45.3020)+EXP((STATES(1) - 73.0800)/19.2310))
      ALGBRC(14) = CONSTS(43) - (CONSTS(43) - CONSTS(44))/(1.00000+CONSTS(45)/STATES(21) ** CONSTS(46))
      ALGBRC(31) = CONSTS(47)/ALGBRC(14)
      ALGBRC(37) =  CONSTS(48)*ALGBRC(14)
      ALGBRC(8) = 0.0140000/(1.00000+EXP(- (STATES(1) - 40.0000)/9.00000))
      ALGBRC(25) =  0.00100000*EXP(- STATES(1)/45.0000)
      ALGBRC(34) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(25))
      ALGBRC(40) = 1.00000/(ALGBRC(8)+ALGBRC(25))
      ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+57.0000)/5.00000))
      ALGBRC(29) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000))
      ALGBRC(35) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000))
      ALGBRC(41) = 1.00000/(ALGBRC(29)+ALGBRC(35))
      ALGBRC(13) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000))
      ALGBRC(30) = 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(36) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(30))
      ALGBRC(42) = 6.65000/(ALGBRC(13)+ALGBRC(30))
      ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+13.5000)/6.00000))
      ALGBRC(17) = TERNRY(STATES(1).EQ.- 35.0000, - 35.0000, TERNRY(STATES(1).EQ.0.00000, 1.00000e-05, STATES(1))
      ALGBRC(33) = ( - 0.0283900*(ALGBRC(17)+35.0000))/(EXP(- (ALGBRC(17)+35.0000)/2.50000) - 1.00000) - ( 0.0849000*ALGBRC(17))/(EXP(- ALGBRC(17)/4.80000) - 1.00000)
      ALGBRC(39) = TERNRY(STATES(1).EQ.5.00000, 5.00001, STATES(1))
      ALGBRC(44) = ( 0.0114300*(ALGBRC(39) - 5.00000))/(EXP((ALGBRC(39) - 5.00000)/2.50000) - 1.00000)
      ALGBRC(46) = 1.00000/(ALGBRC(33)+ALGBRC(44))
      ALGBRC(32) =  CONSTS(1)*CONSTS(13)*(STATES(1) - CONSTS(14))*STATES(6)*STATES(7)
      ALGBRC(15) =  CONSTS(1)*CONSTS(10)*(STATES(1) - CONSTS(9))*STATES(3)*STATES(4)*STATES(5)
      ALGBRC(48) =  CONSTS(1)*0.383300*CONSTS(19)*(STATES(1) - CONSTS(76))*STATES(14) ** 2.00000
      ALGBRC(49) =  CONSTS(1)*0.616700*CONSTS(19)*(STATES(1) - CONSTS(77))*STATES(14) ** 2.00000
      ALGBRC(50) = ALGBRC(48)+ALGBRC(49)
      ALGBRC(51) =  CONSTS(1)*CONSTS(21)*(STATES(1) - CONSTS(22))*STATES(15)*STATES(16)
      ALGBRC(38) =  CONSTS(1)*CONSTS(15)*(STATES(1) - CONSTS(77))*( 0.600000*STATES(9)+ 0.400000*STATES(8))*STATES(10)
      ALGBRC(43) =  CONSTS(1)*CONSTS(16)*(STATES(1) - CONSTS(78))*STATES(11) ** 2.00000
      ALGBRC(45) =  CONSTS(1)*CONSTS(17)*(STATES(1) - CONSTS(77))*STATES(12)*STATES(13)
      ALGBRC(47) =  CONSTS(1)*CONSTS(18)*(STATES(1) - CONSTS(77))*STATES(13)
      ALGBRC(53) = ( CONSTS(1)*CONSTS(26))/( (1.00000+CONSTS(24)/CONSTS(8) ** 1.20000)*(1.00000+CONSTS(25)/CONSTS(5) ** 1.30000)*(1.00000+EXP(- ((STATES(1) - CONSTS(76))+120.000)/30.0000)))
      ALGBRC(55) = EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
      ALGBRC(61) = 1.00000+ (CONSTS(41)/CONSTS(40))*(1.00000+EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))+ (CONSTS(6)/CONSTS(38))*(1.00000+ (CONSTS(6)/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37)))
      ALGBRC(63) = ( (( (CONSTS(6)/CONSTS(38))*CONSTS(6))/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37))*EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(61)
      ALGBRC(62) = ( (CONSTS(41)/CONSTS(40))*EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))/ALGBRC(61)
      ALGBRC(59) = EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
      ALGBRC(64) =  ALGBRC(55)*CONSTS(81)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(62)*ALGBRC(59)*(CONSTS(80)+ALGBRC(55))
      ALGBRC(56) = 1.00000+ (STATES(2)/CONSTS(33))*(1.00000+EXP(( - CONSTS(29)*STATES(1))/CONSTS(79))+CONSTS(5)/CONSTS(36))+ (CONSTS(5)/CONSTS(34))*(1.00000+ (CONSTS(5)/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32)))
      ALGBRC(57) = ( (STATES(2)/CONSTS(33))*EXP(( - CONSTS(29)*STATES(1))/CONSTS(79)))/ALGBRC(56)
      ALGBRC(58) = ( (( (CONSTS(5)/CONSTS(34))*CONSTS(5))/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32))*EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(56)
      ALGBRC(60) =  ALGBRC(59)*CONSTS(80)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(55)*ALGBRC(57)*(CONSTS(81)+ALGBRC(59))
      ALGBRC(65) =  ALGBRC(58)*CONSTS(80)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(57)*ALGBRC(63)*(CONSTS(80)+ALGBRC(55))
      ALGBRC(66) =  ALGBRC(63)*CONSTS(81)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(58)*ALGBRC(62)*(CONSTS(81)+ALGBRC(59))
      ALGBRC(67) = ( CONSTS(1)*CONSTS(28)*( ALGBRC(60)*ALGBRC(62) -  ALGBRC(64)*ALGBRC(57)))/(ALGBRC(64)+ALGBRC(60)+ALGBRC(65)+ALGBRC(66))
      ALGBRC(54) =  CONSTS(1)*CONSTS(27)*(STATES(1) - CONSTS(9))
      ALGBRC(52) =  CONSTS(1)*CONSTS(23)*(STATES(1) - CONSTS(76))
      ALGBRC(70) =  CONSTS(62)*STATES(2)*(1.00000 - STATES(26)) -  CONSTS(67)*STATES(26)
      ALGBRC(68) =  CONSTS(42)*STATES(18)*(STATES(21) - STATES(2))
      ALGBRC(69) = (STATES(2) - STATES(23))/CONSTS(51)
      ALGBRC(72) =  CONSTS(59)*STATES(23)*(1.00000 - STATES(27)) -  CONSTS(64)*STATES(27)
      ALGBRC(71) = CONSTS(53)/(1.00000+CONSTS(54)/STATES(23))
      ALGBRC(73) = (STATES(22) - STATES(21))/CONSTS(52)
      ALGBRC(74) =  CONSTS(61)*STATES(23)*(1.00000 - (STATES(28)+STATES(24))) -  CONSTS(65)*STATES(28)
      ALGBRC(75) =  CONSTS(63)*STATES(21)*(1.00000 - STATES(29)) -  CONSTS(68)*STATES(29)
      ALGBRC(76) =  CONSTS(62)*STATES(23)*(1.00000 - STATES(25)) -  CONSTS(67)*STATES(25)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END