I use my mesh using pointwise ( 7M DoF, example case because of GPU memory limit 10GB. )

I attached my ini file. (How can ini file or msh file directly?)

```
[backend-cuda]
device-id = round-robin
mpi-type = standard
[constants]
;gamma = 1.4
;Pr = 0.72
nu = 1.6667e-5
Uin = 0.9902681
Vin = 0.1391731
w=0
Pc = 1
ac-zeta = 2.5 ;3
[solver]
system = ac-navier-stokes ;incompressible n-s
order = 3 ;order of polynomial solution basis
;anti-alias = flux
[solver-time-integrator]
formulation = dual
scheme = sdirk33 ;backward-euler
pseudo-scheme =rk34 ;euler
tstart = 0.0
tend =14
dt = 0.00005
pseudo-dt = 0.000005
controller = none
pseudo-niters-min = 1
pseudo-niters-max = 1
pseudo-resid-norm = uniform ;l2
pseudo-resid-tol = 1e-8
pseudo-controller = none ;local-pi
;atol = 1e-6
;safety-fact = 0.9
;min-fact = 0.3 ;0.99
;max-fact = 2.5 ;1.005
;pseudo-dt-max-mult = 2.5
[solver-dual-time-integrator-multip]
pseudo-dt-fact = 2.3
;cycle = [(1, 1)]
;cycle = [(2, 1), (1, 1), (0, 2), (1, 1), (2, 1)]
;cycle = [(3, 1), (2, 1), (1, 1), (0, 2), (1, 1), (2, 1), (3, 3)]
cycle = [(3, 1), (2, 2), (1, 4), (0, 8), (1, 4), (2, 2), (3, 1)]
[solver-interfaces]
riemann-solver = rusanov
ldg-beta = 0.5 ;idk
ldg-tau = 0.1 ;idk
[solver-interfaces-line]
flux-pts = gauss-legendre
quad-deg = 10
quad-pts = gauss-legendre
[solver-interfaces-tri]
flux-pts = williams-shunn
quad-deg = 10
quad-pts = williams-shunn
[solver-interfaces-quad]
flux-pts = gauss-legendre
quad-deg = 10
quad-pts = gauss-legendre
[solver-elements-tri]
soln-pts = williams-shunn
[solver-elements-quad]
soln-pts = gauss-legendre
[solver-elements-hex]
soln-pts = gauss-legendre
quad-deg = 10
quad-pts = gauss-legendre
[solver-elements-tet]
soln-pts = shunn-ham
quad-deg = 10
quad-pts = shunn-ham
[solver-elements-pri]
soln-pts = williams-shunn~gauss-legendre
quad-deg = 10
quad-pts = williams-shunn~gauss-legendre
[solver-elements-pyr]
soln-pts = gauss-legendre
quad-deg = 10
quad-pts = witherden-vincent
[soln-plugin-writer]
dt-out = 0.02
basedir = .
basename = {t:.4f}
[soln-plugin-nancheck]
nsteps = 10
[soln-plugin-pseudostats]
flushsteps = 10
file = residual.csv
header = true
[soln-plugin-fluidforce-upper]
nsteps = 10
file = upper_forces.csv
header = true
morigin = (0.0, 0.0, 0.0)
[soln-plugin-fluidforce-lower]
nsteps = 10
file = lower_forces.csv
header = true
morigin = (0.0, 0.0, 0.0)
[soln-bcs-inlet]
type = ac-char-riem-inv
;rho = 0.0
u = Uin
v = Vin
w = 0
p = Pc
ac-zeta = 2.5 ;3
[soln-bcs-outlet]
type = ac-char-riem-inv
;rho = 0.0
u = Uin
v = Vin
w = 0
p = Pc
ac-zeta = 2.5 ;3
[soln-bcs-upper]
type = no-slp-wall
[soln-bcs-lower]
type = no-slp-wall
[soln-ics]
;rho = rhoc
u = Uin
v = Vin
w = 0
p = Pc
rho = 1.0
```

Also for question3, do you mean span-wise average can be get in Paraview (not pyfr)?