C
C There are a total of 10 entries in the algebraic variable array.
C There are a total of 4 entries in each of the rate and state variable arrays.
C There are a total of 21 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C CONSTS(1) is Cm in component membrane (femtoF).
C STATES(1) is V in component membrane (millivolt).
C ALGBRC(5) is ICa in component Ca_current (femtoA).
C ALGBRC(7) is IK in component rapid_K_current (femtoA).
C ALGBRC(10) is Il in component leak_current (femtoA).
C ALGBRC(8) is Is1 in component slow_K_current (femtoA).
C ALGBRC(9) is Is2 in component very_slow_K_current (femtoA).
C CONSTS(2) is Vm in component Ca_current (millivolt).
C CONSTS(3) is VCa in component Ca_current (millivolt).
C CONSTS(4) is gCa in component Ca_current (picoS).
C ALGBRC(1) is minf in component Ca_current (dimensionless).
C CONSTS(5) is sm in component Ca_current (millivolt).
C CONSTS(6) is VK in component rapid_K_current (millivolt).
C CONSTS(7) is gK in component rapid_K_current (picoS).
C STATES(2) is n in component rapid_K_current (dimensionless).
C CONSTS(8) is lambda in component rapid_K_current (dimensionless).
C CONSTS(9) is tnbar in component rapid_K_current (dimensionless).
C CONSTS(10) is Vn in component rapid_K_current (millivolt).
C CONSTS(11) is sn in component rapid_K_current (millivolt).
C ALGBRC(6) is taun in component rapid_K_current (dimensionless).
C ALGBRC(2) is ninf in component rapid_K_current (dimensionless).
C CONSTS(12) is gs1 in component slow_K_current (picoS).
C STATES(3) is s1 in component slow_K_current (dimensionless).
C ALGBRC(3) is s1inf in component slow_K_current (dimensionless).
C CONSTS(13) is Vs1 in component slow_K_current (millivolt).
C CONSTS(14) is ss1 in component slow_K_current (millivolt).
C CONSTS(15) is taus1 in component slow_K_current (dimensionless).
C CONSTS(16) is Vs2 in component very_slow_K_current (millivolt).
C STATES(4) is s2 in component very_slow_K_current (dimensionless).
C ALGBRC(4) is s2inf in component very_slow_K_current (dimensionless).
C CONSTS(17) is ss2 in component very_slow_K_current (millivolt).
C CONSTS(18) is gs2 in component very_slow_K_current (picoS).
C CONSTS(19) is taus2 in component very_slow_K_current (dimensionless).
C CONSTS(20) is gl in component leak_current (picoS).
C CONSTS(21) is Vl in component leak_current (millivolt).
C RATES(1) is d/dt V in component membrane (millivolt).
C RATES(2) is d/dt n in component rapid_K_current (dimensionless).
C RATES(3) is d/dt s1 in component slow_K_current (dimensionless).
C RATES(4) is d/dt s2 in component very_slow_K_current (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 4524
      STATES(1) = -43
      CONSTS(2) = -22
      CONSTS(3) = 100
      CONSTS(4) = 280
      CONSTS(5) = 7.5
      CONSTS(6) = -80
      CONSTS(7) = 1300
      STATES(2) = 0.03
      CONSTS(8) = 1.1
      CONSTS(9) = 9.09
      CONSTS(10) = -9
      CONSTS(11) = 10
      CONSTS(12) = 20
      STATES(3) = 0.1
      CONSTS(13) = -40
      CONSTS(14) = 0.5
      CONSTS(15) = 1000
      CONSTS(16) = -42
      STATES(4) = 0.434
      CONSTS(17) = 0.4
      CONSTS(18) = 32
      CONSTS(19) = 120000
      CONSTS(20) = 25
      CONSTS(21) = -40
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/(1.00000+EXP((CONSTS(13) - STATES(1))/CONSTS(14)))
      RATES(3) = (ALGBRC(3) - STATES(3))/( CONSTS(15)*1.00000)
      ALGBRC(4) = 1.00000/(1.00000+EXP((CONSTS(16) - STATES(1))/CONSTS(17)))
      RATES(4) = (ALGBRC(4) - STATES(4))/( CONSTS(19)*1.00000)
      ALGBRC(6) = CONSTS(9)/(1.00000+EXP((STATES(1) - CONSTS(10))/CONSTS(11)))
      ALGBRC(2) = 1.00000/(1.00000+EXP((CONSTS(10) - STATES(1))/CONSTS(11)))
      RATES(2) = ( CONSTS(8)*(ALGBRC(2) - STATES(2)))/( ALGBRC(6)*1.00000)
      ALGBRC(1) = 1.00000/(1.00000+EXP((CONSTS(2) - STATES(1))/CONSTS(5)))
      ALGBRC(5) =  CONSTS(4)*ALGBRC(1)*(STATES(1) - CONSTS(3))
      ALGBRC(7) =  CONSTS(7)*STATES(2)*(STATES(1) - CONSTS(6))
      ALGBRC(10) =  CONSTS(20)*(STATES(1) - CONSTS(21))
      ALGBRC(8) =  CONSTS(12)*STATES(3)*(STATES(1) - CONSTS(6))
      ALGBRC(9) =  CONSTS(18)*STATES(4)*(STATES(1) - CONSTS(6))
      RATES(1) = - (ALGBRC(5)+ALGBRC(7)+ALGBRC(10)+ALGBRC(8)+ALGBRC(9))/CONSTS(1)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/(1.00000+EXP((CONSTS(13) - STATES(1))/CONSTS(14)))
      ALGBRC(4) = 1.00000/(1.00000+EXP((CONSTS(16) - STATES(1))/CONSTS(17)))
      ALGBRC(6) = CONSTS(9)/(1.00000+EXP((STATES(1) - CONSTS(10))/CONSTS(11)))
      ALGBRC(2) = 1.00000/(1.00000+EXP((CONSTS(10) - STATES(1))/CONSTS(11)))
      ALGBRC(1) = 1.00000/(1.00000+EXP((CONSTS(2) - STATES(1))/CONSTS(5)))
      ALGBRC(5) =  CONSTS(4)*ALGBRC(1)*(STATES(1) - CONSTS(3))
      ALGBRC(7) =  CONSTS(7)*STATES(2)*(STATES(1) - CONSTS(6))
      ALGBRC(10) =  CONSTS(20)*(STATES(1) - CONSTS(21))
      ALGBRC(8) =  CONSTS(12)*STATES(3)*(STATES(1) - CONSTS(6))
      ALGBRC(9) =  CONSTS(18)*STATES(4)*(STATES(1) - CONSTS(6))
      RETURN
      END