C
C There are a total of 56 entries in the algebraic variable array.
C There are a total of 22 entries in each of the rate and state variable arrays.
C There are a total of 132 entries in the constant variable array.
C
C
C VOI is time in component Time (time_units).
C CONSTS(1) is T in component Environment (Temperature_units).
C CONSTS(2) is T_exp in component Environment (Temperature_units).
C CONSTS(3) is F in component Environment (F_units).
C CONSTS(4) is R in component Environment (R_units).
C CONSTS(5) is Q10Ca in component Environment (dimensionless).
C CONSTS(6) is Q10K in component Environment (dimensionless).
C CONSTS(7) is Q10Na in component Environment (dimensionless).
C CONSTS(8) is Ca_o in component Environment (millimolar).
C CONSTS(9) is Na_o in component Environment (millimolar).
C CONSTS(10) is K_o in component Environment (millimolar).
C CONSTS(11) is Cl_o in component Environment (millimolar).
C CONSTS(105) is T_correction_Na in component Environment (dimensionless).
C CONSTS(106) is T_correction_K in component Environment (dimensionless).
C CONSTS(107) is T_correction_Ca in component Environment (dimensionless).
C CONSTS(108) is T_correction_BK in component Environment (conductance_units).
C CONSTS(109) is FoRT in component Environment (Inverse_Voltage_units).
C CONSTS(110) is RToF in component Environment (voltage_units).
C CONSTS(12) is Cm in component ICC_Membrane (capacitance_units).
C CONSTS(13) is Vol in component ICC_Membrane (volume_units).
C CONSTS(14) is P_cyto in component ICC_Membrane (dimensionless).
C CONSTS(111) is V_cyto in component ICC_Membrane (volume_units).
C CONSTS(15) is fc in component ICC_Membrane (dimensionless).
C STATES(1) is Vm in component ICC_Membrane (voltage_units).
C STATES(2) is Ca_i in component ICC_Membrane (millimolar).
C CONSTS(16) is Na_i in component ICC_Membrane (millimolar).
C CONSTS(17) is K_i in component ICC_Membrane (millimolar).
C CONSTS(18) is Cl_i in component ICC_Membrane (millimolar).
C ALGBRC(29) is J_leak in component PU_unit (millimolar_per_second).
C ALGBRC(32) is I_Na in component I_Na (current_units).
C ALGBRC(14) is I_Ltype in component I_Ltype (current_units).
C ALGBRC(17) is I_VDDR in component I_VDDR (current_units).
C ALGBRC(30) is I_kv11 in component I_kv11 (current_units).
C ALGBRC(22) is I_BK in component I_BK (current_units).
C ALGBRC(27) is I_ERG in component I_ERG (current_units).
C ALGBRC(18) is I_CaCl in component I_CaCl (current_units).
C ALGBRC(34) is I_NSCC in component I_NSCC (current_units).
C ALGBRC(25) is I_bk in component I_bk (current_units).
C ALGBRC(15) is J_PMCA in component J_PMCA (millimolar_per_second).
C ALGBRC(1) is d_inf_Ltype in component d_Ltype (dimensionless).
C CONSTS(112) is tau_d_Ltype in component d_Ltype (time_units).
C STATES(3) is d_Ltype in component d_Ltype (dimensionless).
C ALGBRC(2) is f_inf_Ltype in component f_Ltype (dimensionless).
C CONSTS(113) is tau_f_Ltype in component f_Ltype (time_units).
C STATES(4) is f_Ltype in component f_Ltype (dimensionless).
C ALGBRC(3) is f_ca_inf_Ltype in component f_ca_Ltype (dimensionless).
C CONSTS(114) is tau_f_ca_Ltype in component f_ca_Ltype (time_units).
C STATES(5) is f_ca_Ltype in component f_ca_Ltype (dimensionless).
C ALGBRC(13) is E_Ca in component I_Ltype (voltage_units).
C CONSTS(19) is G_max_Ltype in component I_Ltype (conductance_units).
C CONSTS(20) is J_max_PMCA in component J_PMCA (millimolar_per_second).
C ALGBRC(4) is d_inf_VDDR in component d_VDDR (dimensionless).
C CONSTS(115) is tau_d_VDDR in component d_VDDR (time_units).
C STATES(6) is d_VDDR in component d_VDDR (dimensionless).
C ALGBRC(5) is f_inf_VDDR in component f_VDDR (dimensionless).
C CONSTS(116) is tau_f_VDDR in component f_VDDR (time_units).
C STATES(7) is f_VDDR in component f_VDDR (dimensionless).
C ALGBRC(16) is E_Ca in component I_VDDR (voltage_units).
C CONSTS(21) is G_max_VDDR in component I_VDDR (conductance_units).
C ALGBRC(6) is d_inf_CaCl in component d_CaCl (dimensionless).
C CONSTS(22) is tau_d_CaCl in component d_CaCl (time_units).
C STATES(8) is d_CaCl in component d_CaCl (dimensionless).
C CONSTS(117) is E_Cl in component I_CaCl (voltage_units).
C CONSTS(23) is G_max_CaCl in component I_CaCl (conductance_units).
C ALGBRC(20) is d_BK in component d_BK (dimensionless).
C CONSTS(118) is E_K in component I_BK (voltage_units).
C CONSTS(24) is G_max_BK in component I_BK (conductance_units).
C CONSTS(119) is E_K in component I_bk (voltage_units).
C CONSTS(25) is G_max_bk in component I_bk (conductance_units).
C ALGBRC(7) is d_inf_ERG in component d_ERG (dimensionless).
C CONSTS(120) is tau_d_ERG in component d_ERG (time_units).
C STATES(9) is d_ERG in component d_ERG (dimensionless).
C CONSTS(121) is E_K in component I_ERG (voltage_units).
C CONSTS(26) is G_max_ERG in component I_ERG (conductance_units).
C ALGBRC(8) is d_inf_kv11 in component d_kv11 (dimensionless).
C CONSTS(122) is tau_d_kv11 in component d_kv11 (time_units).
C STATES(10) is d_kv11 in component d_kv11 (dimensionless).
C ALGBRC(9) is f_inf_kv11 in component f_kv11 (dimensionless).
C CONSTS(123) is tau_f_kv11 in component f_kv11 (time_units).
C STATES(11) is f_kv11 in component f_kv11 (dimensionless).
C CONSTS(124) is E_K in component I_kv11 (voltage_units).
C CONSTS(27) is G_max_kv11 in component I_kv11 (conductance_units).
C ALGBRC(10) is d_inf_Na in component d_Na (dimensionless).
C CONSTS(125) is tau_d_Na in component d_Na (time_units).
C STATES(12) is d_Na in component d_Na (dimensionless).
C ALGBRC(11) is f_inf_Na in component f_Na (dimensionless).
C CONSTS(126) is tau_f_Na in component f_Na (time_units).
C STATES(13) is f_Na in component f_Na (dimensionless).
C CONSTS(127) is E_Na in component I_Na (voltage_units).
C CONSTS(28) is G_max_Na in component I_Na (conductance_units).
C STATES(14) is Ca_PU in component PU_unit (millimolar).
C ALGBRC(12) is d_inf_NSCC in component d_NSCC (dimensionless).
C CONSTS(29) is tau_d_NSCC in component d_NSCC (time_units).
C STATES(15) is d_NSCC in component d_NSCC (dimensionless).
C CONSTS(128) is E_NSCC in component I_NSCC (voltage_units).
C CONSTS(30) is G_max_NSCC in component I_NSCC (conductance_units).
C CONSTS(31) is NaPerm_o_Kperm in component I_NSCC (dimensionless).
C CONSTS(32) is P_PU in component PU_unit (dimensionless).
C CONSTS(33) is P_mito in component PU_unit (dimensionless).
C CONSTS(34) is P_ER in component PU_unit (dimensionless).
C CONSTS(129) is V_MITO in component PU_unit (volume_units).
C CONSTS(131) is V_ER in component PU_unit (volume_units).
C CONSTS(130) is V_PU in component PU_unit (volume_units).
C CONSTS(35) is fe in component PU_unit (dimensionless).
C CONSTS(36) is fm in component PU_unit (dimensionless).
C STATES(16) is Ca_m in component PU_unit (millimolar).
C STATES(17) is Ca_ER in component PU_unit (millimolar).
C STATES(18) is ADP_m in component PU_unit (millimolar).
C STATES(19) is ADP_i in component PU_unit (millimolar).
C STATES(20) is NADH_m in component PU_unit (millimolar).
C STATES(21) is h in component PU_unit (dimensionless).
C CONSTS(37) is IP3 in component PU_unit (millimolar).
C STATES(22) is deltaPsi in component PU_unit (voltage_units).
C CONSTS(38) is deltapH in component PU_unit (dimensionless).
C CONSTS(39) is Cmito in component PU_unit (millifarads).
C CONSTS(40) is K_res in component PU_unit (dimensionless).
C CONSTS(41) is r1 in component PU_unit (dimensionless).
C CONSTS(42) is r2 in component PU_unit (dimensionless).
C CONSTS(43) is r3 in component PU_unit (dimensionless).
C CONSTS(44) is ra in component PU_unit (rate_constants_units).
C CONSTS(45) is rb in component PU_unit (rate_constants_units).
C CONSTS(46) is rc1 in component PU_unit (rate_constants_units).
C CONSTS(47) is rc2 in component PU_unit (rate_constants_units).
C CONSTS(48) is deltaPsi_B in component PU_unit (voltage_units).
C CONSTS(49) is g in component PU_unit (dimensionless).
C CONSTS(50) is K_F1 in component PU_unit (millimolar).
C CONSTS(51) is Pi_m in component PU_unit (millimolar).
C CONSTS(52) is p1 in component PU_unit (dimensionless).
C CONSTS(53) is p2 in component PU_unit (dimensionless).
C CONSTS(54) is p3 in component PU_unit (dimensionless).
C CONSTS(55) is pa in component PU_unit (rate_constants_units).
C CONSTS(56) is pb in component PU_unit (rate_constants_units).
C CONSTS(57) is pc1 in component PU_unit (rate_constants_units).
C CONSTS(58) is pc2 in component PU_unit (rate_constants_units).
C CONSTS(59) is frac in component PU_unit (dimensionless).
C CONSTS(60) is K_act in component PU_unit (millimolar).
C CONSTS(61) is na in component PU_unit (dimensionless).
C CONSTS(62) is deltaPsi_star in component PU_unit (voltage_units).
C CONSTS(63) is K_Na in component PU_unit (millimolar).
C CONSTS(64) is K_Ca in component PU_unit (millimolar).
C CONSTS(65) is K_trans in component PU_unit (millimolar).
C CONSTS(66) is L in component PU_unit (dimensionless).
C CONSTS(67) is b in component PU_unit (dimensionless).
C CONSTS(68) is beta_max in component PU_unit (rate_constants_units).
C CONSTS(69) is beta1 in component PU_unit (per_millimolar).
C CONSTS(70) is beta2 in component PU_unit (per_millimolar).
C CONSTS(71) is beta3 in component PU_unit (per_millimolar).
C CONSTS(72) is beta4 in component PU_unit (per_millimolar).
C CONSTS(73) is beta5 in component PU_unit (per_millimolar).
C CONSTS(74) is beta6 in component PU_unit (per_millimolar).
C CONSTS(75) is beta7 in component PU_unit (per_millimolar).
C CONSTS(76) is KCa_PDH in component PU_unit (millimolar).
C CONSTS(77) is u1 in component PU_unit (dimensionless).
C CONSTS(78) is u2 in component PU_unit (dimensionless).
C CONSTS(79) is n in component PU_unit (dimensionless).
C CONSTS(80) is K_Glc in component PU_unit (millimolar).
C CONSTS(81) is nhyd in component PU_unit (dimensionless).
C CONSTS(82) is K_hyd in component PU_unit (rate_constants_units).
C CONSTS(83) is J_ERleak in component PU_unit (rate_constants_units).
C CONSTS(84) is Jmax_IP3 in component PU_unit (rate_constants_units).
C CONSTS(85) is d_IP3 in component PU_unit (millimolar).
C CONSTS(86) is d_ACT in component PU_unit (millimolar).
C CONSTS(87) is d_INH in component PU_unit (millimolar).
C CONSTS(88) is tauh in component PU_unit (time_units).
C CONSTS(89) is Jmax_serca in component PU_unit (millimolar_per_second).
C CONSTS(90) is k_serca in component PU_unit (millimolar).
C CONSTS(91) is conc in component PU_unit (millimolar).
C CONSTS(92) is Jmax_uni in component PU_unit (rate_constants_units).
C CONSTS(93) is Jmax_NaCa in component PU_unit (millimolar_per_second).
C CONSTS(94) is J_max_leak in component PU_unit (rate_constants_units).
C CONSTS(95) is rho_res in component PU_unit (millimolar).
C CONSTS(96) is rho_F1 in component PU_unit (millimolar).
C CONSTS(97) is g_H in component PU_unit (millimolar_per_second_per_millivolt).
C CONSTS(98) is J_red_basal in component PU_unit (millimolar_per_second).
C CONSTS(99) is Jmax_ANT in component PU_unit (millimolar_per_second).
C CONSTS(100) is J_hyd_max in component PU_unit (millimolar_per_second).
C CONSTS(101) is Glc in component PU_unit (millimolar).
C CONSTS(102) is total_NAD_m in component PU_unit (millimolar).
C CONSTS(103) is total_ANP_m in component PU_unit (millimolar).
C CONSTS(104) is total_ANP_i in component PU_unit (millimolar).
C ALGBRC(35) is NAD_m in component PU_unit (millimolar).
C ALGBRC(39) is ATP_m in component PU_unit (millimolar).
C ALGBRC(40) is ADP_mfree in component PU_unit (millimolar).
C ALGBRC(44) is ADP3_m in component PU_unit (millimolar).
C ALGBRC(45) is ATP4_m in component PU_unit (millimolar).
C ALGBRC(46) is ATP_i in component PU_unit (millimolar).
C ALGBRC(52) is ADP_ifree in component PU_unit (millimolar).
C ALGBRC(54) is ADP3_i in component PU_unit (millimolar).
C ALGBRC(53) is MgADP_i in component PU_unit (millimolar).
C ALGBRC(55) is ATP4_i in component PU_unit (millimolar).
C ALGBRC(19) is J_ERout in component PU_unit (millimolar_per_second).
C ALGBRC(21) is J_SERCA in component PU_unit (millimolar_per_second).
C ALGBRC(23) is MWC in component PU_unit (millimolar).
C ALGBRC(24) is J_uni in component PU_unit (millimolar_per_second).
C ALGBRC(26) is J_NaCa in component PU_unit (millimolar_per_second).
C ALGBRC(49) is J_red in component PU_unit (millimolar_per_second).
C ALGBRC(48) is J_pTCA in component PU_unit (millimolar_per_second).
C ALGBRC(41) is A_F1 in component PU_unit (voltage_units).
C ALGBRC(42) is J_pF1 in component PU_unit (millimolar_per_second).
C ALGBRC(43) is J_HF1 in component PU_unit (millimolar_per_second).
C ALGBRC(36) is A_res in component PU_unit (voltage_units).
C ALGBRC(37) is J_o in component PU_unit (millimolar_per_second).
C ALGBRC(47) is J_glyTotal in component PU_unit (millimolar_per_second).
C ALGBRC(28) is f_PDHa in component PU_unit (dimensionless).
C ALGBRC(38) is J_Hres in component PU_unit (millimolar_per_second).
C ALGBRC(56) is J_ANT in component PU_unit (millimolar_per_second).
C ALGBRC(31) is PMF in component PU_unit (voltage_units).
C ALGBRC(33) is J_Hleak in component PU_unit (millimolar_per_second).
C ALGBRC(50) is J_pGly in component PU_unit (millimolar_per_second).
C CONSTS(132) is J_hydSS in component PU_unit (millimolar_per_second).
C ALGBRC(51) is J_hyd in component PU_unit (millimolar_per_second).
C RATES(1) is d/dt Vm in component ICC_Membrane (voltage_units).
C RATES(2) is d/dt Ca_i in component ICC_Membrane (millimolar).
C RATES(3) is d/dt d_Ltype in component d_Ltype (dimensionless).
C RATES(4) is d/dt f_Ltype in component f_Ltype (dimensionless).
C RATES(5) is d/dt f_ca_Ltype in component f_ca_Ltype (dimensionless).
C RATES(6) is d/dt d_VDDR in component d_VDDR (dimensionless).
C RATES(7) is d/dt f_VDDR in component f_VDDR (dimensionless).
C RATES(8) is d/dt d_CaCl in component d_CaCl (dimensionless).
C RATES(9) is d/dt d_ERG in component d_ERG (dimensionless).
C RATES(10) is d/dt d_kv11 in component d_kv11 (dimensionless).
C RATES(11) is d/dt f_kv11 in component f_kv11 (dimensionless).
C RATES(12) is d/dt d_Na in component d_Na (dimensionless).
C RATES(13) is d/dt f_Na in component f_Na (dimensionless).
C RATES(15) is d/dt d_NSCC in component d_NSCC (dimensionless).
C RATES(20) is d/dt NADH_m in component PU_unit (millimolar).
C RATES(18) is d/dt ADP_m in component PU_unit (millimolar).
C RATES(19) is d/dt ADP_i in component PU_unit (millimolar).
C RATES(14) is d/dt Ca_PU in component PU_unit (millimolar).
C RATES(16) is d/dt Ca_m in component PU_unit (millimolar).
C RATES(17) is d/dt Ca_ER in component PU_unit (millimolar).
C RATES(22) is d/dt deltaPsi in component PU_unit (voltage_units).
C RATES(21) is d/dt h in component PU_unit (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 310
      CONSTS(2) = 297
      CONSTS(3) = 96.4846
      CONSTS(4) = 8.3144
      CONSTS(5) = 2.1
      CONSTS(6) = 1.5
      CONSTS(7) = 2.45
      CONSTS(8) = 2.5
      CONSTS(9) = 137
      CONSTS(10) = 7
      CONSTS(11) = 134
      CONSTS(12) = 0.025
      CONSTS(13) = 1e-12
      CONSTS(14) = 0.7
      CONSTS(15) = 0.01
      STATES(1) = -67
      STATES(2) = 0.00000993087
      CONSTS(16) = 30
      CONSTS(17) = 120
      CONSTS(18) = 88
      STATES(3) = 0
      STATES(4) = 1
      STATES(5) = 1
      CONSTS(19) = 2
      CONSTS(20) = 0.088464
      STATES(6) = 0
      STATES(7) = 1
      CONSTS(21) = 3
      CONSTS(22) = 0.03
      STATES(8) = 0
      CONSTS(23) = 10.1
      CONSTS(24) = 23
      CONSTS(25) = 0.15
      STATES(9) = 0
      CONSTS(26) = 2.5
      STATES(10) = 0
      STATES(11) = 1
      CONSTS(27) = 6.3
      STATES(12) = 0
      STATES(13) = 1
      CONSTS(28) = 20
      STATES(14) = 0.0000902
      CONSTS(29) = 0.35
      STATES(15) = 0
      CONSTS(30) = 12.15
      CONSTS(31) = 1.056075
      CONSTS(32) = 0.001
      CONSTS(33) = 0.12871
      CONSTS(34) = 0.1
      CONSTS(35) = 0.01
      CONSTS(36) = 0.0003
      STATES(16) = 0.000136
      STATES(17) = 0.007299
      STATES(18) = 2.60093454
      STATES(19) = 0.0077282
      STATES(20) = 0.101476
      STATES(21) = 0.9397
      CONSTS(37) = 0.0006
      STATES(22) = 164.000044
      CONSTS(38) = -0.4
      CONSTS(39) = 0.006995
      CONSTS(40) = 1.35e18
      CONSTS(41) = 2.077e-18
      CONSTS(42) = 1.728e-9
      CONSTS(43) = 1.059e-26
      CONSTS(44) = 6.394e-10
      CONSTS(45) = 1.762e-13
      CONSTS(46) = 2.656e-19
      CONSTS(47) = 8.632e-27
      CONSTS(48) = 50
      CONSTS(49) = 0.85
      CONSTS(50) = 1.71e9
      CONSTS(51) = 20
      CONSTS(52) = 1.346e-8
      CONSTS(53) = 7.739e-7
      CONSTS(54) = 6.65e-15
      CONSTS(55) = 1.656e-5
      CONSTS(56) = 3.373e-7
      CONSTS(57) = 9.651e-14
      CONSTS(58) = 4.845e-19
      CONSTS(59) = 0.5
      CONSTS(60) = 0.00038
      CONSTS(61) = 2.8
      CONSTS(62) = 91
      CONSTS(63) = 9.4
      CONSTS(64) = 0.003
      CONSTS(65) = 0.006
      CONSTS(66) = 50
      CONSTS(67) = 0.5
      CONSTS(68) = 2.055
      CONSTS(69) = 1.66
      CONSTS(70) = 0.0249
      CONSTS(71) = 4
      CONSTS(72) = 2.83
      CONSTS(73) = 1.3
      CONSTS(74) = 2.66
      CONSTS(75) = 0.16
      CONSTS(76) = 0.00005
      CONSTS(77) = 15
      CONSTS(78) = 1.1
      CONSTS(79) = 2
      CONSTS(80) = 8.7
      CONSTS(81) = 2.7
      CONSTS(82) = 0.05125
      CONSTS(83) = 1.666667
      CONSTS(84) = 50000
      CONSTS(85) = 0.00025
      CONSTS(86) = 0.001
      CONSTS(87) = 0.0014
      CONSTS(88) = 4
      CONSTS(89) = 1.8333
      CONSTS(90) = 0.00042
      CONSTS(91) = 0.001
      CONSTS(92) = 5000
      CONSTS(93) = 0.05
      CONSTS(94) = 0.01
      CONSTS(95) = 0.4
      CONSTS(96) = 0.7
      CONSTS(97) = 0.0033333
      CONSTS(98) = 0.3333
      CONSTS(99) = 15
      CONSTS(100) = 0.037625
      CONSTS(101) = 1
      CONSTS(102) = 8
      CONSTS(103) = 12
      CONSTS(104) = 2
      CONSTS(105) = CONSTS(7) ** (CONSTS(1) - CONSTS(2))/10.0000
      CONSTS(106) = CONSTS(6) ** (CONSTS(1) - CONSTS(2))/10.0000
      CONSTS(107) = CONSTS(5) ** (CONSTS(1) - CONSTS(2))/10.0000
      CONSTS(108) =  1.10000*(CONSTS(1) - CONSTS(2))
      CONSTS(109) = CONSTS(3)/( CONSTS(4)*CONSTS(1))
      CONSTS(110) = ( CONSTS(4)*CONSTS(1))/CONSTS(3)
      CONSTS(111) =  CONSTS(13)*CONSTS(14)
      CONSTS(112) =  CONSTS(107)*0.00100000
      CONSTS(113) =  CONSTS(107)*0.0860000
      CONSTS(114) =  CONSTS(107)*0.00200000
      CONSTS(115) =  CONSTS(107)*0.00600000
      CONSTS(116) =  CONSTS(107)*0.0400000
      CONSTS(117) =  CONSTS(110)*log(CONSTS(18)/CONSTS(11))
      CONSTS(118) =  CONSTS(110)*log(CONSTS(10)/CONSTS(17))
      CONSTS(119) =  CONSTS(110)*log(CONSTS(10)/CONSTS(17))
      CONSTS(120) =  CONSTS(106)*0.00300000
      CONSTS(121) =  CONSTS(110)*log(CONSTS(10)/CONSTS(17))
      CONSTS(122) =  CONSTS(106)*0.00500000
      CONSTS(123) =  CONSTS(106)*0.00500000
      CONSTS(124) =  CONSTS(110)*log(CONSTS(10)/CONSTS(17))
      CONSTS(125) =  CONSTS(105)*0.00300000
      CONSTS(126) =  CONSTS(105)*0.00160000
      CONSTS(127) =  CONSTS(110)*log(CONSTS(9)/CONSTS(16))
      CONSTS(128) =  CONSTS(110)*log((CONSTS(10)+ CONSTS(9)*CONSTS(31))/(CONSTS(17)+ CONSTS(16)*CONSTS(31)))
      CONSTS(129) =  CONSTS(13)*CONSTS(33)
      CONSTS(130) =  CONSTS(13)*CONSTS(32)
      CONSTS(131) =  CONSTS(13)*CONSTS(34)
      CONSTS(132) = CONSTS(100)/(1.00000+CONSTS(80)/CONSTS(101) ** CONSTS(81))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(21) = ( 1.00000*(CONSTS(87) -  STATES(21)*(STATES(14)+CONSTS(87))))/CONSTS(88)
      ALGBRC(1) = 1.00000/(1.00000+EXP((STATES(1)+17.0000)/- 4.30000))
      RATES(3) = (ALGBRC(1) - STATES(3))/CONSTS(112)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+43.0000)/8.90000))
      RATES(4) = (ALGBRC(2) - STATES(4))/CONSTS(113)
      ALGBRC(3) = 1.00000 - 1.00000/(1.00000+EXP(((STATES(2) - 0.000100000) - 0.000214000)/- 1.31000e-05))
      RATES(5) = (ALGBRC(3) - STATES(5))/CONSTS(114)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+26.0000)/- 6.00000))
      RATES(6) = (ALGBRC(4) - STATES(6))/CONSTS(115)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+66.0000)/6.00000))
      RATES(7) = (ALGBRC(5) - STATES(7))/CONSTS(116)
      ALGBRC(6) = 1.00000/(1.00000+0.000140000/STATES(2) ** 3.00000)
      RATES(8) = (ALGBRC(6) - STATES(8))/CONSTS(22)
      ALGBRC(7) = 0.200000+0.800000/(1.00000+EXP((STATES(1)+20.0000)/- 1.80000))
      RATES(9) = (ALGBRC(7) - STATES(9))/CONSTS(120)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+25.0000)/- 7.70000))
      RATES(10) = (ALGBRC(8) - STATES(10))/CONSTS(122)
      ALGBRC(9) = 0.500000+0.500000/(1.00000+EXP((STATES(1)+44.8000)/4.40000))
      RATES(11) = (ALGBRC(9) - STATES(11))/CONSTS(123)
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+47.0000)/- 4.80000))
      RATES(12) = (ALGBRC(10) - STATES(12))/CONSTS(125)
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+78.0000)/7.00000))
      RATES(13) = (ALGBRC(11) - STATES(13))/CONSTS(126)
      ALGBRC(12) = 1.00000/(1.00000+7.45000e-05/STATES(14) ** - 85.0000)
      RATES(15) = (ALGBRC(12) - STATES(15))/CONSTS(29)
      ALGBRC(19) =  ( CONSTS(84)*CONSTS(37)/(CONSTS(37)+CONSTS(85)) ** 3.00000*STATES(14)/(STATES(14)+CONSTS(86)) ** 3.00000*STATES(21) ** 3.00000+CONSTS(83))*(STATES(17) - STATES(14))
      ALGBRC(21) = ( CONSTS(89)*STATES(14) ** 2.00000)/(CONSTS(90) ** 2.00000+STATES(14) ** 2.00000)
      RATES(17) =  CONSTS(35)*(ALGBRC(21) - ALGBRC(19))
      ALGBRC(23) = ( (( CONSTS(91)*STATES(14))/CONSTS(65))*1.00000+STATES(14)/CONSTS(65) ** 3.00000)/(1.00000+STATES(14)/CONSTS(65) ** 4.00000+CONSTS(66)/1.00000+STATES(14)/CONSTS(60) ** CONSTS(61))
      ALGBRC(24) = ( CONSTS(92)*(ALGBRC(23) -  STATES(16)*EXP( - 2.00000*CONSTS(109)*(STATES(22) - CONSTS(62))))*2.00000*CONSTS(109)*(STATES(22) - CONSTS(62)))/(1.00000 - EXP( - 2.00000*CONSTS(109)*(STATES(22) - CONSTS(62))))
      ALGBRC(26) = ( CONSTS(93)*EXP( CONSTS(67)*CONSTS(109)*(STATES(22) - CONSTS(62))))/( (1.00000+CONSTS(63)/CONSTS(16) ** CONSTS(79))*(1.00000+CONSTS(64)/STATES(16)))
      RATES(16) =  CONSTS(36)*(ALGBRC(24) - ALGBRC(26))
      ALGBRC(29) =  CONSTS(94)*(STATES(14) - STATES(2))
      ALGBRC(13) =  0.500000*CONSTS(110)*log(CONSTS(8)/STATES(2))
      ALGBRC(14) =  CONSTS(19)*STATES(4)*STATES(3)*STATES(5)*(STATES(1) - ALGBRC(13))
      ALGBRC(16) =  0.500000*CONSTS(110)*log(CONSTS(8)/STATES(2))
      ALGBRC(17) =  CONSTS(21)*STATES(7)*STATES(6)*(STATES(1) - ALGBRC(16))
      ALGBRC(15) = ( CONSTS(20)*1.00000)/(1.00000+0.000298000/STATES(2))
      RATES(2) =  CONSTS(15)*(( - 1.00000*ALGBRC(14)+ - 1.00000*ALGBRC(17))/( 2.00000*1.00000e+06*1.00000e+06*CONSTS(3)*CONSTS(111))+ALGBRC(29)+ - 1.00000*ALGBRC(15))
      RATES(14) =  CONSTS(15)*(( (ALGBRC(26) - ALGBRC(24))*CONSTS(129))/CONSTS(130)+( (ALGBRC(19) - ALGBRC(21))*CONSTS(131))/CONSTS(130)+( - 1.00000*ALGBRC(29)*CONSTS(111))/CONSTS(130))
      ALGBRC(32) =  CONSTS(28)*STATES(13)*STATES(12)*(STATES(1) - CONSTS(127))
      ALGBRC(30) =  CONSTS(27)*STATES(11)*STATES(10)*(STATES(1) - CONSTS(124))
      ALGBRC(20) = 1.00000/(1.00000+EXP(STATES(1)/- 17.0000 -  2.00000*log(STATES(2)/0.00100000)))
      ALGBRC(22) =  (CONSTS(24)+CONSTS(108))*ALGBRC(20)*(STATES(1) - CONSTS(118))
      ALGBRC(27) =  CONSTS(26)*STATES(9)*(STATES(1) - CONSTS(121))
      ALGBRC(18) =  CONSTS(23)*STATES(8)*(STATES(1) - CONSTS(117))
      ALGBRC(34) =  CONSTS(30)*STATES(15)*(STATES(1) - CONSTS(128))
      ALGBRC(25) =  CONSTS(25)*(STATES(1) - CONSTS(119))
      RATES(1) =  (( - 1.00000*1.00000)/CONSTS(12))*(ALGBRC(32)+ALGBRC(14)+ALGBRC(17)+ALGBRC(30)+ALGBRC(27)+ALGBRC(22)+ALGBRC(18)+ALGBRC(34)+ALGBRC(25)+ ALGBRC(15)*2.00000*1.00000e+06*1.00000e+06*CONSTS(3)*CONSTS(111))
      ALGBRC(46) = CONSTS(104) - STATES(19)
      ALGBRC(47) = ( CONSTS(68)*(1.00000+ CONSTS(69)*CONSTS(101))*CONSTS(70)*CONSTS(101)*ALGBRC(46))/(1.00000+ CONSTS(71)*ALGBRC(46)+ (1.00000+ CONSTS(72)*ALGBRC(46))*CONSTS(73)*CONSTS(101)+ (1.00000+ CONSTS(74)*ALGBRC(46))*CONSTS(75)*CONSTS(101))
      ALGBRC(28) = 1.00000/(1.00000+ CONSTS(78)*(1.00000+CONSTS(77)/1.00000+STATES(16)/CONSTS(76) ** 2.00000))
      ALGBRC(49) = CONSTS(98)+ 6.39440*ALGBRC(28)*ALGBRC(47)
      ALGBRC(35) = CONSTS(102) - STATES(20)
      ALGBRC(36) =  CONSTS(110)*log(( CONSTS(40)* STATES(20) ** (1.0 / 2))/ ALGBRC(35) ** (1.0 / 2))
      ALGBRC(37) = ( CONSTS(95)*0.500000*( ( CONSTS(44)*10.0000 **  6.00000*CONSTS(38)+ CONSTS(46)*EXP( 6.00000*CONSTS(48)*CONSTS(109)))*EXP( ALGBRC(36)*CONSTS(109))+ - 1.00000*CONSTS(44)*EXP( CONSTS(49)*6.00000*CONSTS(109)*STATES(22))+ CONSTS(47)*EXP( CONSTS(109)*ALGBRC(36))*EXP( CONSTS(109)*STATES(22)*6.00000*CONSTS(49))))/( (1.00000+ CONSTS(41)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( CONSTS(109)*CONSTS(48)*6.00000)+ (CONSTS(42)+ CONSTS(43)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( CONSTS(109)*STATES(22)*CONSTS(49)*6.00000))
      RATES(20) = ALGBRC(49) - ALGBRC(37)
      ALGBRC(48) = CONSTS(98)/3.00000+ 0.840000*ALGBRC(28)*ALGBRC(47)
      ALGBRC(39) = CONSTS(103) - STATES(18)
      ALGBRC(40) =  0.800000*STATES(18)
      ALGBRC(41) =  CONSTS(110)*log(( CONSTS(50)*ALGBRC(39))/( ALGBRC(40)*CONSTS(51)))
      ALGBRC(42) = ( CONSTS(96)*- 1.00000*( ( CONSTS(55)*10.0000 **  3.00000*CONSTS(38)+ CONSTS(57)*EXP( 3.00000*CONSTS(109)*CONSTS(48)))*EXP( CONSTS(109)*ALGBRC(41))+ - 1.00000*CONSTS(55)*EXP( 3.00000*CONSTS(109)*STATES(22))+ CONSTS(58)*EXP( CONSTS(109)*ALGBRC(41))*EXP( 3.00000*CONSTS(109)*STATES(22))))/( (1.00000+ CONSTS(52)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*CONSTS(48))+ (CONSTS(53)+ CONSTS(54)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*STATES(22)))
      ALGBRC(44) =  0.450000*ALGBRC(40)
      ALGBRC(45) =  0.0500000*ALGBRC(39)
      ALGBRC(52) =  0.300000*STATES(19)
      ALGBRC(54) =  0.450000*ALGBRC(52)
      ALGBRC(55) =  0.0500000*ALGBRC(46)
      ALGBRC(56) = ( CONSTS(99)*(1.00000 -  (( ALGBRC(55)*ALGBRC(44))/( ALGBRC(54)*ALGBRC(45)))*EXP( - 1.00000*CONSTS(109)*STATES(22))))/( (1.00000+ (ALGBRC(55)/ALGBRC(54))*EXP( - 1.00000*CONSTS(59)*CONSTS(109)*STATES(22)))*(1.00000+ALGBRC(44)/ALGBRC(45)))
      RATES(18) = ALGBRC(56)+ - 1.00000*ALGBRC(48)+ - 1.00000*ALGBRC(42)
      ALGBRC(50) =  0.150000*ALGBRC(47)
      ALGBRC(51) =  CONSTS(82)*ALGBRC(46)+CONSTS(132)
      RATES(19) = ( - 1.00000*ALGBRC(56)*CONSTS(129))/CONSTS(111)+ALGBRC(51)+ - 1.00000*ALGBRC(50)
      ALGBRC(43) = ( - 1.00000*CONSTS(96)*3.00000*( CONSTS(55)*10.0000 **  3.00000*CONSTS(38)*EXP( CONSTS(109)*ALGBRC(41))+ CONSTS(56)*10.0000 **  3.00000*CONSTS(38)+ - 1.00000*(CONSTS(55)+CONSTS(56))*EXP( 3.00000*CONSTS(109)*STATES(22))))/( (1.00000+ CONSTS(52)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*CONSTS(48))+ (CONSTS(53)+ CONSTS(54)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*STATES(22)))
      ALGBRC(38) = ( CONSTS(95)*3.96600*( CONSTS(44)*10.0000 **  6.00000*CONSTS(38)*EXP( CONSTS(109)*ALGBRC(36))+ CONSTS(45)*10.0000 **  6.00000*CONSTS(38)+ - 1.00000*(CONSTS(44)+CONSTS(45))*EXP( CONSTS(49)*CONSTS(109)*STATES(22)*6.00000)))/( (1.00000+ CONSTS(41)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( 6.00000*CONSTS(109)*CONSTS(48))+ (CONSTS(42)+ CONSTS(43)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( CONSTS(49)*6.00000*CONSTS(109)*STATES(22)))
      ALGBRC(31) = STATES(22) -  2.30300*CONSTS(110)*CONSTS(38)
      ALGBRC(33) =  CONSTS(97)*ALGBRC(31)
      RATES(22) =  (( - 1.00000*CONSTS(3)*CONSTS(129)*1.00000e+06*1.00000)/CONSTS(39))*(ALGBRC(33)+ - 1.00000*ALGBRC(38)+ALGBRC(56)+ALGBRC(43)+ 2.00000*ALGBRC(24))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = 1.00000/(1.00000+EXP((STATES(1)+17.0000)/- 4.30000))
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+43.0000)/8.90000))
      ALGBRC(3) = 1.00000 - 1.00000/(1.00000+EXP(((STATES(2) - 0.000100000) - 0.000214000)/- 1.31000e-05))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+26.0000)/- 6.00000))
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+66.0000)/6.00000))
      ALGBRC(6) = 1.00000/(1.00000+0.000140000/STATES(2) ** 3.00000)
      ALGBRC(7) = 0.200000+0.800000/(1.00000+EXP((STATES(1)+20.0000)/- 1.80000))
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+25.0000)/- 7.70000))
      ALGBRC(9) = 0.500000+0.500000/(1.00000+EXP((STATES(1)+44.8000)/4.40000))
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+47.0000)/- 4.80000))
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+78.0000)/7.00000))
      ALGBRC(12) = 1.00000/(1.00000+7.45000e-05/STATES(14) ** - 85.0000)
      ALGBRC(19) =  ( CONSTS(84)*CONSTS(37)/(CONSTS(37)+CONSTS(85)) ** 3.00000*STATES(14)/(STATES(14)+CONSTS(86)) ** 3.00000*STATES(21) ** 3.00000+CONSTS(83))*(STATES(17) - STATES(14))
      ALGBRC(21) = ( CONSTS(89)*STATES(14) ** 2.00000)/(CONSTS(90) ** 2.00000+STATES(14) ** 2.00000)
      ALGBRC(23) = ( (( CONSTS(91)*STATES(14))/CONSTS(65))*1.00000+STATES(14)/CONSTS(65) ** 3.00000)/(1.00000+STATES(14)/CONSTS(65) ** 4.00000+CONSTS(66)/1.00000+STATES(14)/CONSTS(60) ** CONSTS(61))
      ALGBRC(24) = ( CONSTS(92)*(ALGBRC(23) -  STATES(16)*EXP( - 2.00000*CONSTS(109)*(STATES(22) - CONSTS(62))))*2.00000*CONSTS(109)*(STATES(22) - CONSTS(62)))/(1.00000 - EXP( - 2.00000*CONSTS(109)*(STATES(22) - CONSTS(62))))
      ALGBRC(26) = ( CONSTS(93)*EXP( CONSTS(67)*CONSTS(109)*(STATES(22) - CONSTS(62))))/( (1.00000+CONSTS(63)/CONSTS(16) ** CONSTS(79))*(1.00000+CONSTS(64)/STATES(16)))
      ALGBRC(29) =  CONSTS(94)*(STATES(14) - STATES(2))
      ALGBRC(13) =  0.500000*CONSTS(110)*log(CONSTS(8)/STATES(2))
      ALGBRC(14) =  CONSTS(19)*STATES(4)*STATES(3)*STATES(5)*(STATES(1) - ALGBRC(13))
      ALGBRC(16) =  0.500000*CONSTS(110)*log(CONSTS(8)/STATES(2))
      ALGBRC(17) =  CONSTS(21)*STATES(7)*STATES(6)*(STATES(1) - ALGBRC(16))
      ALGBRC(15) = ( CONSTS(20)*1.00000)/(1.00000+0.000298000/STATES(2))
      ALGBRC(32) =  CONSTS(28)*STATES(13)*STATES(12)*(STATES(1) - CONSTS(127))
      ALGBRC(30) =  CONSTS(27)*STATES(11)*STATES(10)*(STATES(1) - CONSTS(124))
      ALGBRC(20) = 1.00000/(1.00000+EXP(STATES(1)/- 17.0000 -  2.00000*log(STATES(2)/0.00100000)))
      ALGBRC(22) =  (CONSTS(24)+CONSTS(108))*ALGBRC(20)*(STATES(1) - CONSTS(118))
      ALGBRC(27) =  CONSTS(26)*STATES(9)*(STATES(1) - CONSTS(121))
      ALGBRC(18) =  CONSTS(23)*STATES(8)*(STATES(1) - CONSTS(117))
      ALGBRC(34) =  CONSTS(30)*STATES(15)*(STATES(1) - CONSTS(128))
      ALGBRC(25) =  CONSTS(25)*(STATES(1) - CONSTS(119))
      ALGBRC(46) = CONSTS(104) - STATES(19)
      ALGBRC(47) = ( CONSTS(68)*(1.00000+ CONSTS(69)*CONSTS(101))*CONSTS(70)*CONSTS(101)*ALGBRC(46))/(1.00000+ CONSTS(71)*ALGBRC(46)+ (1.00000+ CONSTS(72)*ALGBRC(46))*CONSTS(73)*CONSTS(101)+ (1.00000+ CONSTS(74)*ALGBRC(46))*CONSTS(75)*CONSTS(101))
      ALGBRC(28) = 1.00000/(1.00000+ CONSTS(78)*(1.00000+CONSTS(77)/1.00000+STATES(16)/CONSTS(76) ** 2.00000))
      ALGBRC(49) = CONSTS(98)+ 6.39440*ALGBRC(28)*ALGBRC(47)
      ALGBRC(35) = CONSTS(102) - STATES(20)
      ALGBRC(36) =  CONSTS(110)*log(( CONSTS(40)* STATES(20) ** (1.0 / 2))/ ALGBRC(35) ** (1.0 / 2))
      ALGBRC(37) = ( CONSTS(95)*0.500000*( ( CONSTS(44)*10.0000 **  6.00000*CONSTS(38)+ CONSTS(46)*EXP( 6.00000*CONSTS(48)*CONSTS(109)))*EXP( ALGBRC(36)*CONSTS(109))+ - 1.00000*CONSTS(44)*EXP( CONSTS(49)*6.00000*CONSTS(109)*STATES(22))+ CONSTS(47)*EXP( CONSTS(109)*ALGBRC(36))*EXP( CONSTS(109)*STATES(22)*6.00000*CONSTS(49))))/( (1.00000+ CONSTS(41)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( CONSTS(109)*CONSTS(48)*6.00000)+ (CONSTS(42)+ CONSTS(43)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( CONSTS(109)*STATES(22)*CONSTS(49)*6.00000))
      ALGBRC(48) = CONSTS(98)/3.00000+ 0.840000*ALGBRC(28)*ALGBRC(47)
      ALGBRC(39) = CONSTS(103) - STATES(18)
      ALGBRC(40) =  0.800000*STATES(18)
      ALGBRC(41) =  CONSTS(110)*log(( CONSTS(50)*ALGBRC(39))/( ALGBRC(40)*CONSTS(51)))
      ALGBRC(42) = ( CONSTS(96)*- 1.00000*( ( CONSTS(55)*10.0000 **  3.00000*CONSTS(38)+ CONSTS(57)*EXP( 3.00000*CONSTS(109)*CONSTS(48)))*EXP( CONSTS(109)*ALGBRC(41))+ - 1.00000*CONSTS(55)*EXP( 3.00000*CONSTS(109)*STATES(22))+ CONSTS(58)*EXP( CONSTS(109)*ALGBRC(41))*EXP( 3.00000*CONSTS(109)*STATES(22))))/( (1.00000+ CONSTS(52)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*CONSTS(48))+ (CONSTS(53)+ CONSTS(54)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*STATES(22)))
      ALGBRC(44) =  0.450000*ALGBRC(40)
      ALGBRC(45) =  0.0500000*ALGBRC(39)
      ALGBRC(52) =  0.300000*STATES(19)
      ALGBRC(54) =  0.450000*ALGBRC(52)
      ALGBRC(55) =  0.0500000*ALGBRC(46)
      ALGBRC(56) = ( CONSTS(99)*(1.00000 -  (( ALGBRC(55)*ALGBRC(44))/( ALGBRC(54)*ALGBRC(45)))*EXP( - 1.00000*CONSTS(109)*STATES(22))))/( (1.00000+ (ALGBRC(55)/ALGBRC(54))*EXP( - 1.00000*CONSTS(59)*CONSTS(109)*STATES(22)))*(1.00000+ALGBRC(44)/ALGBRC(45)))
      ALGBRC(50) =  0.150000*ALGBRC(47)
      ALGBRC(51) =  CONSTS(82)*ALGBRC(46)+CONSTS(132)
      ALGBRC(43) = ( - 1.00000*CONSTS(96)*3.00000*( CONSTS(55)*10.0000 **  3.00000*CONSTS(38)*EXP( CONSTS(109)*ALGBRC(41))+ CONSTS(56)*10.0000 **  3.00000*CONSTS(38)+ - 1.00000*(CONSTS(55)+CONSTS(56))*EXP( 3.00000*CONSTS(109)*STATES(22))))/( (1.00000+ CONSTS(52)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*CONSTS(48))+ (CONSTS(53)+ CONSTS(54)*EXP( CONSTS(109)*ALGBRC(41)))*EXP( 3.00000*CONSTS(109)*STATES(22)))
      ALGBRC(38) = ( CONSTS(95)*3.96600*( CONSTS(44)*10.0000 **  6.00000*CONSTS(38)*EXP( CONSTS(109)*ALGBRC(36))+ CONSTS(45)*10.0000 **  6.00000*CONSTS(38)+ - 1.00000*(CONSTS(44)+CONSTS(45))*EXP( CONSTS(49)*CONSTS(109)*STATES(22)*6.00000)))/( (1.00000+ CONSTS(41)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( 6.00000*CONSTS(109)*CONSTS(48))+ (CONSTS(42)+ CONSTS(43)*EXP( CONSTS(109)*ALGBRC(36)))*EXP( CONSTS(49)*6.00000*CONSTS(109)*STATES(22)))
      ALGBRC(31) = STATES(22) -  2.30300*CONSTS(110)*CONSTS(38)
      ALGBRC(33) =  CONSTS(97)*ALGBRC(31)
      ALGBRC(53) =  0.550000*ALGBRC(52)
      RETURN
      END