subroutine xso_Neumann use modglobal, only : ie, ihc, rk3step, dt, dxi, nsv use modfields, only :sv0, svm real rk3coef integer n, m rk3coef = dt/(4.-dble(rk3step)) do n = 1, nsv do m = 1, ihc sv0(ie + m, :, :, n) = sv0(ie, :, :, n) svm(ie + m, :, :, n) = svm(ie, :, :, n) end do end do end subroutine xso_Neumann