EA.tem

begintemplate sEA		// create a new template object
public soma, kl, ampapost, gabaapost, gababpost, PYlist, TClist, REgabaalist, REgabablist

objectvar ampapost, gabaapost, gababpost, PYlist, TClist, REgabaalist, REgabablist

create soma // soma[1]			// one compartment of about 29000 um2
soma {
  nseg = 1
  diam = 96
  L = 96
  cm = 1
}

objectvar kl

proc init() { local v_potassium, v_sodium

objectvar kl
kl = new kleak()


  v_potassium = -100		// potassium reversal potential 
  v_sodium = 50			// sodium reversal potential 

  soma {
	diam = 96		// geometry 
	L = 96			// so that area is about 29000 um2
	nseg = 1
	Ra = 100

	insert pas		// leak current 
	e_pas = -70		// original from Rinzel
	//e_pas = -100
	g_pas = 1e-5		//1e-5

	kl.loc(0.5)		// K-leak
	Erev_kleak = v_potassium
	kl.gmax = 0.004		// (uS)
				// conversion: x(uS) = x(mS/cm2)*29000e-8*1e3
				//		     = x(mS/cm2) * 0.29



     insert inak2005  //Nav1.1
	ekf = v_potassium
	enat = v_sodium
		//dont change these! - did parameter optimization, bigger Na gives spindles that don't stop  
	gnatbar_inak2005 = 0.15		//0.15
	gkfbar_inak2005  = 0.06

	insert it		// T-current 
	cai = 2.4e-4 
	cao = 2 
	eca = 120 
	gcabar_it = 0.002 
	shift_it = 2

	insert iar		// h-current
	//eh = -40		// reversal
	eh = -40
	nca_iar = 4		// nb of binding sites for Ca++ on protein
	k2_iar = 0.0004		// decay of Ca++ binding on protein
	cac_iar = 0.002		// half-activation of Ca++ binding
	nexp_iar = 1		// nb of binding sites on Ih channel
	k4_iar = 0.001		// decay of protein binding on Ih channel
	Pc_iar = 0.01		// half-activation of binding on Ih channel
	ginc_iar = 2		// augm of conductance of bound Ih
	//ginc_iar = 1.25 		// augm of conductance of bound Ih
	ghbar_iar = 2e-5	// low Ih for slow oscillations


	insert cad		// calcium decay
	depth_cad = 1
	taur_cad = 5
	cainf_cad = 2.4e-4
	kt_cad = 0		// no pump

	PYlist = new List()
	TClist = new List()
	REgabaalist = new List()
	REgabablist = new List()

	gababpost = new List()

  }
}
endtemplate sEA