create dend1
create dend2
create SpineStalk
create SpineHead
create Soma


dend1 {L = 630 // lambda=1261 for 1 um Dia
 diam = 1 
 nseg=315
 insert pas g_pas=1/10000 e_pas=-70 Ra=150 cm=1}

dend2 {L = 630 diam = 1
nseg=315
insert pas g_pas=1/10000 e_pas=-70 Ra=150 cm=1}

SpineStalk {L=1 diam= 0.1
nseg=10
insert pas g_pas=1/10000 e_pas=-70 Ra=150 cm=1}

SpineHead {L=0.7 diam=1.0
nseg=1
insert canrgc eca=55
insert calrgc eca=45
insert pas g_pas=1/10000 e_pas=-70 Ra=150 cm=1}

Soma {L=25 diam=25
nseg=3
insert pas g_pas=1/10000 e_pas=-70 Ra=150 cm=1}

connect Soma(1),dend1(0)
connect dend1(1), dend2(0)
connect SpineStalk(0), dend1(1)
connect SpineHead(0), SpineStalk(1)

forall{
// nseg=int(L/10)+1
// insert pas
// g_pas=1/10000
// Ra=110
global_ra=150
e_pas=-70
// eca=55
// eca=45
}





/*------------------------------------*/
// SPINEHEAD STMULATION
 access SpineHead
 objectvar salpha // IChead

// IChead= new IClamp(.5)
// IChead.del=10
// IChead.dur=10
// IChead.amp=.01

 salpha = new AlphaSynapse(0.50)
 salpha.onset = 10
 salpha.tau = 10
 salpha.gmax = 0.002

// ------------------------------------
// DENDRITE STIMULATION
 access dend1
 objectvar  dalpha // ddalpha,ICdend

// ICdend = new IClamp(1)

 dalpha = new AlphaSynapse(1)
 dalpha.onset = 10
 dalpha.tau = 10
 dalpha.gmax = 0.002

// ddalpha =  new AlphaSynapse(0.99)
// ddalpha.onset = 10
// ddalpha.tau = 10
// ddalpha.gmax = 0.01

// ------------------------------------
// objectvar vbox1, tempobj1
// proc show_alphas(){
      //  vbox1 = new VBox()
      // vbox1.intercept(1)

      //  xpanel("")
     //   xvalue("Spine Head  gmax","salpha.gmax", 1, "")
     //   xvalue("Spine Base (a) gmax", "dalpha.gmax", 1, "")
     //   xvalue("Dendrite (b) gmax", "ddalpha.gmax", 1, "")
     //   xvalue("Ca++ channel density","CABAR", 1,"set_cabar()") 

      //  xpanel()
      //  vbox1.intercept(0)
       // vbox1.map
//}

// access soma
access SpineHead

tstop=100
vlow=-80
vhigh=40
dt=0.01
Vinit=-70

 // load_proc ("nrnmainmenu")
 // Load_file("spine.ses")