| ! 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 |