C     USER INPUT FOR ADAPTIVE MESH CONSTRAINT
C
      SUBROUTINE UMESHMOTION(UREF,ULOCAL,NODE,NNDOF,
     $     LNODETYPE,ALOCAL,NDIM,TIME,DTIME,PNEWDT,
     $     KSTEP,KINC,KMESHSWEEP,JMATYP,JGVBLOCK,LSMOOTH)
C
      include 'ABA_PARAM.INC'
C
C     USER DEFINED DIMENSION STATEMENTS
C
      CHARACTER*80 PARTNAME
      DIMENSION ARRAY(1000),JPOS(15)
      DIMENSION ULOCAL(*)
      DIMENSION JGVBLOCK(*),JMATYP(*)
      DIMENSION ALOCAL(NDIM,*)
      Parameter (NELEMMAX=100)
      DIMENSION JELEMLIST(NELEMMAX),JELEMTYPE(NELEMMAX)
      DIMENSION TIME(2)
C
C     The dimensions of the variables ARRAY and JARRAY
C     must be set equal to or greater than 15
C
      NELEMS = NELEMMAX
      JTYP = 0
      CALL GETNODETOELEMCONN(NODE,NELEMS,
     $     JELEMLIST,JELEMTYPE,JRCD,JGVBLOCK)
      CALL GETVRMAVGATNODE(NODE,JTYP,'FLVEL',ARRAY,JRCD,
     $     JELEMLIST,NELEMS,JMATYP,JGVBLOCK)
      FVMAG = SQRT(ARRAY(1)**2+ARRAY(2)**2)
C
C     Ablation proportional to the local fluid velocity
      ULOCAL(NDIM) = ULOCAL(NDIM) - FVMAG * UREF
C
      RETURN
      END