anchor.diff:  
           if (nsd .eq. 3) then
           where (ibits(iBC,14,3) .eq. 7)
  -          BC(:,elas3) = BCtmp(:,20) * BCtmp(:,17)
  -          BC(:,elas4) = BCtmp(:,20) * BCtmp(:,18)
  -          BC(:,elas5) = BCtmp(:,20) * BCtmp(:,19)
  +          BC(:,elas3) = BCtmp(:,20) * BCtmp(:,17) * Delt(1)
  +          BC(:,elas4) = BCtmp(:,20) * BCtmp(:,18) * Delt(1)
  +          BC(:,elas5) = BCtmp(:,20) * BCtmp(:,19) * Delt(1)
           endwhere
           endif
         endif ! end if iALE = 2
  
positive.code:           BCtmp(:, 23) = BCtmp(:, 17) * BCtmp(:, 23)
            BCtmp(:, 24) = BCtmp(:, 17) * BCtmp(:, 24)
       &                 - BCtmp(:, 21) * BCtmp(:, 20)
            BC(:,elas5)  = BCtmp(:, 24) / BCtmp(:, 23)
            BC(:,elas6)  = BCtmp(:, 22) / BCtmp(:, 23)
          endwhere
          endif
  
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) 
  