Sznajd Model

; Sznajd Model
; version 1.0
; 21/10/14
;
; First attempt in NetLogo, by Cicero Juliao

patches-own
[
  szpin  ; state A = +1,  state B = -1
]


globals
[
  magnetization  ; global magnetization of the lattice
  pa timestep sizeLattice
  
]

to setup
  ca
  ask patches
  [ifelse random-float 1.00 <= 0.50 [ set szpin 1] [set szpin -1]  ; random initial configuration
  ]
  
  ask patches
  [
  ifelse szpin = 1 [set magnetization magnetization + 1]  [set magnetization magnetization - 1]
  
  ifelse szpin = 1 [set pcolor blue] [set pcolor red]
  ]
  
  set sizeLattice max-pxcor + max-pxcor + 1
end 

to go
  sznajdModel     ; it implement the model
  calculateMagn   ; it calculates global magnetization
  doGraphics      ; it plot the graphs
  set timestep timestep + 1
end 

to calculateMagn
  set magnetization 0
  ask patches
  [ifelse szpin = 1 [ set magnetization magnetization + 1]
    [set magnetization magnetization - 1]
  ]
  set magnetization magnetization / sizeLattice
end 

to doGraphics
  set-current-plot "Magnetization"
  plotxy timestep magnetization
end 

to sznajdModel
  set pa random-pxcor          ; coordinate x
  ask patch (pa + 2) 0 [set szpin ([szpin] of (patch pa 0))]
  ask patch (pa - 1) 0 [set szpin ([szpin] of (patch (pa + 1) 0))]

  ask patches
  [
    ifelse szpin = 1 [set pcolor blue] [set pcolor red]
  ]
end