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

import VrampExp 
import imp
imp.reload(VrampExp)


Vinit=[-160]
paraKeys=['Clamp_parameters/Cai','Clamp_parameters/Nai','Clamp_parameters/Ki',]
paras=[8.97e-5,4,140]
outKeys=['outputs/b_inf','outputs/g_inf','outputs/tau_b','outputs/tau_g','outputs/d_inf',
         'outputs/f_inf','outputs/tau_d','outputs/tau_f2','outputs/m_inf','outputs/h_inf',
         'outputs/tau_m','outputs/tau_h','outputs/q_inf','outputs/r_inf','outputs/tau_q',
         'outputs/tau_r1','outputs/tau_r2','outputs/p_inf','outputs/k_inf','outputs/tau_p',
         'outputs/tau_k1','outputs/tau_k2','outputs/s_inf','outputs/x_inf','outputs/tau_s',
         'outputs/tau_x','outputs/y_inf','outputs/tau_y','outputs/n_inf','outputs/c_inf',
         'outputs/tau_c1','outputs/tau_c2',]

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

for i, sedversion in enumerate(sedversions):
    simfile='../Voltage_ramp_experiment_%s.sedml' % sedversion
    simFign='Vramp_%s' % sedversion
    VrampExp.VrampExp(simFign,simfile,Vinit,paraKeys,paras,outKeys)