Reading the statistics

[dineshkumarreddyg36]

Can someone tell me how to load the statistic files that are being written on giving "initstat" some value?

[mathrack]

Hello,

https://github.com/fschuch/xcompact3d_toolbox/
or

#! /usr/bin/env python3
#
import numpy as np
import argparse
#
parser = argparse.ArgumentParser(description='Print min / max and average in space.')
parser.add_argument('files', nargs='+', help='list of .bin files')
args = parser.parse_args()
#
# Incompact3d grid
#
nx = 128
ny = 65
nz = 64
f = open("yp.dat")
yp = np.fromfile(f, sep='\n')
f.close()
#
for binfile in args.files:
    # Read 3D binary file
    f = open(binfile, 'rb')
    data = np.fromfile(f) #, '<f4')
    f.close()
    # Print the file name and min / max
    print("File : " + binfile)
    print("   Min : " + str(np.min(data)))
    print("   Max : " + str(np.max(data)))
    # Map data in a 3D array
    array = np.reshape(data, (nx,ny,nz), order='F')
    # Compute spatial average and write a text file
    f = open(binfile+'.txt', 'w')
    for j in range(ny):
        val=np.sum(array[:,j,:])/np.float(nx*nz)
        f.write(str('{:.6e}'.format(yp[j])) + " " + str('{:.6e}'.format(val)) + '\n')
    f.close()

[pbartholomew08]

It should be writing out to separate files like ux, etc. However it appears to have a bug

[pbartholomew08]

#166

[pbartholomew08]

A fix is now available on pull request https://github.com/xcompact3d/Incompact3d/tree/166-turbulence-statistics-are-not-written-out-properly

Using this you should be able to treat the eddy viscosity files in the same way as ux, uy, etc. Although for visualisation you will need to write snapshot.xdmf files (you can adapt those written for the main flow variables for this).

[pbartholomew08]

This has now been merged, if you update to latest you should be able to output turbulent viscosity data per timestep

[dineshkumarreddyg36]

It's working now.
Thank you