anchor.diff:        m%udot_TP    = (/u%TPMesh%TranslationVel( :,1), u%TPMesh%RotationVel(:,1)/)
         m%udotdot_TP = (/u%TPMesh%TranslationAcc( :,1), u%TPMesh%RotationAcc(:,1)/)
         ! Inputs on interior nodes:
  -      CALL ConstructUFL( u, p, m, m%UFL )
  +      CALL GetExtForceOnInternalDOF( u, p, m, m%UFL )
   
         !________________________________________
         ! Set motion outputs on y%Y2mesh
  
positive.code:    nColumns = MembersCol
  endif
  DO I = 1, p%NMembers
     CALL ReadAry( UnIn, SDInputFile, Dummy_IntAry, nColumns, 'Members', 'Member number and connectivity ', ErrStat2,ErrMsg2, UnEc); if(Failed()) return
     Init%Members(I,1:nColumns) = Dummy_IntAry(1:nColumns)
  ENDDO   
  IF (Check( p%NMembers < 1 , 'NMembers must be > 0')) return
  
negative.code: !> Construct force vector on internal DOF (L) from the values on the input mesh 
  !! First, the full vector of external forces is built on the non-reduced DOF
  !! Then, the vector is reduced using the Tred matrix
  SUBROUTINE GetExtForceOnInternalDOF( u, p, m, UFL )
     type(SD_InputType),     intent(in   )  :: u ! Inputs
     type(SD_ParameterType), intent(in   )  :: p ! Parameters
     type(SD_MiscVarType),   intent(inout)  :: m ! Misc, for storage optimization of Fext and Fext_red
  