46 integer :: nxf,nyf,nzf,n,ix,jy,kz,ixs,nzfmax,nmax
47 real :: field(0:nxmax-1,0:nymax-1,nzfmax,nmax),xshiftaux(0:nxmax-1)
58 if (nxshift.ne.0)
then
60 if (ix.ge.nxshift)
then
65 xshiftaux(ixs)=field(ix,jy,kz,n)
68 field(ix,jy,kz,n)=xshiftaux(ix)
75 field(nxf,jy,kz,n)=field(0,jy,kz,n)
subroutine shift_field(field, nxf, nyf, nzfmax, nzf, nmax, n)