exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Write forces dat file

Browse source
This commit is contained in:
Hrvoje Jasak 2016-08-09 13:25:49 +01:00
parent 3143f456cc
commit 520b77015d
1 changed files with 6 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
bin/plotForces.py
View file

@ -35,7 +35,6 @@ for line in lines:
mvy.append(float(match.group(12))) mvy.append(float(match.group(12)))
mvz.append(float(match.group(13))) mvz.append(float(match.group(13)))
# combine pressure and viscous forces and moments # combine pressure and viscous forces and moments
fx = fpx fx = fpx
fy = fpy fy = fpy
@ -53,6 +52,12 @@ for i in range(1,len(t)):
my[i] += mvy[i] my[i] += mvy[i]
mz[i] += mvz[i] mz[i] += mvz[i]
# write clean data file
outfile=open('forces.dat','w')
for data in zip(t,fx,fy,fz):
outfile.write(' '.join([str(d) for d in data])+'\n')
# plot forces # plot forces
import pylab import pylab
pylab.xlabel('iteration') pylab.xlabel('iteration')