anchor.diff:              ! Tell our nice users what is about to happen that may take a while:
            CALL WrScr ( ' Calculating second order difference-frequency force using the full quadratic transfer function.' )
   
  -         CALL DiffQTF_InitCalc( InitInp, p, ThisBodyNum, DiffQTFData, DiffQTFForce, ErrMsgTmp, ErrStatTmp )
  +         CALL DiffQTF_InitCalc( InitInp, p, DiffQTFData, DiffQTFForce, ErrMsgTmp, ErrStatTmp )
            CALL SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName )
            IF ( ErrStat >= AbortErrLev ) THEN
               CALL CleanUp
  
positive.code:             ! 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 SetErrStat( ErrStatTmp, ErrMsgTmp, ErrStat, ErrMsg, RoutineName )
           IF ( ErrStat >= AbortErrLev ) THEN
              CALL CleanUp
  
negative.code: 
              ! Set rotation
              ! NOTE: RotateZMatrixT is the rotation from local to global.
           RotateZMatrixT(:,1) = (/  cos(InitInp%PtfmRefztRot(IBody)), -sin(InitInp%PtfmRefztRot(IBody)) /)
           RotateZMatrixT(:,2) = (/  sin(InitInp%PtfmRefztRot(IBody)),  cos(InitInp%PtfmRefztRot(IBody)) /)
  
              ! Loop through all the frequencies
           DO J=1,InitInp%NStepWave2
  