load_proc("nrnmainmenu")


// --------------------------------------------------------------
// redefine some things in stdrun.hoc
// --------------------------------------------------------------

tstop = 50
steps_per_ms = 40
dt = 0.025
secondorder = 1

// --------------------------------------------------------------
// cell geometry
// --------------------------------------------------------------

create soma 
soma { nseg = 1 L = 30 diam = 30 }
create soma1 
soma1 { nseg = 1 L = 30 diam = 30 }

//RNG
objref rnd
soma1 rnd = new RNG(0.5)
rnd.seed = 2365877 



// --------------------------------------------------------------
// passive & active membrane 
// --------------------------------------------------------------

ra        = 200
global_ra = ra
rm        = 40000
c_m       = 0.75
v_init    = -70
celsius = 23

Ek = -90
Ena = 60
            

gna_dend = 30
gna_node = 30000
gk = 500
frac_gk_apical = 0
gk_node = 0


proc install_channels() {
  
  // passive
  forall {
    insert pas
    Ra = ra
    cm = c_m 
    g_pas = 1/rm
    e_pas = v_init
  }

  // Na+ channels
/*
  forall insert na3
  soma.gbar_na3 = gna_dend

  forall if (ismembrane("na3")) ena = Ena     // per S&S ena = 90
*/
  // K+ channels

  soma insert kd3
  soma.gbar_kd3 = gk


  soma1 insert sk
  soma1.gbar_sk = gk
  soma1 { D_sk = int(gbar_sk/gama_sk) } 

  forall if (ismembrane("kd3")) ek = Ek      // per S&S ek = -100
  forall if (ismembrane("sk")) ek = Ek         // per S&S ek = -100
}

install_channels()


// --------------------------------------------------------------
// stimulus
// --------------------------------------------------------------

objectvar st
soma  st = new IClamp(0.5)
st.amp = 0.2
st.del = 3
st.dur = 40

objectvar st1
soma1  st1 = new IClamp(0.5)
st1.amp = 0.2
st1.del = 3
st1.dur = 40


// --------------------------------------------------------------
// create useful graphs & panels
// --------------------------------------------------------------

access soma
nrnmainmenu()
nrncontrolmenu()
newPlotV()
graphItem.addvar("soma1.v(.5)",7,1)
init()
print n_kd3,soma1.n_sk
//run()