ccwes_waveamplitude Subroutine

public pure subroutine ccwes_waveamplitude(omega, domega, hs, wave_amplitude)

中国沿海波能谱模型计算频段有义波幅 China Coastal Wave Energy Spectrum

Arguments

Type IntentOptional Attributes Name
real(kind=rk), intent(in) :: omega

波浪角频率, rad/s
real(kind=rk), intent(in) :: domega

频率间隔, rad/s
real(kind=rk), intent(in) :: hs

有义波高, m
real(kind=rk), intent(out) :: wave_amplitude

频段有义波幅, m

Contents

Source Code

ccwes_waveamplitude

Source Code

    pure subroutine ccwes_waveamplitude(omega, domega, hs, wave_amplitude)
        real(kind=rk), intent(in) :: omega  !! 波浪角频率, rad/s
        real(kind=rk), intent(in) :: domega !! 频率间隔, rad/s
        real(kind=rk), intent(in) :: hs     !! 有义波高, m
        real(kind=rk), intent(out) :: wave_amplitude !! 频段有义波幅, m
        real(kind=rk) :: s

        associate (u => 6.28_rk*sqrt(hs))
            s = 0.74_rk*exp(-(g/(u*omega))**2)/omega**5
            wave_amplitude = sqrt(2*s*domega)  ! A = sqrt(2*S_\omega*d\omega)
        end associate

    end subroutine ccwes_waveamplitude