# To reproduce the data needed for clamping experiments in associated original paper,

import IVFigPlot
import imp
imp.reload(IVFigPlot)

sedversions=['0','1','2','3']

simFign='FigA1F'
V_clamp=[-60, 120, 120, -50, 80, 10]
paraKeys=['Clamp_parameters/Cai']
paras=[9.43e-5]
outKeys=['outputs/I_CaL']
current= r'$I_{CaL}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)

simFign='FigA2F'
V_clamp=[-80, 100, 100, -70, 80, 10]
paraKeys=['Clamp_parameters/Cai']
paras=[9.43e-5]
outKeys=['outputs/I_CaT']
current= r'$I_{CaT}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)

simFign='FigA3F'
V_clamp=[-90, 50, 50, -50, 80, 10]
paraKeys=['Clamp_parameters/Nai']
paras=[4]
outKeys=['outputs/I_Na']
current= r'$I_{Na}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)

simFign='FigA4F'
V_clamp=[-80, 8000, 8000, -60, 80, 10]
paraKeys=['Clamp_parameters/Ki']
paras=[140]
outKeys=['outputs/I_K1']
current= r'$I_{K1}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)

simFign='FigA5F'
V_clamp=[-80, 8000, 8000, -60, 80, 10]
paraKeys=['Clamp_parameters/Ki']
paras=[140]
outKeys=['outputs/I_K2']
current= r'$I_{K2}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)


simFign='FigA6F'
V_clamp=[-80, 150, 150, -80, 80, 10]
paraKeys=['Clamp_parameters/Ki']
paras=[140]
outKeys=['outputs/I_KA']
current= r'$I_{KA}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)

simFign='FigA7D'
V_clamp=[-10, 3000, 3000, -20, -130, -10]
paraKeys=['Clamp_parameters/Nai','Clamp_parameters/Ki']
paras=[4,140]
outKeys=['outputs/I_h']
current= r'$I_{h}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)

simFign='FigA8E'
V_clamp=[0, 100, 100, -60, 50, 10]
paraKeys=['Clamp_parameters/Cai']
paras=[9.43e-5]
outKeys=['outputs/I_Cl']
current= r'$I_{Cl}$'
IVFigPlot.IVFigPlot(simFign,sedversions,current,outKeys)