anchor.diff:                       ! First get the wave amplitude -- must be reconstructed from the WaveElevC0 array.  First index is the real (1) or
                        ! imaginary (2) part.  Divide by NStepWave2 to remove the built in normalization in WaveElevC0.  Note that the phase
                        ! shift associated with the (x,y) location is accounted for by the WaveElevxyPrime0 variable.
  -!                  WaveElevC_n = CMPLX( InitInp%WaveElevC0(1,n), InitInp%WaveElevC0(2,n)) / InitInp%NStepWave2
                     WaveElevC_n =  WaveElevC0Norm(n)
    
   
  
positive.code:       !! The equation is given by:
        !!
        !! \f$ B_{nm}^-(z, \omega_n, \omega_m, \theta_n, \theta_m) =\frac{g^2}{\omega_n \omega_m}
        !!             \cdot    \frac{1}{4} \frac{\cosh\left[k_{nm}^-(h+z)\right]}{\cosh\left[k_{nm}^-(h)\right]}
        !!                      \frac{D_{nm}^-}{\omega_n - \omega_m}   \f$
        !!
        FUNCTION TransFuncB_minus(n,m,k_n,k_m,z)
  
negative.code:                   D_minus     =  TransFuncD_minus(n,m,k_n,k_m,R_n,R_m)
  
                       !> Calculate the value of 
                       !!    \f$ L^-_{nm} = \frac{1}{8} \left[ 
                       !!             \frac{D^-_{nm} - |\vec{k}_n| |\vec{k}_m| \cos(\theta_n - \theta_m) - R_n R_m}{\sqrt{R_n R_m}}
                       !!          +  (R_n+R_m) \right] \f$
                       !!
                       !!    The value of \f$ D^-_{nm} \f$ is found from by the ::TransFuncD_minus routine.
  
                    L_minus  =  (( D_minus - k_n * k_m * COS(D2R*InitInp%WaveDirArr(n) - D2R*InitInp%WaveDirArr(m)) - R_n * R_m )/SQRT( R_n * R_m ) + R_n + R_m) / 8.0_ReKi !4.0_ReKi
  
                       ! Calculate the terms \f$ n,m \f$ necessary for calculations
  
  