Nematodes

;Every tick is a generation, where each organism reproduces asexually, and its offspring inherit its fertility

breed[males a-male]
breed[hermaphrodites a-hermaphrodites]

globals
[
  ;fertility is used in a slider and can determine probability of turtles

  min-fertility
  max-fertility
  max-fertility-range
  degree-of-fertility
  fertility-rate
]

turtles-own
[
  energy    ;energy for each kind of worm
  male?     ; if true, worm born is male
  hermaphrodite? ; if true, the worm born is hermaphrodite
  fertility ; the whole number of fertility
  fertility-remainder ; the fractional part
]

to setup
  clear-all
   setup-worms
   set min-fertility 1
   set max-fertility-range 3
   reset-ticks
   set fertility-rate(10 ^ log-fertility-rate)
   set degree-of-fertility 0.1
   setup-worms
end 

to setup-worms
  clear-patches
  clear-turtles
  clear-all-plots
  clear-ticks


  ;sets population of two distinct kind of worms randomly across viewer


  create-males carrying-capacity-for-males

  [
    let temp-fertility min-fertility
    setxy random-xcor random-ycor         ; make and set random turtle locations
    set energy 10
    set color blue
    set size 2
    set male? false
    set hermaphrodite? false
    set shape "caterpillar"
    set fertility 1
    set fertility-remainder 0


    ; each worm has a 50% chance of starting out male.

     if (random-float 100 < 50)
    [
      set male? true
      set hermaphrodite? false

    ]

  ]

  create-hermaphrodites carrying-capacity-for-hermaphrodites ;same as initial number
  [
    let temp-fertility min-fertility
    setxy random-xcor random-ycor         ; make and set random turtle locations
    set energy 10
    set color red
    set size 2
    set male? false
    set hermaphrodite? false
    set shape "caterpillar"
    set fertility 1
    set fertility-remainder 0
     ; each worm has a 50% chance of starting out hermaphrodite

     if (random-float 100 < 50)
    [
      set male? false
      set hermaphrodite? true

    ]




  ]


  reset-ticks
end 

to go
  reproduce
  cease
  cease2

  tick
end 

;preference for hermaphrodites on and off switch
;don't mix with the breed function and variable on Netlogo, this means breed with males or not

to breed-males
  ask hermaphrodites
  [ if pcolor = blue [ set pcolor red]] ; asking turtles and turtle procedure
  ifelse preference-for-males?
      [set label males]
      [set label ""]
end 

to assign-color ; turtle procedure
  if male?
    [set color blue]
  if hermaphrodite?
    [set color red]
end 

to wiggle  ;; turtle procedure
  rt random 40
  fd random 40
  if not can-move? 1 [ rt 180 ]
end 


; reproduction for the nematodes

to reproduce
ask Hermaphrodites
  [

    ifelse (random-float 100) < (100 )
      [ hatch  1 [ wiggle ]]
      [ hatch   1 [wiggle ]]
    set energy 10
  ]
end 


; kill worms in excess of carrying capacity. This limits the number of worms on viewer.
; note that red and blue worms have equal probability of dying

to cease2   ;When there are too many hermaphrodites in the area
  let num-turtles count turtles
  if num-turtles <= carrying-capacity-for-hermaphrodites
    [ stop ]
  let chance-to-die (num-turtles - carrying-capacity-for-hermaphrodites) / num-turtles
  ask turtles
  [
    if random-float 1.0 < chance-to-die
      [ die ]
  ]
end 

to cease     ; when there are too many males in the area
  let num-turtles count turtles
  if num-turtles <= carrying-capacity-for-males
    [ stop ]
  let chance-to-die (num-turtles - carrying-capacity-for-males) / num-turtles
  ask turtles
  [
    if random-float 1.0 < chance-to-die
      [ die ]
  ]
end 


;worms can die from overpopulation or insufficient amount of energy

to die-by-energy
  if energy <=  0 [ die ] ; When energy dips at or below zero, the worm dies
end 
;;problem with males reproducing, males cannot reproduce on their own... Only hermaphrodites