---trafficmanagementmendoza.nlogo

NETLOGO

   set list-Period50 lput (total-delay / total-finished-cars) list-Period50
  ]  
  
 
  set model "Reservation" 
  set granularity 1 
  set period 30
  set beta 0 
  set alpha 0.1
  set spawn-probability-N  .001
  set spawn-probability-S  .001 
  set spawn-probability-E  .0 
  set spawn-probability-W  .0 
  set ticks-by-second 50   
  ask turtles [die]     
  set list-reservation []
  while [spawn-probability-N <= .02]
  [
    start-variables 
    ask turtles [die]
     while [clock <= 20000] [
          go                
     ] 
   set spawn-probability-N spawn-probability-N + .001
   set spawn-probability-S  spawn-probability-N   
   set list-reservation lput (total-delay / total-finished-cars) list-reservation 
  ]
  do-plotting3 
end

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;
to experiment4 
  set lanes  1
  set plotting false
  set model "Reservation" 
  set granularity 1 
  set period 30
  set beta 0 
  set alpha 0.1
  set spawn-probability-N  .001
  set spawn-probability-S  .001 
  set spawn-probability-E  .000 
  set spawn-probability-W  .000  
    set ticks-by-second 50
  setup 
  while [spawn-probability-N <= .02]
  [
    start-variables 
    ask turtles [die]
     while [clock <= 20000] [   
        go                     
     ] 
    do-plotting2a 
   set spawn-probability-N spawn-probability-N + .001
   set spawn-probability-S  spawn-probability-N    
  ]
  
  set granularity 2  
  set period 30 
  set beta 0
  set alpha 0.1
  set spawn-probability-N  .001
  set spawn-probability-S  .001 
  set spawn-probability-E  .001 
  set spawn-probability-W  .001 
  set ticks-by-second 50
  ask turtles [die]     
  while [spawn-probability-N <= .02]
  [
    start-variables 
    ask turtles [die]
     while [clock <= 20000] [
        go        
     ]  
    do-plotting2a
   set spawn-probability-N spawn-probability-N + .001
   set spawn-probability-S  spawn-probability-N   
  ]  
end

;;;;;;;;;;;;;;;;;;;;;;;
;;PLOTTING PROCEDURES;;
;;;;;;;;;;;;;;;;;;;;;;;
to do-plotting1
  set-current-plot "Average Delay"  
  set-current-plot-pen "default" 
  ifelse total-finished-cars > 0 [ plot total-delay / total-finished-cars]
  [plot 0]  
  set-current-plot "Maximum Delay" 
  plot maximum-delay 
end

to do-plotting 
  set-current-plot "Figure a"  
  set-current-plot-pen "Theoretical" 
  plot expected-average-delay
  set-current-plot-pen "Empirical"
  plot total-delay / total-finished-cars

  set-current-plot "Figure b" 
  set-current-plot-pen "Theoretical" 
  plot expected-maximum-delay
  set-current-plot-pen "Empirical"
  plot maximum-delay 
end


to do-plotting2a
  set-current-plot "Figure c"  
  set-current-plot-pen "Granu. 1" 
  plot total-delay / total-finished-cars
end

to do-plotting2b
  set-current-plot "Figure c"  
  set-current-plot-pen "Granu. 2" 
  plot total-delay / total-finished-cars
end

to do-plotting3
  set-current-plot "Figure d"  
  foreach list-period10 [
    set-current-plot-pen "Period 10" 
    plot ?  
  
    set-current-plot-pen "Period 30" 
    plot first list-period30
    set list-period30 but-first list-period30
  
    set-current-plot-pen "Period 50" 
    plot first list-period50
    set list-period50 but-first list-period50
  
    set-current-plot-pen "Reservation" 
    plot first list-reservation  
    set list-reservation but-first list-reservation
  ]
end
@#$#@#$#@
GRAPHICS-WINDOW
176
87
587
519
200
200
1.0
0
10
1
1
1
0

CC-WINDOW
5
541
1050
636
Command Center

BUTTON
8
10
107
43
Setup
setup
NIL
1
T
OBSERVER
T
NIL

BUTTON
9
46
92
79
Go
go
T
1
T
OBSERVER
T
NIL

MONITOR
602
10
659
59
NIL
clock
3
1

PLOT
592
162
815
313
Figure a
Period or Alpha
Average Delay
0.0
0.0
0.0
10.0
true
true
PENS
"Theoretical" 1.0 0 -16776961 true
"Empirical" 1.0 0 -65536 true

BUTTON
110
11
176
44
go once
go
NIL
1
T
OBSERVER
T
NIL

CHOICE
10
97
148
142
model
model
"Overpass" "Trafic Light" "Reservation"
2

SLIDER
10
143
173
176
lanes
lanes
1
6
3
1
1
NIL

SLIDER
8
176
173
209
spawn-probability-N
spawn-probability-N
0.0
1.0
0.015
0.0010
1
NIL

MONITOR
603
111
670
160
Total Delay
total-delay\n
3
1

SLIDER
12
316
170
349
period
period
10
200
30
1
1
NIL

SLIDER
12
350
171
383
alpha
alpha
0
1
0.45
0.05
1
NIL

SLIDER
12
384
171
417
beta
beta
0
1
0.1
0.1
1
NIL

MONITOR
795
10
890
59
Finished Cars
total-finished-cars
3
1

MONITOR
671
111
768
160
Average Delay
total-delay / total-finished-cars 
3
1

BUTTON
13
461
90
494
Experiment 1
experiment1
NIL
1
T
OBSERVER
T
NIL

PLOT
817
162
1041
312
Figure b
Period or Alpha
Maximum Delay
0.0
0.0
0.0
10.0
true
true
PENS
"Theoretical" 1.0 0 -16776961 true
"Empirical" 1.0 0 -65536 true

BUTTON
96
462
167
495
Experiment 2
experiment2
NIL
1
T
OBSERVER
T
NIL

SLIDER
12
419
171
452
granularity
granularity
1
50
2
1
1
NIL

MONITOR
768
60
870
109
Expec. Max. Delay
expected-maximum-delay
3
1

MONITOR
671
61
768
110
Expec. Avg. Delay
expected-average-delay
3
1

MONITOR
770
110
844
159
Max. Delay
maximum-delay
3
1

MONITOR
671
10
725
59
NIL
total-cars
3
1

MONITOR
727
10
794