~/datastore/nrn-5.8/i686/bin/nrnivmodl

"/home/sterratt/datastore/nrn-5.8/i686/bin/nocmodl" psd

Translating psd.mod into psd.c

nocmodl: diffeq.y:194: difeq_yyerror: Assertion `0' failed.

make: *** [psd.lo] Aborted

It seems to be a problem with the DERIVATIVE block (no error when I delete it), but I can't get any further than that. (I've commented out a fair amount to simplify things). Any help would be much appreciated.

David.

psd.mod:

Code: Select all

```
COMMENT
psd.mod "Postsynaptic density" based on Rubin et al (2005) "Calcium
Time Course as a signal for STDSP".
ENDCOMMENT
NEURON {
POINT_PROCESS psd
: RANGE P, V, A, B, D, W
RANGE V
USEION ca READ cai
GLOBAL pHC, pHN, aHC, aHN, thetav, sigmav, thetad, sigmad, thetab, sigmab, taup, kp, taua, kd, tauv, alphav, taud, alphad, taub, alphab, cp, cd, tauw, p, alphaw, d, betaw
}
UNITS {
(molar) = (/liter)
(mM) = (millimolar)
(uM) = (micromolar)
}
PARAMETER {
pHC = 4 (uM)
pHN = 4
aHC = 0.6 (uM)
aHN = 3
thetav = 2 (uM)
sigmav = -0.05
thetad = 2.6
sigmad = -0.01
thetab = 0.55
sigmab = -0.02
taup = 500 (ms)
kp = -0.1
taua = 5 (ms)
kd = -0.002
tauv = 10 (ms)
alphav = 1.0
taud = 250 (ms)
alphad = 1.0
taub = 40 (ms)
alphab = 5.0
cp = 5
cd = 4
tauw = 500 (ms)
p = 0.3
alphaw = 0.8
d = 0.01
betaw = 0.6
}
ASSIGNED {
cai (mM)
}
: STATE { P V A B D W }
STATE { V }
BREAKPOINT {
SOLVE state METHOD cnexp
}
DERIVATIVE state {
: P' = (10 * hill( cai,pHC, pHN) - cp*A*P )/taup
V' = (alphav * sigm(cai,thetav,sigmav) - V)/tauv
: A' = ( hill( cai,aHC, aHN) - A )/taua
: B' = (alphab * sigmoid(A, thetab,sigmab) - B - cd*B*V)/taub
: D' = (alphad * sigmoid(D, thetad,sigmad) - D )/taud
: W' = (alphaw * sigmoid(P, p, kp) - betaw * sigmoid(D, d, kd) - W)/tauw
}
INITIAL {
: P = 0
V = 0
: A = 0
: B = 0
: D = 0
: W = 1
}
FUNCTION hill(x (mM), HC (uM), HN) {
hill = 1 / ( (x/((0.001)*HC))^HN + 1 )
}
FUNCTION sigm(x, theta, sigma) {
sigm = 1 / ( 1 + exp((x-theta)/sigma) )
}
```