anchor.diff:                          Omega1 = (J + K) * InitInp%WaveDOmega        ! the mth frequency -- \mu^- + n = m
                           Omega2 = K * InitInp%WaveDOmega              ! the nth frequency
   
  -
                              ! Find the Wave amplitudes 1 and 2
                           aWaveElevC1 = CMPLX( InitInp%WaveElevC0(1,J+K), InitInp%WaveElevC0(2,J+K), SiKi)  / InitInp%NStepWave2
                           aWaveElevC2 = CMPLX( InitInp%WaveElevC0(1,K),   InitInp%WaveElevC0(2,K),   SiKi)  / InitInp%NStepWave2
   
                              ! Set the (omega1,omega2,beta1,beta2) point we are looking for.
                           Coord4 = (/ REAL(Omega1,SiKi), REAL(Omega2,SiKi), InitInp%WaveDirArr(J+K), InitInp%WaveDirArr(K) /)
  +
  +                           ! Apply local Z rotation to heading angle (degrees) to put wave direction into the local (rotated) body frame
  +                        Coord4(3) = Coord4(3) - RotateZdegOffset
  +                        Coord4(4) = Coord4(4) - RotateZdegOffset
   
                              ! get the interpolated value for F(omega1,omega2,beta1,beta2)  --> QTF_Value
                           CALL WAMIT_Interp4D_Cplx( Coord4, TmpData4D, DiffQTFData%Data4D%WvFreq1, DiffQTFData%Data4D%WvFreq2, &
  
positive.code:                            RETURN
                          ENDIF
  
                             ! Calculate this value and add it to what we have so far.
                          TmpHMinusC = TmpHMinusC + aWaveElevC1 * CONJG(aWaveElevC2) * QTF_Value
  
  
negative.code: 
           ! Wave information and QTF temporary
        COMPLEX(SiKi)                                      :: QTF_Value            !< Temporary complex number for QTF
        COMPLEX(SiKi), ALLOCATABLE                         :: TmpComplexArr(:,:)   !< Temporary complex array for frequency domain of one complete load component
        COMPLEX(SiKi)                                      :: TmpHMinusC           !< Temporary variable for holding the current value of \f$ H^- \f$
        COMPLEX(SiKi)                                      :: aWaveElevC1          !< Wave elevation of the first  frequency component.  NStepWave2 factor removed.
        COMPLEX(SiKi)                                      :: aWaveElevC2          !< Wave elevation of the second frequency component.  NStepWave2 factor removed.
  