anchor.diff:  
   
                  ! Calculation of B_minus
  -            TransFuncB_minus  =  InitInp%Gravity*InitInp%Gravity / ( 4.0_ReKi * Omega_n * Omega_m ) &
  -                                 * COSHNumOvrCOSHDen(k_nm, InitInp%WtrDpth, z)  * D_minus / ( Omega_n - Omega_m )
  +            TransFuncB_minus  =  InitInp%Gravity*InitInp%Gravity / ( 8.0_ReKi * Omega_n * Omega_m ) &          
  +                                 * COSHNumOvrCOSHDen(k_nm, InitInp%WtrDpth, z)  * D_minus / ( Omega_n - Omega_m )       ! 4.0_ReKi
   
   
            ENDIF
  
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:                   ! 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 =  WaveElevC0Norm(n)
   
                   !> ### Calculate the array of \f$ K^+(\omega_n) \f$ for the first term of the velocity, acceleration, and pressure. ###
  