anchor.diff:        if (iALE .eq. 2) then 
   c
   c----------------Re-organize Mesh Elastic Velocities-----------------
  +c.... convert mesh velocity to mesh displacement
   c
           elas3 = ndof + 2
           elas4 = ndof + 3
  
positive.code:      &                 - BCtmp(:, 19) * BCtmp(:, 22)
            BCtmp(:, 20) = BCtmp(:, 23) * BCtmp(:, 20)
       &                 - BCtmp(:, 19) * BCtmp(:, 23)
            BC(:,elas3)  = BCtmp(:, 20) / BCtmp(:, 18)
            BC(:,elas4)  = BCtmp(:, 17) / BCtmp(:, 18)
  c
            BCtmp(:, 21) = BCtmp(:, 18) * BCtmp(:, 21)
  
negative.code:      &                 - BCtmp(:, 18) * BCtmp(:, 23)
            BCtmp(:, 20) = BCtmp(:, 22) * BCtmp(:, 20)
       &                 - BCtmp(:, 18) * BCtmp(:, 23)
            BC(:,elas3)  = BCtmp(:, 20) / BCtmp(:, 17) * Delt(1)
            BC(:,elas4)  = BCtmp(:, 19) / BCtmp(:, 17)
  c
            BCtmp(:, 22) = BCtmp(:, 17) * BCtmp(:, 22) 
  