Hi
In prosou.f there is the code where the source terms are added when you add some sources in your model, for example.
There is a part for the mass and a part for the dynamic equation.
In my opinion, you could add, in this subroutine, a term to FU (x axis) or FV (y axis) or both (if necessary) to force the flow at the location you want.
You could look in this subroutine for the source termes of wind or other physical proceses to better understand the method
As the source terms could be applied to each node, you could only limit the area by selecting only few nodes
Kind regards