anchor.diff:       &                 - BCtmp(:, 19) * BCtmp(:, 22)
             BCtmp(:, 20) = BCtmp(:, 23) * BCtmp(:, 20)
        &                 - BCtmp(:, 19) * BCtmp(:, 24)
  -          BC(:,elas3)  = BCtmp(:, 20) / BCtmp(:, 17)
  +          BC(:,elas3)  = BCtmp(:, 20) / BCtmp(:, 17) * Delt(1)
             BC(:,elas4)  = BCtmp(:, 18) / BCtmp(:, 17)
   c
             BCtmp(:, 22) = BCtmp(:, 17) * BCtmp(:, 22)
  
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: c
          where (ibits(iBC,14,3) .eq. 2)
            tmp     = BCtmp(:,17)**2 + BCtmp(:,18)**2 + BCtmp(:,19)**2
            BC(:,elas3) = tmp * BCtmp(:,20) / BCtmp(:,18) * Delt(1)
            BC(:,elas4) =       BCtmp(:,17) / BCtmp(:,18)
            BC(:,elas5) =       BCtmp(:,19) / BCtmp(:,18)
          endwhere
  