Thanks for reply.
Your first option is, in my opinion, a little bit redundant because require a double pass: first run the simulation with DFLUX subroutine, writing a text file containing association between node id and flux value; second execute a python script that read that text and update...