;*********************************************
;Author: Micah Brown
;Model: Maker/Moocher Cells

;*********************************************


breed [ cells cell ]
breed [ theo theophylline-molecule ]

turtles-own [
  energy 
  lifespan
  time-since-reproduction
]

;SETUP

to setup
  ca
  setup-background
  setup-turtles ;turtles set up in random location, ensuring no overlaps. 
  reset-ticks
end 

to setup-background
  ask patches [
    set pcolor white
  
  ]
end 


;------------------------------

to setup-turtles

  
  set-default-shape cells "circle 2"
  set-default-shape theo "lightning"

  
  create-cells (initial-amount-makers + initial-amount-moochers) [
    set color black
    setxy  random-xcor random-ycor
    set size 10
    ifelse who < initial-amount-makers
      [set color red
       set time-since-reproduction random maker-reproduction-downtime]
      [set color cyan
       set time-since-reproduction random moocher-reproduction-downtime]
    set energy initial-energy
    set lifespan initial-lifespan
     
     
  ]
  
  ;overlaps handled by killing overlaping cell :/
  ask turtles [
    if (count turtles in-radius (size + 1) > 1) [
      die
    ]
  ]
  
  create-theo initial-theophylline [
    set color black
    setxy random-xcor random-ycor
    set size 3
    set lifespan initial-lifespan / ( 2 / 3 )
    
  
  ]
end 

;------------------------------

to go
  ;stops when everyone has died
  if (count cells = 0)
    [stop]
  
  move-turtles
  degredation
  reproduce 
  tick
end 

;------------------------------

to move-turtles
  
  ask theo [
    
    let nearest-cell min-one-of other cells [distance myself]  
    let nearest-cell-distance distance nearest-cell
    
    ;test whether theophylline is in "membrane" 4
    if ( nearest-cell-distance <= 6 ) and ( nearest-cell-distance) >= 4 [ 
      ifelse count cells in-cone 6.0 180 < 1
        [if random-float 1 > probability-diffusion-out ;if facing (and moving) away from cell
          [ face nearest-cell
            fd 2 
          ] 
        ]
        [if random-float 1 > probability-diffusion-in ;else must be facing (and moving) towards a cell
          [ face nearest-cell 
            set heading heading + 180
            fd 1.5 
          ] 
        ]
    ]
    random-movement-theo
    set lifespan lifespan - 1
  ] 

  ask cells [
    random-movement-cell
    if count cells in-radius 10.1 > 1 [  ;ensures cells don't "overlap"
      rt 180 
      fd 1
    ]
      
    if color = red [ ;only makers produce theophylline
      if random-float 1 <= theo-production / 2 [

        hatch-theo 1 [
          set size 3 
          set color black
          set lifespan initial-lifespan / ( 2 / 3)  ;we give theophylline a lifespan equal to (1.5) of a cell's.
          setxy ( xcor + (random-float (4 * sqrt 2) ) - (2 * sqrt 2) ) (ycor + (random-float (4 * sqrt 2) ) - (2 * sqrt 2) )
        ]
      ]    
    ]
    if count theo in-radius 6 > 0
      [ set energy energy + log ( count theo in-radius 6 ) 2 + 1 ] ;energy gained from having theophylline is scaled logarithmically 
    
    if energy >= 100 [ ;maximum energy is 100
      set energy 100 
    ]
   
   
    set energy energy - .5
    set lifespan lifespan - 1
    set time-since-reproduction time-since-reproduction + 1
  ] 
end 

;--------------------------------------

to degredation
  
  ask cells [
    if energy <= 0 [
      if random 10 = 1 [die] ;randomness gives more variation
    ]
  ]
  ask turtles [
    if lifespan <= 0 [
      if random 10 = 1 [die]
    ]
  ]
end 

;-----------------------------

;local randomness preserved in order to have theophylline inside of cell move the same way the cell moves.

to random-movement-cell
  with-local-randomness [
    ifelse random 2 = 1
      [right random 30 forward .5]
      [right random -30 forward .5]
  ]
  ask other theo in-radius 6 [     ; moves its theophylline contained in cell in the same direction that it just moved
    ifelse random 2 = 1
      [right random 30 forward .5]
      [right random -30 forward .5]]
end   

;------------------------------

to random-movement-theo
  
    ifelse random 2 = 1
      [right random 30 forward 1 ] 
      [right random -30 forward 1 ]
end  

;--------------------------------

to reproduce
  ;this is an approximate for how many cells will fit into the space given the coordinates that the user sets.
  let carrying-capacity (max-pxcor * max-pxcor) / 50
  if (count cells > carrying-capacity) [ 
    stop
  ]
  
  
  ask cells [
    if ( (color = red) and (energy >= maker-reproduction-energy-threshold) and time-since-reproduction >= maker-reproduction-downtime ) or
       ( (color = cyan) and (energy >= moocher-reproduction-energy-threshold) and (time-since-reproduction >= moocher-reproduction-downtime ) ) [
      set time-since-reproduction 0
      let child-color color
      let half-energy energy / 2
      set energy half-energy
      let num-theo count theo in-radius 6
      
      ;kills half of the theophylline in the cell (but "adds" them to the new cell, later)
      repeat num-theo * 0.5 [
        ask one-of theo in-radius 6 [die]    
      ]
        
      hatch-cells 1 [ 
        set color child-color
        set lifespan initial-lifespan
        set energy half-energy
        set time-since-reproduction 0
        while [count cells in-radius (size + 1) > 1] [ ;attempts to set child cell close to parent cell
          setxy xcor + ((random count cells) - (count cells) / 2 ) ycor + ((random count cells) - (count cells) / 2 )
        ] 
        hatch-theo ( num-theo / 2 ) [
          set size 3
          set color black
        ]
      ] 
    ]
  ]  
end 

;-----------------------------------