Readme for the model associated with the paper

Traub RD, Contreras D, Cunningham MO, Murray H, Lebeau FE, Roopun A, 
Bibbig A, Wilent WB, Higley M, Whittington MA.
A single-column thalamocortical network model exhibiting gamma 
oscillations, sleep spindles and epileptogenic bursts.
J Neurophysiol. 2004 Nov 3; [Epub ahead of print]

This is a conversion of the original model code to the g77 compiler.
To run the model first build it by typing

make

on a cluster with mpi installed.  Then you can start up the simulation
with:

mpirun -machinefile hostnames.txt -np 14 ./groucho 

where hostnames.txt holds the names of the nodes that you would like to use.
We use the supplied file job1 and just type
./job1&

Here are some notes that were supplied with the original model that
will not work with this g77 version:

1st I compile all the integration subroutines, with compile_integration,
which uses ifc.

Then type "make", which compiles the main program (with mpifwhatever) 
and links everything to give an executable. Then "qsub job", which 
ships the lot to the master node of the Linux cluster; that in turn 
loads the stuff onto 14 cpu's.

groucho.f = main program
compile_integration = does what it says, compiles integration subroutines
makefile creates an executable for "job"
job = what gets submitted, with command "qsub job"
