www.neuron.yale.edu

Code: Select all

START = 0 // first time at which field is calculated
INTERVAL = 1000 // between successive calculations of field
objref fih
fih = new FInitializeHandler("initcalc()")
proc initcalc() {
  cvode.event(START, "docalc()")
// print "launched event to trigger field calc at ", START
}
proc docalc() {
// print "calculating field at t = ", t
  system("command_that_computes_field")
  cvode.event(t + INTERVAL, "docalc()")
// print "launched event to trigger next calc"
}

Code: Select all

strdef cmdstr
proc docalc() {
// print "calculating field at t = ", t
  sprint(cmdstr, "calcfield field%d", t) // assumes program name is calcfield
    // and argument is name of output file
  system(cmdstr)
  cvode.event(t + INTERVAL, "docalc()")
// print "launched event to trigger next calc"
}

Code: Select all

def pysave_imembrane(): #Hines' pybs.mod method calls this procedure at every time step
  
  for gid in cells:
    ... #This code saves all the currents for all the cells to binary files
 
     #Want to call executable  "a.out" every 1000 ms  
     if ((index1 % 1000) == 0 and index1>0 and gid==0):
       os.system("/home/cat/neuron/eaps/a.out")