---rgw.nlogo

NETLOGO

globals [field stappen lijst]
turtles-own [ energy age]
breeds [rabbits]

to setup
  ca
  grow-grass-and-weeds
  set-default-shape rabbits "rabbit"
  create-custom-rabbits number [
    set color white
    randomize-position
    set age 0
    set energy random 10  ;start with a random amt. of energy
  ]
  setup-plot
  do-plot
  set field (count patches)
 end

to go
 just-add
 ;refill
 ;grow-grass-and-weeds
 move
  do-plot
 if not any? turtles
 [ stop ]
end

to test
locals [gras onkruid]
ca
set field (count patches)
set lijst []
set gras []
set onkruid []
repeat 10 [
  set stappen 0
  while [((count (patches with [pcolor = black])) / field) > 0.1] 
              [ ;refill
                just-add
                ;grow-grass-and-weeds 
                set stappen (stappen + 1)]
  show stappen
  set lijst (fput stappen lijst)
  set gras (fput (count patches with [pcolor = green ]) gras)
  set onkruid (fput (count patches with [pcolor = violet]) onkruid) 
  cp
 ]
 show mean lijst
 show gras
 show onkruid
end

to een-test
locals [g v]
set g []
set v []
repeat 50 [  
  cp
  ;refill
  ;just-add
  grow-grass-and-weeds
  set g (fput (count patches with [pcolor = green ]) g)
  set v (fput (count patches with [pcolor = violet]) v) 
  ]
show word "grass: " (mean g)
show word "weed: " (mean v)
end

to move
  ask rabbits 
  [ step
    eat-grass
    eat-weeds
    set age (age + 1)
    reproduce
    death 
    ]
end

;ondanks dat het gras over de weeds heen groeit, is bij 90% begroeing de 
;verdeling tussen gras en onkruid 1-1
; en ook hier duurt het +- 76 tellen
; bij een leeg veld wordt er bijna 15 en iets meer dan 14 van resp gras & weeds toegevoegd
to grow-grass-and-weeds
  ask patches [
    if pcolor = black [
      if ( random-float 1000 ) < weeds-grow-rate
        [ set pcolor violet ]
      if ( random-float 1000 ) < grass-grow-rate
        [ set pcolor green ]
  ] 
  ]
end

; refill levert een 1-1 verdeling op van grass en weeds. het veld raakt vol in ongeveer 76 tellen
; bij een leeg veld wordt er inderdaad 16 en 13,5 van resp gras & weeds toegevoegd
to refill
  ask patches [
     if ( random-float 1000 ) < weeds-grow-rate
        [ set pcolor violet ]
      if ( random-float 1000 ) < grass-grow-rate
        [ set pcolor green ]
  ] 
end

;zoalng de refill kleiner is dan 1/10e gaat het goed
to just-add
locals [n]
set n (count patches with [pcolor = black])
ifelse n > 28
  [
    ask random-n-of 15 (patches with [pcolor = black]) [set pcolor green]
    ask random-n-of 13 (patches with [pcolor = black]) [set pcolor violet]
  ]
  [ask random-n-of round (n * (15 / 28)) (patches with [pcolor = black]) [set pcolor green]
   ask random-n-of round (n * (13 / 28)) (patches with [pcolor = black]) [set pcolor violet]   ]
end


to step   ;moving takes some energy
  rt random-float 50
  lt random-float 50
  fd 1
  set energy ( energy - 0.5 )
end

to eat-grass       ;;gain "grass-energy" by eating grass
  if ( pcolor = green )
  [ stamp black
    set energy ( energy + grass-energy ) ]
end

to eat-weeds     ;;gain "weed-energy" by eating weeds
  if ( pcolor = violet )
  [ stamp black
    set energy ( energy + weed-energy )]
end

to reproduce     ;;give birth to a new rabbit, but it takes lots of energy
  if ( energy > birth-threshold )
    [ set energy ( energy / 2 )
      hatch 1 [ fd 1 set age 0] ]
end

to death     ;;die if you run out of energy
  if ( energy < 0 )
    [ plot-age age die ]
end

to plot-age [n]
set-current-plot "age"
set-current-plot-pen "age"
plot n
end

to do-plot
  set-current-plot "Populations"
  set-current-plot-pen "grass"
  plot ( count patches with [pcolor = green] )
  set-current-plot-pen "rabbits"
  plot count turtles
  set-current-plot-pen "weeds"
  plot ( count patches with [pcolor = violet] )
end

to randomize-position
  setxy random-float screen-size-x
        random-float screen-size-y
end

to setup-plot ;; set up plotting
  set-current-plot "Populations"
  set-plot-y-range 0 number
end


; *** NetLogo Model Copyright Notice ***
;
; This model was created as part of the project:
; PARTICIPATORY SIMULATIONS: NETWORK-BASED DESIGN FOR SYSTEMS LEARNING IN
; CLASSROOMS.  The project gratefully acknowledges the support of the
; National Science Foundation (REPP program) -- grant number REC #9814682.
;
; Copyright 2001 by Uri Wilensky.  Updated 2002.  All rights reserved.
;
; Permission to use, modify or redistribute this model is hereby granted,
; provided that both of the following requirements are followed:
; a) this copyright notice is included.
; b) this model will not be redistributed for profit without permission
;    from Uri Wilensky.
; Contact Uri Wilensky for appropriate licenses for redistribution for
; profit.
;
; To refer to this model in academic publications, please use:
; Wilensky, U. (2001).  NetLogo Rabbits Grass Weeds model.
; http://ccl.northwestern.edu/netlogo/models/RabbitsGrassWeeds.
; Center for Connected Learning and Computer-Based Modeling,
; Northwestern University, Evanston, IL.
;
; In other publications, please use:
; Copyright 1998 by Uri Wilensky.  All rights reserved.  See
; http://ccl.northwestern.edu/netlogo/models/RabbitsGrassWeeds
; for terms of use.
;
; *** End of NetLogo Model Copyright Notice ***
@#$#@#$#@
GRAPHICS-WINDOW
296
13
740
478
15
15
14.0
1
10
1
1
1
0

CC-WINDOW
5
525
749
620
Command Center

BUTTON
22
26
77
59
setup
setup
NIL
1
T
OBSERVER
T
NIL

BUTTON
87
26
142
59
go
go
T
1
T
OBSERVER
T
NIL

SLIDER
3
108
141
141
grass-grow-rate
grass-grow-rate
0.0
20.0
20.0
1.0
1
NIL

SLIDER
3
141
150
174
weeds-grow-rate
weeds-grow-rate
0.0
20.0
20.0
1.0
1
NIL

SLIDER
143
108
274
141
grass-energy
grass-energy
0.0
10.0
5.0
0.5
1
NIL

SLIDER
152
141
274
174
weed-energy
weed-energy
0.0
10.0
5.0
0.5
1
NIL

SLIDER
3
76
112
109
number
number
0.0
500.0
90.0
1.0
1
NIL

SLIDER
114
76
274
109
birth-threshold
birth-threshold
0.0
20.0
15.0
1.0
1
NIL

PLOT
3
176
274
372
Populations
Time
Pop
0.0
100.0
0.0
111.0
true
true
PENS
"grass" 1.0 0 -11352576 true
"rabbits" 1.0 0 -65536 true
"weeds" 1.0 0 -8716033 true

MONITOR
51
385
140
434
count rabbits
count rabbits
1
1

BUTTON
149
25
212
58
NIL
test
NIL
1
T
OBSERVER
T