TITLE Mod file for component: Component(id=MN_na type=ionChannelHH)

COMMENT

    This NEURON file has been generated by org.neuroml.export (see https://github.com/NeuroML/org.neuroml.export)
         org.neuroml.export  v1.4.1
         org.neuroml.model   v1.4.1
         jLEMS               v0.9.7.3

ENDCOMMENT

NEURON {
    SUFFIX aIN_na
    USEION na READ ena WRITE ina
    RANGE gmax, ina
}

UNITS {
    (mA) = (milliamp)
    (mV) = (millivolt)
    (S) = (siemens)
}

PARAMETER {
    gmax = 0  (S/cm2)
}

ASSIGNED {
    v (mV)
    ena (mV)
    ina (mA/cm2)
    na_minf
    na_hinf
    na_mtau (ms)
    na_htau (ms)
    rate_h
    rate_m
    
}

STATE {
    h 
    m
}

INITIAL {
    ena = 50.0
    rates()
    h = na_hinf
    m = na_minf
    
}

BREAKPOINT {
    SOLVE states METHOD cnexp
    ina = gmax * m*m*m * h * (v - ena)
}

DERIVATIVE states {
    rates()
    m' =  ( na_minf  -  m ) /  na_mtau
    h' =  ( na_hinf - h ) / na_htau
}

UNITSOFF
PROCEDURE rates() {LOCAL alpha_m, beta_m, alpha_h, beta_h
    alpha_m = 8.67/(0.5 + exp((-13.01+v)/(-18.56)))
    beta_m =  5.73/(1.0 + exp((2.99+v)/9.69))
    alpha_h = 0.04/exp((15.8+v)/26.0)
    beta_h =  4.08/(0.001 + exp((-19.09+v)/-10.21))

    na_mtau = 1/( alpha_m + beta_m)
    na_htau = 1/( alpha_h + beta_h)
    na_minf = alpha_m * na_mtau
    na_hinf = alpha_h * na_htau
    
}
UNITSON