! Retorna, donada una matriu d'spins S d'un sistema de mida LxL,
! la magnetització del sistema.
real*8 function magne(S, L)
  implicit none
  integer*4, intent(in) :: L
  integer*2, intent(out) :: S(1:128,1:128)
  integer*4 :: i, j
  real*8 :: mag
  mag = 0d0

  do i = 1, L
    do j = 1, L
      mag = mag + S(i, j)
    enddo
  enddo

  magne = mag
  return
end function magne