anchor.diff:              ! Create the node on the mesh
            CALL MeshPositionNode (u%Mesh, ThisBodyNum, XYZloc, ErrStatTmp, ErrMsgTmp, orientation )
            CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName)
  - 
  +
               ! Create the mesh element
            CALL MeshConstructElement (  u%Mesh, ELEMENT_POINT, ErrStatTmp, ErrMsgTmp, ThisBodyNum )
            CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName)
         ENDDO
  - 
  +
         CALL MeshCommit ( u%Mesh, ErrStatTmp, ErrMsgTmp )
         CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName)
         IF ( ErrStat >= AbortErrLev ) THEN
  
positive.code: ThisBodyNum=1
           !> If the MnDrift method will be used, call the subroutine to calculate the force time series
        IF ( p%MnDriftF ) THEN
           
              ! Tell our nice users what is about to happen that may take a while:
           CALL WrScr ( ' Calculating second order mean drift force.' )
  
  
negative.code:          Data3D%WvFreq1( 1 )                 = 0.0_SiKi
           Data3D%WvFreq1( 2 )                 = MAX( TmpWvFreq1(1) - 10.0_SiKi*EPSILON(0.0_SiKi), 0.0_SiKi )  ! make sure the Frequencies are still monotonically increasing
        ENDIF
  
           ! add the two points for step-change after last entered frequency
        Data3D%WvFreq1( Data3D%NumWvFreq1-1 )     = Data3D%WvFreq1(Data3D%NumWvFreq1-2) + 10.0_SiKi*EPSILON(0.0_SiKi)
        Data3D%WvFreq1( Data3D%NumWvFreq1   )     = HUGE(1.0_SiKi)/5 ! floating overflow occurs later with arithmetic so I divided by a small constant
  