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)