anchor.diff:            dy = ( y2 - y1 )/Init%NDiv
             dz = ( z2 - z1 )/Init%NDiv
             
  -          d1 = TempPropsB(Prop1, 5)
  -          t1 = TempPropsB(Prop1, 6)
  -
  -          d2 = TempPropsB(Prop2, 5)
  -          t2 = TempPropsB(Prop2, 6)
  +          d1 = TempProps(Prop1, 5)
  +          t1 = TempProps(Prop1, 6)
  +
  +          d2 = TempProps(Prop2, 5)
  +          t2 = TempProps(Prop2, 6)
             
             dd = ( d2 - d1 )/Init%NDiv
             dt = ( t2 - t1 )/Init%NDiv
  
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), TempPropsB(Prop1, 4), d1+dd, t1+dt, TempPropsB)           
                 kelem = kelem + 1
                 CALL SetNewElem(kelem, Node1, knode, eType, Prop1, kprop, p)  
                 nprop = kprop
  
negative.code:        !  We don't know how many properties will be needed, so allocated to size MaxNProp
         MaxNProp   = Init%NPropSetsB + Init%NElem*NNE ! Maximum possible number of property sets (temp): This is property set per element node, for all elements (bjj, added Init%NPropSets to account for possibility of entering many unused prop sets)
         CALL AllocAry(TempMembers, p%NMembers,    MembersCol , 'TempMembers', ErrStat2, ErrMsg2); if(Failed()) return
         CALL AllocAry(TempProps,  MaxNProp,      PropSetsBCol,'TempProps',  ErrStat2, ErrMsg2); if(Failed()) return
         TempProps = -9999.
         TempMembers                      = p%Elems(1:p%NMembers,:)
         TempProps(1:Init%NPropSetsB, :) = Init%PropSetsB   
  
         kelem = 0
         knode = Init%NJoints
         kprop = Init%NPropSetsB
         DO I = 1, p%NMembers !the first p%NMembers rows of p%Elems contain the element information
            ! Member data
            Node1 = TempMembers(I, 2)
  