argdict = {'maxi':60,'mini':20,'freq':4} # this is probably overkill
CONTOG=0#1.000000e+00
ALL = 0 # 1 for Fig 4C6

p_init = dict({})

p_init['bothamp'] = 0

p_init['nprint'] = 11
p_init['taperkv4'] = 1

p_init['ampa_ratio']=0

p_init['cicrtog'] = 0

p_init['tausk'] = 15

p_init['von'] = 0
p_init['voff']=-80

p_init['constrict']=0.5
p_init['dsk'] = 0.05
p_init['dcicr'] = 0.01
p_init['dense'] = 3

p_init['DISABLE']=1

p_init['distscale']=1
p_init['knmda'] = 4

p_init['girk'] = 0

p_init['kmatp']=1

p_init['glutscale']=1000

p_init['sigamp'] = 100e-3 # peak current
p_init['sigslope'] = 1000 #ms

p_init['gnmda'] =0e-6 #nmda
p_init['cmg'] = 1.2 #1.2 for NMDA, 0 for 'AMPA' or light activated channels
p_init['cafrac']=0

p_init['ecl'] = -60 # Chloride reversal for GABAa

p_init['gabaamp']=0
p_init['gababamp']=0
p_init['chrhodamp']=0
#3640
p_init['regamp']=0#5e-3 # this is producing inputs where it should not
p_init['regtype'] = 'gabaa'
p_init['freq'] = 50
p_init['number']=0

p_init['nmdaamp']=0

p_init['won']=50
p_init['woff']=2950
p_init['fon'] =100
p_init['foff']=100

p_init['musc'] = 0

p_init['deexcite']=0
p_init['inhibit']=0
p_init['noise']=1
p_init['pausedur'] = 500
p_init['ncells']=1
p_init['syntype']='gabaa'

p_init['kass'] = 0 # shift in kv4 slope
p_init['kasv'] = 10 # shift in activation/inactivation Kv4


p_init['both']=1
p_init['kchip']=0.56*(1-CONTOG)+0.05*CONTOG#0.5*(1-CONTOG)# Kv4 with slow kinetics
p_init['fchip']=0.44*(1-CONTOG)+0.95*CONTOG#0.4+0.5*CONTOG # Kv4 with fast kinetics (for viral deletion the total conductance may change)

p_init['fixkv4'] = 0
p_init['qslow'] = 0
p_init['qfast'] =0

p_init['tfast'] = 25

p_init['taukv4'] = 1.5 # scale factor for max inactivation time constant of slow kv4 (1 = 250 ms) depriciated
#

p_init['fast'] = 0
p_init['na_shift'] = 0.0 # mshift positive decreases excitability
p_init['na_hshift'] = 0#-5.0 # hshift
p_init['na_sshift'] = 0.0

p_init['aiscor'] = 0
p_init['allcor'] = 0
p_init['ra'] = 1# 1 = 100 ohm-cm
# making the cell smaller works?

p_init['gkatp'] = 0
p_init['gbk'] = 1

p_init['bks'] = 1

p_init['meta'] = 0.1
p_init['kadp'] = 0.0005
p_init['skatp'] = 1

p_init['pfcal'] = 0.6

p_init['glna'] = 2.75+2.25*CONTOG
p_init['glk'] = 3+3*CONTOG # 1 = 10 uS/cm2
p_init['kmark'] = 0 # use markov potassium (post manuscript code)

p_init['gcan'] =8
p_init['capump'] = 0.3# CA pump
p_init['gka'] = 0.67 # combination of peak conductance of 1, qhalf = -65 (5 kasv) creates peak currents in real morphologies comparable to Jochen's data
p_init['gcal'] = 0.4 # L type
p_init['coup'] = 0

#p_init[16] = -0
p_init['ais_na'] = 200 #sodium conductance in AIS (x1000e-6)
p_init['ais_k'] = 15 # kdr conductance in AIS (x500ex6
#p_init[12] = 0
#p_init[14] = 5 #3
p_init['gcat'] = 1.0*CONTOG # CaV3
#0.5
catbuff = 0.0 # 1 means CaT is not coupled to SK
taukv4 = 1 # 25 ms inactivation time constant
capval =1 #
ashift = -10
SCALE =1





p_init['chop'] = 1
p_init['dst'] = 300  #'distal' distance

p_init['rtype']='nexp' # poisson is actually recursive poisson.  Don't use it


p_init['catfix'] = 0

phi_h = 1.5 # speed scale factor for CaT
p_init['cm'] = 1
p_init['scale'] =1

p_init['proxyadp'] = 0
p_init['ratio'] = 1 # fraction of NaV that is slow inactivating

p_init['gnahh'] =15.0*(1+CONTOG) #poorly named
p_init['nafrac'] = 0.0


  #  0.6 - 0.5*(CONTOG+ALL)#-0.5 # 0.6 0.1
p_init['ghcn'] = 0.5*(CONTOG+ALL)#+0.3 # 0, 0.5
p_init['mbar'] = -1 # fix HCN activation

p_init['gkca'] =(0.5+0.5*(CONTOG+ALL))#0.4+0.6*CONTOG#+0.3 # 0.2, 1
p_init['dca'] = 0.2  #reduced mobility allows for larger calcium swings
p_init['local'] = 0 # 1 or 0 for use of microdomain for SK


p_init['control'] = 0
p_init['all'] = 0


p_init['nspeed'] = 0.2

p_init['gkhh'] =2.5+2.5*CONTOG+ALL*1.5#1.5 4
p_init['nshift']=0#-5*CONTOG #
p_init['nslope'] = 12

p_init['tcicr'] = 0

p_init['size'] =500*(1+1*CONTOG) # conventional size lowered a bit so that failure remains at 100 pA ramp

p_init['slow'] =(1-CONTOG)*0.67+(2*CONTOG)# 1 = 40 1/s 3 (120/s) matches Ding data, 2 (80 /s) matches NMDA data.

p_init['slow_ais'] = 0.1

# stimulus parameter
p_init['nkatp']=3
p_init['idel'] = 2000 # 
p_init['idur'] = 4000 # width of V-ramp
p_init['basal'] = 0#(-25e-3) # initial hyperpolarizing holding step for bidirectional ramp

p_init['iamp'] =0#(100.0e-3)  #translates to 50, 100 pA for 300, 600 um change this and in dep_fun so that the basal is the special case
p_init['offset'] = 200
p_init['i2amp'] =0
p_init['i2dur']=10

p_init['gkm'] = 0.75
p_init['gkerg'] =0 

p_init['catdist'] = 50
p_init['tmhalf'] = -52
p_init['thhalf'] = -81
p_init['tphi'] = 1.5
p_init['catbuff']=0
#p_init[26] = .0