Hello,
Yes, those subroutine are the right one. However, there are not yet available in 3D.
In 3D, I suggest you use the subroutine SOURCE and loop over the NSCE. Then you know what are the ins and the outs from the list of coordinates in your steering file and your know the node number these correspond to with ISCE.
You can therefore say something like
if I.EQ.1 then
Q1 = FLOWFUNCTION( H(ISCE(I)), H(ISCE(I+1)) )
QSCE(1) = Q1
elseif I.EQ.2 then
QSCE(2) = -Q1
endif
Of course you have to define your FLOWFUNCTION ...
Hope this helps,
Sébastien