pIKdrs = numpy.array(IKdrs)
pIKCa = numpy.array(IKCa)
pIm = numpy.array(Im)
pIl = numpy.array(Il)
pINa = numpy.array(INa)
pICaT = numpy.array(ICaT)
pICaL = numpy.array(ICaL)
pIh = numpy.array(Ih)
r0 = [
pIKa[TStart:TEnd], pIKdrf[TStart:TEnd], pIKdrs[TStart:TEnd],
pIKCa[TStart:TEnd], pIm[TStart:TEnd], pIl[TStart:TEnd], pINa[TStart:TEnd],
pICaT[TStart:TEnd], pICaL[TStart:TEnd], pIh[TStart:TEnd]
]
fig0 = plotCurrentscape(vvec0[TStart:TEnd], r0)
# fig0.savefig('PLOTfiles/' + cellname + '_Currentscape_' + modelname + '_StartTime_' + str(spike0+starttime) + currsize + '.pdf',dpi=500)
fig0.savefig('PLOTfiles/' + cellname + '_Currentscape_' + modelname +
'_StartTime_' + currsize + '_' + str(0) + '.png',
dpi=500)
fig0.clf()
plt.close(fig0)
startcount = 0
for starttime in starttimes:
h.f(Nsyns[0], freqs[0], spike0 + starttime) # Sets synaptic inputs
TStart = int(spike0 * 10) - int(ISIBase * 15)
TEnd = int(spike0 * 10) + int(ISIBase * 25)
vvec0 = numpy.array(v_vec_d)
pIKdrs = numpy.array(IKdrs)
pIKCa = numpy.array(IKCa)
pIm = numpy.array(Im)
pIl = numpy.array(Il)
pINa = numpy.array(INa)
pICaT = numpy.array(ICaT)
pICaL = numpy.array(ICaL)
pIh = numpy.array(Ih)
r0 = [
pIKa[TStart:], pIKdrf[TStart:], pIKdrs[TStart:],
pIKCa[TStart:], pIm[TStart:], pIl[TStart:], pINa[TStart:],
pICaT[TStart:], pICaL[TStart:], pIh[TStart:]
]
fig0 = plotCurrentscape(vvec0[TStart:], r0)
fig0.savefig('PLOTfiles/' + cellname + '_Currentscape_' +
modelname + '_Frequency_' + str(rate) +
'Hz_SynapseNum_' + str(synnum) + '.pdf',
dpi=500)
fig0.savefig('PLOTfiles/' + cellname + '_Currentscape_' +
modelname + '_Frequency_' + str(rate) +
'Hz_SynapseNum_' + str(synnum) + '.png',
dpi=500)
fig0.clf()
plt.close(fig0)
vvec = numpy.array(v_vec)
tvec = numpy.array(t_vec)
numpy.save(
'Output/' + cellname + modelname + '_NumberOfSyns_' + str(synnum) +
tvec0 = numpy.array(t_vec)
pIKa = numpy.array(IKa)
pIKdrf = numpy.array(IKdrf)
pIKdrs = numpy.array(IKdrs)
pIKCa = numpy.array(IKCa)
pIm = numpy.array(Im)
pIl = numpy.array(Il)
pINa = numpy.array(INa)
pICaT = numpy.array(ICaT)
pICaL = numpy.array(ICaL)
pIh = numpy.array(Ih)
r0 = [pIKa,pIKdrf,pIKdrs,pIKCa,pIm,pIl,pINa,pICaT,pICaL,pIh]
fig0 = plotCurrentscape(vvec0, r0)
fig0.savefig('PLOTfiles/' + cellname + '_Currentscape0_' + modelname + '.pdf',dpi=500)
fig0.savefig('PLOTfiles/' + cellname + '_Currentscape0_' + modelname + '.png',dpi=500)
fig0.clf()
plt.close(fig0)
h.ic_step.amp = -0.09 + (TTX_Hold-DNQX_Hold)
h.run()
vvec01 = numpy.array(v_vec)
vvec010 = numpy.array(v_vec0)
tvec01 = numpy.array(t_vec)
pIKa = numpy.array(IKa)
pIKdrf = numpy.array(IKdrf)
pIKdrs = numpy.array(IKdrs)
pIKCa = numpy.array(IKCa)