SUBROUTINE MPC(UE,A,JDOF,MDOF,N,JTYPE,X,U,UINIT,MAXDOF, + LMPC,KSTEP,KINC,TIME,NT,NF,TEMP,FIELD,LTRAN,TRAN) C INCLUDE 'ABA_PARAM.INC' C DIMENSION A(N),JDOF(N),X(6,N),U(MAXDOF,N),UINIT(MAXDOF,N), + TIME(2),TEMP(NT,N),FIELD(NF,NT,N),LTRAN(N),TRAN(3,3,N) C LMPC=1 IF (KSTEP.GE.5) THEN LMPC=0 end if C IF (JTYPE.EQ.1) THEN C C---------- Set X displacement of node a = node b JDOF(1) = 1 JDOF(2) = 1 C A(1)=1. A(2)=-1. UE=U(1,2) END IF C IF (JTYPE.EQ.2) THEN C C---------- Set Y displacement of node a = node b JDOF(1) = 2 JDOF(2) = 2 C A(1)= 1. A(2)=-1. UE=U(2,2) END IF C RETURN END