stim1.del = 100000
stim2.amp=0
xopen("./_reset.hoc")
nmda_mg(1)
sprint(input_file,"./stim_files2/stimxout_s0_prep000_DANO.txt")
init()

saveplots("./data3/data_pre_spine0.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_pre_spine0_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_pre_spine0_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_pre_spine0_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_pre_spine0_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_pre_spine0_caTblock.txt")
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_pre_spine0_mgfree.txt")


xopen("./_reset.hoc")
nmda_mg(1)
sprint(input_file,"./stim_files2/stimxout_s1_prep000_DANO.txt")
init()

saveplots("./data3/data_pre_spine1.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_pre_spine1_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_pre_spine1_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_pre_spine1_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_pre_spine1_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_pre_spine1_caTblock.txt")
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_pre_spine1_mgfree.txt")



////add current

stim1.del = 0
xopen("./_reset.hoc")
nmda_mg(1)
sprint(input_file,"./stim_files2/stimxout_s0_prep000_DANO.txt")
init()

saveplots("./data3/data_c_pre_spine0.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_c_pre_spine0_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_c_pre_spine0_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_c_pre_spine0_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_c_pre_spine0_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_c_pre_spine0_caTblock.txt")
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_c_pre_spine0_mgfree.txt")


xopen("./_reset.hoc")
nmda_mg(1)
sprint(input_file,"./stim_files2/stimxout_s1_prep000_DANO.txt")
init()

saveplots("./data3/data_c_pre_spine1.txt")
xopen("./_reset.hoc")
chnmdawt(0)
saveplots("./data3/data_c_pre_spine1_nmdablock.txt")
xopen("./_reset.hoc")
chmgluwt(0)
saveplots("./data3/data_c_pre_spine1_mglublock.txt")
xopen("./_reset.hoc")
champawt(0)
saveplots("./data3/data_c_pre_spine1_ampablock.txt")
xopen("./_reset.hoc")
set_caL(0)
set_caL13(0)
saveplots("./data3/data_c_pre_spine1_caLblock.txt")
xopen("./_reset.hoc")
set_cat(0)
saveplots("./data3/data_c_pre_spine1_caTblock.txt")
xopen("./_reset.hoc")
nmda_mg(0)
saveplots("./data3/data_c_pre_spine1_mgfree.txt")