//genesis
float dt = 20e-6
str disk = "results/"
int amp = {1 / dt}
str a = (amp)
int i, k, l
str el, name
float t
/*********************************************************************
** Simple Granule cell model script (#1)
** Carl Piaf BBF 1994
*********************************************************************/
str filename = (disk) @ "G_soma_GABA_alpha3ms"
/* always include these default definitions! */
include defaults
str cellpath = "/Granule"
/* Purkinje cell constants */
include Gran_const.g
/* special scripts to create the prototypes */
include Gran_chan_tab.g
//include Gran_synchan2_biexpGABA.g
include Gran_synchan.g
include Gran_comp.g
/* Set the clocks */
for (i = 0; i <= 7; i = i + 1)
setclock {i} {dt}
end
setclock 8 2.0e-5
setclock 9 1
/* To ensure that all subsequent elements are made in the library */
if (! {exists /library/granule})
create neutral /library/granule
end
ce /library/granule
/* These make the prototypes of channels and compartments that can be
** invoked in .p files */
make_Granule_chans
make_Granule_syns
make_Granule_comps
/*
call Gran_InNa TABSAVE tabInNa37.data
call Gran_KDr TABSAVE tabKDr37.data
call Gran_KA TABSAVE tabKA37.data
call Gran_CaHVA TABSAVE tabCaHVA37.data
call Gran_H TABSAVE tabH37.data
call Moczyd_KC TABSAVE tabKCa37.data
*/
//make_Vmgraph
/* create the model and set up the run cell mode */
// read cell data from .p file
readcell Gran1M0.p {cellpath} // Gran1M0.p {cellpath}
create neutral /library/granule/soma/mf_presyn
disable /library/granule/soma/mf_presyn
setfield /library/granule/soma/mf_presyn z 0
// Comment out whichever one to switch it off
addmsg /library/granule/soma/mf_presyn /Granule/soma/GABAA ACTIVATION z
//addmsg /library/granule/soma/mf_presyn /Granule/soma/GABAB ACTIVATION z
//addmsg /library/granule/soma/mf_presyn /Granule/soma/mf_NMDA ACTIVATION z
//addmsg /library/granule/soma/mf_presyn /Granule/soma/mf_AMPA ACTIVATION z
setfield /Granule/soma/GABAA gmax {{getfield /Granule/soma/GABAA gmax} * 10} // = * 45 / 1.5 / 3
setfield /Granule/soma/GABAA tau1 {{getfield /Granule/soma/GABAA tau2} / 3}
setfield /Granule/soma/GABAA tau2 {{getfield /Granule/soma/GABAA tau2} / 3}
/*
create randomspike /random
setfield ^ rate 100 abs_refract 0.01
create spikegen /random/spike
setfield ^ thresh 0.5
addmsg /random /random/spike INPUT state
addmsg /random/spike {cellpath}/soma/mf_NMDA SPIKE
setfield {cellpath}/soma/mf_NMDA synapse[0].weight {dt} \
synapse[0].delay 0.0
*/
//read_hines -vm -method 11 Gran1M0.p {cellpath}
/* Create the output element */
create asc_file /output/plot_out
//create disk_out /output/plot_out
useclock /output/plot_out 8
enable /output
enable /output/plot_out
ce {cellpath}
// setup the hines solver
echo preparing hines solver...
create hsolve solve
ce solve
// if this is set then reset will NOT change Vm in Hines
setfield . path "../##[][TYPE=compartment]" comptmode 1 chanmode 4
/* Vm_reset 1 !! geeft foutboodschap cannot find field 'Vm_reset'
on '/Granule/solve' */
call . SETUP
setmethod 11
// setmethod 0
// addmsg {cellpath}/soma /output/plot_out SAVE Vm
/*
el = ({findsolvefield {cellpath}/solve soma/H Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#2
el = ({findsolvefield {cellpath}/solve soma/Ca_pool Ca})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#3
*/
el = ({findsolvefield {cellpath}/solve soma Vm}) //#4
addmsg {cellpath}/solve /output/plot_out SAVE {el}
/*
el = ({findsolvefield {cellpath}/solve soma/Moczyd_KC Gk})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#5
el = ({findsolvefield {cellpath}/solve soma/KDr Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#6
el = ({findsolvefield {cellpath}/solve soma/InNa Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#7
el = ({findsolvefield {cellpath}/solve soma/KA Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#8
el = ({findsolvefield {cellpath}/solve soma/CaHVA Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#9
el = ({findsolvefield {cellpath}/solve soma/mf_AMPA Gk})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#10
el = ({findsolvefield {cellpath}/solve soma/mf_AMPA Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#10
el = ({findsolvefield {cellpath}/solve soma/mf_NMDA Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#0
el = ({findsolvefield {cellpath}/solve soma/Mg_BLOCK Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#0
el = ({findsolvefield {cellpath}/solve soma/Mg_BLOCK Gk})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //#0
*/
el = ({findsolvefield {cellpath}/solve soma/GABAA Gk})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //
el = ({findsolvefield {cellpath}/solve soma/GABAA Ik})
addmsg {cellpath}/solve /output/plot_out SAVE {el} //
setfield /output/plot_out filename {filename} initialize 1 leave_open 1 \
flush 1
echo Output to {filename}
//check
reset
// Synaptic stimulation protocol
setfield {cellpath}/soma inject 20.0e-12
call {cellpath}/solve HPUT {cellpath}/soma
step 0.3 -time
/*
setfield {cellpath}/soma Em -0.065
setfield {cellpath}/soma Vm -0.0625
setfield {cellpath}/soma initVm -0.0625
call {cellpath}/solve HPUT {cellpath}/soma
*/
// step 0.05 -time
setfield /library/granule/soma/mf_presyn z {amp}
step 1
setfield /library/granule/soma/mf_presyn z 0
step 100e-3 -time
step 0.3 -t
// Current injection protocol
/*
step 0.1 -time
setfield {cellpath}/soma inject 20.0e-12
call {cellpath}/solve HPUT {cellpath}/soma
step 0.5 -time
call {cellpath}/solve HSAVE {cellpath}/soma
setfield {cellpath}/soma inject 0.0e-12
call {cellpath}/solve HPUT {cellpath}/soma
step 0.2 -time
setfield {cellpath}/soma inject 0.0e-12
call {cellpath}/solve HPUT {cellpath}/soma
step 1.0 -time
*/