program p2e1
  implicit none
  integer*4, PARAMETER :: L = 32
  integer*2 :: S(1:L, 1:L)
  real*8 :: MAGNE, ENERG
  integer*4 :: PBC(0:L+1), I

  PBC(0) = L
  PBC(L + 1) = 1
  do I = 1, L
    PBC(I) = I
  enddo

  call WRITECONFIG(S, L)
  print *, "Magnetització:"
  print *, MAGNE(S, L)
  print *, ""
  print *, "Energia:"
  print *, ENERG(S, L, PBC)
endprogram p2e1