anchor.diff:              ! Tell our nice users what is about to happen that may take a while:
            CALL WrScr ( ' Calculating second order sum-frequency force using the full quadratic transfer function.' )
   
  -         CALL SumQTF_InitCalc( InitInp, p, ThisBodyNum, SumQTFData, SumQTFForce, ErrMsgTmp, ErrStatTmp )
  +         CALL SumQTF_InitCalc( InitInp, p, SumQTFData, SumQTFForce, ErrMsgTmp, ErrStatTmp )
            CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName )
            IF ( ErrStat >= AbortErrLev ) THEN
               CALL CleanUp
  
positive.code:    !! Since the frequency range of the QTF has not yet been checked, we will do that now.
     !!
     !-------------------------------------------------------------------------------------------------------------------------------
     SUBROUTINE MnDrift_InitCalc( InitInp, p, ThisBodyNum, MnDriftData, MnDriftForce, ErrMsg, ErrStat )
  
        IMPLICIT NONE
  
        TYPE(WAMIT2_InitInputType),         INTENT(IN   )  :: InitInp              !< Input data for initialization routine
        TYPE(WAMIT2_ParameterType),         INTENT(IN   )  :: p                    !< Parameters
        INTEGER(IntKi),                     INTENT(IN   )  :: ThisBodyNum          !< Which body are we on
        TYPE(W2_DiffData_Type),             INTENT(INOUT)  :: MnDriftData          !< Data storage for the MnDrift method.  Set to INOUT in case we need to convert 4D to 3D
        REAL(SiKi),  ALLOCATABLE,           INTENT(  OUT)  :: MnDriftForce(:)      !< Force data.  Index 1 is the force component.  Constant for all time.
        CHARACTER(*),                       INTENT(  OUT)  :: ErrMsg
  
negative.code:             ! Tell our nice users what is about to happen that may take a while:
           CALL WrScr ( " Calculating second order difference-frequency force using the Newman's approximation." )
  
           CALL NewmanApp_InitCalc( InitInp, p, NewmanAppData, NewmanAppForce, ErrMsgTmp, ErrStatTmp )
           CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName )
           IF ( ErrStat >= AbortErrLev ) THEN
              CALL CleanUp
  