This is a part of my program
I try to change dendrites' parameters by the modification of fiber's diameter through a window.
Code: Select all
/création d'une procédure, pour récupérer le diamètre et modifier
load_file("cell.ses")
if (CellBuild[0].continuous==0) {
CellBuild[0].continuous=1
}
//definition parametre du grand axone
celsius=37 //Temperature (°c)
axonnodes=21 //Nombre de noeud (SI)
paranodes1=40 //Nombre de paranode à gauche (SI)
paranodes2=40 //Nombre de paranode à droite (SI)
axoninter=120 //Nombre d'internode (SI)
axontotal=221 //Nombre total (SI)
D_droite=5.7 //Diametre de la fibre (um)
paralength1=3 //Longueur de myelin entre les paranode (um)
nodelength=1.0 //Longueur d'un noeud (um)
space_p1=0.002 //largeur de myelin entre les paranodes (um)
space_p2=0.004 //largeur des paranodes (um)
space_i=0.004 //largeur des internodes (um)
rhoa=0.7e6 //Resistivite (ohm/um)
mycm=0.1 //Capacite du Myelin (uF/cm²)
mygm=0.001 //Conductivite du Myelin (S/cm²)
//case selon la taille de la fibre
proc affectation_droite(){
if (D_droite==5.7){
g=0.605
axonD=3.4 //Diametre des internodes (um)
nodeD=1.9 //Diametre des noeuds (um)
paraD1=1.9 //Diametre du myelin (um)
paraD2=3.4 //Diametre des paranodes (um)
deltax=500 //Longueur separant 2 nodes (um)
paralength2=35 //Longueur des paranodes (um)
nl=80} //Nombre de "myelin lemella" (SI)
if (D_droite==7.3){
g=0.630
axonD=4.6 //Diametre des internodes (um)
nodeD=2.4 //Diametre des noeuds (um)
paraD1=2.4 //Diametre du myelin (um)
paraD2=4.6 //Diametre des paranodes (um)
deltax=750 //Longueur separant 2 nodes (um)
paralength2=38 //Longueur des paranodes (um)
nl=100} //Nombre de "myelin lemella" (SI)
if (D_droite==8.7){
g=0.661
axonD=5.8 //Diametre des internodes (um)
nodeD=2.8 //Diametre des noeuds (um)
paraD1=2.8 //Diametre du myelin (um)
paraD2=5.8 //Diametre des paranodes (um)
deltax=1000 //Longueur separant 2 nodes (um)
paralength2=40 //Longueur des paranodes (um)
nl=110} //Nombre de "myelin lemella" (SI)
if (D_droite==10.0){
g=0.690
axonD=6.9 //Diametre des internodes (um)
nodeD=3.3 //Diametre des noeuds (um)
paraD1=3.3 //Diametre du myelin (um)
paraD2=6.9 //Diametre des paranodes (um)
deltax=1150 //Longueur separant 2 nodes (um)
paralength2=46 //Longueur des paranodes (um)
nl=120} //Nombre de "myelin lemella" (SI)
if (D_droite==11.5){
g=0.700
axonD=8.1 //Diametre des internodes (um)
nodeD=3.7 //Diametre des noeuds (um)
paraD1=3.7 //Diametre du myelin um)
paraD2=8.1 //Diametre des paranodes (um)
deltax=1250 //Longueur separant 2 nodes (um)
paralength2=50 //Longueur des paranodes (um)
nl=130} //Nombre de "myelin lemella" (SI)
if (D_droite==12.8){
g=0.719
axonD=9.2 //Diametre des internodes (um)
nodeD=4.2 //Diametre des noeuds (um)
paraD1=4.2 //Diametre du myelin (um)
paraD2=9.2 //Diametre des paranodes (um)
deltax=1350 //Longueur separant 2 nodes (um)
paralength2=54 //Longueur des paranodes (um)
nl=135} //Nombre de "myelin lemella" (SI)
if (D_droite==14.0){
g=0.739
axonD=10.4 //Diametre des internodes (um)
nodeD=4.7 //Diametre des noeuds (um)
paraD1=4.7 //Diametre du myelin (um)
paraD2=10.4 //Diametre des paranodes (um)
deltax=1400 //Longueur separant 2 nodes (um)
paralength2=56 //Longueur des paranodes (um)
nl=140} //Nombre de "myelin lemella" (SI)
if (D_droite==15.0){
g=0.767
axonD=11.5 //Diametre des internodes (um)
nodeD=5.0 //Diametre des noeuds (um)
paraD1=5.0 //Diametre du myelin (um)
paraD2=11.5 //Diametre des paranodes (um)
deltax=1450 //Longueur separant 2 nodes (um)
paralength2=58 //Longueur des paranodes (um)
nl=145} //Nombre de "myelin lemella" (SI)
if (D_droite==16.0){
g=0.791
axonD=12.7 //Diametre des internodes (um)
nodeD=5.5 //Diametre des noeuds (um)
paraD1=5.5 //Diametre du myelin (um)
paraD2=12.7 //Diametre des paranodes (um)
deltax=1500 //Longueur separant 2 nodes (um)
paralength2=60 //Longueur des paranodes (um)
nl=150} //Nombre de "myelin lemella" (SI)
//resistivite de l'axone (ohm.cm)
Rpn0=(rhoa*.01)/(PI*((((nodeD/2)+space_p1)^2)-((nodeD/2)^2)))
}
CellBuild[0].continuous=0
proc setD_droite(){
D_droite = $1
affectation_droite()
for (i=0;i<=6;i=i+2) dend[i]{
nseg=100 //nombre de segments (SI)
diam=nodeD //diametre (um)
L=100*nodelength //longueur (um)
Ra=rhoa/10000 //resistivite (ohm)
cm=2 //capacite du noeud (uF/cm²)
insert hh //ajout parametre Hodgkin-H
insert extracellular //ajout parametre Extra-cellul
xraxial=Rpn0 //resistivite longitudinal
xg=1e10 //conductance
xc=0} //capacite
}
//titre du tableau
xpanel("Section Fibres Droite", 0)
//case pour regler le retard
xvalue("Diametre Fibre Droite (um)", "D_droite", 1, "setD_droite(D_droite)", 0, 0)
//position de la fenetre
xpanel(38,790)
CellBuild[0].continuous=1
But i don't want this decrease and i want to obtain "the nervous conduction speed".
Can you help me.