/*----------------------------------------------------------------
Geometry file for thalamic reticular neuron
- 3 compartments
- 1 soma, 1 proximal dendrite, 1 distal dendrite
Obtained from merging the cell from rat nRt, stained by Huguenard & Prince,
Stanford. The cell was traced in the Salk Institute by A. Destexhe.
Starting from a 800 compartment cell, the geometry was reduced
into an equivalent geometry of 10 compartments, by the method
of Bush & Sejnowski (J Neurosci Methods 46: 159-166, 1993).
In a second step, the 10-compartment geometry was merged into 3 compartments
using the same method.
Structure of the reduced compartmental model:
1--0 1--0 1--0
(dend1) (soma)
Length and diameters:
---------------------
Soma : L=34.546 diam=14.075 S=1527.55
dend1[0] : L=103.24 diam=5.56 S=1803.32
dend1[1] : L=190.69 diam=3.06 S=3636.475
=> total dendritic surface = 5439.794 um2
The reconstructed cell was 13302.9 um2
=> reduction factor for dendrites is 2.445
Input resistance:
-----------------
dendritic correction Ri (Mohms)
1 400
2.445 214
3 180
3.5 158
4 143
For more details, see:
Destexhe, A., Contreras, D., Steriade, M., Sejnowski, T.J. and Huguenard,
J.R. In vivo, in vitro and computational analysis of dendritic calcium
currents in thalamic reticular neurons. J. Neurosci. 16: 169-185, 1996.
See also:
http://www.cnl.salk.edu/~alain
http://cns.fmed.ulaval.ca
----------------------------------------------------------------*/
print "----------------------------------------------"
print "Simplified reticular cell with 3 compartments"
print "----------------------------------------------"
/* create compartments (sections) */
create soma[1], \
dend1[2]
/* connect compartments */
soma connect dend1[0](0),1
dend1[0] connect dend1[1](0),1
/* setup 3d geometry : */
soma { /* define soma first */
nseg = 1
pt3dclear()
pt3dadd(0, 0, 0, 14.0754) /* x,y,z,diam */
pt3dadd(0, 34.546, 0, 14.0754)
}
dend1[0] { /* beginning dendrites */
nseg = 1
pt3dclear()
pt3dadd(0, 34.546, 0, 5.56)
pt3dadd(0, 137.786, 0, 5.56)
}
dend1[1] { /* beginning dendrites */
nseg = 1
pt3dclear()
pt3dadd(0, 137.786, 0, 3.06)
pt3dadd(0, 328.476, 0, 3.06)
}