* simulation from:
Lytton WW; Contreras D; Destexhe A; Steriade M. Dynamic interactions
determine partial thalamic quiescence in a computer network model of
spike-and-wave seizures. J Neurophysiol, 77:1679-1696, 1997.
Fig 7 on p 1688  (jnphys77:1679.pdf)

Simulation of figure 7 may take several minutes. fig7.gif shows the output of
the simulation.

A phase plane analysis of a two cell interaction between a
thalamocortical neuron (TC) and a thalamic reticularis neuron (RE).
Figure legend : Analysis of spontaneous transition from mutual to
quiescent mode with bursting viewed in the RE T-channel m/h plane.
  A. RE neuron voltage trace (A1) and TC voltage trace (A2) shows 8
bursts of mutual activity followed by switch into the quiescent
state.  Color for each cycle give the key for the phase planes
below.  Interburst intervals are depicted in panel Ba.  Alternating
bursts are shown in panels Bb and Bc.  Initial oscillation appears
to be approximately period-2 with burst alternating with an
interburst depolarization.
  B. Apparent approximate period-2 oscillation illustrated to be
period-4: four cycles occur before the cycle repeats.  In each phase
plane, abscissa gives the value of mT (T channel activation),
ordinate the value of hT (T channel inactivation).  
Trajectories run counterclockwise.  Interbursts alternated
between large cycles (solid line type) and small cycles (dashed line type)
  Burst Bb
showed a progression of gradually increasing hT, while burst Bc
was more coherent and showed no regular shift.  The inner loop in Ba
gradually collapsed to terminate the oscillation.  The subsequent
interburst (grey dashed) was similar in extent to large cycle.
  C. Inner and outer interburst cycles seen in the voltage-calcium plane
corresponded to Ba-inner and Ba-outer respectively.  The collapse of the inner
cycle is not as clear as in Ba, however.  Trajectories run clockwise.

* COMPILE and RUN:
  nrnivmodl
RUN: 
  $CPU/special init.hoc -
RE-RUN:
  simsim() allows you to simulate a simulation by graphing out the
    saved data slowly so that you can see the phase planes unfold in time
    (can be slowed down by setting disp_delay to a large number)

* notes
1) graphics box comes up with graphs set to "New View" to allow closer
inspection of phase planes
2) use 'vecpanel()' call to bring up interactive tool for generating
other graphics
3) original figure only showed data from about t=600 rather than from beginning

* hoc list: from ~/nrniv/tcre2 geom.hoc v5, network.hoc v39, params.hoc v108
 run.hoc v105

init.hoc    master hoc file
geom.hoc    defines single compartment cells and column
network.hoc creates cells, connectivity and stimulation
params.hoc  sets parameters
run.hoc     additional initialization routines, graphing

boxes.hoc   make large arrays of graphs. not used here due to
specialized layout
declist.hoc list handling routines
decvec.hoc  vector handling routines
default.hoc flags and values
grvec.hoc   main vector graphic routines
local.hoc   additional initialization routines
nrnoc.hoc   main initialization routines
simctrl.hoc simulation control package using RCS

* mod list
HH_traub.mod            SUFFIX hh2
IAHP_destexhe.mod:      SUFFIX iahp
ICAN_destexhe.mod:      SUFFIX ican
IT2_huguenard.mod:      SUFFIX it2
IT_wang.mod             SUFFIX it
Ih_old.mod:             SUFFIX iar
calciumpump_destexhe.mod SUFFIX cad
gen.mod:                POINT_PROCESS gen
nmda.mod:               POINT_PROCESS NMDA
passiv.mod:             SUFFIX Pass
presyn.mod:             POINT_PROCESS PRESYN
pulse.mod:              POINT_PROCESS PULSE
AMPA.mod                POINT_PROCESS AMPA
NMDA.mod                POINT_PROCESS NMDA
GABAB1.mod              POINT_PROCESS GABAB1
GABALOW.mo              POINT_PROCESS GABALOW
rand.mod:               SUFFIX nothing

* mechanisms and processes used
Pass     leak current 
hh2      Hodgin-Huxley INa and IK 
it       TC IT current 
iar      Ih current 
it2      reticular IT current 
iahp     slow IAHP 
ican     slow ICAN 
cad      calcium pump 

gen      random spike generator
AMPA
NMDA
IClamp

Changelog
---------
2022-05: Updated MOD files to contain valid C++ and be compatible with the
         upcoming versions 8.2 and 9.0 of NEURON.