I'm trying to make a mechanism to introduce a 'temperature clamp' in my cell model. As a first approach, I wrote a NMODL code similar to the Example 3 of the Hines and Carnevale NMODL manual.
Code: Select all
TITLE Mechanism for a temperature clamp generating square T pulse.
NEURON {
SUFFIX TClamp
RANGE del, dur, temp1, temp2
}
PARAMETER {
del = 100 (ms)
dur = 300 (ms)
temp1 = 35 (degC)
temp2 = 20 (degC)
}
ASSIGNED {
celsius (degC)
}
BREAKPOINT {
at_time (del)
at_time (del+dur)
if (t < del +dur && t > del) {celsius = temp2} else {celsius = temp1}
}
celsius is a special NEURON variable that should not be assigned a value in a model description file
Can I change the temperature of the system as a function of simulation time?
Thank you in advance.