Matriz Agroecológica

patches-own [
  tipo
  nuevos_organismos

  prob_mortalidad
  prob_movimiento

  grupo

  vals
  estado
  tiempo_desde_quemado
]

turtles-own [
  vivo?
]

globals [
  num_grupos
  p             ;; parámetro para crear paisajes con el ruido de perlin
]

to setup
;  clear-all
  clear-globals
  clear-ticks
  clear-turtles
  clear-patches
  clear-drawing
  clear-output
  inicializar_graficas

  inicializar_paisaje
  inicializar_poblacion

  colorear_cuadros

  reset-ticks
end 

to go


  ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  ;;;;;;;;;; MODELO FAHRIG ;;;;;;;;;;
  ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  ask patches [
    sprout nuevos_organismos [
      inicialiar_agente
    ]
    actualizar_variables_celda
  ]
  ask turtles with [ vivo? = false ] [ die ]
  ask turtles [
    foreach shuffle [ "mortalidad" "movimiento" "reproduccion" ] [ procedimiento -> run procedimiento ]
  ]
  mortalidad_sobreabundancia
  ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

  if MODIFICAR_MATRIZ? [
    actualizar_matriz_paisaje
  ]

  if random-float 1.0 < PROB_INCENDIOS [
    let origen one-of patches
    iniciar_incendio [pxcor] of origen [pycor] of origen
  ]

  extender_incendio
  ask patches with [ estado = "quemandose" ][ ask turtles-here [ die ]]
  colorear_cuadros


  tick
end 

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;; MODELO FAHRIG ;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

to mortalidad
  if random-float 1.00 < prob_mortalidad [
    set vivo? false
  ]
end 

to movimiento
  if random-float 1.00 < prob_movimiento  [
    if FORMA_MOVIMIENTO = "saltos (fahrig)" [
      set heading random 360
      jump one-of [ 1 2 3 4 ]
    ]
    if FORMA_MOVIMIENTO = "caminata aleatoria 360" [
      set heading random 360
      fd 1
    ]
    if FORMA_MOVIMIENTO = "caminata aleatoria 90" [
      rt random-float 90
      lt random-float 90
      jump 1
    ]
    if FORMA_MOVIMIENTO = "caminata aleatoria 4 vecinos" [
      set heading one-of [ 0 90 180 270 ]
      fd 1
    ]
    if FORMA_MOVIMIENTO = "movimiento browniano" [
      set heading random 360
      fd random-normal 0 1
    ]
    if FORMA_MOVIMIENTO = "vuelo de levy" [
      set heading random 360
      fd .5 * (random-float 1) ^ (-1 / 2.5)
    ]
  ]
end 

to reproduccion
  if tipo = "habitat" and random-float 1.00 < PROB_REPRODUCCION [
    set nuevos_organismos nuevos_organismos + 1
  ]
end 

to mortalidad_sobreabundancia
  ask patches [
    if count turtles-here > CAPACIDAD_MAX_CELDA [
      ask n-of ( count turtles-here - CAPACIDAD_MAX_CELDA ) turtles-here [ die ]
    ]
  ]
end 

to actualizar_variables_celda
  set nuevos_organismos 0

  if tipo = "habitat" [
    set prob_mortalidad PROB_MORTALIDAD_HABITAT
    set prob_movimiento PROB_MOVIMIENTO_HABITAT
  ]

  if tipo = "matriz 1" [
    set prob_mortalidad PROB_MORTALIDAD_MATRIZ_1
    set prob_movimiento PROB_MOVIMIENTO_MATRIZ_1
  ]

  if tipo = "matriz 2" [
    set prob_mortalidad PROB_MORTALIDAD_MATRIZ_2
    set prob_movimiento PROB_MOVIMIENTO_MATRIZ_2
  ]
end 
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

to inicializar_paisaje
  ask patches [
    set grupo -999
    set estado "NA"
    set tiempo_desde_quemado -999
  ]
  ;; se crean los parches de hábitat de acuerdo al algoritmo seleccionado
  if PAISAJE = "solo hábitat" [ inicializar_paisaje_solo_habitat ]
  if PAISAJE = "cuadricula regular" [ inicializar_paisaje_cuadricula_regular ]
  if PAISAJE = "fahrig fragmentado" [ inicializar_paisaje_fragmentado ]
  if PAISAJE = "ruido de perlin" [ inicializar_paisaje_ruido_perlin ]
  if PAISAJE = "bordes" [ inicializar_paisaje_bordes ]

   ;el resto del pasiaje se crea de forma aleatoria
 ask patches with [ tipo != "habitat" ][
   ifelse random-float 1.0 > PROP_MATRIZ_2 [ set tipo "matriz 1" ][ set tipo "matriz 2"]
 ]
  ask patches with [ tipo != "habitat" ][ set tipo "matriz 1" ]
  actualizar_matriz_paisaje

   if PAISAJE = "Zaachila" [ inicializar_paisaje_zaa]

  ask patches [ actualizar_variables_celda ]
  formar_vecindades
end 

to actualizar_matriz_paisaje
  while [ count patches with [ tipo = "matriz 2" ] > round (PROP_MATRIZ_2 * count patches with [ tipo != "habitat" ]) ] [
    ;; hay más celdas de matríz 2 que la que debería
    ask one-of patches with [ tipo = "matriz 2" ] [
      set tipo "matriz 1"
    ]
  ]
  while [ count patches with [ tipo = "matriz 2" ] < round (PROP_MATRIZ_2 * count patches with [ tipo != "habitat" ]) ] [
    ;; hay más celdas de matríz 2 que la que debería
    ask one-of patches with [ tipo = "matriz 1" ] [
      set tipo "matriz 2"
    ]
  ]
end 

to inicializar_paisaje_ruido_perlin
  ask patches [ set vals [] ]
  foreach range OCTAVAS [
    set p shuffle range 256
    ask patches [
      set vals lput (perlin_noise (0.3 * pxcor) (0.3 * pycor)) vals
    ]
  ]
  ask patches [
    if (sum vals > UMBRAL) [ set tipo "habitat" ]
  ]
end 

to inicializar_paisaje_fragmentado
  while [ count patches with [ tipo = "habitat" ] != round (COVERTURA * count patches ) ] [
    ask one-of patches [
      if tipo != "habitat" [
        if random-float 1.00 < FRAGMENTACION or any? neighbors4 with [ tipo = "habitat" ][
          set tipo "habitat"
        ]
      ]
    ]
  ]
end 

to inicializar_paisaje_solo_habitat
  ask patches [ set tipo "habitat" ]
end 

to inicializar_paisaje_bordes

  ; Asignar el tipo "habitat" a las primeras ocho columnas a la izquierda
  ask patches with [pxcor < -10] [
  set tipo "habitat"
]

  ; Asignar el tipo "habitat" a las ultimas ocho columnas
  ask patches with [pxcor > 10] [
  set tipo "habitat"
]
end 

to inicializar_paisaje_cuadricula_regular
  generar_cuadricula_regular_habitat
end 

to generar_cuadricula_regular_habitat
  let tamanio_cuadro LONG_CUADRO
  let tamanio_espacio LONG_ESPACIO

  let contador_espacio_y 0
  let contador_cuadro_y 0

  foreach sort remove-duplicates [pycor] of patches [
    cy ->
    set contador_espacio_y contador_espacio_y + 1
    if contador_espacio_y > tamanio_espacio [

      let contador_espacio_x 0
      let contador_cuadro_x 0

      foreach sort remove-duplicates [pxcor] of patches [
        cx ->
        set contador_espacio_x contador_espacio_x + 1
        if contador_espacio_x > tamanio_espacio [
          ask patch cx cy [
            set tipo "habitat"
          ]
          set contador_cuadro_x contador_cuadro_x + 1
          if contador_cuadro_x > tamanio_cuadro [
            set contador_espacio_x 0
            set contador_cuadro_x 0
          ]
        ]
      ]
      set contador_cuadro_y contador_cuadro_y + 1
      if contador_cuadro_y > tamanio_cuadro [
        set contador_cuadro_y 0
        set contador_espacio_y 0
      ]
    ]
  ]
end 

to  inicializar_paisaje_zaa

  let data [[-9 -7 "matriz 2"] [-13 2 "habitat"] [-11 -7 "matriz 2"] [9 -11 "matriz 2"] [12 16 "matriz 1"] [-1 8 "matriz 2"] [-16 2 "habitat"] [6 16 "matriz 2"] [-3 -4 "matriz 1"] [2 6 "matriz 1"] [-9 8 "habitat"] [-12 14 "habitat"] [-11 -4 "matriz 1"] [0 14 "matriz 2"] [10 10 "matriz 2"] [-6 9 "habitat"] [11 1 "matriz 1"] [11 13 "matriz 2"] [-5 -8 "matriz 1"] [4 -10 "matriz 2"] [0 9 "matriz 2"] [8 11 "matriz 1"] [15 -6 "matriz 2"] [9 -7 "matriz 2"] [-12 0 "matriz 2"] [-11 9 "habitat"] [-14 4 "habitat"] [9 -10 "matriz 2"] [-15 -12 "matriz 2"] [-5 -6 "matriz 1"] [9 -12 "matriz 2"] [14 -13 "habitat"] [6 15 "matriz 2"] [4 -3 "matriz 1"] [-10 9 "habitat"] [-15 -16 "matriz 1"] [-7 8 "habitat"] [-7 10 "matriz 2"] [-11 5 "habitat"] [-2 6 "matriz 1"] [11 0 "matriz 1"] [-15 10 "habitat"] [0 -6 "matriz 1"] [14 -11 "habitat"] [-3 7 "matriz 2"] [-9 9 "habitat"] [5 -14 "matriz 2"] [0 -10 "matriz 2"] [-14 -4 "matriz 2"] [2 16 "matriz 2"] [-10 1 "habitat"] [-11 -9 "matriz 2"] [-14 13 "habitat"] [-12 4 "habitat"] [-14 -3 "matriz 2"] [7 -4 "matriz 1"] [3 -5 "matriz 1"] [-13 15 "habitat"] [7 -16 "matriz 2"] [-13 -11 "matriz 2"] [6 -15 "matriz 2"] [16 3 "matriz 1"] [7 8 "matriz 2"] [-8 15 "habitat"] [-2 13 "matriz 2"] [15 -9 "habitat"] [-11 -14 "matriz 1"] [-14 -16 "matriz 1"] [-7 -11 "matriz 1"] [-16 -13 "matriz 2"] [-11 15 "habitat"] [-9 4 "habitat"] [-16 -8 "matriz 2"] [1 -1 "matriz 1"] [8 6 "matriz 2"] [-10 11 "habitat"] [11 -3 "matriz 2"] [0 15 "matriz 2"] [3 13 "matriz 2"] [-3 14 "matriz 2"] [-16 0 "habitat"] [6 1 "matriz 1"] [-8 -7 "matriz 1"] [-3 11 "matriz 2"] [14 -5 "matriz 2"] [1 4 "matriz 2"] [-2 5 "matriz 2"] [4 -5 "matriz 1"] [-2 -1 "matriz 2"] [14 -12 "habitat"] [-8 5 "matriz 1"] [-5 16 "habitat"] [3 16 "matriz 2"] [7 9 "matriz 2"] [-6 16 "habitat"] [2 -2 "matriz 1"] [-9 -1 "matriz 2"] [12 8 "matriz 2"] [14 11 "matriz 2"] [11 -2 "matriz 1"] [-16 15 "habitat"] [12 6 "matriz 1"] [7 -5 "matriz 2"] [11 5 "matriz 1"] [16 0 "matriz 2"] [-1 -15 "matriz 2"] [9 -14 "habitat"] [13 15 "matriz 2"] [8 -3 "matriz 2"] [7 -3 "matriz 1"] [-8 -13 "matriz 1"] [10 11 "matriz 2"] [-9 -14 "matriz 1"] [11 4 "matriz 1"] [-7 -7 "matriz 1"] [13 16 "matriz 1"] [-8 1 "matriz 2"] [11 -13 "habitat"] [-2 8 "matriz 2"] [11 10 "matriz 2"] [-1 15 "matriz 1"] [-7 -5 "matriz 2"] [-5 15 "habitat"] [-6 3 "matriz 2"] [9 0 "matriz 1"] [-11 -12 "matriz 1"] [-16 10 "habitat"] [4 -12 "matriz 2"] [3 -13 "matriz 2"] [-2 2 "matriz 2"] [4 10 "matriz 2"] [7 14 "matriz 2"] [15 -4 "matriz 1"] [8 12 "matriz 2"] [6 -9 "matriz 2"] [16 5 "matriz 1"] [9 -6 "matriz 2"] [16 -3 "matriz 2"] [-16 7 "habitat"] [7 15 "matriz 2"] [-16 6 "habitat"] [9 -8 "matriz 2"] [14 -9 "habitat"] [-9 15 "habitat"] [12 -4 "matriz 2"] [11 -9 "habitat"] [-10 -5 "matriz 2"] [13 -14 "habitat"] [3 2 "matriz 2"] [-7 -2 "matriz 1"] [-5 0 "matriz 2"] [16 -4 "matriz 2"] [-16 16 "habitat"] [2 -1 "matriz 1"] [-4 -6 "matriz 1"] [-8 -3 "matriz 2"] [-5 -7 "matriz 1"] [6 4 "matriz 2"] [-4 -15 "matriz 1"] [-14 11 "habitat"] [10 -4 "matriz 2"] [15 4 "matriz 1"] [2 8 "matriz 2"] [-9 1 "habitat"] [-10 -3 "matriz 2"] [-6 -13 "matriz 1"] [-2 -5 "matriz 1"] [-7 -4 "matriz 2"] [-10 -11 "matriz 1"] [-4 3 "matriz 2"] [16 -15 "matriz 1"] [-7 5 "habitat"] [-13 -3 "matriz 1"] [-14 5 "habitat"] [-10 10 "habitat"] [-3 -1 "matriz 2"] [5 -9 "matriz 2"] [-6 11 "habitat"] [-6 13 "matriz 2"] [12 7 "matriz 2"] [1 -9 "matriz 1"] [3 -15 "matriz 1"] [-5 -9 "matriz 1"] [7 -13 "matriz 2"] [-8 11 "habitat"] [10 -10 "matriz 1"] [2 0 "matriz 2"] [5 4 "matriz 1"] [-5 3 "matriz 2"] [-1 5 "matriz 2"] [5 -10 "matriz 2"] [7 -6 "matriz 2"] [8 3 "matriz 1"] [0 0 "matriz 2"] [-13 8 "habitat"] [-4 -11 "matriz 1"] [-11 -15 "matriz 1"] [-9 14 "habitat"] [-5 8 "matriz 2"] [16 14 "matriz 2"] [10 8 "matriz 2"] [13 10 "matriz 2"] [16 4 "matriz 1"] [-2 -12 "matriz 2"] [16 11 "matriz 1"] [1 -4 "matriz 1"] [7 -1 "matriz 1"] [-10 4 "habitat"] [-12 -8 "matriz 2"] [10 -9 "matriz 2"] [13 -12 "habitat"] [15 12 "matriz 2"] [-3 5 "matriz 2"] [6 -10 "matriz 1"] [5 -2 "matriz 1"] [1 -3 "matriz 1"] [15 15 "matriz 2"] [-6 14 "habitat"] [-11 2 "matriz 2"] [13 1 "matriz 1"] [8 -2 "matriz 1"] [5 -8 "matriz 2"] [-15 7 "habitat"] [13 13 "matriz 1"] [-10 16 "habitat"] [12 -3 "matriz 2"] [5 15 "matriz 2"] [-12 -11 "matriz 2"] [-16 14 "habitat"] [-14 14 "habitat"] [-9 -4 "matriz 1"] [10 -3 "matriz 2"] [14 16 "matriz 2"] [11 12 "matriz 2"] [3 -12 "matriz 1"] [-7 -12 "matriz 1"] [8 2 "matriz 1"] [-9 -15 "matriz 1"] [13 -9 "habitat"] [10 -6 "matriz 2"] [1 -2 "matriz 1"] [-2 0 "matriz 2"] [-8 -1 "matriz 2"] [8 -5 "matriz 2"] [14 6 "matriz 1"] [-5 4 "matriz 2"] [14 12 "matriz 2"] [2 2 "matriz 2"] [-14 -7 "matriz 2"] [15 -14 "habitat"] [0 12 "matriz 1"] [2 13 "matriz 2"] [15 -13 "habitat"] [-14 7 "habitat"] [-8 13 "habitat"] [5 -3 "matriz 1"] [12 -2 "matriz 2"] [5 9 "matriz 2"] [9 7 "matriz 2"] [-1 14 "matriz 2"] [8 -8 "matriz 1"] [-14 9 "habitat"] [-12 -2 "matriz 1"] [3 6 "matriz 2"] [-14 -11 "matriz 1"] [11 -14 "habitat"] [12 2 "matriz 1"] [4 -8 "matriz 1"] [15 14 "matriz 2"] [15 -11 "habitat"] [-13 1 "matriz 2"] [6 -16 "matriz 2"] [-5 -15 "matriz 2"] [3 -9 "matriz 2"] [-6 -12 "matriz 1"] [-13 4 "habitat"] [-3 -12 "matriz 2"] [14 -3 "matriz 1"] [11 11 "matriz 2"] [-10 3 "habitat"] [-5 -3 "matriz 2"] [11 -6 "matriz 2"] [-15 8 "habitat"] [-12 5 "habitat"] [4 3 "matriz 2"] [15 -5 "matriz 2"] [9 -9 "matriz 2"] [-9 10 "habitat"] [-14 0 "matriz 2"] [3 -8 "matriz 2"] [-9 0 "matriz 2"] [12 -9 "habitat"] [5 7 "matriz 2"] [5 -4 "matriz 1"] [1 14 "matriz 2"] [-5 11 "matriz 2"] [4 7 "matriz 2"] [-2 -15 "matriz 2"] [-1 -2 "matriz 1"] [-6 -15 "matriz 2"] [8 16 "matriz 2"] [-4 11 "matriz 2"] [-7 4 "habitat"] [13 11 "matriz 2"] [-7 0 "matriz 1"] [-11 1 "habitat"] [-16 8 "habitat"] [-6 -11 "matriz 1"] [-4 -13 "matriz 2"] [2 -4 "matriz 1"] [14 3 "matriz 1"] [-2 -3 "matriz 2"] [-7 -14 "matriz 1"] [-15 14 "habitat"] [-1 -14 "matriz 1"] [10 -14 "habitat"] [9 8 "matriz 2"] [-5 -10 "matriz 1"] [-13 -15 "matriz 1"] [8 -13 "matriz 2"] [16 -7 "habitat"] [-14 -1 "matriz 1"] [6 0 "matriz 1"] [8 5 "matriz 2"] [14 -1 "matriz 1"] [-14 -6 "matriz 1"] [-11 4 "habitat"] [14 0 "matriz 2"] [6 -2 "matriz 1"] [-11 14 "habitat"] [12 9 "matriz 2"] [16 10 "matriz 2"] [-10 -1 "matriz 2"] [8 -6 "matriz 2"] [-1 3 "matriz 2"] [-1 -6 "matriz 1"] [-13 -6 "matriz 2"] [-1 16 "matriz 2"] [2 3 "matriz 1"] [2 -5 "matriz 1"] [4 12 "matriz 2"] [-1 1 "matriz 2"] [5 -12 "matriz 2"] [15 13 "matriz 1"] [-13 10 "habitat"] [-6 -6 "matriz 1"] [-8 14 "habitat"] [8 7 "matriz 2"] [-11 0 "matriz 2"] [-9 -12 "matriz 1"] [1 6 "matriz 2"] [-3 3 "matriz 1"] [-10 8 "habitat"] [7 -14 "matriz 2"] [-16 1 "habitat"] [12 0 "matriz 1"] [-4 7 "matriz 2"] [11 2 "matriz 1"] [12 5 "matriz 1"] [-15 16 "habitat"] [-15 -3 "matriz 2"] [-5 -13 "matriz 1"] [-7 3 "matriz 1"] [8 8 "matriz 2"] [-10 -14 "matriz 1"] [-12 -16 "matriz 1"] [-16 -9 "matriz 1"] [1 1 "matriz 2"] [-15 15 "habitat"] [16 -13 "habitat"] [6 5 "matriz 2"] [4 16 "matriz 2"] [2 5 "matriz 2"] [-4 -7 "matriz 1"] [0 8 "matriz 2"] [-7 -9 "matriz 1"] [11 15 "matriz 1"] [11 -8 "habitat"] [2 -7 "matriz 1"] [-11 16 "habitat"] [-10 -13 "matriz 1"] [-6 -14 "matriz 2"] [-14 -10 "matriz 2"] [4 5 "matriz 2"] [-3 -16 "matriz 2"] [0 2 "matriz 2"] [0 7 "matriz 2"] [-13 -7 "matriz 2"] [2 -11 "matriz 1"] [6 -8 "matriz 2"] [-12 3 "habitat"] [-8 -15 "matriz 1"] [2 9 "matriz 1"] [-4 -1 "matriz 2"] [-12 7 "habitat"] [-4 -8 "matriz 1"] [-5 -5 "matriz 1"] [14 -14 "habitat"] [-10 -10 "matriz 1"] [0 -13 "matriz 2"] [15 7 "matriz 1"] [4 -2 "matriz 1"] [-16 -2 "matriz 2"] [-8 -14 "matriz 1"] [-12 9 "habitat"] [3 7 "matriz 1"] [-8 0 "matriz 2"] [4 8 "matriz 2"] [11 16 "matriz 2"] [-3 -13 "matriz 2"] [3 14 "matriz 2"] [-3 12 "matriz 2"] [-13 -8 "matriz 2"] [-15 -6 "matriz 2"] [-13 7 "habitat"] [1 -14 "matriz 2"] [0 5 "matriz 2"] [-14 -12 "matriz 2"] [-11 -8 "matriz 2"] [4 -9 "matriz 2"] [9 6 "matriz 2"] [-13 9 "habitat"] [13 -3 "matriz 2"] [12 -14 "habitat"] [-1 7 "matriz 2"] [13 -15 "matriz 2"] [9 13 "matriz 2"] [-7 14 "habitat"] [-6 -1 "matriz 2"] [8 13 "matriz 2"] [15 1 "matriz 1"] [-1 -13 "matriz 2"] [-8 9 "habitat"] [8 -14 "matriz 2"] [-8 -12 "matriz 1"] [-9 -11 "matriz 1"] [9 16 "matriz 2"] [1 -15 "matriz 2"] [-15 4 "habitat"] [3 1 "matriz 2"] [11 -11 "habitat"] [9 -15 "habitat"] [0 -1 "matriz 1"] [5 -16 "matriz 2"] [2 -16 "matriz 2"] [13 -7 "habitat"] [-3 -7 "matriz 1"] [-7 -6 "matriz 1"] [-13 13 "habitat"] [-12 -12 "matriz 2"] [5 5 "matriz 2"] [12 -15 "habitat"] [13 7 "matriz 1"] [-2 3 "matriz 2"] [-13 -14 "matriz 1"] [11 -4 "matriz 2"] [8 0 "matriz 1"] [0 -4 "matriz 1"] [-15 -15 "matriz 1"] [10 -8 "matriz 1"] [-4 -4 "matriz 1"] [-8 6 "habitat"] [1 -6 "matriz 1"] [-1 0 "matriz 2"] [7 16 "matriz 2"] [-2 -4 "matriz 1"] [9 15 "matriz 2"] [15 -7 "habitat"] [-13 -5 "matriz 2"] [12 -1 "matriz 2"] [-12 -1 "matriz 2"] [0 16 "matriz 2"] [6 7 "matriz 2"] [12 -5 "matriz 2"] [-9 -16 "matriz 1"] [16 -12 "habitat"] [14 -2 "matriz 2"] [-3 4 "matriz 1"] [3 15 "matriz 2"] [10 -13 "habitat"] [-8 -16 "matriz 1"] [10 2 "matriz 1"] [-4 4 "matriz 2"] [-11 -16 "matriz 1"] [-2 -14 "matriz 1"] [7 -12 "matriz 2"] [-11 -2 "matriz 1"] [5 -15 "matriz 2"] [5 12 "matriz 1"] [15 2 "matriz 1"] [16 8 "matriz 1"] [-14 1 "matriz 2"] [0 11 "matriz 2"] [9 -3 "matriz 2"] [14 9 "matriz 2"] [4 0 "matriz 1"] [12 -6 "matriz 2"] [-15 12 "habitat"] [-5 13 "matriz 2"] [-9 12 "habitat"] [-3 2 "matriz 1"] [-13 -2 "matriz 2"] [8 10 "matriz 2"] [16 2 "matriz 1"] [-10 15 "habitat"] [1 -10 "matriz 1"] [1 13 "matriz 2"] [-7 13 "habitat"] [-1 10 "matriz 2"] [-5 6 "matriz 1"] [-5 -12 "matriz 1"] [-10 -9 "matriz 2"] [-16 -15 "matriz 1"] [2 -12 "matriz 2"] [-13 3 "habitat"] [-15 -4 "matriz 2"] [6 11 "matriz 2"] [6 -12 "matriz 2"] [-2 16 "matriz 2"] [10 6 "matriz 2"] [2 4 "matriz 2"] [-12 -7 "matriz 2"] [-4 -3 "matriz 2"] [13 5 "matriz 1"] [-16 -6 "matriz 2"] [-2 -16 "matriz 2"] [14 2 "matriz 1"] [12 12 "matriz 2"] [11 3 "matriz 1"] [10 -1 "matriz 1"] [6 10 "matriz 2"] [-7 12 "habitat"] [13 -13 "habitat"] [5 8 "matriz 2"] [-1 11 "matriz 1"] [-12 10 "habitat"] [10 7 "matriz 2"] [7 5 "matriz 2"] [7 1 "matriz 1"] [-3 -9 "matriz 1"] [-5 -11 "matriz 1"] [6 -7 "matriz 2"] [-1 -7 "matriz 1"] [-6 -16 "matriz 2"] [-7 -8 "matriz 1"] [5 3 "matriz 2"] [-9 16 "habitat"] [12 3 "matriz 1"] [-3 10 "matriz 2"] [13 -2 "matriz 1"] [-12 8 "habitat"] [-3 13 "matriz 2"] [-8 -4 "matriz 2"] [-1 -16 "matriz 2"] [8 15 "matriz 2"] [-15 -10 "matriz 2"] [15 8 "matriz 1"] [-12 6 "habitat"] [1 8 "matriz 1"] [-4 -10 "matriz 1"] [-10 -8 "matriz 2"] [1 16 "matriz 2"] [8 -1 "matriz 1"] [16 -6 "matriz 1"] [4 -1 "matriz 1"] [8 -9 "matriz 1"] [-5 12 "matriz 1"] [-15 11 "habitat"] [-10 6 "habitat"] [6 -11 "matriz 2"] [15 -15 "matriz 1"] [7 10 "matriz 2"] [13 -6 "matriz 2"] [11 -16 "habitat"] [2 -6 "matriz 1"] [16 16 "matriz 1"] [7 2 "matriz 1"] [-16 4 "habitat"] [-10 5 "habitat"] [-12 -10 "matriz 2"] [-3 0 "matriz 2"] [-3 9 "matriz 2"] [-10 0 "matriz 2"] [-5 2 "matriz 2"] [-2 -7 "matriz 1"] [16 -10 "habitat"] [10 -15 "habitat"] [-10 -6 "matriz 2"] [5 2 "matriz 2"] [-7 1 "matriz 2"] [11 8 "matriz 2"] [10 14 "matriz 2"] [-8 10 "habitat"] [-6 6 "matriz 2"] [-12 -9 "matriz 1"] [15 5 "matriz 1"] [-7 -1 "matriz 2"] [-9 13 "habitat"] [1 -12 "matriz 2"] [-16 9 "habitat"] [-16 13 "habitat"] [16 -5 "matriz 2"] [-2 -6 "matriz 1"] [4 1 "matriz 2"] [13 14 "matriz 2"] [11 -15 "habitat"] [3 -6 "matriz 1"] [-11 6 "habitat"] [-8 2 "matriz 2"] [1 15 "matriz 1"] [16 -2 "matriz 2"] [3 3 "matriz 2"] [5 11 "matriz 1"] [-11 7 "habitat"] [-11 10 "habitat"] [-9 -5 "matriz 2"] [8 -15 "habitat"] [-7 -10 "matriz 1"] [7 0 "matriz 1"] [0 -8 "matriz 1"] [11 -1 "matriz 2"] [-10 -16 "matriz 1"] [-1 -3 "matriz 1"] [-7 7 "habitat"] [2 1 "matriz 2"] [9 12 "matriz 2"] [-11 12 "habitat"] [13 6 "matriz 1"] [-10 2 "habitat"] [-6 -7 "matriz 1"] [-13 -4 "matriz 2"] [-11 13 "habitat"] [0 -15 "matriz 2"] [-16 -1 "matriz 2"] [-13 12 "habitat"] [4 -11 "matriz 2"] [-14 -14 "matriz 1"] [13 8 "matriz 1"] [3 -3 "matriz 1"] [5 14 "matriz 2"] [12 1 "matriz 1"] [-16 -12 "matriz 1"] [-9 -2 "matriz 2"] [-5 -2 "matriz 2"] [-15 6 "habitat"] [1 3 "matriz 2"] [-14 8 "habitat"] [6 9 "matriz 2"] [-15 -14 "matriz 2"] [4 4 "matriz 1"] [5 -6 "matriz 2"] [6 2 "matriz 1"] [-16 -7 "matriz 2"] [-7 -13 "matriz 1"] [10 16 "matriz 2"] [9 -4 "matriz 2"] [-2 -9 "matriz 1"] [-9 -8 "matriz 2"] [-6 -2 "matriz 2"] [14 -7 "habitat"] [-4 6 "matriz 2"] [-6 -8 "matriz 1"] [-8 -9 "matriz 1"] [16 -11 "habitat"] [4 -14 "matriz 2"] [-13 -10 "matriz 2"] [11 -7 "habitat"] [10 9 "matriz 2"] [0 -3 "matriz 1"] [-3 -15 "matriz 1"] [-12 2 "habitat"] [7 -2 "matriz 1"] [3 -2 "matriz 1"] [-7 -16 "matriz 1"] [-10 13 "habitat"] [-12 15 "habitat"] [14 -15 "matriz 1"] [0 -2 "matriz 1"] [2 10 "matriz 2"] [9 14 "matriz 1"] [13 9 "matriz 2"] [-6 5 "matriz 1"] [4 2 "matriz 2"] [3 -14 "matriz 2"] [-12 -14 "matriz 1"] [3 -11 "matriz 1"] [-2 12 "matriz 2"] [1 0 "matriz 2"] [12 14 "matriz 2"] [4 9 "matriz 2"] [8 -16 "habitat"] [-14 -9 "matriz 1"] [1 10 "matriz 2"] [15 -10 "habitat"] [-3 -14 "matriz 2"] [13 -1 "matriz 2"] [-12 -15 "matriz 1"] [7 3 "matriz 1"] [15 11 "matriz 2"] [-8 12 "habitat"] [-9 7 "habitat"] [-16 -10 "matriz 2"] [16 6 "matriz 1"] [-14 12 "habitat"] [-16 -5 "matriz 2"] [3 0 "matriz 2"] [-14 3 "habitat"] [-1 -8 "matriz 1"] [-11 -13 "matriz 1"] [-1 -1 "matriz 2"] [0 13 "matriz 2"] [13 -5 "matriz 2"] [-16 -11 "matriz 2"] [2 7 "matriz 2"] [2 -3 "matriz 1"] [-6 2 "matriz 1"] [10 12 "matriz 1"] [15 3 "matriz 1"] [-10 -15 "matriz 1"] [7 -7 "matriz 2"] [16 -1 "matriz 2"] [10 -11 "matriz 2"] [-15 -9 "matriz 2"] [7 -15 "matriz 2"] [12 -11 "habitat"] [-2 7 "matriz 2"] [-4 15 "matriz 2"] [12 -7 "habitat"] [2 12 "matriz 2"] [11 7 "matriz 2"] [7 7 "matriz 2"] [12 13 "matriz 1"] [-3 6 "matriz 2"] [-4 -14 "matriz 1"] [5 -13 "matriz 2"] [2 -9 "matriz 1"] [-16 -14 "matriz 2"] [-14 16 "habitat"] [9 9 "matriz 2"] [15 -8 "habitat"] [-9 3 "habitat"] [14 -4 "matriz 2"] [1 7 "matriz 2"] [10 -2 "matriz 2"] [-14 -15 "matriz 1"] [-6 0 "matriz 2"] [-15 -2 "matriz 2"] [-13 11 "habitat"] [-2 1 "matriz 2"] [-7 9 "habitat"] [-15 -7 "matriz 2"] [15 -12 "habitat"] [13 -8 "habitat"] [14 -16 "matriz 2"] [-6 15 "habitat"] [3 -16 "matriz 2"] [-4 0 "matriz 1"] [4 -15 "matriz 2"] [-4 8 "matriz 2"] [7 12 "matriz 2"] [-1 9 "matriz 2"] [5 -5 "matriz 1"] [10 -5 "matriz 2"] [11 6 "matriz 2"] [-13 -12 "matriz 1"] [-15 9 "habitat"] [7 11 "matriz 2"] [-15 13 "habitat"] [-8 -2 "matriz 2"] [1 11 "matriz 2"] [6 3 "matriz 2"] [6 6 "matriz 2"] [-5 1 "matriz 2"] [6 14 "matriz 2"] [8 9 "matriz 2"] [9 -1 "matriz 1"] [0 -11 "matriz 2"] [-3 -2 "matriz 2"] [4 -13 "matriz 2"] [5 -11 "matriz 2"] [-9 -13 "matriz 1"] [-11 -6 "matriz 1"] [-12 -13 "matriz 1"] [2 -8 "matriz 2"] [4 -4 "matriz 1"] [-5 7 "matriz 1"] [3 -7 "matriz 2"] [5 0 "matriz 1"] [-16 11 "habitat"] [9 5 "matriz 2"] [-6 -4 "matriz 1"] [-10 7 "habitat"] [15 -3 "matriz 2"] [0 -14 "matriz 2"] [-2 4 "matriz 2"] [-1 2 "matriz 2"] [-3 -3 "matriz 2"] [6 -5 "matriz 2"] [-6 -3 "matriz 2"] [-9 -10 "matriz 1"] [9 11 "matriz 2"] [0 -16 "matriz 2"] [-15 5 "habitat"] [-6 12 "habitat"] [-4 -5 "matriz 1"] [-2 -10 "matriz 1"] [0 -7 "matriz 1"] [-9 11 "habitat"] [-10 12 "habitat"] [-8 7 "habitat"] [-3 16 "matriz 1"] [-7 -15 "matriz 1"] [14 15 "matriz 2"] [-16 12 "habitat"] [4 6 "matriz 2"] [-3 -10 "matriz 1"] [4 -6 "matriz 1"] [1 5 "matriz 2"] [-16 -16 "matriz 1"] [-8 3 "matriz 2"] [-4 16 "matriz 2"] [-13 16 "habitat"] [1 -16 "matriz 2"] [-14 2 "habitat"] [16 1 "matriz 2"] [14 8 "matriz 1"] [-1 -4 "matriz 1"] [10 0 "matriz 1"] [6 -3 "matriz 1"] [-5 -14 "matriz 2"] [-3 15 "matriz 2"] [-5 9 "matriz 2"] [-9 6 "habitat"] [-13 5 "habitat"] [-12 1 "matriz 2"] [2 -13 "matriz 2"] [-10 -7 "matriz 2"] [3 10 "matriz 2"] [-15 1 "habitat"] [-13 -1 "matriz 2"] [-16 -4 "matriz 2"] [1 -8 "matriz 1"] [4 15 "matriz 2"] [-9 -9 "matriz 1"] [-13 0 "matriz 2"] [-4 -16 "matriz 1"] [11 14 "matriz 2"] [-4 12 "matriz 2"] [12 -12 "habitat"] [14 -8 "habitat"] [9 -13 "matriz 2"] [-4 13 "matriz 2"] [7 6 "matriz 2"] [0 4 "matriz 1"] [7 -8 "matriz 1"] [8 -11 "matriz 2"] [-9 2 "habitat"] [13 -10 "habitat"] [8 -7 "matriz 2"] [-4 -12 "matriz 1"] [12 -8 "habitat"] [12 10 "matriz 2"] [-13 6 "habitat"] [1 -11 "matriz 2"] [8 -10 "matriz 2"] [16 7 "matriz 1"] [-5 10 "matriz 1"] [10 4 "matriz 1"] [0 -9 "matriz 1"] [-1 -10 "matriz 1"] [6 13 "matriz 1"] [6 -13 "matriz 2"] [2 11 "matriz 2"] [3 12 "matriz 2"] [-6 -5 "matriz 1"] [-13 -9 "matriz 2"] [-10 -4 "matriz 2"] [6 8 "matriz 1"] [3 4 "matriz 2"] [-11 11 "habitat"] [-14 6 "habitat"] [-3 -5 "matriz 1"] [9 -16 "habitat"] [6 12 "matriz 2"] [3 -4 "matriz 1"] [-15 -11 "matriz 2"] [5 1 "matriz 1"] [-5 -4 "matriz 2"] [-12 -4 "matriz 2"] [-5 14 "matriz 2"] [5 6 "matriz 1"] [-2 -8 "matriz 1"] [12 -16 "matriz 1"] [-4 10 "matriz 2"] [-11 -10 "matriz 1"] [10 -16 "habitat"] [13 12 "matriz 2"] [2 14 "matriz 2"] [-4 1 "matriz 2"] [-15 0 "matriz 2"] [2 -10 "matriz 1"] [5 -7 "matriz 1"] [-5 -16 "matriz 1"] [-8 -8 "matriz 1"] [-7 16 "habitat"] [-11 -11 "matriz 1"] [-6 -9 "matriz 1"] [4 -16 "matriz 2"] [16 -8 "habitat"] [3 -10 "matriz 2"] [3 -1 "matriz 1"] [14 13 "matriz 2"] [-4 14 "matriz 2"] [16 -14 "habitat"] [-3 1 "matriz 1"] [-9 -6 "matriz 2"] [4 13 "matriz 2"] [-14 -5 "matriz 2"] [9 10 "matriz 1"] [16 -9 "habitat"] [1 12 "matriz 2"] [9 -5 "matriz 1"] [6 -4 "matriz 1"] [15 9 "matriz 2"] [11 -12 "habitat"] [9 1 "matriz 1"] [8 -4 "matriz 2"] [-3 -8 "matriz 1"] [10 -7 "matriz 2"] [11 -10 "habitat"] [7 4 "matriz 2"] [-2 14 "matriz 2"] [2 15 "matriz 2"] [-11 -5 "matriz 2"] [-8 -5 "matriz 2"] [-11 -1 "matriz 1"] [0 -5 "matriz 1"] [14 14 "matriz 1"] [14 1 "matriz 1"] [-15 -5 "matriz 1"] [16 9 "matriz 2"] [16 -16 "matriz 1"] [8 1 "matriz 1"] [-14 15 "habitat"] [-14 -8 "matriz 2"] [-7 11 "habitat"] [16 12 "matriz 2"] [-3 8 "matriz 2"] [-10 -12 "matriz 1"] [12 15 "matriz 2"] [5 10 "matriz 2"] [-2 -11 "matriz 2"] [7 -10 "matriz 2"] [15 10 "matriz 2"] [4 11 "matriz 2"] [-13 14 "habitat"] [16 15 "matriz 1"] [-7 -3 "matriz 2"] [-2 10 "matriz 2"] [0 -12 "matriz 2"] [-10 -2 "matriz 2"] [10 -12 "habitat"] [-16 5 "habitat"] [-15 2 "habitat"] [13 -16 "matriz 2"] [-8 -10 "matriz 1"] [13 2 "matriz 1"] [-12 -5 "matriz 1"] [-1 12 "matriz 2"] [9 3 "matriz 1"] [9 2 "matriz 1"] [1 -13 "matriz 1"] [0 1 "matriz 2"] [15 16 "matriz 2"] [3 11 "matriz 2"] [15 -1 "matriz 2"] [-2 15 "matriz 2"] [2 -15 "matriz 2"] [14 -6 "matriz 1"] [5 -1 "matriz 1"] [5 16 "matriz 2"] [-8 4 "matriz 1"] [-2 11 "matriz 2"] [-5 -1 "matriz 2"] [-12 11 "habitat"] [-3 -11 "matriz 1"] [13 3 "matriz 1"] [13 4 "matriz 1"] [10 5 "matriz 2"] [1 -7 "matriz 1"] [-11 3 "habitat"] [-6 -10 "matriz 1"] [-13 -16 "matriz 1"] [8 14 "matriz 2"] [-8 -11 "matriz 1"] [0 10 "matriz 1"] [-6 8 "habitat"] [-2 -13 "matriz 1"] [14 4 "matriz 1"] [-12 13 "habitat"] [15 -16 "matriz 1"] [-1 13 "matriz 1"] [-14 10 "habitat"] [-12 12 "habitat"] [13 -11 "habitat"] [-7 6 "matriz 2"] [0 3 "matriz 2"] [6 -6 "matriz 2"] [7 -9 "matriz 2"] [-1 -5 "matriz 1"] [-14 -13 "matriz 2"] [15 -2 "matriz 2"] [-12 -3 "matriz 2"] [1 9 "matriz 2"] [-15 -1 "matriz 2"] [12 -13 "habitat"] [12 -10 "habitat"] [-1 -12 "matriz 2"] [-1 4 "matriz 2"] [14 5 "matriz 1"] [0 6 "matriz 2"] [-8 16 "habitat"] [-13 -13 "matriz 2"] [-4 9 "matriz 2"] [5 13 "matriz 2"] [-12 16 "habitat"] [14 7 "matriz 1"] [-9 5 "habitat"] [1 2 "matriz 2"] [16 13 "matriz 1"] [14 10 "matriz 1"] [12 11 "matriz 2"] [-7 2 "matriz 2"] [13 0 "matriz 2"] [-6 4 "matriz 1"] [6 -14 "matriz 1"] [-6 1 "matriz 2"] [-15 -8 "matriz 2"] [15 6 "matriz 1"] [-1 -9 "matriz 1"] [-16 -3 "matriz 2"] [8 -12 "matriz 2"] [9 4 "matriz 1"] [-10 14 "habitat"] [-4 -2 "matriz 2"] [6 -1 "matriz 1"] [3 5 "matriz 1"] [-9 -3 "matriz 2"] [12 4 "matriz 1"] [10 1 "matriz 1"] [-5 5 "matriz 2"] [-11 8 "habitat"] [-12 -6 "matriz 2"] [3 8 "matriz 1"] [10 3 "matriz 1"] [-4 2 "matriz 2"] [-6 7 "matriz 2"] [10 13 "matriz 2"] [3 9 "matriz 1"] [11 -5 "matriz 2"] [4 -7 "matriz 1"] [9 -2 "matriz 2"] [8 4 "matriz 2"] [-4 -9 "matriz 1"] [-2 -2 "matriz 2"] [-6 10 "matriz 2"] [15 0 "matriz 2"] [-14 -2 "matriz 1"] [1 -5 "matriz 1"] [-2 9 "matriz 2"] [-1 6 "matriz 2"] [-15 -13 "matriz 1"] [4 14 "matriz 1"] [-11 -3 "matriz 2"] [-8 -6 "matriz 2"] [-15 3 "habitat"] [-1 -11 "matriz 1"] [-16 3 "habitat"] [-8 8 "habitat"] [2 -14 "matriz 2"] [-7 15 "habitat"] [10 15 "matriz 2"] [-3 -6 "matriz 1"] [-4 5 "matriz 2"] [13 -4 "matriz 2"] [11 9 "matriz 2"] [14 -10 "habitat"] [7 -11 "matriz 2"] [7 13 "matriz 2"]]

 foreach data [ [c] ->

    ask patch (item 0 c) (item 1 c) [set tipo (item 2 c) ]]
end 

to formar_vecindades
  let num_grupo 1
  ask patches [
    if (tipo = "habitat" and grupo = -999) [
      formar_vecindad num_grupo
      set num_grupo num_grupo + 1
    ]
  ]
  set num_grupos max [grupo] of patches
end 

to formar_vecindad [ num_grupo ]
  set grupo num_grupo
  ask neighbors4 with [ tipo = "habitat" and grupo = -999 ][
    formar_vecindad num_grupo
  ]
end 

to colorear_cuadros
  if COLOREAR_POR = "tipo" [
    ask patches [
      if tipo = "habitat"  [ set pcolor green ]
      if tipo = "matriz 1" [ set pcolor brown + 1 ]
      if tipo = "matriz 2" [ set pcolor brown  ]
    ]
  ]
  if COLOREAR_POR = "presencia individuos en parche" [
    ask patches [
      if tipo != "habitat" [ set pcolor black + 1 ]
      if tipo = "habitat" [
        ifelse any? turtles with [ grupo = [grupo] of myself] [ set pcolor gray + 2][ set pcolor black ]
      ]
    ]
  ]
  if COLOREAR_POR = "densidad individuos en parche" [
    ask patches [
      if tipo != "habitat" [ set pcolor black + 1 ]
      if tipo = "habitat" [
        set pcolor scale-color magenta (count turtles with [ grupo = [grupo] of myself])  0 (count turtles)
      ]
    ]
  ]
  ask patches [
    let escala 1
    if estado = "quemandose" [ set pcolor red + (tiempo_desde_quemado * escala / TIEMPO_QUEMADO) ]
  ]
end 

to inicializar_poblacion
  create-turtles TAMANIO_POBLACION_INICIAL [
    move-to one-of patches
    inicialiar_agente
  ]
end 

to inicialiar_agente
  set color red
  set vivo? true
end 

to-report fade [t]
  report t * t * t * ( t * ( t * 6 - 15 ) + 10 )
end 

to-report lerp [ a b t ]
  report a + t * (b - a)
end 

to-report grad [ hash x y]
  let vectores  [[1 0] [-1 0] [0 1] [0 -1]]
  let vec item (hash mod 4) vectores
  report (item 0 vec) * x + (item 1 vec) * y
end 

to-report perlin_noise [ x y ]
  let xi (floor x) mod 255
  let yi (floor y) mod 255
  let xf x - floor x
  let yf y - floor y

  let u fade xf
  let v fade yf

  let aa (grad (item (((item ((xi    ) mod 255) p) + yi    ) mod 255) p)  xf       yf     )
  let ab (grad (item (((item ((xi    ) mod 255) p) + yi + 1) mod 255) p)  xf      (yf - 1))
  let ba (grad (item (((item ((xi + 1) mod 255) p) + yi    ) mod 255) p) (xf - 1)  yf     )
  let bb (grad (item (((item ((xi + 1) mod 255) p) + yi + 1) mod 255) p) (xf - 1) (yf - 1))

  let x1 lerp aa ba u
  let x2 lerp ab bb u

  report lerp x1 x2 v
end 

to incendiar
  if mouse-down? [ iniciar_incendio mouse-xcor mouse-ycor ]
end 

to iniciar_incendio [ x y ]
  ask patch x y [
    set estado "quemandose"
    set tiempo_desde_quemado 0
    ask turtles-here [ die ]
  ]
end 

to extender_incendio
  ask patches with [ estado = "quemandose" and tiempo_desde_quemado = 1][
    ask neighbors4 with [ estado = "NA" and tipo = [tipo] of myself ] [
      if random-float 1.0 < PROB_QUEMARSE [
        set estado "quemandose"
        set tiempo_desde_quemado 0
      ]
    ]
  ]
  ask patches with [ estado = "quemandose" ][
    set tiempo_desde_quemado tiempo_desde_quemado + 1
    if tiempo_desde_quemado > TIEMPO_QUEMADO [ set estado "NA" set tiempo_desde_quemado -999 ]
  ]
end 

to EDITAR_PAISAJE_MANUALMENTE
  if mouse-down?
  [ ask patch mouse-xcor mouse-ycor [ set tipo TIPO_A_AGREGAR ]
    ask patches [ set grupo -999 ]

    colorear_cuadros
    display
    formar_vecindades
  ]
end 

to INICIALIZAR_GRAFICAS_EXPLORACION
  let nombres_graficas [ "exploracion poblacion" "exploracion parches ocupados" ]
   foreach nombres_graficas [
    nombre ->
    set-current-plot nombre
    clear-plot
  ]
end 

to inicializar_graficas
  let nombres_graficas [ "tamaño población" "fracción parches hábitat ocupados" ]
  foreach nombres_graficas [
    nombre ->
    set-current-plot nombre
    clear-plot
  ]
end 

to EXPLORAR_AUTOMATICAMENTE
  set PROP_MATRIZ_2 1.0
  repeat (round 1 / intervalo_explorar) [
    repeat num_repeticiones [
      setup
      repeat tiempo_final_exploracion [ go ]
    ]
    set PROP_MATRIZ_2 PROP_MATRIZ_2 - intervalo_explorar
  ]
end 

to-report fraccion_parches_ocupados
  ifelse any? turtles with [grupo = -999]
  [ report ((length remove-duplicates [grupo] of turtles) - 1) / num_grupos ]
  [ report ((length remove-duplicates [grupo] of turtles)) / num_grupos ]
end