I am attempting to run the model Meuth et al 2005 (https://senselab.med.yale.edu/ModelDB/S ... del=121600). However, when I initially run the model (version 4 for instance) the plot only displays data for the first time step and nothing else. On the MacOS I also get the following notifications in Terminal:
Jun 26 10:04:04 dh130141./domain removed/.ca nrniv[2109] <Error>: kCGErrorInvalidConnection: CGSGetWindowTags: Invalid connection
Jun 26 10:04:04 dh130141./domain removed/.ca nrniv[2109] <Error>: kCGErrorFailure: Set a breakpoint @ CGErrorBreakpoint() to catch errors as they are logged.
/domain removed/ - inserted myself for privacy; my university hosting site was displayed here.
I found that I can overcome this error by running the simulation once, getting the errors, then reassigning the variable celsius in Terminal, and running the simulation again. After doing so the simulation will run completely and the errors will not occur again.
This is a fairly annoying and time consuming fix, and I feel that I must be doing something wrong or there is an error in code. Is anyone able to give me an indication on what is going on?
Thank you.
PS. To save time, here is the code I am referencing (I've added a few stimuli) from Meuth et al 2005:
load_file("nrngui.hoc")
v_init = -70
tstop = 1000
create soma
access soma
{L = 30 diam = 20}
insert leak
insert HH
insert iA
insert iC
insert iH
insert iT
insert iNaP
insert iL
celsius = 35
nao = 145
{ki = 135 ko = 3.1}
{cli = 7 clo = 120}
pk_leak=8e-6
gkbar_iA=0.0055
gkbar_iC=0.00018
ghbar_iH=0.0006
eh= -40
pcabar_iT=0.00045
g_iNaP=17e-5
objref stimulus
stimulus = new IClamp(0.5)
stimulus.amp = 0.1
stimulus.del = 200
stimulus.dur = 30
objref stimulus2
stimulus2 = new IClamp(0.5)
stimulus2.amp = 0.2
stimulus2.del = 400
stimulus2.dur = 30
objref stimulus3
stimulus3 = new IClamp(0.5)
stimulus3.amp = 0.3
stimulus3.del = 600
stimulus3.dur = 30
psection()