anchor.diff:            Init%MemberNodes(I,           1) = Node1
             Init%MemberNodes(I, Init%NDiv+1) = Node2
             
  -          IF  ( ( .not. EqualRealNos(TempPropsBeam(Prop1, 2),TempPropsBeam(Prop2, 2) ) ) &
  -           .OR. ( .not. EqualRealNos(TempPropsBeam(Prop1, 3),TempPropsBeam(Prop2, 3) ) ) &
  -           .OR. ( .not. EqualRealNos(TempPropsBeam(Prop1, 4),TempPropsBeam(Prop2, 4) ) ) )  THEN
  +          IF  ( ( .not. EqualRealNos(TempPropsBeams(Prop1, 2),TempPropsBeams(Prop2, 2) ) ) &
  +           .OR. ( .not. EqualRealNos(TempPropsBeams(Prop1, 3),TempPropsBeams(Prop2, 3) ) ) &
  +           .OR. ( .not. EqualRealNos(TempPropsBeams(Prop1, 4),TempPropsBeams(Prop2, 4) ) ) )  THEN
             
                CALL Fatal(' Material E,G and rho in a member must be the same')
                RETURN
  
positive.code:       ! Column 6: member type
        p%Elems(iMem, iMType) = Init%Members(iMem, iMType) ! 
     END DO !iMem, loop through members
     
     ! Initialize TempMembers
     TempMembers = p%Elems(1:p%NMembers,:)
     
     ! Initialize Temp property set, first user defined sets
     TempPropsBeam = 0
     TempPropsBeam(1:Init%NPropSets, :) = Init%PropSets   
     print*,'TempPropsBeam(:,1:5)'
     DO I=1,Init%NPropSets
        print*,TempPropsBeam(I,1:5)
     enddo
  
  
     TempPropsCable = 0
     TempPropsCable(1:Init%NCablePropSets, :) = Init%CablePropSets
     print*,'TempPropsCable(:,:)'
     DO I=1,Init%NCablePropSets
        print*,TempPropsCable(I,1:4)
     enddo
  
     TempPropsRigid = 0
     TempPropsRigid(1:Init%NRigidPropSets, :) = Init%RigidPropSets
     print*,'TempPropsRigid(:,:)'
     DO I=1,Init%NRigidPropSets
        print*,TempPropsRigid(I,:)
  
negative.code:                ! create a new property set 
                 ! k, E, G, rho, d, t, Init
                 kprop = kprop + 1
                 CALL SetNewProp(kprop, TempPropsBeams(Prop1, 2), TempPropsBeams(Prop1, 3), TempPropsBeams(Prop1, 4), d1+dd, t1+dt, TempPropsBeams)           
                 kelem = kelem + 1
                 CALL SetNewElem(kelem, Node1, knode, eType, Prop1, kprop, p)  
                 nprop = kprop              
  