CoCalc -- fbed.F90

GitHub Repository: ElmerCSC/elmerfem
Path: blob/devel/elmerice/examples/Test_MshGlacierSynthetic/fbed.F90
³²⁰⁶ views

1
FUNCTION fbed(x,y) 
2
USE types
3
IMPLICIT NONE
4
REAL(KIND=dp) :: x, y, fbed
5
REAL(KIND=dp) :: Lxy, fsurf
6

7
Lxy = 5.0e3_dp
8

9
fbed = fsurf(x,y) -1000.0d0*(1.0 - 0.5d0*SIN(2.0*Pi*x/Lxy)*SIN(2.0*Pi*y/Lxy))
10

11
END FUNCTION fbed
12

13