proc celldef() {
  topol()
  subsets()
  geom()
  biophys()
}

create soma, dend, axon

proc topol() { local i
  // connect dend(0), soma(1)
  basic_shape()
}
proc basic_shape() {
  soma {pt3dclear() pt3dadd(0, 0, 0, 1) pt3dadd(15, 0, 0, 1)}
  dend {pt3dclear() pt3dadd(15, 0, 0, 1) pt3dadd(515, 0, 0, 1)}

}

objref all
proc subsets() {
  all = new SectionList()
    soma all.append()
    dend all.append()
}

proc geom() {
  forsec all {  }
  soma {  L = 20  diam = 20  nseg = 1}
  dend {  L = 500  diam = 1  nseg = 99}
}

proc biophys() {
  forsec all {
    cm = 1
    Ra = 35.4
  }

  dend {
     insert VDpas
        gmax_VDpas = 0.0047  
        gmin_VDpas = 0.00060  
        e50_VDpas = -30
        s_VDpas = -6
        f_VDpas = 0.096 
    }
}

access soma

celldef()