| ! { dg-do compile } |
| |
| module yemdyn |
| implicit none |
| integer, parameter :: jpim = selected_int_kind(9) |
| integer, parameter :: jprb = selected_real_kind(13,300) |
| real(kind=jprb) :: elx |
| real(kind=jprb), allocatable :: xkcoef(:) |
| integer(kind=jpim),allocatable :: ncpln(:), npne(:) |
| end module yemdyn |
| |
| subroutine suedyn |
| |
| use yemdyn |
| |
| implicit none |
| |
| integer(kind=jpim) :: jm, jn |
| real(kind=jprb) :: zjm, zjn, zxxx |
| |
| jn=0 |
| do jm=0,ncpln(jn) |
| zjm=real(jm,jprb) / elx |
| xkcoef(npne(jn)+jm) = - zxxx*(zjm**2)**0.5_jprb |
| end do |
| |
| end subroutine suedyn |
| |