xopen("./_reset.hoc")
nmda_mg(1)
sprint(input_file,"./stimxout_jns_sqwave_noinput.dat")
init()

stim4.del=200
stim4a.del=210
stim4b.del=220

saveplots("./data3/data_post3_02.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_post3_02_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_post3_02_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_post3_02_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_post3_02_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_post3_02_caTblock.txt")
/*
xopen("./_reset.hoc")
set_kas(0)
set_kasd(0)
saveplots("./data3/data_post3_02_K12block.txt")
xopen("./_reset.hoc")
set_kaf(0)
set_kafd(0)
saveplots("./data3/data_post3_02_K42block.txt")
xopen("./_reset.hoc")
set_naf(0)
set_nap(0)
saveplots("./data3/data_post3_02_Naproxblock.txt")
xopen("./_reset.hoc")
set_nafd(0)
set_napd(0)
saveplots("./data3/data_post3_02_Nadistblock.txt")
*/
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_post3_02_mgfree.txt")




xopen("./_reset.hoc")
nmda_mg(1)
stim4.del=200
stim4a.del=10000
stim4b.del=10000

saveplots("./data3/data_post1_02.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_post1_02_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_post1_02_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_post1_02_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_post1_02_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_post1_02_caTblock.txt")
xopen("./_reset.hoc")
set_kas(0)
set_kasd(0)
saveplots("./data3/data_post1_02_K12block.txt")
xopen("./_reset.hoc")
set_kaf(0)
set_kafd(0)
saveplots("./data3/data_post1_02_K42block.txt")
xopen("./_reset.hoc")
set_naf(0)
set_nap(0)
saveplots("./data3/data_post1_02_Naproxblock.txt")
xopen("./_reset.hoc")
set_nafd(0)
set_napd(0)
saveplots("./data3/data_post1_02_Nadistblock.txt")
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_post1_02_mgfree.txt")


////add current

stim1.del = 0
nmda_mg(1)
xopen("./_reset.hoc")
stim4.del=10000
stim4a.del=10000
stim4b.del=10000
stim4u.del=200
stim4ua.del=210
stim4ub.del=220

saveplots("./data3/data_c_post3_02.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_c_post3_02_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_c_post3_02_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_c_post3_02_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_c_post3_02_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_c_post3_02_caTblock.txt")
/*
xopen("./_reset.hoc")
set_kas(0)
set_kasd(0)
saveplots("./data3/data_c_post3_02_K12block.txt")
xopen("./_reset.hoc")
set_kaf(0)
set_kafd(0)
saveplots("./data3/data_c_post3_02_K42block.txt")
xopen("./_reset.hoc")
set_naf(0)
set_nap(0)
saveplots("./data3/data_c_post3_02_Naproxblock.txt")
xopen("./_reset.hoc")
set_nafd(0)
set_napd(0)
saveplots("./data3/data_c_post3_02_Nadistblock.txt")
*/
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_c_post3_02_mgfree.txt")


xopen("./_reset.hoc")
nmda_mg(1)
stim4u.del=200
stim4ua.del=10000
stim4ub.del=10000

saveplots("./data3/data_c_post1_02.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_c_post1_02_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_c_post1_02_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_c_post1_02_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_c_post1_02_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_c_post1_02_caTblock.txt")
/*
xopen("./_reset.hoc")
set_kas(0)
set_kasd(0)
saveplots("./data3/data_c_post1_02_K12block.txt")
xopen("./_reset.hoc")
set_kaf(0)
set_kafd(0)
saveplots("./data3/data_c_post1_02_K42block.txt")
xopen("./_reset.hoc")
set_naf(0)
set_nap(0)
saveplots("./data3/data_c_post1_02_Naproxblock.txt")
xopen("./_reset.hoc")
set_nafd(0)
set_napd(0)
saveplots("./data3/data_c_post1_02_Nadistblock.txt")
*/
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_c_post1_02_mgfree.txt")