anchor.diff:              IVDOF = VTL2G(K + NEQ*(J-1) + NEQ*NNDE*(M-1))
               VT1W(K,J)=VT1E(IVDOF)*DIR
               VT2W(K,J)=VT2E(IVDOF)*DIR
  +C MWF DEBUG
  +            WRITE(6,*)'EL_VEL_PREP IDVE.EQ.2 M= ',M,' J= ',J,
  +     &           ' IEM= ',IEM,' XW(',K,',',J,')= ',XW(K,J),
  +     &           ' IVDOF= ',IVDOF,' VT1W(',K,',',J,')= ',VT1W(K,J),
  +     &           ' VT2W(',K,',',J,')= ',VT2W(K,J)
  +C MWF DEBUG            
             ENDDO
           ENDDO
         ELSE IF (IDVE.EQ.3) THEN
  
positive.code: C ===== FOR THE CASE THAT MPT(IPT)=0, I.E., THE VERY FIRST
  C     TRACKING: LOOP OVER ALL CONNECTED ELEMENTS
  C
  C MWF DEBUG
  C        WRITE(LU_O,*)' ENTERING NODE TRACKING STEP IPT= ', IPT
          DO 150 I=NLRL(NP)+1,NLRL(NP+1)
            M=LRL(I)
  C ELENOD IN GENERAL SELECTS THE NUMBER OF NODES PER ELEMENT USING IE
  
negative.code:       DOUBLE PRECISION DEQ,TS,DTS,T,SDT,DT0
        DOUBLE PRECISION DIR
  C     
        INTEGER I,IPT,NODE,NP,M,K
        INTEGER IDSDT,I1,I2,I3,M2,M3,N1,N2,N3,J1,J2,J3
  C
  C
  