quad10/fenomens/lab/p4/pbc.f90

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