anchor.diff:        RETURN
      ENDIF
   
  -   ! ---
  +   ! --- Reindexing JointsID and MembersID into Nodes and Elems arrays
  +   ! NOTE: need NNode and NElem 
      CALL SD_ReIndex_CreateNodesAndElems(Init,p, ErrStat2, ErrMsg2);  if(Failed()) return
      
  -   ! Initialize TempMembers
  -   TempMembers = p%Elems(1:p%NMembers,:)
  -   
  -   ! Initialize Temp property set, first user defined sets
  -   TempPropsBeams = 0
  -   TempPropsCable = 0
  -   TempPropsRigid = 0
  -   TempPropsBeams(1:Init%NPropSets, :)      = Init%PropSets   
  -   TempPropsCable(1:Init%NCablePropSets, :) = Init%CablePropSets
  -   TempPropsRigid(1:Init%NRigidPropSets, :) = Init%RigidPropSets
  -   print*,'TempPropsBeams(:,1:5)'
  -   DO I=1,Init%NPropSets
  -      print*,TempPropsBeams(I,1:5)
  -   enddo
  -   print*,'TempPropsCable(:,:)'
  -   DO I=1,Init%NCablePropSets
  -      print*,TempPropsCable(I,1:4)
  -   enddo
  -
  -   print*,'TempPropsRigid(:,:)'
  -   DO I=1,Init%NRigidPropSets
  -      print*,TempPropsRigid(I,:)
  -   enddo
   
      
      ! Initialize boundary constraint vector
  
positive.code:              
               IF ( CreateNewProp ) THEN   
                    ! create a new property set 
                    ! k, E, G, rho, d, t, Init
                    
                    kprop = kprop + 1
                    CALL SetNewProp(kprop, TempPropsBeams(Prop1, 2), TempPropsBeams(Prop1, 3), Init%PropSets(Prop1, 4), d1 + J*dd, t1 + J*dt,  TempPropsBeams)           
                    kelem = kelem + 1
  
negative.code:    CALL AllocAry(Init%BCs,        6*p%NReact,    2,          'Init%BCs',        ErrStat2, ErrMsg2); CALL SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,'SD_Discrt') !!!! RRD: THIS MAY NEED TO CHANGE IF NOT ALL NODES ARE RESTRAINED
     CALL AllocAry(Init%IntFc,      6*Init%NInterf,2,          'Init%IntFc',      ErrStat2, ErrMsg2); CALL SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,'SD_Discrt')
     
     CALL AllocAry(TempReacts,      p%NReact,      ReactCol,   'TempReacts',      ErrStat2, ErrMsg2); CALL SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,'SD_Discrt')
  
     IF ( ErrStat >= AbortErrLev ) THEN
  