! Funció que imposa les condicions de contorn periòdiques, enviant l'índex N al
! valor corresponent en l'interval [1, L].
!
! Prerequisit: N està a l'interval [0, L + 1].
integer*4 function PBC(N, L)
  integer*4, intent(in) :: N, L
  if (N == 0) then
    PBC = L
  else if (N == L + 1) then
    PBC = 1
  else
    PBC = N
  endif
  return
endfunction