PARAMETER(nxblock=nx/nxp, nyblock=ny/nyp, nxlocal=nxblock+1, nylocal=nyblock+1) |
REAL u(0:nxlocal,0:nylocal), |
unew(0:nxlocal,0:nylocal), |
f(0:nxlocal,0:nylocal) |
REAL coord(1:2), dims(1:2), periods(1,2) |
dims(1) = nxp; dims(2) = nyp; |
periods(1) = .false.; |
periods(2) = .false. |
reorder = .true. |
ndim = 2 |