REAL, ALLOCATEBLE :: u(:,:)
CALL RANDOM_SEED
CALL RANDOM_NUMBER(u)

REAL r(0:nx, 0:ny)
error_value = SQRT (SUM(r*r))

REAL u(0:nx, 0:ny)
u(0:nx:jmp2, jmp:ny:jmp2) = &
u(0:nx:jmp2, jmp:ny:jmp2) + &
0.5*CSHIFT(u(0:nx:jmp2, jmp2:ny:jmp2),1,2)