C
C There are a total of 31 entries in the algebraic variable array.
C There are a total of 10 entries in each of the rate and state variable arrays.
C There are a total of 13 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 C in component membrane (microF_per_cm2).
C ALGBRC(1) is i_Na in component fast_sodium_current (microA_per_cm2).
C ALGBRC(21) is i_si in component secondary_inward_current (microA_per_cm2).
C ALGBRC(23) is i_K2 in component pacemaker_potassium_current (microA_per_cm2).
C ALGBRC(25) is i_x1 in component plateau_potassium_current1 (microA_per_cm2).
C ALGBRC(27) is i_x2 in component plateau_potassium_current2 (microA_per_cm2).
C ALGBRC(28) is i_qr in component transient_chloride_current (microA_per_cm2).
C ALGBRC(29) is i_K1 in component time_independent_outward_current (microA_per_cm2).
C ALGBRC(30) is i_Na_b in component sodium_background_current (microA_per_cm2).
C ALGBRC(31) is i_Cl_b in component chloride_background_current (microA_per_cm2).
C CONSTS(2) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(3) is g_Na in component fast_sodium_current (milliS_per_cm2).
C STATES(2) is m in component fast_sodium_current_m_gate (dimensionless).
C STATES(3) is h in component fast_sodium_current_h_gate (dimensionless).
C ALGBRC(2) is alpha_m in component fast_sodium_current_m_gate (per_millisecond).
C ALGBRC(12) is beta_m in component fast_sodium_current_m_gate (per_millisecond).
C ALGBRC(3) is alpha_h in component fast_sodium_current_h_gate (per_millisecond).
C ALGBRC(13) is beta_h in component fast_sodium_current_h_gate (per_millisecond).
C CONSTS(4) is g_si in component secondary_inward_current (milliS_per_cm2).
C CONSTS(5) is g_si_ in component secondary_inward_current (milliS_per_cm2).
C CONSTS(6) is E_si in component secondary_inward_current (millivolt).
C STATES(4) is d in component secondary_inward_current_d_gate (dimensionless).
C STATES(5) is f in component secondary_inward_current_f_gate (dimensionless).
C ALGBRC(11) is d1 in component secondary_inward_current_d1_gate (dimensionless).
C ALGBRC(4) is alpha_d in component secondary_inward_current_d_gate (per_millisecond).
C ALGBRC(14) is beta_d in component secondary_inward_current_d_gate (per_millisecond).
C ALGBRC(5) is alpha_f in component secondary_inward_current_f_gate (per_millisecond).
C ALGBRC(15) is beta_f in component secondary_inward_current_f_gate (per_millisecond).
C ALGBRC(22) is I_K2 in component pacemaker_potassium_current (microA_per_cm2).
C CONSTS(7) is E_K in component pacemaker_potassium_current (millivolt).
C STATES(6) is s in component pacemaker_potassium_current_s_gate (dimensionless).
C ALGBRC(6) is alpha_s in component pacemaker_potassium_current_s_gate (per_millisecond).
C ALGBRC(16) is beta_s in component pacemaker_potassium_current_s_gate (per_millisecond).
C CONSTS(8) is E_s in component pacemaker_potassium_current_s_gate (millivolt).
C ALGBRC(24) is I_x1 in component plateau_potassium_current1 (microA_per_cm2).
C STATES(7) is x1 in component plateau_potassium_current1_x1_gate (dimensionless).
C ALGBRC(7) is alpha_x1 in component plateau_potassium_current1_x1_gate (per_millisecond).
C ALGBRC(17) is beta_x1 in component plateau_potassium_current1_x1_gate (per_millisecond).
C ALGBRC(26) is I_x2 in component plateau_potassium_current2 (microA_per_cm2).
C STATES(8) is x2 in component plateau_potassium_current2_x2_gate (dimensionless).
C ALGBRC(8) is alpha_x2 in component plateau_potassium_current2_x2_gate (per_millisecond).
C ALGBRC(18) is beta_x2 in component plateau_potassium_current2_x2_gate (per_millisecond).
C CONSTS(9) is E_Cl in component transient_chloride_current (millivolt).
C CONSTS(10) is g_qr in component transient_chloride_current (milliS_per_cm2).
C STATES(9) is q in component transient_chloride_current_q_gate (dimensionless).
C STATES(10) is r in component transient_chloride_current_r_gate (dimensionless).
C ALGBRC(9) is alpha_q in component transient_chloride_current_q_gate (per_millisecond).
C ALGBRC(19) is beta_q in component transient_chloride_current_q_gate (per_millisecond).
C ALGBRC(10) is alpha_r in component transient_chloride_current_r_gate (per_millisecond).
C ALGBRC(20) is beta_r in component transient_chloride_current_r_gate (per_millisecond).
C CONSTS(11) is E_K1 in component time_independent_outward_current (millivolt).
C CONSTS(12) is g_Nab in component sodium_background_current (milliS_per_cm2).
C CONSTS(13) is g_Clb in component chloride_background_current (milliS_per_cm2).
C RATES(1) is d/dt V in component membrane (millivolt).
C RATES(2) is d/dt m in component fast_sodium_current_m_gate (dimensionless).
C RATES(3) is d/dt h in component fast_sodium_current_h_gate (dimensionless).
C RATES(4) is d/dt d in component secondary_inward_current_d_gate (dimensionless).
C RATES(5) is d/dt f in component secondary_inward_current_f_gate (dimensionless).
C RATES(6) is d/dt s in component pacemaker_potassium_current_s_gate (dimensionless).
C RATES(7) is d/dt x1 in component plateau_potassium_current1_x1_gate (dimensionless).
C RATES(8) is d/dt x2 in component plateau_potassium_current2_x2_gate (dimensionless).
C RATES(9) is d/dt q in component transient_chloride_current_q_gate (dimensionless).
C RATES(10) is d/dt r in component transient_chloride_current_r_gate (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -78.041367
      CONSTS(1) = 10
      CONSTS(2) = 40
      CONSTS(3) = 150
      STATES(2) = 0.02566853
      STATES(3) = 0.78656359
      CONSTS(4) = 0.8
      CONSTS(5) = 0.04
      CONSTS(6) = 70
      STATES(4) = 0.00293135
      STATES(5) = 0.80873917
      CONSTS(7) = -110
      STATES(6) = 0.75473994
      CONSTS(8) = -52
      STATES(7) = 0.02030289
      STATES(8) = 0.0176854
      CONSTS(9) = -70
      CONSTS(10) = 2.5
      STATES(9) = 3.11285794
      STATES(10) = 0.13500116
      CONSTS(11) = -30
      CONSTS(12) = 0.105
      CONSTS(13) = 0.01
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = ( 1.00000*(STATES(1)+47.0000))/(1.00000 - EXP(- (STATES(1)+47.0000)/10.0000))
      ALGBRC(12) =  40.0000*EXP( - 0.0560000*(STATES(1)+72.0000))
      RATES(2) =  ALGBRC(2)*(1.00000 - STATES(2)) -  ALGBRC(12)*STATES(2)
      ALGBRC(3) =  0.00850000*EXP( - 0.184000*(STATES(1)+71.0000))
      ALGBRC(13) = 2.50000/(EXP( - 0.0820000*(STATES(1)+10.0000))+1.00000)
      RATES(3) =  ALGBRC(3)*(1.00000 - STATES(3)) -  ALGBRC(13)*STATES(3)
      ALGBRC(4) = ( 0.00200000*(STATES(1)+40.0000))/(1.00000 - EXP( - 0.100000*(STATES(1)+40.0000)))
      ALGBRC(14) =  0.0200000*EXP( - 0.0888000*(STATES(1)+40.0000))
      RATES(4) =  ALGBRC(4)*(1.00000 - STATES(4)) -  ALGBRC(14)*STATES(4)
      ALGBRC(5) =  0.000987000*EXP( - 0.0400000*(STATES(1)+60.0000))
      ALGBRC(15) = 0.0200000/(EXP( - 0.0870000*(STATES(1)+26.0000))+1.00000)
      RATES(5) =  ALGBRC(5)*(1.00000 - STATES(5)) -  ALGBRC(15)*STATES(5)
      ALGBRC(6) = ( 0.00100000*(STATES(1) - CONSTS(8)))/(1.00000 - EXP( - 0.200000*(STATES(1) - CONSTS(8))))
      ALGBRC(16) =  5.00000e-05*EXP( - 0.0670000*(STATES(1) - CONSTS(8)))
      RATES(6) =  ALGBRC(6)*(1.00000 - STATES(6)) -  ALGBRC(16)*STATES(6)
      ALGBRC(7) = ( 0.000500000*EXP((STATES(1)+50.0000)/12.1000))/(1.00000+EXP((STATES(1)+50.0000)/17.5000))
      ALGBRC(17) = ( 0.00130000*EXP(- (STATES(1)+20.0000)/16.6700))/(1.00000+EXP(- (STATES(1)+20.0000)/25.0000))
      RATES(7) =  ALGBRC(7)*(1.00000 - STATES(7)) -  ALGBRC(17)*STATES(7)
      ALGBRC(8) = ( 0.000127000*1.00000)/(1.00000+EXP(- (STATES(1)+19.0000)/5.00000))
      ALGBRC(18) = ( 0.000300000*EXP(- (STATES(1)+20.0000)/16.6700))/(1.00000+EXP(- (STATES(1)+20.0000)/25.0000))
      RATES(8) =  ALGBRC(8)*(1.00000 - STATES(8)) -  ALGBRC(18)*STATES(8)
      ALGBRC(9) = ( 0.00800000*STATES(1))/(1.00000 - EXP( - 0.100000*STATES(1)))
      ALGBRC(19) =  0.0800000*EXP( - 0.0888000*STATES(1))
      RATES(9) =  ALGBRC(9)*(1.00000 - STATES(9)) -  ALGBRC(19)*STATES(9)
      ALGBRC(10) =  0.000180000*EXP( - 0.0400000*(STATES(1)+80.0000))
      ALGBRC(20) = 0.0200000/(EXP( - 0.0870000*(STATES(1)+26.0000))+1.00000)
      RATES(10) =  ALGBRC(10)*(1.00000 - STATES(10)) -  ALGBRC(20)*STATES(10)
      ALGBRC(1) =  CONSTS(3)*STATES(2) ** 3.00000*STATES(3)*(STATES(1) - CONSTS(2))
      ALGBRC(11) = 1.00000/(1.00000+EXP( - 0.150000*(STATES(1)+40.0000)))
      ALGBRC(21) =  CONSTS(4)*STATES(4)*STATES(5)*(STATES(1) - CONSTS(6))+ CONSTS(5)*ALGBRC(11)*(STATES(1) - CONSTS(6))
      ALGBRC(22) = ( 2.80000*(EXP((STATES(1) - CONSTS(7))/25.0000) - 1.00000))/(EXP((STATES(1)+60.0000)/12.5000)+EXP((STATES(1)+60.0000)/25.0000))
      ALGBRC(23) =  ALGBRC(22)*STATES(6)
      ALGBRC(24) = ( 1.20000*(EXP((STATES(1)+95.0000)/25.0000) - 1.00000))/EXP((STATES(1)+45.0000)/25.0000)
      ALGBRC(25) =  STATES(7)*ALGBRC(24)
      ALGBRC(26) = 25.0000+ 1.00000*0.385000*STATES(1)
      ALGBRC(27) =  STATES(8)*ALGBRC(26)
      ALGBRC(28) =  CONSTS(10)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(9))
      ALGBRC(29) = ALGBRC(22)/2.80000+( 0.200000*(STATES(1) - CONSTS(11)))/(1.00000 - EXP(- (STATES(1) - CONSTS(11))/25.0000))
      ALGBRC(30) =  CONSTS(12)*(STATES(1) - CONSTS(2))
      ALGBRC(31) =  CONSTS(13)*(STATES(1) - CONSTS(9))
      RATES(1) = - (ALGBRC(1)+ALGBRC(21)+ALGBRC(23)+ALGBRC(25)+ALGBRC(27)+ALGBRC(28)+ALGBRC(29)+ALGBRC(30)+ALGBRC(31))/CONSTS(1)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = ( 1.00000*(STATES(1)+47.0000))/(1.00000 - EXP(- (STATES(1)+47.0000)/10.0000))
      ALGBRC(12) =  40.0000*EXP( - 0.0560000*(STATES(1)+72.0000))
      ALGBRC(3) =  0.00850000*EXP( - 0.184000*(STATES(1)+71.0000))
      ALGBRC(13) = 2.50000/(EXP( - 0.0820000*(STATES(1)+10.0000))+1.00000)
      ALGBRC(4) = ( 0.00200000*(STATES(1)+40.0000))/(1.00000 - EXP( - 0.100000*(STATES(1)+40.0000)))
      ALGBRC(14) =  0.0200000*EXP( - 0.0888000*(STATES(1)+40.0000))
      ALGBRC(5) =  0.000987000*EXP( - 0.0400000*(STATES(1)+60.0000))
      ALGBRC(15) = 0.0200000/(EXP( - 0.0870000*(STATES(1)+26.0000))+1.00000)
      ALGBRC(6) = ( 0.00100000*(STATES(1) - CONSTS(8)))/(1.00000 - EXP( - 0.200000*(STATES(1) - CONSTS(8))))
      ALGBRC(16) =  5.00000e-05*EXP( - 0.0670000*(STATES(1) - CONSTS(8)))
      ALGBRC(7) = ( 0.000500000*EXP((STATES(1)+50.0000)/12.1000))/(1.00000+EXP((STATES(1)+50.0000)/17.5000))
      ALGBRC(17) = ( 0.00130000*EXP(- (STATES(1)+20.0000)/16.6700))/(1.00000+EXP(- (STATES(1)+20.0000)/25.0000))
      ALGBRC(8) = ( 0.000127000*1.00000)/(1.00000+EXP(- (STATES(1)+19.0000)/5.00000))
      ALGBRC(18) = ( 0.000300000*EXP(- (STATES(1)+20.0000)/16.6700))/(1.00000+EXP(- (STATES(1)+20.0000)/25.0000))
      ALGBRC(9) = ( 0.00800000*STATES(1))/(1.00000 - EXP( - 0.100000*STATES(1)))
      ALGBRC(19) =  0.0800000*EXP( - 0.0888000*STATES(1))
      ALGBRC(10) =  0.000180000*EXP( - 0.0400000*(STATES(1)+80.0000))
      ALGBRC(20) = 0.0200000/(EXP( - 0.0870000*(STATES(1)+26.0000))+1.00000)
      ALGBRC(1) =  CONSTS(3)*STATES(2) ** 3.00000*STATES(3)*(STATES(1) - CONSTS(2))
      ALGBRC(11) = 1.00000/(1.00000+EXP( - 0.150000*(STATES(1)+40.0000)))
      ALGBRC(21) =  CONSTS(4)*STATES(4)*STATES(5)*(STATES(1) - CONSTS(6))+ CONSTS(5)*ALGBRC(11)*(STATES(1) - CONSTS(6))
      ALGBRC(22) = ( 2.80000*(EXP((STATES(1) - CONSTS(7))/25.0000) - 1.00000))/(EXP((STATES(1)+60.0000)/12.5000)+EXP((STATES(1)+60.0000)/25.0000))
      ALGBRC(23) =  ALGBRC(22)*STATES(6)
      ALGBRC(24) = ( 1.20000*(EXP((STATES(1)+95.0000)/25.0000) - 1.00000))/EXP((STATES(1)+45.0000)/25.0000)
      ALGBRC(25) =  STATES(7)*ALGBRC(24)
      ALGBRC(26) = 25.0000+ 1.00000*0.385000*STATES(1)
      ALGBRC(27) =  STATES(8)*ALGBRC(26)
      ALGBRC(28) =  CONSTS(10)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(9))
      ALGBRC(29) = ALGBRC(22)/2.80000+( 0.200000*(STATES(1) - CONSTS(11)))/(1.00000 - EXP(- (STATES(1) - CONSTS(11))/25.0000))
      ALGBRC(30) =  CONSTS(12)*(STATES(1) - CONSTS(2))
      ALGBRC(31) =  CONSTS(13)*(STATES(1) - CONSTS(9))
      RETURN
      END