1FUNCTION fbed(x,y) 2USE types 3IMPLICIT NONE 4REAL(KIND=dp) :: x, y, fbed 5REAL(KIND=dp) :: Lxy, fsurf 6 7Lxy = 100.0 8 9fbed = 10.0*(COS(6.0*Pi*x/Lxy)*COS(6.0*Pi*y/Lxy) -1.0) 10 11END FUNCTION fbed 12 13