Hello,
The bed slope is not a variable that is kept across the program, it is built when necessary.
Two options:
1) You need the largest local slope (a scalar value) : there is an example of call to subroutine slope.f in subroutine propag.f, in this case the result is given for every triangle (the slope is constant on every triangle, the bottom topography ZF being linear).
2) You need a vector with gradient in x and y : this is a standard call to subroutine vector, with 'GRADF X' and 'GRADF Y'. To understand how VECTOR works see the guide for programming in the Telemac system. The result of VECTOR is an integral and to get nodal values one must divide by the integral of test functions, or multiply by UNSV2D%R(I) for a point I, UNSV2D being the inverse of the integral of test functions.
This is it, to start with,
With best regards,
Jean-Michel Hervouet