REAL u(0:nc, 0:ny), fact, avg
u = fact * (u - avg)

REAL r(0:nx, 0:ny), u(0:nx, 0:ny), f(0:nx, 0:ny)
r(0:nx:jmp, 0:ny:jmp)=-u(0:nx:jmp, 0:ny:jmp) + &
牋0.25*(CSHIFT(u(0:nx:jmp,0:ny:jmp),1,1) + &
牋CSHIFT(u(0:nx:jmp,0:ny:jmp), -1,1) + &
牋CSHIFT(u(0:nx:jmp,0:ny:jmp), 1,2) + &
牋CSHIFT(u(0:nx:jmp,0:ny:jmp), -1,2) + &
f(0:nx:jmp, 0:ny:jmp)

REAL r(:,:), u(:,:), f(:,:)
r = - u+0.25*(CSHIFT(u,1,1)+CSHIFT(U,-1,1)+&
燙SHIFT(u,1,2)+CSHIFT(u,-1,2)) + f