Yu_Thomas_Gas_Station

globals [
  grid-x-b
  grid-y-b
  roads
  homes
  gasstation 
  gas-time
  ]

turtles-own [ 
  car
  energy
  home-x
  home-y
  dest-x
  dest-y
  fdest-x
  fdest-y
  direction
  intersection
  gas-refuel-time
  ]

patches-own [ 
  count-neighbors
  count-homes
  gas-level
  gas-price
  ]

to setup
  clear-all
  
  setup-globals
  setup-patches
  setup-homes
  setup-gasstation
  setup-dest
  
  set-default-shape turtles "car"
  
  ask patches [
    if pcolor = green [
      sprout 1
       ]             ;; sprout a car on each home
    if pcolor = grey [
      set gas-level 2000
      set gas-price random-normal 3.5 1
      ]  ;; set oil level of each gas station
    ]
  
  ask patches [
    if pcolor = grey [
    set plabel gas-level
    set plabel-color black
    ]
  ]
  
  ask turtles [
    set energy 300
    set car turtles
    set home-x [pxcor] of min-one-of (patches with [pcolor = green]) [distance myself]
    set home-y [pycor] of min-one-of (patches with [pcolor = green]) [distance myself]
  
    set car turtles
    set dest-x [pxcor] of one-of (patches with [pcolor = orange])
    set dest-y [pycor] of one-of (patches with [pcolor = orange])
    
    move-to min-one-of (roads with [not any? turtles-on self]) [distance myself] 
    ]
end 

to setup-globals
  set grid-x-b world-width / grid-size
  set grid-y-b world-width / grid-size
end 

to setup-patches
  ask patches [
  set pcolor brown]
  
  set roads patches with [
    ( floor (( pxcor + max-pxcor - floor (grid-x-b - 1)) mod grid-x-b) = 0 ) or
    (floor((pycor + max-pycor) mod grid-y-b) = 0 )
    ]
  ask roads [ set pcolor white ]
end 

to setup-homes
  ask patches [ 
    set count-neighbors count neighbors4 with [pcolor = white] 
    ]
  
  ask patches [
    if pcolor = brown and count-neighbors > 0 [
      if random-float 10 > (10 - home-prob) [
        set pcolor green ] 
      ] 
    ]
  
  set homes patches with [
    pcolor = green 
    ]
end 

to setup-gasstation
  ask patches [ 
    set count-neighbors count neighbors4 with [pcolor = white] ]
  
  ask patches [
    if pcolor = brown and count-neighbors = 1 [
      if random-float 10 > (10 - gas-prob) [
        set pcolor grey ] 
      ] 
    ]
  
  set gasstation patches with [
    pcolor = grey 
    ]
end 

to setup-dest
  ask patches [
    set count-neighbors count neighbors4 with [pcolor = white] 
    set count-homes count neighbors with [pcolor = green ] 
    ]
  
  ask patches [
    if pcolor = brown and pcolor != grey and pcolor != green and count-neighbors > 0 and count-homes = 0 [
      if random-float 10 > (10 - dest-prob) [
        set pcolor orange ] 
      ] 
    ]
end 

to setup-car-home
end 

;to update-gas-price
;  set label " gas-price "
;end

to go
  ask turtles [
    facexy (dest-x)
           (dest-y) 
    
    ifelse energy < 40 [
      ifelse gas-searching-mode = "Mode 1"
      [ move-to-gas-mode1 ] [
        if gas-searching-mode = "Mode 2"
        [ move-to-gas-mode2 ] ]
      ]
      [ 
      ifelse distancexy dest-x dest-y < 1    ;; determine whether the car has reached its destination     
        [ ifelse distancexy home-x home-y < 1 and distancexy home-x home-y = distancexy dest-x dest-y
          [ set dest-x [pxcor] of one-of (patches with [pcolor = orange])                ;; if it is a home, it assignes a new destination
            set dest-y [pycor] of one-of (patches with [pcolor = orange]) ]
          [ set dest-x home-x                ;; if it is a home, it assignes a new destination
            set dest-y home-y ] 
        ]  
        [ move-easy ]  
      ]
  ] ;; if it is an actual destination, it proceed to move closer to it
  
  ask patches [ 
    check-gas-level1
    update-label 
  ]
  
  tick
end      

to move-easy
  facexy (dest-x)
         (dest-y)
         
  if [pcolor] of patch-ahead 1 != white
    [ rt 45 ]
  if [pcolor] of patch-ahead 1 != white
    [ lt 90 ]          
  fd 1
  
  set energy energy - 1   
end 

to move-to-gas-mode1
  
  set fdest-x [pxcor] of min-one-of (patches with [pcolor = grey]) [distance myself]
  set fdest-y [pycor] of min-one-of (patches with [pcolor = grey]) [distance myself]
  
  facexy (fdest-x)
         (fdest-y)
         
;  if [pcolor] of patch-ahead 1 != white
;    [ rt 45 ]
;  if [pcolor] of patch-ahead 1 != white
;    [ lt 90 ]   
        
  fd 1
  
  set energy energy - 1
  
  if distancexy fdest-x fdest-y < 1
    [ set gas-level (gas-level - (300 - energy)) 
      set energy 300 
    ]
end 

to move-to-gas-mode2
  if energy < 40 [
    let cheapest-gas count patches in-radius (10) with [pcolor = gray] 
    ifelse cheapest-gas < 2
      [ set fdest-x [pxcor] of min-one-of (patches with [pcolor = grey]) [distance myself]
        set fdest-y [pycor] of min-one-of (patches with [pcolor = grey]) [distance myself] ]
      [ set fdest-x [pxcor] of min-one-of (patches in-radius 10 with [pcolor = grey]) [gas-price]
        set fdest-y [pycor] of min-one-of (patches in-radius 10 with [pcolor = grey]) [gas-price] ]
  ]
  
  facexy (fdest-x)
         (fdest-y)

;  if [pcolor] of patch-ahead 1 != white
;    [ rt 45 ]
;  if [pcolor] of patch-ahead 1 != white
;    [ lt 90 ]  
  
  fd 1
  
  set energy energy - 1
  
  if distancexy fdest-x fdest-y < 1
    [ set gas-level (gas-level - (300 - energy)) 
      set energy 300 
    ]
end 

to check-gas-level1
  if gas-level < 300 [
    set gas-level 2000    ;; gas is replenished
    set gas-time ticks
    set gas-price ( gas-price + (ticks - 1000) / 1000 )]  ;; formula needs adjustment
end          

to update-label
  if pcolor = grey [
    set plabel gas-level ]
end