program p1e3
  implicit none
  integer*4, PARAMETER :: L = 64
  integer*2 :: S(1:L, 1:L)
  real*8 :: MAGNE

  call WRITECONFIG(S, L)
  print *, "Magnetització:"
  print *, MAGNE(S, L)
endprogram p1e3