[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode
C        USER INPUT FOR ADAPTIVE MESH CONSTRAINT
           SUBROUTINE UMESHMOTION(UREF,ULOCAL,NODE,NNDOF,LNODETYPE,ALOCAL,NDIM,TIME,DTIME,PNEWDT,KSTEP,KINC,KMESHSWEEP,JMATYP,JGVBLOCK,LSMOOTH)
          INCLUDE 'ABA_PARAM.INC'
C      USER DEFINED DIMENSION STATEMENTS
        DIMENSION ULOCAL(*)   
        DIMENSION ALOCAL(NDIM,*)
        DIMENSION JGVBLOCK(*),JMATYP(*)
        DIMENSION  TIME(2)        
        DIMENSION  ARRAY(15)
        DIMENSION  JELEMLIST(10),JELEMTYPE(10)
        double precision k,PEEQ
C      
        NELEMS=10
        k=1.0
C      
       CALL GETNODETOELEMCONN(NODE,NELEMS,JELEMLIST,JELEMTYPE,JRCD,JGVBLOCK)
C     
       CALL GETVRMAVGATNODE(NODE,JTYP,'PE',ARRAY,JRCD,JELEMLIST,NELEMS,JMATYP,JGVBLOCK)
C     
       PEEQ=ARRAY(7)
C     
       if(PEEQ.GT.0.0) then
          ULOCAL(NDIM)=-k*PEEQ
       end if
       RETURN
       END