anchor.diff:               
                IF ( CreateNewProp ) THEN   
                     ! create a new property set 
  -                  ! k, E, G, rho, d, t, Init
  -                  
  +                  ! 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
  
positive.code:                   nprop = kprop                
               ELSE
                    kelem = kelem + 1
                    CALL SetNewElem(kelem, knode-1, knode, eType, nprop, nprop, p)         
                     
               ENDIF
            ENDDO
            
            ! the element connect to Node2
            kelem = kelem + 1
            CALL SetNewElem(kelem, knode, Node2, eType, nprop, Prop2, p)                
  
         ENDDO ! loop over all members
  
      ELSE ! NDiv = 1
  
         Init%MemberNodes(1:p%NMembers, 1:2) = p%Elems(1:Init%NElem, 2:3)   
  
      ENDIF ! if NDiv is greater than 1
  
      ! set the props in Init
      Init%NProp = kprop
      CALL AllocAry(Init%Props, Init%NProp, PropSetsCol,  'Init%Props', ErrStat2, ErrMsg2); CALL SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,'SD_Discrt')
      IF (ErrStat >= AbortErrLev ) THEN
         CALL CleanUp_Discrt()
         RETURN
      ENDIF
      !Init%Props(1:kprop, 1:Init%PropSetsCol) = TempProps
      Init%Props = TempPropsBeams(1:Init%NProp, :)  !!RRD fixed it on 1/23/14 to account for NDIV=1
      print*,'Init%Props'
      print*, Init%Props
  
  
negative.code:       RETURN
     ENDIF
     
     ! Total number of element   
     Init%NElem = p%NMembers*Init%NDiv                          
     ! Total number of nodes - Depends on division and number of nodes per element
     Init%NNode = Init%NJoints + ( Init%NDiv - 1 )*p%NMembers 
     Init%NNode = Init%NNode + (NNE - 2)*Init%NElem
     !bjj: replaced with max value instead of NNE: Init%MembersCol = Init%MembersCol + (NNE - 2) 
     MaxNProp   = Init%NPropSets + 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)
     
     ! check the number of interior modes
     IF ( p%Nmodes > 6*(Init%NNode - Init%NInterf - p%NReact) ) THEN
  