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) DIMENSION ULOCAL(*) DIMENSION JGVBLOCK(*),JMATYP(*) DIMENSION ALOCAL(NDIM,*) PARAMETER (NELEMMAX=100) DIMENSION JELEMLIST(NELEMMAX),JELEMTYPE(NELEMMAX) PARAMETER (CHARLENGTH = 5.0D-3,ELINC = 0.1D0) DIMENSION TIME(2) C C The dimensions of the variables ARRAY and JARRAY C must be set equal to or greater than 15 C LOCNUM = 0 JRCD = 0 PARTNAME = ' ' PEEQ = 0.0D0 JTYP = 0 CALL GETPARTINFO(NODE,0,PARTNAME,LOCNUM,JRCD) NELEMS = NELEMMAX CALL GETNODETOELEMCONN(NODE,NELEMS, $ JELEMLIST,JELEMTYPE,JRCD,JGVBLOCK) CALL GETVRMAVGATNODE(NODE,JTYP,'PE',ARRAY,JRCD, $ JELEMLIST,NELEMS,JMATYP,JGVBLOCK) PEEQ = ARRAY(7) IF (PEEQ.GE.0.02D0) THEN ULOCAL(NDIM) = ULOCAL(NDIM) - PEEQ * UREF END IF RETURN END