C
C There are a total of 29 entries in the algebraic variable array.
C There are a total of 9 entries in each of the rate and state variable arrays.
C There are a total of 45 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C CONSTS(1) is RT in component parameters (joule_per_mole).
C STATES(1) is V_m in component membrane (millivolt).
C CONSTS(2) is C_m in component membrane (picofarad).
C ALGBRC(17) is I_Na in component I_Na (picoampere).
C ALGBRC(19) is I_T in component I_T (picoampere).
C ALGBRC(21) is I_KCa in component I_KCa (picoampere).
C ALGBRC(26) is I_Ca in component I_Ca (picoampere).
C ALGBRC(27) is I_L in component I_L (picoampere).
C ALGBRC(29) is I_K in component I_K (picoampere).
C ALGBRC(14) is I_stim in component I_stim (picoampere).
C CONSTS(3) is I_mag in component I_stim (picoampere).
C CONSTS(4) is t_0 in component I_stim (millisecond).
C CONSTS(5) is t_1 in component I_stim (millisecond).
C ALGBRC(4) is H_0 in component I_stim (dimensionless).
C ALGBRC(10) is H_1 in component I_stim (dimensionless).
C CONSTS(6) is g_Na in component I_Na (nanosiemens).
C CONSTS(7) is E_Na in component I_Na (millivolt).
C STATES(2) is m in component m (dimensionless).
C STATES(3) is h in component h (dimensionless).
C CONSTS(8) is E_m in component m (millivolt).
C CONSTS(9) is V_alpha_m in component m (millivolt).
C CONSTS(10) is V_beta_m in component m (millivolt).
C ALGBRC(1) is alpha_m in component m (per_millisecond).
C ALGBRC(7) is beta_m in component m (per_millisecond).
C CONSTS(11) is alpha_m_max in component m (per_millivolt_per_millisecond).
C CONSTS(12) is beta_m_max in component m (per_millisecond).
C CONSTS(13) is E_h in component h (millivolt).
C CONSTS(14) is V_alpha_h in component h (millivolt).
C CONSTS(15) is V_beta_h in component h (millivolt).
C ALGBRC(2) is alpha_h in component h (per_millisecond).
C ALGBRC(8) is beta_h in component h (per_millisecond).
C CONSTS(16) is alpha_h_max in component h (per_millisecond).
C CONSTS(17) is beta_h_max in component h (per_millisecond).
C CONSTS(18) is R_s in component I_T (per_nanosiemens).
C CONSTS(19) is C_t in component I_T (picofarad).
C STATES(4) is V_T in component V_T (millivolt).
C CONSTS(20) is g_KCa in component I_KCa (nanosiemens).
C CONSTS(21) is E_K in component I_K (millivolt).
C STATES(5) is o in component o (dimensionless).
C ALGBRC(20) is w in component w (dimensionless).
C STATES(6) is c in component calcium_handling (dimensionless).
C ALGBRC(16) is o_oinf in component o (dimensionless).
C ALGBRC(3) is alpha_Vm in component o (per_millisecond).
C ALGBRC(9) is beta_Vm in component o (per_millisecond).
C ALGBRC(13) is tau_Vm in component o (millisecond).
C CONSTS(22) is d_1 in component o (dimensionless).
C CONSTS(23) is d_2 in component o (dimensionless).
C CONSTS(24) is k_1 in component o (dimensionless).
C CONSTS(25) is k_2 in component o (dimensionless).
C CONSTS(26) is a_bar in component o (per_millisecond).
C CONSTS(27) is b_bar in component o (per_millisecond).
C CONSTS(28) is kd in component w (dimensionless).
C STATES(7) is cer in component calcium_handling (dimensionless).
C ALGBRC(28) is j_mem in component calcium_handling (picoampere).
C ALGBRC(22) is j_leak in component calcium_handling (picoampere).
C ALGBRC(23) is j_serca in component calcium_handling (picoampere).
C ALGBRC(24) is j_er in component calcium_handling (picoampere).
C CONSTS(29) is alpha in component calcium_handling (dimensionless).
C CONSTS(30) is k_pmca in component calcium_handling (dimensionless).
C CONSTS(31) is k_serca in component calcium_handling (picoampere).
C CONSTS(32) is p_leak in component calcium_handling (picoampere).
C CONSTS(33) is f_er in component calcium_handling (per_picoampere_per_millisecond).
C CONSTS(34) is f_cyt in component calcium_handling (per_picoampere_per_millisecond).
C CONSTS(35) is v_cytver in component calcium_handling (dimensionless).
C STATES(8) is d in component d (dimensionless).
C ALGBRC(25) is g_Ca in component I_Ca (nanosiemens).
C CONSTS(36) is g_Ca0 in component I_Ca (nanosiemens_per_millivolt).
C CONSTS(37) is E_Ca in component I_Ca (millivolt).
C ALGBRC(15) is alpha_d in component d (per_millisecond).
C ALGBRC(18) is beta_d in component d (per_millisecond).
C ALGBRC(11) is tau_d in component d (millisecond).
C ALGBRC(5) is d_infinity in component d (dimensionless).
C CONSTS(38) is E_L in component I_L (millivolt).
C CONSTS(39) is g_Lmax in component I_L (nanosiemens).
C STATES(9) is n in component n (dimensionless).
C CONSTS(40) is g_K in component I_K (nanosiemens).
C CONSTS(41) is E_n in component n (millivolt).
C CONSTS(42) is V_alpha_n in component n (millivolt).
C CONSTS(43) is V_beta_n in component n (millivolt).
C ALGBRC(6) is alpha_n in component n (per_millisecond).
C ALGBRC(12) is beta_n in component n (per_millisecond).
C CONSTS(44) is alpha_n_max in component n (per_millivolt_per_millisecond).
C CONSTS(45) is beta_n_max in component n (per_millisecond).
C RATES(1) is d/dt V_m in component membrane (millivolt).
C RATES(2) is d/dt m in component m (dimensionless).
C RATES(3) is d/dt h in component h (dimensionless).
C RATES(4) is d/dt V_T in component V_T (millivolt).
C RATES(5) is d/dt o in component o (dimensionless).
C RATES(6) is d/dt c in component calcium_handling (dimensionless).
C RATES(7) is d/dt cer in component calcium_handling (dimensionless).
C RATES(8) is d/dt d in component d (dimensionless).
C RATES(9) is d/dt n in component n (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 2577
      STATES(1) = -75
      CONSTS(2) = 0.009
      CONSTS(3) = 2
      CONSTS(4) = 5
      CONSTS(5) = 6
      CONSTS(6) = 0.45
      CONSTS(7) = 50
      STATES(2) = 0
      STATES(3) = 1
      CONSTS(8) = -42
      CONSTS(9) = 10
      CONSTS(10) = 18
      CONSTS(11) = 0.208
      CONSTS(12) = 2.081
      CONSTS(13) = -41
      CONSTS(14) = 14.7
      CONSTS(15) = 7.6
      CONSTS(16) = 0.0156
      CONSTS(17) = 3.382
      CONSTS(18) = 15
      CONSTS(19) = 0.04
      STATES(4) = -70
      CONSTS(20) = 0.5
      CONSTS(21) = -70
      STATES(5) = 0.15
      STATES(6) = 0.15
      CONSTS(22) = 0.84
      CONSTS(23) = 1
      CONSTS(24) = 0.18
      CONSTS(25) = 0.011
      CONSTS(26) = 0.48
      CONSTS(27) = 0.28
      CONSTS(28) = 0.18
      STATES(7) = 200
      CONSTS(29) = 4.5e-6
      CONSTS(30) = 0.2
      CONSTS(31) = 0.4
      CONSTS(32) = 0.0005
      CONSTS(33) = 0.01
      CONSTS(34) = 0.01
      CONSTS(35) = 5
      STATES(8) = 0
      CONSTS(36) = 0.05
      CONSTS(37) = 50
      CONSTS(38) = -75
      CONSTS(39) = 0.0024
      STATES(9) = 0
      CONSTS(40) = 0.415
      CONSTS(41) = -40
      CONSTS(42) = 7
      CONSTS(43) = 40
      CONSTS(44) = 0.0229
      CONSTS(45) = 0.09616
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(4) = (STATES(1) - STATES(4))/( CONSTS(18)*CONSTS(19))
      ALGBRC(1) = ( CONSTS(11)*(STATES(1) - CONSTS(8)))/(1.00000 - EXP((CONSTS(8) - STATES(1))/CONSTS(9)))
      ALGBRC(7) =  CONSTS(12)*EXP((CONSTS(8) - STATES(1))/CONSTS(10))
      RATES(2) =  ALGBRC(1)*(1.00000 - STATES(2)) -  ALGBRC(7)*STATES(2)
      ALGBRC(2) =  CONSTS(16)*EXP((STATES(1) - CONSTS(13))/CONSTS(14))
      ALGBRC(8) = CONSTS(17)/(1.00000+EXP((CONSTS(13) - STATES(1))/CONSTS(15)))
      RATES(3) =  ALGBRC(2)*(1.00000 - STATES(3)) -  ALGBRC(8)*STATES(3)
      ALGBRC(6) = ( CONSTS(44)*(STATES(1) - CONSTS(41)))/(1.00000 - EXP((CONSTS(41) - STATES(1))/CONSTS(42)))
      ALGBRC(12) =  CONSTS(45)*EXP((CONSTS(41) - STATES(1))/CONSTS(43))
      RATES(9) =  ALGBRC(6)*(1.00000 - STATES(9)) -  ALGBRC(12)*STATES(9)
      ALGBRC(3) = CONSTS(26)/(1.00000+ CONSTS(24)*EXP(( - 2.00000*CONSTS(22)*96.4850*STATES(1))/( CONSTS(1)*STATES(6))))
      ALGBRC(9) = CONSTS(27)/(1.00000+STATES(6)/( CONSTS(25)*EXP(( - 2.00000*CONSTS(23)*96.4850*STATES(1))/CONSTS(1))))
      ALGBRC(13) = 1.00000/(ALGBRC(3)+ALGBRC(9))
      ALGBRC(16) =  ALGBRC(3)*ALGBRC(13)
      RATES(5) = (ALGBRC(16) - STATES(5))/ALGBRC(13)
      ALGBRC(11) = ( 80.0000*1.00000)/cosh( - 0.0310000*(STATES(1)+37.1000))
      ALGBRC(5) = 1.00000/(1.00000+EXP((- 24.6000 - STATES(1))/11.3000))
      ALGBRC(15) = ALGBRC(5)/ALGBRC(11)
      ALGBRC(18) = (1.00000 - ALGBRC(5))/ALGBRC(11)
      RATES(8) =  ALGBRC(15)*(1.00000 - STATES(8)) -  ALGBRC(18)*STATES(8)
      ALGBRC(22) =  CONSTS(32)*(STATES(7) - STATES(6))
      ALGBRC(23) =  CONSTS(31)*STATES(6)
      ALGBRC(24) = ALGBRC(22) - ALGBRC(23)
      RATES(7) =  - CONSTS(33)*CONSTS(35)*ALGBRC(24)
      ALGBRC(25) = ( - CONSTS(36)*STATES(1))/(EXP( 0.117000*STATES(1)) - 1.00000)
      ALGBRC(26) =  ALGBRC(25)*STATES(8) ** 2.00000*(STATES(1) - CONSTS(37))
      ALGBRC(28) =  - CONSTS(29)*ALGBRC(26)*CONSTS(30)*STATES(6)
      RATES(6) =  CONSTS(34)*(ALGBRC(28)+ALGBRC(24))
      ALGBRC(17) =  CONSTS(6)*STATES(2) ** 3.00000*STATES(3)*(STATES(1) - CONSTS(7))
      ALGBRC(19) = (STATES(1) - STATES(4))/CONSTS(18)
      ALGBRC(20) = STATES(6) ** 5.00000/(STATES(6) ** 5.00000+CONSTS(28) ** 5.00000)
      ALGBRC(21) =  CONSTS(20)*STATES(5)*ALGBRC(20)*(STATES(1) - CONSTS(21))
      ALGBRC(27) =  CONSTS(39)*(STATES(1) - CONSTS(38))
      ALGBRC(29) =  CONSTS(40)*STATES(9) ** 4.00000*(STATES(1) - CONSTS(21))
      ALGBRC(4) = TERNRY(VOI.LT.CONSTS(4), 0.00000, 1.00000)
      ALGBRC(10) = TERNRY(VOI.LT.CONSTS(5), 0.00000, 1.00000)
      ALGBRC(14) =  CONSTS(3)*(ALGBRC(4) - ALGBRC(10))
      RATES(1) = (ALGBRC(14) - (ALGBRC(17)+ALGBRC(26)+ALGBRC(29)+ALGBRC(27)+ALGBRC(19)+ALGBRC(21)))/CONSTS(2)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = ( CONSTS(11)*(STATES(1) - CONSTS(8)))/(1.00000 - EXP((CONSTS(8) - STATES(1))/CONSTS(9)))
      ALGBRC(7) =  CONSTS(12)*EXP((CONSTS(8) - STATES(1))/CONSTS(10))
      ALGBRC(2) =  CONSTS(16)*EXP((STATES(1) - CONSTS(13))/CONSTS(14))
      ALGBRC(8) = CONSTS(17)/(1.00000+EXP((CONSTS(13) - STATES(1))/CONSTS(15)))
      ALGBRC(6) = ( CONSTS(44)*(STATES(1) - CONSTS(41)))/(1.00000 - EXP((CONSTS(41) - STATES(1))/CONSTS(42)))
      ALGBRC(12) =  CONSTS(45)*EXP((CONSTS(41) - STATES(1))/CONSTS(43))
      ALGBRC(3) = CONSTS(26)/(1.00000+ CONSTS(24)*EXP(( - 2.00000*CONSTS(22)*96.4850*STATES(1))/( CONSTS(1)*STATES(6))))
      ALGBRC(9) = CONSTS(27)/(1.00000+STATES(6)/( CONSTS(25)*EXP(( - 2.00000*CONSTS(23)*96.4850*STATES(1))/CONSTS(1))))
      ALGBRC(13) = 1.00000/(ALGBRC(3)+ALGBRC(9))
      ALGBRC(16) =  ALGBRC(3)*ALGBRC(13)
      ALGBRC(11) = ( 80.0000*1.00000)/cosh( - 0.0310000*(STATES(1)+37.1000))
      ALGBRC(5) = 1.00000/(1.00000+EXP((- 24.6000 - STATES(1))/11.3000))
      ALGBRC(15) = ALGBRC(5)/ALGBRC(11)
      ALGBRC(18) = (1.00000 - ALGBRC(5))/ALGBRC(11)
      ALGBRC(22) =  CONSTS(32)*(STATES(7) - STATES(6))
      ALGBRC(23) =  CONSTS(31)*STATES(6)
      ALGBRC(24) = ALGBRC(22) - ALGBRC(23)
      ALGBRC(25) = ( - CONSTS(36)*STATES(1))/(EXP( 0.117000*STATES(1)) - 1.00000)
      ALGBRC(26) =  ALGBRC(25)*STATES(8) ** 2.00000*(STATES(1) - CONSTS(37))
      ALGBRC(28) =  - CONSTS(29)*ALGBRC(26)*CONSTS(30)*STATES(6)
      ALGBRC(17) =  CONSTS(6)*STATES(2) ** 3.00000*STATES(3)*(STATES(1) - CONSTS(7))
      ALGBRC(19) = (STATES(1) - STATES(4))/CONSTS(18)
      ALGBRC(20) = STATES(6) ** 5.00000/(STATES(6) ** 5.00000+CONSTS(28) ** 5.00000)
      ALGBRC(21) =  CONSTS(20)*STATES(5)*ALGBRC(20)*(STATES(1) - CONSTS(21))
      ALGBRC(27) =  CONSTS(39)*(STATES(1) - CONSTS(38))
      ALGBRC(29) =  CONSTS(40)*STATES(9) ** 4.00000*(STATES(1) - CONSTS(21))
      ALGBRC(4) = TERNRY(VOI.LT.CONSTS(4), 0.00000, 1.00000)
      ALGBRC(10) = TERNRY(VOI.LT.CONSTS(5), 0.00000, 1.00000)
      ALGBRC(14) =  CONSTS(3)*(ALGBRC(4) - ALGBRC(10))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END