anchor.diff:  c.... original mov+rot+shrk
   c
         if ( casenumber .eq. 6 ) then
  -        if (lstep .eq. 0) then
  -c        dyn_org    = lstep * 0.04
  -        dyn_org    = 2.625e-3
  +        xtsl     = 1e-4
  +        dyn_org  = 2.625e-3 + lstep * xtsl
  +        rotf     = 180.0
           do i = 1,numnp
             if ( (ibits(iBC(i),14,3) .eq. 7) .and.
        &         (x(i,1) .lt. 6.65e-3) .and. (x(i,1) .gt. -1.4e-3) .and.
        &         ((x(i,2)*x(i,2) + x(i,3)*x(i,3)) .lt. 2.89e-6) ) then
               if( x(i,2) .ge. 0.6e-3 ) then ! top
                 if ( x(i,1) .le. (dyn_org-1.025e-3) ) then ! top tail
  -                disp(i,1) = -0.05 * (x(i,1)-0.5e-3) !+ 0.02
  -     &                    + (x(i,1)-0.5e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,2)-1.125e-3) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)-1.125e-3) !+ 0.01
  -     &                    + (x(i,2)-1.125e-3) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,1)-0.5e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-0.5e-3) + xtsl
  +     &                    + (x(i,1)-0.5e-3) * (cos(pi/rotf) - 1.0)
  +     &                    - (x(i,2)-1.125e-3) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)-1.125e-3) !+ 0.01
  +     &                    + (x(i,2)-1.125e-3) * (cos(pi/rotf) - 1.0)
  +     &                    + (x(i,1)-0.5e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 else if ( x(i,1) .ge. (dyn_org+1.025e-3) ) then ! top head
  -                disp(i,1) = -0.05 * (x(i,1)-4.75e-3) !+ 0.02
  -     &                    + (x(i,1)-4.75e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,2)-1.125e-3) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)-1.125e-3) !+ 0.01
  -     &                    + (x(i,2)-1.125e-3) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,1)-4.75e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-4.75e-3) + xtsl
  +     &                    + (x(i,1)-4.75e-3) * (cos(pi/rotf) - 1.0)
  +     &                    - (x(i,2)-1.125e-3) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)-1.125e-3) !+ 0.01
  +     &                    + (x(i,2)-1.125e-3) * (cos(pi/rotf) - 1.0)
  +     &                    + (x(i,1)-4.75e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 else ! top middle
  -                disp(i,1) = -0.05 * (x(i,1)-2.625e-3) !+ 0.02
  -     &                    + (x(i,1)-2.625e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,2)-1.125e-3) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)-1.125e-3) !+ 0.01
  -     &                    + (x(i,2)-1.125e-3) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,1)-2.625e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-2.625e-3) + xtsl
  +     &                    + (x(i,1)-2.625e-3) * (cos(pi/rotf) - 1.0)
  +     &                    - (x(i,2)-1.125e-3) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)-1.125e-3) !+ 0.01
  +     &                    + (x(i,2)-1.125e-3) * (cos(pi/rotf) - 1.0)
  +     &                    + (x(i,1)-2.625e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 endif ! end top--switch head middle tail
   
               else if( x(i,2) .le. -0.6e-3 ) then ! bottom
                 if ( x(i,1) .le. (dyn_org-1.025e-3) ) then ! bottom tail
  -                disp(i,1) = -0.05 * (x(i,1)-0.5e-3) !+ 0.02
  -     &                    - (x(i,1)-0.5e-3) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,2)+1.125e-3) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)+1.125e-3) !- 0.01
  -     &                    - (x(i,2)+1.125e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,1)-0.5e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-0.5e-3) + xtsl
  +     &                    - (x(i,1)-0.5e-3) * (cos(pi/rotf) - 1.0)
  +     &                    + (x(i,2)+1.125e-3) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)+1.125e-3) !- 0.01
  +     &                    - (x(i,2)+1.125e-3) * (cos(pi/rotf) - 1.0)
  +     &                    - (x(i,1)-0.5e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 else if ( x(i,1) .ge. (dyn_org+1.025e-3) ) then ! bottom head
  -                disp(i,1) = -0.05 * (x(i,1)-4.75e-3) !+ 0.02
  -     &                    - (x(i,1)-4.75e-3) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,2)+1.125e-3) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)+1.125e-3) !- 0.01
  -     &                    - (x(i,2)+1.125e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,1)-4.75e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-4.75e-3) + xtsl
  +     &                    - (x(i,1)-4.75e-3) * (cos(pi/rotf) - 1.0)
  +     &                    + (x(i,2)+1.125e-3) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)+1.125e-3) !- 0.01
  +     &                    - (x(i,2)+1.125e-3) * (cos(pi/rotf) - 1.0)
  +     &                    - (x(i,1)-4.75e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 else ! bottom middle
  -                disp(i,1) = -0.05 * (x(i,1)-2.625e-3) !+ 0.02
  -     &                    - (x(i,1)-2.625e-3) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,2)+1.125e-3) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)+1.125e-3) !- 0.01
  -     &                    - (x(i,2)+1.125e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,1)-2.625e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-2.625e-3) + xtsl
  +     &                    - (x(i,1)-2.625e-3) * (cos(pi/rotf) - 1.0)
  +     &                    + (x(i,2)+1.125e-3) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)+1.125e-3) !- 0.01
  +     &                    - (x(i,2)+1.125e-3) * (cos(pi/rotf) - 1.0)
  +     &                    - (x(i,1)-2.625e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 endif ! end bottom--switch head middle tail
   
               else ! middle
                 if ( x(i,1) .le. (dyn_org-1.025e-3) ) then ! middle tail
  -                disp(i,1) = -0.05 * (x(i,1)-0.5e-3) !+ 0.02
  -     &                    + (x(i,1)-0.5e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,2)) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)) ! + 0.01
  -     &                    + (x(i,2)) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,1)-0.5e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-0.5e-3) + xtsl
  +c     &                    + (x(i,1)-0.5e-3) * (cos(pi/rotf) - 1.0)
  +c     &                    - (x(i,2)) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)) ! + 0.01
  +c     &                    + (x(i,2)) * (cos(pi/rotf) - 1.0)
  +c     &                    + (x(i,1)-0.5e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 else if ( x(i,1) .ge. (dyn_org+1.025e-3) ) then ! middle head
  -                disp(i,1) = -0.05 * (x(i,1)-4.75e-3) !+ 0.02
  -     &                    + (x(i,1)-4.75e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,2)) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)) ! + 0.01
  -     &                    + (x(i,2)) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,1)-4.75e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-4.75e-3) + xtsl
  +c     &                    + (x(i,1)-4.75e-3) * (cos(pi/rotf) - 1.0)
  +c     &                    - (x(i,2)) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)) ! + 0.01
  +c     &                    + (x(i,2)) * (cos(pi/rotf) - 1.0)
  +c     &                    + (x(i,1)-4.75e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 else ! middle middle
  -                disp(i,1) = -0.05 * (x(i,1)-2.625e-3) !+ 0.02
  -     &                    + (x(i,1)-2.625e-3) * (cos(pi/600) - 1.0)
  -     &                    - (x(i,2)) *  sin(pi/600)
  -                disp(i,2) = -0.05 * (x(i,2)) ! + 0.01
  -     &                    + (x(i,2)) * (cos(pi/600) - 1.0)
  -     &                    + (x(i,1)-2.625e-3) *  sin(pi/600)
  +                disp(i,1) = -0.05 * (x(i,1)-2.625e-3) + xtsl
  +c     &                    + (x(i,1)-2.625e-3) * (cos(pi/rotf) - 1.0)
  +c     &                    - (x(i,2)) *  sin(pi/rotf)
  +                disp(i,2) = -0.05 * (x(i,2)) ! + 0.01
  +c     &                    + (x(i,2)) * (cos(pi/rotf) - 1.0)
  +c     &                    + (x(i,1)-2.625e-3) *  sin(pi/rotf)
                   disp(i,3) = -0.05 * x(i,3)
   
                 endif ! end middle--switch head middle tail
  
positive.code:             BC(i,1:3)   = disp(i,1:3)
            endif ! end if inside cylinder
          enddo ! end loop numnp
          endif ! first time step
        endif ! end if case 6
  c
  c.... end test case 6
  
negative.code:       real*8    acc
        real*8    totalForce(3),    objMass
        dimension Forin(3),         Forout(3)
  c
  c.... dynamic origin, x translation, rotation frequence
        real*8    dyn_org,   xtsl,  rotf
  c
        if (elasFDC .gt. 0) then
          casenumber = elasFDC
  