/usr/share/nrn/demo/release/khhchan.mod is in neuron 7.5-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 | TITLE HH k channel channel
: Hodgkin - Huxley k channel
NEURON {
SUFFIX HHk
USEION k READ ek WRITE ik
RANGE gkbar, ik
GLOBAL inf, tau
THREADSAFE
}
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
}
PARAMETER {
gkbar=.036 (mho/cm2) <0,1e9>
ek (mV) : -77 suggested, default set by NEURON
}
STATE {
n
}
ASSIGNED {
v (mV)
celsius (degC) : 16
ik (mA/cm2)
inf
tau (ms)
}
INITIAL {
rate(v*1(/mV))
n = inf
}
BREAKPOINT {
SOLVE states METHOD cnexp
ik = gkbar*n*n*n*n*(v - ek)
}
DERIVATIVE states { : exact when v held constant
rate(v*1(/mV))
n' = (inf - n)/tau
}
UNITSOFF
FUNCTION alp(v(mV)) { LOCAL q10
v = -v - 65
q10 = 3^((celsius - 6.3)/10)
alp = q10 * .01*expM1(v + 10, 10)
}
FUNCTION bet(v(mV)) { LOCAL q10
v = -v - 65
q10 = 3^((celsius - 6.3)/10)
bet = q10 * .125*exp(v/80)
}
FUNCTION expM1(x,y) {
if (fabs(x/y) < 1e-6) {
expM1 = y*(1 - x/y/2)
}else{
expM1 = x/(exp(x/y) - 1)
}
}
PROCEDURE rate(v) {LOCAL a, b :rest = -70
TABLE inf, tau DEPEND celsius FROM -100 TO 100 WITH 200
a = alp(v) b=bet(v)
tau = 1/(a + b)
inf = a/(a + b)
}
UNITSON
|