Vector flow maps using PyTOUGH

Hi All,

Can anyone tell me, how to plot flow vectors using t2listing class?

I can easily plot scalar field from my results using layer_plot method, but do not know what to assign to flow attribute of the layer_plot method.

I have all kinds of output files in my folder, i.e. OUTPUT.dat, OUTPUT_CONNE.csv, OUTPUT_ELEME.csv.

Please look at the code below. I have successfuly created both flux_matrix and grid objects, but I simply do not know what assign to flow attribute. I am interested to show flow vectors, but it can be heat flow as well. Any help is much appreciated.

from t2listing import *
from mulgrids import *
from t2data import *
import numpy as np
import matplotlib.pyplot as plt
import scipy as sc

geo = mulgrid('geom.dat')
bm = geo.block_mapping(geo)
results = t2listing('OUTPUT')

results.last() # moves to the last time step
p = results.element['PRES']

limits = ((0, 3000),(0,5000))

dat = t2data()
dat.grid = t2grid().fromgeo(geo)
fm = dat.grid.flux_matrix(geo, bm)

geo.layer_plot(-4275,p/1.0e5, 'Pressure ', 'bar',  block_names=False, plot_limits = limits, colourmap='coolwarm', flux_matrix=fm, grid = dat.grid, flow = ?)

Best regards,
Maciej