/* load_file ("stdgui.hoc") */
// load_proc ("nrnmainmenu")
chdir ("e:\AmaMartB")
load_file ("martB.nrn")

forall nseg=int(L/3)+1
forall insert pas
forall g_pas=1/70000
forall Ra=110
global_ra=110
forall e_pas=-70
forall insert spike
celsius=22
forall insert cad
forall ek = -75
forall ena = 35
/* Make dendrites 2X soma  */
forall gcabar_spike = 0.0015
forall gkbar_spike = 0.018
forall gabar_spike = 0.054
forall gkcbar_spike = 0.0015
forall gnabar_spike = 0.1 // Same as soma 100% i.e. normal is 0.04
forall taur_cad = 10

/* Make Soma= normal for trigger function */
access soma
gcabar_spike = 0.0015   // normal is 0.0022
gkbar_spike = 0.018   // normal 0.012
gabar_spike = 0.054 // normal  is .036 so 2.22 times
gkcbar_spike = 0.0015  // normal is 0.00005 so 2 times to slow soma down
gnabar_spike = 0.1      // normal is 0.04
taur_cad = 10
g_pas=1/70000           // soma leak present due to electrode?

/* put in stims for on and off */

access dend1[52]
objectvar stimON
stimON=new AlphaSynapse (0.5)
stimON.onset = 250      // ms
stimON.tau = 50         // ms
stimON.gmax = 0.0005    // umho
objectvar stimOFF
stimOFF=new AlphaSynapse (0.5)
stimOFF.onset = 750
stimOFF.tau = 50
stimOFF.gmax = 0.0004

access dend1[63]
objectvar stim2ON
stim2ON=new AlphaSynapse (0.5)
stim2ON.onset = 250     // ms
stim2ON.tau = 50        // ms
stim2ON.gmax = 0.0002   // umho
objectvar stim2OFF
stim2OFF=new AlphaSynapse (0.5)
stim2OFF.onset = 750
stim2OFF.tau = 50
stim2OFF.gmax = 0.0002

access dend1[60]
objectvar stim3ON
stim3ON=new AlphaSynapse (0.5)
stim3ON.onset = 250     // ms
stim3ON.tau = 50        // ms
stim3ON.gmax = 0.0002   // umho
objectvar stim3OFF
stim3OFF=new AlphaSynapse (0.5)
stim3OFF.onset = 750
stim3OFF.tau = 50
stim3OFF.gmax = 0.0002

access dend1[87]
objectvar stim4ON
stim4ON=new AlphaSynapse (0.1)
stim4ON.onset = 250     // ms
stim4ON.tau = 50        // ms
stim4ON.gmax = 0.0002   // umho
objectvar stim4OFF
stim4OFF=new AlphaSynapse (0.1)
stim4OFF.onset = 750
stim4OFF.tau = 50
stim4OFF.gmax = 0.0002

access dend2[13]
objectvar stim5ON
stim5ON=new AlphaSynapse (0.4)
stim5ON.onset = 250     // ms
stim5ON.tau = 50        // ms
stim5ON.gmax = 0.0005   // umho
objectvar stim5OFF
stim5OFF=new AlphaSynapse (0.4)
stim5OFF.onset = 750
stim5OFF.tau = 50
stim5OFF.gmax = 0.0004

access dend1[80]
objectvar stim6ON
stim6ON=new AlphaSynapse (0.5)
stim6ON.onset = 250     // ms
stim6ON.tau = 50        // ms
stim6ON.gmax = 0.0002   // umho
objectvar stim6OFF
stim6OFF=new AlphaSynapse (0.5)
stim6OFF.onset = 750
stim6OFF.tau = 50
stim6OFF.gmax = 0.0002

access dend1[76]
objectvar stim7ON
stim7ON=new AlphaSynapse (0.3)
stim7ON.onset = 250     // ms
stim7ON.tau = 50        // ms
stim7ON.gmax = 0.0002   // umho
objectvar stim7OFF
stim7OFF=new AlphaSynapse (0.3)
stim7OFF.onset = 750
stim7OFF.tau = 50
stim7OFF.gmax = 0.0002

access dend1[73]
objectvar stim8ON
stim8ON=new AlphaSynapse (0.5)
stim8ON.onset = 250    // ms
stim8ON.tau = 50        // ms
stim8ON.gmax = 0.0002   // umho
objectvar stim8OFF
stim8OFF=new AlphaSynapse (0.5)
stim8OFF.onset = 750
stim8OFF.tau = 50
stim8OFF.gmax = 0.0002

access dend1[40]
objectvar stim9ON
stim9ON=new AlphaSynapse (0.2)
stim9ON.onset = 250     // ms
stim9ON.tau = 50        // ms
stim9ON.gmax = 0.0002   // umho
objectvar stim9OFF
stim9OFF=new AlphaSynapse (0.2)
stim9OFF.onset = 750
stim9OFF.tau = 50
stim9OFF.gmax = 0.0002

access dend1[8]
objectvar stim10ON
stim10ON=new AlphaSynapse (0.2)
stim10ON.onset = 250    // ms
stim10ON.tau = 50       // ms
stim10ON.gmax = 0.0002  // umho
objectvar stim10OFF
stim10OFF=new AlphaSynapse (0.2)
stim10OFF.onset = 750
stim10OFF.tau = 50
stim10OFF.gmax = 0.0002

access dend1[11]
objectvar stim11ON
stim11ON=new AlphaSynapse (0.8)
stim11ON.onset = 250    // ms
stim11ON.tau = 50       // ms
stim11ON.gmax = 0.001   // umho
objectvar stim11OFF
stim11OFF=new AlphaSynapse (0.8)
stim11OFF.onset = 750
stim11OFF.tau = 50
stim11OFF.gmax = 0.001

access dend1[46]
objectvar stim12ON
stim12ON=new AlphaSynapse (0.2)
stim12ON.onset = 250    // ms
stim12ON.tau = 50       // ms
stim12ON.gmax = 0.001   // umho
objectvar stim12OFF
stim12OFF=new AlphaSynapse (0.2)
stim12OFF.onset = 750
stim12OFF.tau = 50
stim12OFF.gmax = 0.001

access dend1[38]
objectvar stim13ON
stim13ON=new AlphaSynapse (0.2)
stim13ON.onset = 250    // ms
stim13ON.tau = 50       // ms
stim13ON.gmax = 0.0002  // umho
objectvar stim13OFF
stim13OFF=new AlphaSynapse (0.2)
stim13OFF.onset = 750
stim13OFF.tau = 50
stim13OFF.gmax = 0.0002


dt=0.025
tstop=2000
steps_per_ms =5
v_init=-70
init()
forall depth_cad=diam/2

// xopen ("martB.ses")