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