-
Notifications
You must be signed in to change notification settings - Fork 4
/
mosinit.py
53 lines (46 loc) · 1.27 KB
/
mosinit.py
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
if __name__ == "__main__":
import sys
import os
import string
from neuron import *
h("strdef simname, allfiles, simfiles, output_file, datestr, uname, osname, comment")
h.simname=simname = "mtlhpc"
h.allfiles=allfiles = "geom.hoc pyinit.py geom.py network.py params.py run.py"
h.simfiles=simfiles = "pyinit.py geom.py network.py params.py run.py"
h("runnum=1")
runnum = 1.0
h.datestr=datestr = "11may20"
h.output_file=output_file = "data/11may20.05"
h.uname=uname = "x86_64"
h.osname=osname="linux"
h("templates_loaded=0")
templates_loaded=0
h("xwindows=1.0")
xwindows = 1.0
h.xopen("nrnoc.hoc")
h.xopen("init.hoc")
from pyinit import *
from geom import *
from networkmsj import *
from params import *
from run import *
# setup washin,washout
"""
import run as Run
Run.olmWash = [0, 1]
Run.basWash = [1, 1]
Run.pyrWashA = [1, 1]
Run.pyrWashB = [1, 1]
Run.washinT = 1e3
Run.washoutT = 2e3
Run.fiwash = h.FInitializeHandler(1,Run.setwash)
"""
h.tstop = 5e3
h.run()
net.rasterplot()
net.calc_lfp()
net.pravgrates()
myg = h.Graph()
net.vlfp.plot(myg,h.dt)
myg.exec_menu("View = plot")
myg.exec_menu("New Axis")