Rappaport Model

globals [
  resources        ;; natural resources available
  ritual-threshold ;; threshold to trigger rituals
  pig-count        ;; count of pigs
  human-count      ;; count of humans
  ritual-count     ;; count of rituals performed
  fertility-stage  ;; current fertility stage
  total-rituals    ;; total count of rituals performed
  rituals-per-tick ;; average rituals per tick
]

turtles-own [
  pig?             ;; indicates if the turtle is a pig
]

patches-own [
  resource-level   ;; resource level on each patch
]

to setup
  clear-all
  ;; Initial setup of agents and patches
  set-default-shape turtles "person"

  ;; Create the initial number of humans
  create-turtles init_humans [
    setxy random-xcor random-ycor
    set color black
    set pig? false
  ]
  ;; Create the initial number of pigs
  create-turtles init_pigs [
    setxy random-xcor random-ycor
    set color pink
    set shape "pig"
    set pig? true
  ]
  ;; Assign random resource levels to patches
  ask patches [
    set resource-level random 100
    patch-recolor
  ]
  ;; Calculate the total resources
  set resources sum [resource-level] of patches
  ;; Set the ritual threshold
  set ritual-threshold threshold_resources
  ;; Initial counts of pigs and humans
  set pig-count count turtles with [pig?]
  set human-count count turtles with [not pig?]
  ;; Reset the tick counter
  reset-ticks
  ;; Initialize the total ritual counter
  set total-rituals 0
  ;; Initialize population graphs
  initialize-population-plots
end 

to patch-recolor ;; patch procedure
  ;; Color patches based on the amount of resources they have
  set pcolor (green + 4.9 - (resource-level / 25))
end 

to go
  ;; Execute simulation steps
  if agricultural-cycles [
    switch-fertility-stages
  ]
  grow-resources
  grow-humans
  grow-pigs
  consume-resources
  check-ritual
  ;; Update pig and human counts
  set pig-count count turtles with [pig?]
  set human-count count turtles with [not pig?]
  ;; Recolor the patches
  ask patches [patch-recolor]
  ;; Calculate the average rituals per tick
  if ticks > 0 [
    set rituals-per-tick total-rituals / ticks
  ]
  ;; Advance one tick
  tick
  ;; Update population graphs
  update-population-plots
end 

to grow-resources
  ;; Increase resources on all patches based on fertility stage and growth rate
  ask patches [
    if agricultural-cycles [
      if fertility-stage = "high" [
        set resource-level resource-level + (2 * resource-growth-rate)
      ]
      if fertility-stage = "low" [
        set resource-level resource-level + (1 * resource-growth-rate)
      ]
    ]
    if not agricultural-cycles [
      set resource-level resource-level + resource-growth-rate
    ]
    ;; Ensure resources do not go below zero
    if resource-level < 0 [
      set resource-level 0
    ]
  ]
  ;; Recalculate total resources
  set resources sum [resource-level] of patches
end 

to switch-fertility-stages
  ;; Switch fertility stages every 50 ticks
  if ticks mod 50 = 0 [
    ask patches [
      ifelse fertility-stage = "high" [
        set fertility-stage "low"
      ] [
        set fertility-stage "high"
      ]
    ]
  ]
end 

to grow-humans
  ;; Grow the human population if it's below 100 (arbitrary limit)
  if (count turtles with [not pig?]) < 100 [
    create-turtles 1 [
      setxy random-xcor random-ycor
      set color black
      set pig? false
    ]
  ]
end 

to grow-pigs
  ;; Grow the pig population if it's below 200 (arbitrary limit)
  if (count turtles with [pig?]) < 200 [
    create-turtles 1 [
      setxy random-xcor random-ycor
      set color pink
      set shape "pig"
      set pig? true
    ]
  ]
end 

to consume-resources
  ;; Consume resources by all turtles
  ask turtles [
    if any? patches with [resource-level > 0] [
      let target-patch one-of patches with [resource-level > 0]
      ask target-patch [
        set resource-level max (list 0 (resource-level - 10))
        set resources resources - 10
      ]
    ]
  ]
end 

to check-ritual
  ;; Check if a ritual should be performed
  if resources < ritual-threshold [
    perform-ritual
  ]
end 

to perform-ritual
  ;; Perform a ritual
  let pigs-to-kill count turtles with [pig?] / 2
  set ritual-count ritual-count + 1
  set total-rituals total-rituals + 1
  ;; Sacrifice pigs
  ask n-of pigs-to-kill turtles with [pig?] [
    die
  ]
  ;; Increase resources due to the ritual
  set resources resources + pigs-to-kill
  ;; Check if the human population should decrease
  if resources < ritual-threshold [
    perform-human-death
  ]
end 

to perform-human-death
  ;; Decrease the human population (one-fifth of humans die if resources are below the threshold)
  let humans-to-kill count turtles with [not pig?] / 5
  ask n-of humans-to-kill turtles with [not pig?] [
    die
  ]
end 

to initialize-population-plots
  ;; Initialize population graphs
  set-current-plot "Population"
  clear-plot
  set-current-plot-pen "Humans"
  set-plot-pen-color blue
  set-current-plot-pen "Pigs"
  set-plot-pen-color red
end 

to update-population-plots
  ;; Update population graphs
  set-current-plot "Population"
  set-current-plot-pen "Humans"
  plot count turtles with [not pig?]
  set-current-plot-pen "Pigs"
  plot count turtles with [pig?]
end