wormsim 0.0.1__py3-none-any.whl

wormsim/components/model/AWAR.hoc

@@ -0,0 +1,165 @@
+begintemplate AWAR
+
+public init, topol, basic_shape, subsets, geom, memb
+public x, y, z, position
+
+
+//  Some methods for referencing the cells
+
+public toString, netInfo, set_initial_v, set_initial_ion_properties, synlist
+
+
+// Section: Soma, id: Soma, parent: $s.parent
+create Soma
+public Soma
+// Section: Axon, id: Axon, parent: Soma
+create Axon
+public Axon
+// Section: Dendrite2, id: Dendrite2, parent: Soma
+create Dendrite2
+public Dendrite2
+
+// Group: {name=soma_group, groups=[], sections=[Soma]}
+objref soma_group
+// Group: {name=axon_group, groups=[], sections=[Axon]}
+objref axon_group
+// Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+objref dendrite_group
+// Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+objref all
+
+objref synlist
+
+proc init() {
+    topol()
+    geom_nseg()
+    subsets()
+    subsets_inhomogeneous()
+    biophys()
+    biophys_inhomogeneous()
+    synlist = new List()
+    x = y = z = 0
+
+}
+
+
+
+
+
+proc topol() {
+    basic_shape()
+    connect Axon(0), Soma(1.0)
+    connect Dendrite2(0), Soma(1.0)
+
+}
+
+
+proc basic_shape() {
+
+    // Section: Soma, id: Soma, parent: $s.parent (entry count 0/100)
+    Soma { pt3dclear() }
+    Soma { pt3dadd( -8.2, -265.875, 42.75, 2.8142495 ) }
+    Soma { pt3dadd( -8.2, -263.0607505, 42.75, 2.8142495 ) }
+    // Section: Axon, id: Axon, parent: Soma (entry count 2/100)
+    Axon { pt3dclear() }
+    Axon { pt3dadd( -8.2, -265.875, 42.75, 0.5830952 ) }
+    Axon { pt3dadd( -8.2, -265.15, 42.75, 0.5830952 ) }
+    Axon { pt3dadd( -8.150001, -264.94998, 42.399998, 0.46904156 ) }
+    Axon { pt3dadd( -8.0, -264.7, 39.7, 0.56568545 ) }
+    Axon { pt3dadd( -8.0, -265.4, 30.75, 0.5744563 ) }
+    Axon { pt3dadd( -7.825, -265.59998, 29.800001, 0.5123476 ) }
+    Axon { pt3dadd( -3.25, -267.6, 26.499998, 0.5744563 ) }
+    Axon { pt3dadd( -2.5, -268.95, 27.5, 0.53851646 ) }
+    Axon { pt3dadd( -2.5, -271.25, 30.5, 0.53851646 ) }
+    Axon { pt3dadd( -2.55, -272.825, 35.975, 0.5612486 ) }
+    Axon { pt3dadd( -3.3, -273.05, 37.6, 0.41231057 ) }
+    Axon { pt3dadd( -4.9, -273.15, 38.125, 0.39051247 ) }
+    Axon { pt3dadd( -7.175, -273.42502, 39.825, 0.295804 ) }
+    Axon { pt3dadd( -8.725, -273.775, 42.449997, 0.44158804 ) }
+    Axon { pt3dadd( -9.25, -274.2, 45.5, 0.5 ) }
+    Axon { pt3dadd( -8.7, -274.65, 48.525, 0.4387482 ) }
+    Axon { pt3dadd( -7.2, -275.0, 51.1, 0.48989794 ) }
+    Axon { pt3dadd( -4.8500004, -275.3, 52.824997, 0.48218253 ) }
+    Axon { pt3dadd( -2.2, -275.35, 53.399998, 0.5 ) }
+    // Section: Dendrite2, id: Dendrite2, parent: Soma (entry count 21/100)
+    Dendrite2 { pt3dclear() }
+    Dendrite2 { pt3dadd( -8.2, -265.875, 42.75, 0.53851646 ) }
+    Dendrite2 { pt3dadd( -8.2, -266.59998, 42.75, 0.53851646 ) }
+    Dendrite2 { pt3dadd( -8.2, -266.95, 42.35, 0.50990194 ) }
+    Dendrite2 { pt3dadd( -8.5, -269.4, 42.899998, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -9.700001, -274.1, 45.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -10.5, -280.69998, 47.9, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -11.7, -289.8, 51.1, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -12.2, -293.19998, 52.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -12.2, -296.6, 52.8, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -11.8, -299.8, 53.3, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -10.400001, -304.6, 54.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -9.3, -320.15, 56.0, 0.5744563 ) }
+    Dendrite2 { pt3dadd( -7.2999997, -331.95, 57.200005, 0.5744563 ) }
+    Dendrite2 { pt3dadd( -6.5, -337.15, 57.8, 0.5744563 ) }
+    Dendrite2 { pt3dadd( -5.6000004, -342.8, 58.600002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -5.0, -348.6, 58.899998, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -5.0, -349.44998, 58.999996, 0.5744563 ) }
+
+}
+
+
+proc subsets() { local i
+
+    // Group: {name=soma_group, groups=[], sections=[Soma]}
+    soma_group = new SectionList()
+    Soma soma_group.append()
+
+    // Group: {name=axon_group, groups=[], sections=[Axon]}
+    axon_group = new SectionList()
+    Axon axon_group.append()
+
+    // Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+    dendrite_group = new SectionList()
+    Dendrite2 dendrite_group.append()
+
+    // Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+	all = new SectionList()
+	Soma all.append()
+	Axon all.append()
+	Dendrite2 all.append()
+}
+
+
+proc subsets_inhomogeneous() {
+}
+
+proc biophys() {
+
+}
+
+func H() { // Heaviside function
+    if ($1>=0) return 1
+    if ($1==0) return 0.5
+    return 0
+}
+
+proc biophys_inhomogeneous() {
+}
+
+proc geom_nseg() {
+    // All sections not mentioned here have nseg = 1
+	Soma nseg = 2
+	Axon nseg = 19
+	Dendrite2 nseg = 17
+}
+
+
+proc position() { local i
+    forsec all {
+        for i = 0, n3d()-1 {
+            pt3dchange(i, $1+x3d(i), $2+y3d(i), $3+z3d(i), diam3d(i))
+        }
+    }
+    x = $1  y = $2  z = $3
+}
+
+
+endtemplate AWAR
+
+

wormsim/components/model/AWBL.hoc

@@ -0,0 +1,165 @@
+begintemplate AWBL
+
+public init, topol, basic_shape, subsets, geom, memb
+public x, y, z, position
+
+
+//  Some methods for referencing the cells
+
+public toString, netInfo, set_initial_v, set_initial_ion_properties, synlist
+
+
+// Section: Soma, id: Soma, parent: $s.parent
+create Soma
+public Soma
+// Section: Axon, id: Axon, parent: Soma
+create Axon
+public Axon
+// Section: Dendrite2, id: Dendrite2, parent: Soma
+create Dendrite2
+public Dendrite2
+
+// Group: {name=soma_group, groups=[], sections=[Soma]}
+objref soma_group
+// Group: {name=axon_group, groups=[], sections=[Axon]}
+objref axon_group
+// Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+objref dendrite_group
+// Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+objref all
+
+objref synlist
+
+proc init() {
+    topol()
+    geom_nseg()
+    subsets()
+    subsets_inhomogeneous()
+    biophys()
+    biophys_inhomogeneous()
+    synlist = new List()
+    x = y = z = 0
+
+}
+
+
+
+
+
+proc topol() {
+    basic_shape()
+    connect Axon(0), Soma(1.0)
+    connect Dendrite2(0), Soma(1.0)
+
+}
+
+
+proc basic_shape() {
+
+    // Section: Soma, id: Soma, parent: $s.parent (entry count 0/100)
+    Soma { pt3dclear() }
+    Soma { pt3dadd( 5.5, -266.225, 43.1, 2.523886 ) }
+    Soma { pt3dadd( 5.5, -263.701114, 43.1, 2.523886 ) }
+    // Section: Axon, id: Axon, parent: Soma (entry count 2/100)
+    Axon { pt3dclear() }
+    Axon { pt3dadd( 5.5, -266.225, 43.1, 0.34641016 ) }
+    Axon { pt3dadd( 5.5, -265.59998, 43.1, 0.34641016 ) }
+    Axon { pt3dadd( 5.5, -265.35, 42.649998, 0.37416574 ) }
+    Axon { pt3dadd( 5.55, -265.05002, 39.5, 0.42426404 ) }
+    Axon { pt3dadd( 4.0750003, -265.4, 30.4, 0.5678908 ) }
+    Axon { pt3dadd( 3.5500002, -265.59998, 29.800001, 0.53851646 ) }
+    Axon { pt3dadd( -1.0, -267.6, 26.499998, 0.6 ) }
+    Axon { pt3dadd( -1.8, -268.95, 27.5, 0.53851646 ) }
+    Axon { pt3dadd( -1.8, -271.25, 30.5, 0.53851646 ) }
+    Axon { pt3dadd( -1.75, -272.84998, 36.0, 0.42426404 ) }
+    Axon { pt3dadd( -1.0, -273.05, 37.3, 0.41231057 ) }
+    Axon { pt3dadd( 0.65000004, -273.15, 38.125, 0.40311286 ) }
+    Axon { pt3dadd( 2.8999999, -273.42502, 39.825, 0.43011624 ) }
+    Axon { pt3dadd( 4.475, -273.775, 42.449997, 0.44158804 ) }
+    Axon { pt3dadd( 4.975, -274.2, 45.5, 0.42720017 ) }
+    Axon { pt3dadd( 4.4500003, -274.65, 48.525, 0.4924429 ) }
+    Axon { pt3dadd( 2.8999999, -275.0, 51.1, 0.48989794 ) }
+    Axon { pt3dadd( 0.6, -275.3, 52.824997, 0.5123476 ) }
+    Axon { pt3dadd( -2.1, -275.35, 53.399998, 0.5 ) }
+    // Section: Dendrite2, id: Dendrite2, parent: Soma (entry count 21/100)
+    Dendrite2 { pt3dclear() }
+    Dendrite2 { pt3dadd( 5.5, -266.225, 43.1, 0.53851646 ) }
+    Dendrite2 { pt3dadd( 5.5, -266.85, 43.1, 0.53851646 ) }
+    Dendrite2 { pt3dadd( 5.5, -267.4, 42.899998, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 5.45, -269.7, 43.45, 0.42426404 ) }
+    Dendrite2 { pt3dadd( 5.45, -274.6, 45.050003, 0.42426404 ) }
+    Dendrite2 { pt3dadd( 5.5, -280.3, 47.4, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 7.4, -289.8, 51.1, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 7.9000006, -293.19998, 52.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 7.9000006, -296.6, 52.8, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 7.5000005, -299.8, 53.3, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 6.2, -304.6, 54.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 5.0, -320.15, 56.0, 0.5744563 ) }
+    Dendrite2 { pt3dadd( 3.1, -331.95, 57.200005, 0.5744563 ) }
+    Dendrite2 { pt3dadd( 2.3, -337.15, 57.8, 0.5744563 ) }
+    Dendrite2 { pt3dadd( 1.3000001, -342.8, 58.600002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 0.70000005, -348.6, 58.899998, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 0.70000005, -349.44998, 58.999996, 0.5744563 ) }
+
+}
+
+
+proc subsets() { local i
+
+    // Group: {name=soma_group, groups=[], sections=[Soma]}
+    soma_group = new SectionList()
+    Soma soma_group.append()
+
+    // Group: {name=axon_group, groups=[], sections=[Axon]}
+    axon_group = new SectionList()
+    Axon axon_group.append()
+
+    // Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+    dendrite_group = new SectionList()
+    Dendrite2 dendrite_group.append()
+
+    // Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+	all = new SectionList()
+	Soma all.append()
+	Axon all.append()
+	Dendrite2 all.append()
+}
+
+
+proc subsets_inhomogeneous() {
+}
+
+proc biophys() {
+
+}
+
+func H() { // Heaviside function
+    if ($1>=0) return 1
+    if ($1==0) return 0.5
+    return 0
+}
+
+proc biophys_inhomogeneous() {
+}
+
+proc geom_nseg() {
+    // All sections not mentioned here have nseg = 1
+	Soma nseg = 2
+	Axon nseg = 19
+	Dendrite2 nseg = 17
+}
+
+
+proc position() { local i
+    forsec all {
+        for i = 0, n3d()-1 {
+            pt3dchange(i, $1+x3d(i), $2+y3d(i), $3+z3d(i), diam3d(i))
+        }
+    }
+    x = $1  y = $2  z = $3
+}
+
+
+endtemplate AWBL
+
+

wormsim/components/model/AWBR.hoc

@@ -0,0 +1,165 @@
+begintemplate AWBR
+
+public init, topol, basic_shape, subsets, geom, memb
+public x, y, z, position
+
+
+//  Some methods for referencing the cells
+
+public toString, netInfo, set_initial_v, set_initial_ion_properties, synlist
+
+
+// Section: Soma, id: Soma, parent: $s.parent
+create Soma
+public Soma
+// Section: Axon, id: Axon, parent: Soma
+create Axon
+public Axon
+// Section: Dendrite2, id: Dendrite2, parent: Soma
+create Dendrite2
+public Dendrite2
+
+// Group: {name=soma_group, groups=[], sections=[Soma]}
+objref soma_group
+// Group: {name=axon_group, groups=[], sections=[Axon]}
+objref axon_group
+// Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+objref dendrite_group
+// Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+objref all
+
+objref synlist
+
+proc init() {
+    topol()
+    geom_nseg()
+    subsets()
+    subsets_inhomogeneous()
+    biophys()
+    biophys_inhomogeneous()
+    synlist = new List()
+    x = y = z = 0
+
+}
+
+
+
+
+
+proc topol() {
+    basic_shape()
+    connect Axon(0), Soma(1.0)
+    connect Dendrite2(0), Soma(1.0)
+
+}
+
+
+proc basic_shape() {
+
+    // Section: Soma, id: Soma, parent: $s.parent (entry count 0/100)
+    Soma { pt3dclear() }
+    Soma { pt3dadd( -9.75, -266.225, 43.1, 2.523886 ) }
+    Soma { pt3dadd( -9.75, -263.701114, 43.1, 2.523886 ) }
+    // Section: Axon, id: Axon, parent: Soma (entry count 2/100)
+    Axon { pt3dclear() }
+    Axon { pt3dadd( -9.75, -266.225, 43.1, 0.41231057 ) }
+    Axon { pt3dadd( -9.75, -265.59998, 43.1, 0.41231057 ) }
+    Axon { pt3dadd( -9.75, -265.35, 42.649998, 0.4358899 ) }
+    Axon { pt3dadd( -9.825, -265.05002, 39.5, 0.5408327 ) }
+    Axon { pt3dadd( -8.349999, -265.4, 30.4, 0.53851646 ) }
+    Axon { pt3dadd( -7.825, -265.59998, 29.800001, 0.5123476 ) }
+    Axon { pt3dadd( -3.25, -267.6, 26.499998, 0.5744563 ) }
+    Axon { pt3dadd( -2.5, -268.95, 27.5, 0.53851646 ) }
+    Axon { pt3dadd( -2.5, -271.25, 30.5, 0.53851646 ) }
+    Axon { pt3dadd( -2.525, -272.84998, 36.0, 0.33541018 ) }
+    Axon { pt3dadd( -3.2749999, -273.05, 37.3, 0.4387482 ) }
+    Axon { pt3dadd( -4.9, -273.15, 38.125, 0.39051247 ) }
+    Axon { pt3dadd( -7.175, -273.42502, 39.825, 0.295804 ) }
+    Axon { pt3dadd( -8.725, -273.775, 42.449997, 0.44158804 ) }
+    Axon { pt3dadd( -9.25, -274.2, 45.5, 0.5 ) }
+    Axon { pt3dadd( -8.7, -274.65, 48.525, 0.4387482 ) }
+    Axon { pt3dadd( -7.2, -275.0, 51.1, 0.48989794 ) }
+    Axon { pt3dadd( -4.8500004, -275.3, 52.824997, 0.48218253 ) }
+    Axon { pt3dadd( -2.2, -275.35, 53.399998, 0.5 ) }
+    // Section: Dendrite2, id: Dendrite2, parent: Soma (entry count 21/100)
+    Dendrite2 { pt3dclear() }
+    Dendrite2 { pt3dadd( -9.75, -266.225, 43.1, 0.46904156 ) }
+    Dendrite2 { pt3dadd( -9.75, -266.85, 43.1, 0.46904156 ) }
+    Dendrite2 { pt3dadd( -9.75, -267.4, 42.899998, 0.5 ) }
+    Dendrite2 { pt3dadd( -9.75, -269.7, 43.45, 0.42426404 ) }
+    Dendrite2 { pt3dadd( -9.75, -274.6, 45.050003, 0.42426404 ) }
+    Dendrite2 { pt3dadd( -9.8, -280.3, 47.4, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -11.7, -289.8, 51.1, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -12.2, -293.19998, 52.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -12.2, -296.6, 52.8, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -11.8, -299.8, 53.3, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -10.400001, -304.6, 54.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -9.3, -320.15, 56.0, 0.5744563 ) }
+    Dendrite2 { pt3dadd( -7.2999997, -331.95, 57.200005, 0.5744563 ) }
+    Dendrite2 { pt3dadd( -6.5, -337.15, 57.8, 0.5744563 ) }
+    Dendrite2 { pt3dadd( -5.6000004, -342.8, 58.600002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -5.0, -348.6, 58.899998, 0.56568545 ) }
+    Dendrite2 { pt3dadd( -5.0, -349.44998, 58.999996, 0.5744563 ) }
+
+}
+
+
+proc subsets() { local i
+
+    // Group: {name=soma_group, groups=[], sections=[Soma]}
+    soma_group = new SectionList()
+    Soma soma_group.append()
+
+    // Group: {name=axon_group, groups=[], sections=[Axon]}
+    axon_group = new SectionList()
+    Axon axon_group.append()
+
+    // Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+    dendrite_group = new SectionList()
+    Dendrite2 dendrite_group.append()
+
+    // Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+	all = new SectionList()
+	Soma all.append()
+	Axon all.append()
+	Dendrite2 all.append()
+}
+
+
+proc subsets_inhomogeneous() {
+}
+
+proc biophys() {
+
+}
+
+func H() { // Heaviside function
+    if ($1>=0) return 1
+    if ($1==0) return 0.5
+    return 0
+}
+
+proc biophys_inhomogeneous() {
+}
+
+proc geom_nseg() {
+    // All sections not mentioned here have nseg = 1
+	Soma nseg = 2
+	Axon nseg = 19
+	Dendrite2 nseg = 17
+}
+
+
+proc position() { local i
+    forsec all {
+        for i = 0, n3d()-1 {
+            pt3dchange(i, $1+x3d(i), $2+y3d(i), $3+z3d(i), diam3d(i))
+        }
+    }
+    x = $1  y = $2  z = $3
+}
+
+
+endtemplate AWBR
+
+

wormsim/components/model/AWCL.hoc

@@ -0,0 +1,168 @@
+begintemplate AWCL
+
+public init, topol, basic_shape, subsets, geom, memb
+public x, y, z, position
+
+
+//  Some methods for referencing the cells
+
+public toString, netInfo, set_initial_v, set_initial_ion_properties, synlist
+
+
+// Section: Soma, id: Soma, parent: $s.parent
+create Soma
+public Soma
+// Section: Axon, id: Axon, parent: Soma
+create Axon
+public Axon
+// Section: Dendrite2, id: Dendrite2, parent: Soma
+create Dendrite2
+public Dendrite2
+
+// Group: {name=soma_group, groups=[], sections=[Soma]}
+objref soma_group
+// Group: {name=axon_group, groups=[], sections=[Axon]}
+objref axon_group
+// Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+objref dendrite_group
+// Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+objref all
+
+objref synlist
+
+proc init() {
+    topol()
+    geom_nseg()
+    subsets()
+    subsets_inhomogeneous()
+    biophys()
+    biophys_inhomogeneous()
+    synlist = new List()
+    x = y = z = 0
+
+}
+
+
+
+
+
+proc topol() {
+    basic_shape()
+    connect Axon(0), Soma(1.0)
+    connect Dendrite2(0), Soma(1.0)
+
+}
+
+
+proc basic_shape() {
+
+    // Section: Soma, id: Soma, parent: $s.parent (entry count 0/100)
+    Soma { pt3dclear() }
+    Soma { pt3dadd( 3.8, -267.95, 38.95, 2.7676706 ) }
+    Soma { pt3dadd( 3.8, -265.1823294, 38.95, 2.7676706 ) }
+    // Section: Axon, id: Axon, parent: Soma (entry count 2/100)
+    Axon { pt3dclear() }
+    Axon { pt3dadd( 3.8, -267.95, 38.95, 0.6 ) }
+    Axon { pt3dadd( 3.8, -267.5, 38.100002, 0.6 ) }
+    Axon { pt3dadd( 3.85, -267.25, 37.3, 0.46904156 ) }
+    Axon { pt3dadd( 3.9, -266.35, 34.0, 0.5 ) }
+    Axon { pt3dadd( 3.85, -265.55, 30.25, 0.4358899 ) }
+    Axon { pt3dadd( 2.6499999, -266.1, 28.85, 0.50990194 ) }
+    Axon { pt3dadd( -1.05, -267.575, 26.4, 0.5408327 ) }
+    Axon { pt3dadd( -1.85, -268.30002, 26.499998, 0.41231057 ) }
+    Axon { pt3dadd( -1.85, -271.35, 30.4, 0.46904156 ) }
+    Axon { pt3dadd( -1.85, -272.75, 36.0, 0.42426404 ) }
+    Axon { pt3dadd( -1.1, -273.05, 37.6, 0.41231057 ) }
+    Axon { pt3dadd( 0.65000004, -273.15, 38.125, 0.35000002 ) }
+    Axon { pt3dadd( 2.8999999, -273.42502, 39.825, 0.43011624 ) }
+    Axon { pt3dadd( 4.475, -273.775, 42.449997, 0.44158804 ) }
+    Axon { pt3dadd( 4.975, -274.2, 45.5, 0.42720017 ) }
+    Axon { pt3dadd( 4.4500003, -274.65, 48.525, 0.4924429 ) }
+    Axon { pt3dadd( 2.8999999, -275.0, 51.1, 0.48989794 ) }
+    Axon { pt3dadd( 0.6, -275.3, 52.824997, 0.5123476 ) }
+    Axon { pt3dadd( -2.1, -275.35, 53.399998, 0.5 ) }
+    Axon { pt3dadd( -4.8500004, -275.3, 52.824997, 0.33541018 ) }
+    Axon { pt3dadd( -7.1000004, -275.0, 51.1, 0.34641016 ) }
+    Axon { pt3dadd( -8.65, -274.65, 48.525, 0.42720017 ) }
+    // Section: Dendrite2, id: Dendrite2, parent: Soma (entry count 24/100)
+    Dendrite2 { pt3dclear() }
+    Dendrite2 { pt3dadd( 3.8, -267.95, 38.95, 0.6 ) }
+    Dendrite2 { pt3dadd( 3.8, -268.4, 39.8, 0.6 ) }
+    Dendrite2 { pt3dadd( 3.85, -268.9, 40.149998, 0.46904156 ) }
+    Dendrite2 { pt3dadd( 5.25, -271.8, 41.55, 0.46904156 ) }
+    Dendrite2 { pt3dadd( 6.0, -278.4, 46.9, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 6.0, -282.55, 48.45, 0.50990194 ) }
+    Dendrite2 { pt3dadd( 7.4, -289.05002, 51.050003, 0.50990194 ) }
+    Dendrite2 { pt3dadd( 7.9000006, -293.19998, 52.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 7.9000006, -296.6, 52.8, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 7.5000005, -299.8, 53.3, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 6.2, -304.6, 54.100002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 5.0, -320.15, 56.0, 0.5744563 ) }
+    Dendrite2 { pt3dadd( 3.1, -331.95, 57.200005, 0.5744563 ) }
+    Dendrite2 { pt3dadd( 2.3, -337.15, 57.8, 0.5744563 ) }
+    Dendrite2 { pt3dadd( 1.3000001, -342.8, 58.600002, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 0.70000005, -348.6, 58.899998, 0.56568545 ) }
+    Dendrite2 { pt3dadd( 0.70000005, -349.44998, 58.999996, 0.5744563 ) }
+
+}
+
+
+proc subsets() { local i
+
+    // Group: {name=soma_group, groups=[], sections=[Soma]}
+    soma_group = new SectionList()
+    Soma soma_group.append()
+
+    // Group: {name=axon_group, groups=[], sections=[Axon]}
+    axon_group = new SectionList()
+    Axon axon_group.append()
+
+    // Group: {name=dendrite_group, groups=[], sections=[Dendrite2]}
+    dendrite_group = new SectionList()
+    Dendrite2 dendrite_group.append()
+
+    // Group: {name=all, sections=[Soma, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Axon, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2, Dendrite2]}
+	all = new SectionList()
+	Soma all.append()
+	Axon all.append()
+	Dendrite2 all.append()
+}
+
+
+proc subsets_inhomogeneous() {
+}
+
+proc biophys() {
+
+}
+
+func H() { // Heaviside function
+    if ($1>=0) return 1
+    if ($1==0) return 0.5
+    return 0
+}
+
+proc biophys_inhomogeneous() {
+}
+
+proc geom_nseg() {
+    // All sections not mentioned here have nseg = 1
+	Soma nseg = 2
+	Axon nseg = 22
+	Dendrite2 nseg = 17
+}
+
+
+proc position() { local i
+    forsec all {
+        for i = 0, n3d()-1 {
+            pt3dchange(i, $1+x3d(i), $2+y3d(i), $3+z3d(i), diam3d(i))
+        }
+    }
+    x = $1  y = $2  z = $3
+}
+
+
+endtemplate AWCL
+
+