anchor.diff:        integer(IntKi) :: m
         nDOF_Unconstrained=0
         do i = 1,Init%NNode
  -         if (Init%Nodes(i,iJointType) == idJointCantilever ) then
  +         if (int(Init%Nodes(i,iJointType)) == idJointCantilever ) then
               nDOF_Unconstrained = nDOF_Unconstrained + 6
            else
               m = Init%NodesConnE(i,1) ! Col1: number of elements connected to this joint
  
positive.code:       integer(IntKi) :: nDOF
        integer(IntKi) :: aID, ia ! assembly ID, and index in IRA
        integer(IntKi) :: iNode
        integer(IntKi) :: I
        integer(IntKi) :: nc !< Number of DOF after constraints applied
        integer(IntKi) :: nj
  
        ! --- Misc inits
        nullify(IDOFNew)
  
negative.code:    
     ! total unconstrained degrees of freedom of the system 
     Init%TDOF = nDOF_Unconstrained()
     print*,'nDOF_unconstrained:',Init%TDOF, ' (if all Cantilever, it would be: ',6*Init%NNode,')'
  
     CALL AllocAry( Init%K, Init%TDOF, Init%TDOF , 'Init%K',  ErrStat2, ErrMsg2); if(Failed()) return; ! system stiffness matrix 
     CALL AllocAry( Init%M, Init%TDOF, Init%TDOF , 'Init%M',  ErrStat2, ErrMsg2); if(Failed()) return; ! system mass matrix 
  