subroutine mv4psas(NDIM,a,XX,AX,izs,rzs,dzs),6

      implicit none

      include 'champ.cdk'
      include 'chobs.cdk'

      integer  NDIM, i, izs(1),IR
      integer  NTRAJ1
      real*8   AX(NDIM), XX(NDIM)
      real     dl_v(NDIM)
      real*8   rzs(1), dzs(1),a(1)
      real     ZRR(N,N),ZRRR(N,N)
      real*8   J
      COMPLEX  ZW(IKTX,IKTY,52),SW(IKTX,IKTY,52)
       
       GRFLAG = 0
       OBFLAG = 0
       JXFLAG = 0
        
       NTRAJ1 = 52 

       DO IR =1,NDIM
        dl_v(IR) = XX(IR)/sigmao(IR)
       ENDDO

       CALL HT_4d(dl_v,ZW,SW)

       DO NT = 1,NTRAJ1
       DO IKX=1,IKTX
       DO IKY=1,IKTY
        ZTRA2(IKX,IKY,NT) = ZW(IKX,IKY,NT)
       ENDDO
       ENDDO
       ENDDO

       CALL ADJINC
       CALL SQRTBT(ZO,ZR)
       CALL SQRTB(ZR,ZO)
       CALL LININC(J)

       CALL H_4d(ZTRA2,STRA2,dl_v)
       

       DO IR = 1,NDIM
           AX(IR) =dl_v(IR)/sigmao(IR) + XX(IR)
       ENDDO

      end subroutine