anchor.diff:  
            f        = Inv2Pi*Omega
            fp       = 1/Tp
  -         fpOvrf4  = (fp/f)**4.0
  +         fpOvrf4  = (fp/f)**4
   
   
            ! Compute the normalising factor:
  
positive.code:             Sigma = 0.09
           END IF
  
           Alpha    = EXP( ( -0.5*( ( (f/fp) - 1.0 )/Sigma )**2.0 ) )
  
  
           ! Compute the wave spectrum:
  
negative.code:          ! Compute the JONSWAP wave spectrum, unless Omega is zero, in which case,
           !   return zero:
  
        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
  