anchor.diff:        do iNode=2,3
            p%Elems(iMem,iNode) = FINDLOCI(Init%Joints(:,1), Init%Members(iMem, iNode) ) 
            if (p%Elems(iMem,iNode)<=0) then
  -            CALL Fatal(' MemberID '//TRIM(Num2LStr(Init%Members(iMem,1)))//' has JointID'//TRIM(Num2LStr(iNode-1))//' = '// TRIM(Num2LStr(Init%Members(iMem, iNode)))//' which is not in the joint list !')
  +            CALL Fatal(' MemberID '//TRIM(Num2LStr(Init%Members(iMem,1)))//' has JointID'//TRIM(Num2LStr(iNode-1))//' = '// TRIM(Num2LStr(Init%Members(iMem, iNode)))//' which is not in the joint list!')
               return
            endif
         enddo
  
positive.code:          end if
  
           if (p%Elems(iMem,n)<=0) then
              CALL Fatal(' MemberID '//TRIM(Num2LStr(Init%Members(iMem,1)))//' has PropSetID'//TRIM(Num2LStr(n-3))//' = '//TRIM(Num2LStr(Prop))//' which is not in the'//trim(sType)//' table!')
           endif
        END DO !n, loop through property ids         
        ! Column 6: member type
  
negative.code: 
     ! --- Initialize Nodes
     Init%Nodes = 0   
     do I = 1,Init%NJoints
        Init%Nodes(I, 1) = I                ! JointID replaced by index I
        Init%Nodes(I, 2) = Init%Joints(I, 2)
        Init%Nodes(I, 3) = Init%Joints(I, 3)
        Init%Nodes(I, 4) = Init%Joints(I, 4)
     enddo
  
     ! --- Re-Initialize Reactions, pointing to index instead of JointID
     do I = 1, p%NReact
        JointID=p%Reacts(I,1)
        p%Reacts(I,1) = FINDLOCI(Init%Joints(:,1), JointID ) ! Replace JointID with Index
        if (p%Reacts(I,1)<=0) then
           CALL Fatal('Reaction joint table: line '//TRIM(Num2LStr(I))//' refers to JointID '//trim(Num2LStr(JointID))//' which is not in the joint list!')
           return
        endif
     enddo
  
     ! --- Re-Initialize interface joints, pointing to index instead of JointID
     Init%IntFc = 0
     do I = 1, Init%NInterf
        JointID=Init%Interf(I,1)
        Init%Interf(I,1) = FINDLOCI(Init%Joints(:,1), JointID )
        if (Init%Interf(I,1)<=0) then
           CALL Fatal('Interface joint table: line '//TRIM(Num2LStr(I))//' refers to JointID '//trim(Num2LStr(JointID))//' which is not in the joint list!')
           return
        endif
     enddo
  
     ! Change numbering in concentrated mass matrix
     do I = 1, Init%NCMass
        JointID = Init%CMass(I,1)
        Init%CMass(I,1) = FINDLOCI(Init%Joints(:,1), JointID )
        if (Init%Interf(I,1)<=0) then
           CALL Fatal('Concentrated mass table: line '//TRIM(Num2LStr(I))//' refers to JointID '//trim(Num2LStr(JointID))//' which is not in the joint list!')
           return
        endif
     enddo
  
  
     ! --- Initialize Elems, starting with each member as an element (we'll take NDiv into account later)
     p%Elems = 0
     ! --- Replacing "MemberID"  "JointID", and "PropSetID" by simple index in this tables
     DO iMem = 1, p%NMembers
        ! Column 1  : member index (instead of MemberID)
        p%Elems(iMem,     1)  = iMem
  