1     ! constant Pi
2 !radovan: this is more elegant but is unfortunately not standard
3 !   real(REALK), parameter :: PI = acos(-1.0_REALK)
4     real(REALK), parameter :: PI = 3.14159265358979_REALK
5