Hello,
You have an example in telemac2d, called "cone", that uses an emptied version of propag.f that does nothing, and the velocity field and depth is actually given by the initial conditions done in subroutine condin.f. You can inspire from this: your version of propag will give the data of velocity and depth.
However if you give a velocity field and depth from another program, the continuity equation may be approximately solved and you will certainly have problems of mass conservation and monotonicity.
With best regards,
Jean-Michel Hervouet