anchor.diff:                    ! 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)           
  +                  CALL SetNewProp(kprop, TempProps(Prop1, 2), TempProps(Prop1, 3), Init%PropSetsB(Prop1, 4), d1 + J*dd, t1 + J*dt,  TempProps)           
                     kelem = kelem + 1
                     CALL SetNewElem(kelem, knode-1, knode, eType, nprop, kprop, p)
                     nprop = kprop
  
positive.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) = TempPropsB(1:Init%NPropB, 1:PropSetsBCol)
      endif
  
      ! --- Cables and rigid link properties (these cannot be subdivided, so direct copy of inputs)
  
negative.code:           Init%MemberNodes(I,           1) = Node1
            Init%MemberNodes(I, Init%NDiv+1) = Node2
  
            IF  ( ( .not. EqualRealNos(TempProps(Prop1, 2),TempProps(Prop2, 2) ) ) &
             .OR. ( .not. EqualRealNos(TempProps(Prop1, 3),TempProps(Prop2, 3) ) ) &
             .OR. ( .not. EqualRealNos(TempProps(Prop1, 4),TempProps(Prop2, 4) ) ) )  THEN
            
               CALL Fatal(' Material E,G and rho in a member must be the same')
               RETURN
  