The Seacow & Plastic in Ocean Currents

breed [ particles particle ]  ;; the things that are affected by the field
breed [ vectors vector ]    ;; the vectors that affect the particles
breed [ seacows seacow ]
breed [ seaelefants seaelefant ]


globals
[
  max-modulus  ;; the maximum modulus of all the vectors
  step-size    ;; the amount a particle moves forward
  zufall-head
  zufall-border
  zufall-himmelsrichtung
  FULL
  OK
  count-caught
  money
  seacow-number
  plastic-total
  A
  caught
]

seacows-own [volume volume-status]
vectors-own [modulus]  ;; the length of the vector
particles-own [weight]
seaelefants-own [elefant-volume]

; ##################### SETUP ########################

to setup
  clear-all
  setup-vectors
  add-initial-plastic
  add-seacow
  add-seaelefants
  add-marks
  set plastic-total 73400000
  set caught 0
  reset-ticks
end 


; ##################### GO ########################

to go
  ;if ticks >= 30000 [ stop ]
  let stop? false

  add-river-plastic
  add-plastic-flow
  move-particles

  seacow-volume-check

  ask seacows [
    if volume-status = "OK"
    [move-seacows
      collect-plastic
    ]
    if volume-status = "FULL"
    [go-seaelefant]
  ]

  ask seacows [ fd 4 / 10 ]

  move-seaelefants
  ask seaelefants [ fd elefant-speed / 10 ]

  tick
  ;; if one of the particles was going to wrap around the world, stop.
  if stop?
  [ stop ]
end 

; ##################### SEACOWS ########################

to add-seacow
  create-seacows seacows-number
  [
    setxy -117 -117
    set size 5
    set heading 0
    set color white
    set volume 0
    set seacow-number seacow-number + 1
  ]
end 

to move-seacows

  if search-algo = "random"
    [
      search-random
  ]

  if search-algo = "circuit seaelefant"
  [
    ifelse distance one-of seaelefants > max-distance-seacow and 2000 >= volume
    [face one-of seaelefants]
    [rt random 10
      lt random 10
    ]
  ]
end 

to search-random
  ifelse distancexy 0 0 > max-distance-seacow and volume-status = "OK"
  [facexy 0 0]
  [rt random 10
    lt random 10]
end 

to with-the-flow
  if zufall-head > 1
    [ set heading (atan force-x force-y) ]
  if patch-ahead 1 != nobody
    [ forward stroemung / 20 ]
end 

to seacow-volume-check
  ask seacows [
    if volume < 2000
      [set volume-status "OK"]
    if volume >= 2000
      [set volume-status "FULL"
        set label-color red ]

    set label precision volume 2
  ]
end 

to collect-plastic
  ask seacows
  [
    let parthere particles-here
    let total-weight-collected 0
    ask parthere
    [
      let weight-loss (weight * 0.0004)
      set weight weight - weight-loss
      if weight < 10 [die]

      set total-weight-collected total-weight-collected  + weight-loss

    ]
    set volume volume + total-weight-collected
  ]
end 

to go-seaelefant
  ask seacows with [volume >= 2000]
  [face one-of seaelefants ]
  if distance one-of seaelefants < 1
    [ set caught caught + volume
      set plastic-total plastic-total - volume
      set volume-status "OK"
      set volume 0
      set label-color white
      facexy 0 0
      output-show one-of ["fuckin' plastic" "ahoi!" "plastic kills!" "love ya all" "anybody out there?" "attention. some big waves comin in" "save the planet!" "I am so lonely out here" "Arrrrgh!!!" "May day, whiskey bottle empty" "Plastic-free zone"]
  ]
end 


; #################### SEAELEFANT ######################

to add-seaelefants
  create-seaelefants 1
  [
    setxy 40 0
    set size 5
    set heading 0
    set label elefant-volume
    set color yellow
    set elefant-volume 0
  ]
end 

to move-seaelefants
  ask seaelefants
    [
      ifelse distancexy 0 0 > max-dist-center-elefant
      [facexy 0 0]
        [rt random 10
        lt random 10]
  ]

  ;set A max-distance-to-center-or-seaelefant * max-distance-to-center-or-seaelefant * 3.14
  ;ifelse (count particles with [distance one-of seacows < max-distance-to-center-or-seaelefant ] < min-conc-before-chance / A)
  ;[move-to patch-left-and-ahead 45 1]
  ;[facexy 0 100]
  ; [set heading (atan force-x force-y)
end 


; ##################### ADD PLASTIC ########################

to add-initial-plastic

  create-particles amount-classI
  [
    let xpos random-normal 0 world-width * 0.12
    let ypos random-normal 0 world-height * 0.12

    set xpos max list (min list xpos max-pxcor) min-pxcor
    set ypos max list (min list ypos max-pycor) min-pycor

    setxy xpos ypos
    set size 0.1
    set color 125
    set weight 310                                                             ;kg pro Agent
  ]

  create-particles amount-classII
  [
    let xpos random-normal 0 world-width * 0.12
    let ypos random-normal 0 world-height * 0.12

    set xpos max list (min list xpos max-pxcor) min-pxcor
    set ypos max list (min list ypos max-pycor) min-pycor

    setxy xpos ypos
    set size 0.3
    set color 115
    set weight 860
  ]


  create-particles amount-classIII
  [
    let xpos random-normal 0 world-width * 0.12
    let ypos random-normal 0 world-height * 0.12

    set xpos max list (min list xpos max-pxcor) min-pxcor
    set ypos max list (min list ypos max-pycor) min-pycor

    setxy xpos ypos
    set size 0.8
    set color 105
    set weight 8250
  ]

  create-particles amount-classIV
  [

    let xpos random-normal 0 world-width * 0.12
    let ypos random-normal 0 world-height * 0.12

    set xpos max list (min list xpos max-pxcor) min-pxcor
    set ypos max list (min list ypos max-pycor) min-pycor

    setxy xpos ypos
    set size 1.2
    set color 95
    set weight 1060
  ]
end 

to add-plastic-flow
  create-particles plastic-flow-random
  [
    set zufall-himmelsrichtung random 4
    if zufall-himmelsrichtung = 0
    [setxy 200 0]
    if zufall-himmelsrichtung = 1
    [setxy -200 -0]
    if zufall-himmelsrichtung = 2
    [setxy -0 200]
    if zufall-himmelsrichtung = 3
    [setxy 0 -200]
    set size 1
    set heading 0
    set color yellow
  ]
end 

to add-river-plastic
  create-particles plastic-flow-river
  [
    setxy -200 0
    set size 1
    set heading 0
    set color yellow
  ]
end 

to move-particles
  ask particles
  [
    set zufall-head random 100

    if zufall-head < pull-effect
      [ face patch 0 0
        forward stroemung / 3 ]

    if zufall-head < move-randomness
      [set heading random 360
        forward stroemung / 10 ]

    set heading (atan force-x force-y)
    forward stroemung / 10
  ]
end 


; ##################### WORLD ########################

;; report true if we will wrap around if we move forward by step-size

to-report going-to-wrap?  ;; turtle procedure
  let next-patch patch-ahead step-size
  report next-patch = nobody
end 

to add-marks
  ask patches at-points [ [50 0] [50 1] [50 2] [50 -1] [50 -2] [100 0] [100 1] [100 2] [100 -1] [100 -2] ]
  [set pcolor white]
  ask patches at-points [ [0 0] [0 1] [0 2] [0 -1] [0 -2] [1 0] [2 0] [-1 0] [-2 0] ]
  [set pcolor white]
  ask patches at-points [ [-119 -119] [-119 -118] [-119 -117] [-118 -118] [-117 -119] [-117 -118] [-117 -117]]
  [set pcolor white]
end 

; ##################### CREATE VECTOR FIELD ########################

;; create vectors at regular intervals to see the effect of the force
;; at a particular place.

to setup-vectors
  ask patches
  [
    if (pxcor mod 13 = 0) and (pycor mod 13 = 0)
    [
      sprout-vectors 1
      [ setup-vector ]
    ]
  ]
  set max-modulus (max [modulus] of vectors)                ;; draw vector field
  ask vectors [ show-vector ]
end 

;; calculate the horizontal force where the turtle is located

to-report force-x  ;; turtle procedure
  report ycor
end 

;; calculate the vertical force where the turtle is located

to-report force-y  ;; turtle procedure
  report (- xcor)
end 

;; draw the vector using a turtle to display strength and direction of field

to show-vector  ;; turtle procedure
  set modulus (10 * modulus / max-modulus)
  forward modulus
  set color yellow
end 

;; make the turtle become a vector and initialize the vector's variables

to setup-vector  ;; turtle procedure
  set color blue
  pen-down
  if (force-x != 0) or (force-y != 0)
    [ set heading atan force-x force-y ]
  set modulus distancexy 0 0
end 

; Copyright 2017 Dominik Huter