anchor.diff:           Data4D%WvFreq1( 1 )                 = 0.0_SiKi
            Data4D%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
         Data4D%WvFreq1( Data4D%NumWvFreq1-1 )     = Data4D%WvFreq1(Data4D%NumWvFreq1-2) + 10.0_SiKi*EPSILON(0.0_SiKi)
         Data4D%WvFreq1( Data4D%NumWvFreq1   )     = HUGE(1.0_SiKi)/5 ! floating overflow occurs later with arithmetic so I divided by a small constant
  
positive.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
  
negative.code:       INTEGER(IntKi)                                     :: K                    !< Generic counter
        TYPE(FFT_DataType)                                 :: FFT_Data             !< Temporary array for the FFT module we're using. For the first  term in the equation.
        CHARACTER(*), PARAMETER                            :: RoutineName = 'SumQTF_InitCalc'
  
  
           ! Wave information and QTF temporary
        COMPLEX(SiKi)                                      :: QTF_Value           !< Temporary complex number for QTF
  