C
C There are a total of 13 entries in the algebraic variable array.
C There are a total of 5 entries in each of the rate and state variable arrays.
C There are a total of 20 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C STATES(1) is V_s in component soma_compartment (mV).
C CONSTS(1) is V_Na in component soma_compartment (mV).
C CONSTS(2) is V_K in component soma_compartment (mV).
C ALGBRC(9) is I_Na in component soma_compartment (uA_per_cm2).
C ALGBRC(1) is I_K_DR in component soma_compartment (uA_per_cm2).
C CONSTS(3) is g_K_DR in component soma_compartment (mS_per_cm2).
C CONSTS(4) is g_Na in component soma_compartment (mS_per_cm2).
C CONSTS(5) is g_c in component model_parameters (mS_per_cm2).
C STATES(2) is V_D in component dendritic_compartment (mV).
C CONSTS(6) is C_m in component model_parameters (uF_per_cm2).
C CONSTS(7) is p in component model_parameters (dimensionless).
C STATES(3) is n in component gating_variables (dimensionless).
C STATES(4) is h in component gating_variables (dimensionless).
C ALGBRC(5) is m_infinity in component gating_variables (dimensionless).
C CONSTS(8) is V_L in component dendritic_compartment (mV).
C CONSTS(9) is V_NMDA in component dendritic_compartment (mV).
C ALGBRC(2) is I_L in component dendritic_compartment (uA_per_cm2).
C ALGBRC(13) is I_D in component dendritic_compartment (uA_per_cm2).
C ALGBRC(10) is I_pump in component dendritic_compartment (uA_per_cm2).
C ALGBRC(11) is I_NMDA in component dendritic_compartment (uA_per_cm2).
C ALGBRC(12) is I_Na_NMDA in component dendritic_compartment (uA_per_cm2).
C CONSTS(10) is R_pump in component dendritic_compartment (uA_per_cm2).
C CONSTS(11) is alpha in component dendritic_compartment (mMcm2_per_uAs).
C CONSTS(12) is g_NMDA in component dendritic_compartment (mS_per_cm2).
C CONSTS(13) is g_Na_NMDA in component dendritic_compartment (mS_per_cm2).
C CONSTS(14) is g_L in component dendritic_compartment (mS_per_cm2).
C STATES(5) is Na in component dendritic_compartment (mM).
C CONSTS(15) is Na_eq in component dendritic_compartment (mM).
C CONSTS(16) is K_p in component dendritic_compartment (mM).
C CONSTS(17) is Mg_o in component dendritic_compartment (mM).
C CONSTS(18) is K_Mg in component dendritic_compartment (mM).
C CONSTS(19) is q in component dendritic_compartment (mV).
C ALGBRC(6) is phi_Na in component dendritic_compartment (dimensionless).
C CONSTS(20) is phi_Na_eq in component dendritic_compartment (dimensionless).
C ALGBRC(3) is n_infinity in component gating_variables (dimensionless).
C ALGBRC(4) is h_infinity in component gating_variables (dimensionless).
C ALGBRC(7) is tau_h in component gating_variables (second).
C ALGBRC(8) is tau_n in component gating_variables (second).
C RATES(1) is d/dt V_s in component soma_compartment (mV).
C RATES(2) is d/dt V_D in component dendritic_compartment (mV).
C RATES(5) is d/dt Na in component dendritic_compartment (mM).
C RATES(4) is d/dt h in component gating_variables (dimensionless).
C RATES(3) is d/dt n in component gating_variables (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -64
      CONSTS(1) = 55
      CONSTS(2) = -85
      CONSTS(3) = 3.2
      CONSTS(4) = 3.2
      CONSTS(5) = 0.0
      STATES(2) = -25.0
      CONSTS(6) = 1
      CONSTS(7) = 0.5
      STATES(3) = 0.002
      STATES(4) = 1
      CONSTS(8) = -50
      CONSTS(9) = 0
      CONSTS(10) = 18.0
      CONSTS(11) = 0.5
      CONSTS(12) = 1.25
      CONSTS(13) = 1.0
      CONSTS(14) = 0.18
      STATES(5) = 5.09
      CONSTS(15) = 8
      CONSTS(16) = 15
      CONSTS(17) = 1.4
      CONSTS(18) = 10.0
      CONSTS(19) = 12.5
      CONSTS(20) = CONSTS(15) ** 3.00000/(CONSTS(15) ** 3.00000+CONSTS(16) ** 3.00000)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/8.30000))
      ALGBRC(7) =  0.400000*(1.00000+2.00000/(1.00000+EXP((STATES(1)+25.0000)/5.00000)))
      RATES(4) = (ALGBRC(4) - STATES(4))/ALGBRC(7)
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+31.0000)/5.30000))
      ALGBRC(8) = (0.800000+1.60000/(1.00000+EXP( 0.100000*(STATES(1)+25.0000))))/(1.00000+EXP( - 0.100000*(STATES(1)+70.0000)))
      RATES(3) = (ALGBRC(3) - STATES(3))/ALGBRC(8)
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+35.0000)/6.20000))
      ALGBRC(9) =  CONSTS(4)*STATES(4)*(STATES(1) - CONSTS(1))*ALGBRC(5) ** 3.00000
      ALGBRC(1) =  CONSTS(3)*(STATES(1) - CONSTS(2))*STATES(3) ** 2.00000
      RATES(1) = - ( 1000.00*(ALGBRC(1)+ALGBRC(9)+ (CONSTS(5)/CONSTS(7))*(STATES(1) - STATES(2))))/CONSTS(6)
      ALGBRC(2) =  CONSTS(14)*(STATES(2) - CONSTS(8))
      ALGBRC(6) = STATES(5) ** 3.00000/(STATES(5) ** 3.00000+CONSTS(16) ** 3.00000)
      ALGBRC(10) =  CONSTS(10)*(ALGBRC(6) - CONSTS(20))
      ALGBRC(11) =  (CONSTS(12)/(1.00000+ (CONSTS(17)/CONSTS(18))*EXP(- (STATES(2)/CONSTS(19)))))*(STATES(2) - CONSTS(9))
      RATES(2) = - ( 1000.00*(ALGBRC(11)+ALGBRC(10)+ALGBRC(2)+ (CONSTS(5)/(1.00000 - CONSTS(7)))*(STATES(2) - STATES(1))))/CONSTS(6)
      ALGBRC(12) =  (CONSTS(13)/(1.00000+ (CONSTS(17)/CONSTS(18))*EXP(- (STATES(2)/CONSTS(19)))))*(STATES(2) - CONSTS(1))
      RATES(5) =  CONSTS(11)*(- ALGBRC(12) -  3.00000*ALGBRC(10))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/8.30000))
      ALGBRC(7) =  0.400000*(1.00000+2.00000/(1.00000+EXP((STATES(1)+25.0000)/5.00000)))
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+31.0000)/5.30000))
      ALGBRC(8) = (0.800000+1.60000/(1.00000+EXP( 0.100000*(STATES(1)+25.0000))))/(1.00000+EXP( - 0.100000*(STATES(1)+70.0000)))
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+35.0000)/6.20000))
      ALGBRC(9) =  CONSTS(4)*STATES(4)*(STATES(1) - CONSTS(1))*ALGBRC(5) ** 3.00000
      ALGBRC(1) =  CONSTS(3)*(STATES(1) - CONSTS(2))*STATES(3) ** 2.00000
      ALGBRC(2) =  CONSTS(14)*(STATES(2) - CONSTS(8))
      ALGBRC(6) = STATES(5) ** 3.00000/(STATES(5) ** 3.00000+CONSTS(16) ** 3.00000)
      ALGBRC(10) =  CONSTS(10)*(ALGBRC(6) - CONSTS(20))
      ALGBRC(11) =  (CONSTS(12)/(1.00000+ (CONSTS(17)/CONSTS(18))*EXP(- (STATES(2)/CONSTS(19)))))*(STATES(2) - CONSTS(9))
      ALGBRC(12) =  (CONSTS(13)/(1.00000+ (CONSTS(17)/CONSTS(18))*EXP(- (STATES(2)/CONSTS(19)))))*(STATES(2) - CONSTS(1))
      ALGBRC(13) = ALGBRC(11)+ALGBRC(10)+ALGBRC(2)
      RETURN
      END