1! { dg-do compile }
2
3module yemdyn
4   implicit none
5   integer, parameter :: jpim = selected_int_kind(9)
6   integer, parameter :: jprb = selected_real_kind(13,300)
7   real(kind=jprb) :: elx
8   real(kind=jprb), allocatable :: xkcoef(:)
9   integer(kind=jpim),allocatable :: ncpln(:), npne(:)
10end module yemdyn
11
12subroutine suedyn
13
14   use yemdyn
15
16   implicit none
17
18   integer(kind=jpim) :: jm, jn
19   real(kind=jprb) :: zjm, zjn, zxxx
20
21   jn=0
22   do jm=0,ncpln(jn)
23      zjm=real(jm,jprb) / elx
24      xkcoef(npne(jn)+jm) = - zxxx*(zjm**2)**0.5_jprb
25   end do
26
27end subroutine suedyn
28
29