Hi,
I am trying to get a text file containing the time history of quantities related to a single ball (eg force, coordinates,velocity).
What is the best way of doing this?
For resultant force on the ball and ball coordinates I have used the following:
variable z equal z[1]
variable fz equal fz[1]
fix extra all print 1 "$z ${fz}" file Myfile.txt
I do not think it is possible to do something similar with omegax,omegay, omegaz so is there any other way of doing this? I checked and the variable.cpp would need extending to allow the user to specify omega's too. I think this is a simple extension so Christoph, I will try to do this (if there is no other way of getting omegas) and send you one or two extra code lines so that this maybe could be added to the new version?
Thanks for all the help,
George
ckloss | Thu, 09/02/2010 - 18:27
Re: outputting angular velocity of a single grain
Hi George,
yes it should be simple to add that.
However, if you do need it only for post-processing, there is a possibility to do this already:
+ Use dump custom to dump omega(x)/(y)/(z)
In pizza.py type:
+ d=dump("dump.file")
+ omegaxx = d.atom(1,"omegax")
+ file = open('out.omegax','w')
+ print >> file, str(omegaxx.replace('[',' ').replace('[',' ')
+ file.close()
You can then plot it directly in octave:
+ load out.omegax
+ plot(out.omegax)