Regulation of a slow STG rhythm (Nadim et al 1998)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:3511
Frequency regulation of a slow rhythm by a fast periodic input. Nadim, F., Manor, Y., Nusbaum, M. P., Marder, E. (1998) J. Neurosci. 18: 5053-5067
Reference:
1 . Nadim F, Manor Y, Nusbaum MP, Marder E (1998) Frequency regulation of a slow rhythm by a fast periodic input. J Neurosci 18:5053-67 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Stomatogastric ganglion;
Cell Type(s): Stomatogastric Ganglion (STG) Modulatory commissural neuron 1 (MCN1); Stomatogastric Ganglion (STG) interneuron 1 (Int1); Stomatogastric Ganglion (STG) Lateral Gastric (LG) cell;
Channel(s): I Na,t; I K;
Gap Junctions: Gap junctions;
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Temporal Pattern Generation; Invertebrate;
Implementer(s): Nadim, Farzan [Farzan at andromeda.Rutgers.edu];
Search NeuronDB for information about:  I Na,t; I K;
load_file("nrngui.hoc")

strdef file,command
double int1values[10]
double lgvalues[10]
double mcn1values[10]
double mcn1nvalues[3]

/***********************************************************/
pyloric=1000
tmax=0
max_time=10000

while(tmax < max_time){
	tmax = tmax + pyloric
}
dt=.05
secondorder=2
plot_it_too=1
xopen("constants.oc")
vertical_displacement=40
one_out_of_this_many=1/dt
coupling_method=FIXED
/*coupling_method=VDEP*/

xopen("mkcells.oc")
xopen("mkstim.oc")
xopen("mksyn.oc")
xopen("modify.oc")
xopen("init.oc")

proc myrun() {
if(plot_it_too){
	axis(0,tmax,5,-80,-20+1*vertical_displacement,10)
	axis()
	graph() 
	graph("INT1s.v","setcolor(RED)") 
	graph("LGs.v+vertical_displacement","setcolor(BLACK)")
	graph("MCN1n[MCN1_SEG-1].v+2*vertical_displacement","setcolor(BLUE)") 
}

t = 0
j = 0
/* sprint(file,"out.%d",pyloric) */
sprint(file,"out")
wopen(file)

initialize()

while (t < tmax) {

	if(j==one_out_of_this_many){  	

/*		fprint("%f %f %f %f %f %f %f %f %f\n",\
					t, MCN1n[MCN1_SEG-1].v,LGs.v,\
					INT1s.v,mcn1_lg.g*100000,MCN1a.v,\
					LGa.v,INT1a.v,LGn[0].v)
*/
		fprint("%.2f %.1f %.1f %.1f %.1f %.2f\n",\
					t,LGs.v,LGa.v,INT1s.v, MCN1a.v,mcn1_lg.g*100000)
		j=0
	}else{
		j=j+1}


	if(plot_it_too) {
		graph(t)
	}
	fadvance()
}
	if(plot_it_too) {
		graphmode(-1)
	}

	print_init()
wopen()
}

//run()

/***********************************************************/

// minimal modifications so that run is managed by stdrun.hoc
tstop = tmax

proc init() {
	max_time=tstop
	tmax = 0
	while(tmax < max_time){
		tmax = tmax + pyloric
	}
	// from mksyn.oc: needed if manually change tstop
	/*  AB -> INT1 */
	fsyn(tmax/pyloric)
	for(i=0;i<tmax/pyloric;i=i+1) {
		INT1n[PASSIVE_SEG-1] fsyn(i,0.5,pyloric*i,80,0.0018,-70)
	}
	initialize()
}

load_file("all.ses")


Loading data, please wait...