---pyramidal.nlogo

来自「NETLOGO」· NLOGO 代码 · 共 604 行 · 第 1/2 页

breeds [ nmda                                    ; these are the agents in the model
cal                                              ; cal are calcium ions, turtles are magnesium ions
glu1
glu2                                              ; glu are the glutamate molecules
ampa                                           
camk
]
globals [frequency stopper
 response perm ]

to setup
ca
crt 6 
ask turtles [neuron]                             ; drawing the pyramidal cell
channels                                         ; setting up the nmda receptor-channels
axons
glut1
glut2
calcium                                          ; these are the calcum ions
cam
set perm 0
set frequency 0
end

to move1

ifelse pcolor-of patch-here = blue
 [stop][ifelse distance patch 16 25 < distance patch 13 40 
[set heading towards-nowrap patch 16 25 fd 1]
[set heading towards-nowrap patch 13 40 fd 1]]


end

to move2
ifelse pcolor-of patch-here = blue
 [stop][ifelse distance patch -16 25 < distance patch -13 40 
[set heading towards-nowrap patch -16 25 fd 1]
[set heading towards-nowrap patch -13 40 fd 1]]
end

to move3

ifelse pcolor-of patch-here != blue
  [set heading towards-nowrap patch -13 40 fd 1]
    [ifelse any turtles-here with [color = white]
        [set heading random -160 jump random 10] 
        [set heading random 160 jump 10  bind stop]
    ] 
end
to bind
ask camk [move4]
die
end

to move4
if pcolor-of patch-here != blue
[
ifelse distance patch -16 25 < distance patch 13 40 
[set heading towards-nowrap patch -16 25 jump random (distance patch -16 25)
 if pcolor-of patch-here = blue [per die]]
[set heading towards-nowrap patch 13 40 jump random (distance patch 13 40)
if pcolor-of patch-here = blue [per die]]
]
end

to per
if perm < (strength - 1) [set perm (perm + .1)  if perm >= strength [set perm (strength - 1)]]
end

to freq
if frequency < strength [set frequency (frequency + count turtles-here with [color = blue]) 
ask turtles-here with [color = blue] [die] plot response]
end

to action
if frequency >= (strength - perm) [ ask turtles with [who < 3 and who > 0] [ jump 10] ask cal [move3]
                ask turtle 5  
                [ht set shape "box" set size 5
                setxy 0 0
                set heading 180 st
                repeat 200 [set color random 5000 fd .25
                                ] ht
 

                    ]  
                  repeat 5 [plot (-2 * perm)] plot response reuptake]
end

to reuptake
reset-timer 
set frequency 0
set stopper 0
glut1
glut2 
cam
ask turtles with [who < 3 and who > 0] [rt 180 repeat 10 [fd 1 rt 40 lt 40] rt 180]
calcium
end


to axons
ask turtle 4 [ hide-turtle]
ask turtle 5 [hide-turtle]
ask turtle 4 
[setxy 25 35 
paint
set heading 0 jump 5 paint set heading 180 jump 10 paint 
repeat 3 [set heading 45 jump 10 paint set heading 180 jump 3 paint]
setxy -16 24]
ask turtle 5
[setxy -25 35 paint
set heading 0 jump 5 paint set heading 180 jump 10 paint
repeat 3[ set heading (360 - 45) jump 10 paint set heading 180 jump 3 paint]

setxy 16 24 
]
end

to paint
ask patch-here [ set pcolor red]
ask neighbors [set pcolor red ask neighbors 
[set pcolor red ask neighbors [ set pcolor red]]]
end

to neuron
set color red
draw-neuron
end

to channels

ask turtles with [who < 4] [ht set color white set size 3 set shape "box"]
ask turtle 0 [set color cyan setxy -50 -35
set heading 45 
repeat 70[fd 1 lt .45]
ask patch-here [set pcolor blue]
ask neighbors [set pcolor green] set heading random -180]
ask turtle 1 [setxy -50 -35
set heading 45 
repeat 85[fd 1 lt .45]
ask patch-here [set pcolor blue]
ask neighbors [set pcolor green]set heading random -180]
ask turtle 2 [
setxy 50 -35
set heading (360 - 45 )
repeat 70[fd 1 rt .45]
ask patch-here [set pcolor blue]
ask neighbors [set pcolor green]set heading random 180]
ask turtle 3[set color cyan
setxy 50 -35
set heading (360 - 45 )
repeat 85[fd 1 rt .45]
ask patch-here [set pcolor blue]
ask neighbors [set pcolor green]set heading random 180]
ask turtles [show-turtle]
end

to cam
create-camk (strength - (count camk))
ask camk [
set size 1
set color yellow
if xcor = 0 and ycor = 0 [ht setxy 0 32
jump random 15 st]

]
end

to glut1
create-glu1 ((strength / 2)- (count glu1))
ask glu1 [
set shape "circle"
set color blue
set size 3
hide-turtle
setxy 27 35
  set heading random 360 
  jump random 7 show-turtle] 
end

to glut2
create-glu2 ((strength - count glu1) - (count glu2))
ask glu2
[set shape "circle"
set color blue
set size 3
setxy -27 35
set heading random 360
jump random 7 show-turtle]

end

to calcium
create-cal (20 - count cal)
ask cal [

set color green
set size 2
hide-turtle
jump 10
ifelse xcor > 0 
[ setxy -18 30
  set heading random 360 
  jump random 5 show-turtle] 
[setxy -18 40
set heading random 360
jump random 5 show-turtle]
]
end

to draw-neuron
hide-turtle
setxy -50 -35

set heading 45 
repeat 84 [fd 1 lt .45 stamp red]
set heading 0
repeat 10 [ stamp red fd 1]

setxy 50 -35
set heading (360 - 45)

repeat 84 [fd 1 rt .45 stamp red]
set heading 0
repeat 10 [stamp red fd 1]

setxy -50 -43 
set heading 55
repeat 22[fd 1 rt .40 stamp red]
set heading 90
repeat 30 [fd 1 rt 2.5 stamp red]

setxy 50 -43 
set heading (360 - 55)
repeat 22 [fd 1 lt .40 stamp red]
set heading 270
repeat 30 [fd 1 lt 2.5 stamp red]

end

to go
set stopper 1
if perm >= strength [set perm (strength - 1)]
ask glu1 [move1]
ask glu2 [move2]
ask cal  [move3]
ask patches with [pcolor = blue] [freq]
action
if stopper < 1 [stop]
end
@#$#@#$#@
GRAPHICS-WINDOW
600
10
1013
449
50
50
3.9901
1
10
0
1

CC-WINDOW
650
456
975
543
Command Center

BUTTON
35
50
102
83
setup
setup
NIL
1
T
OBSERVER

BUTTON
271
50
334
83
Go
Go
T
1
T
OBSERVER

SLIDER
101
82
273
115
strength