anchor.diff:  
      SUBROUTINE CleanUp_Discrt()
         ! deallocate temp matrices
  -      IF (ALLOCATED(TempPropsB))  DEALLOCATE(TempPropsB)
  +      IF (ALLOCATED(TempProps))   DEALLOCATE(TempProps)
         IF (ALLOCATED(TempMembers)) DEALLOCATE(TempMembers)
      END SUBROUTINE CleanUp_Discrt
   
  
positive.code:                   ! create a new property set 
                    ! k, E, G, rho, d, t, Init                
                    kprop = kprop + 1
                    CALL SetNewProp(kprop, TempPropsB(Prop1, 2), TempPropsB(Prop1, 3), Init%PropSetsB(Prop1, 4), d1 + J*dd, t1 + J*dt,  TempPropsB)           
                    kelem = kelem + 1
                    CALL SetNewElem(kelem, knode-1, knode, eType, nprop, kprop, p)
                    nprop = kprop
  
negative.code:     if (Init%NDiv==1) then
         Init%PropsB(1:Init%NPropB, 1:PropSetsBCol) = Init%PropSetsB(1:Init%NPropB, 1:PropSetsBCol)
      else if (Init%NDiv>1) then
         Init%PropsB(1:Init%NPropB, 1:PropSetsBCol) = TempProps(1:Init%NPropB, 1:PropSetsBCol)
      endif
  
      ! --- Cables and rigid link properties (these cannot be subdivided, so direct copy of inputs)
  