anchor.diff:           ! Compute the JONSWAP wave spectrum, unless Omega is zero, in which case,
            !   return zero:
   
  -      IF ( Omega == 0.0 )  THEN  ! When .TRUE., the formulation below is ill-conditioned; thus, the known value of zero is returned.
  +      IF ( EqualRealNos(Omega, 0.0_ReKi) )  THEN  ! When .TRUE., the formulation below is ill-conditioned; thus, the known value of zero is returned.
   
   
            JONSWAP  = 0.0
  
positive.code:             Sigma = 0.09
           END IF
  
           Alpha    = EXP( ( -0.5*( ( (f/fp) - 1.0 )/Sigma )**2.0 ) )
  
  
           ! Compute the wave spectrum:
  
negative.code: 
           f        = Inv2Pi*Omega
           fp       = 1/Tp
           fpOvrf4  = (fp/f)**4
  
  
           ! Compute the normalising factor:
  