1FUNCTION fsurf(x,y) 2USE types 3IMPLICIT NONE 4REAL(KIND=dp) :: x, y, fsurf 5 6fsurf = - x * TAN(0.5d0*Pi/180.0d0) 7 8END FUNCTION fsurf 9 10