提交: 74c66bd615992ab78b2f0d66fa170824b0d86a75 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修订版: 74c66bd615992ab78b2f0d66fa170824b0d86a75 (tree)
• 时间: 2007-03-10 00:14:57
• 作者: iselllo
• Commiter: iselllo

Log Message

I added the code postprocess2.py, which is able to work out the average
axial velocity along the whole tube as a function of r (radial
coordinate) on a single snapshot of the velocity field.

更改概述

• add: Python-codes/postprocess2.py

差异

diff -r 08a73f562b8e -r 74c66bd61599 Python-codes/postprocess2.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Python-codes/postprocess2.py	Fri Mar 09 15:14:57 2007 +0000
@@ -0,0 +1,50 @@
+#! /usr/bin/env python
+from scipy import *
+import pylab  # used to read the .csv file
+
+q3arr=pylab.load("vpro001") # loading
+
+print 'the dimension of q3arr is', shape(q3arr)
+
+n1=65 # theta grid
+n2=65 # r grid
+n3=49 # z grid
+
+my_mean=zeros(n2)
+
+q3arr=reshape(q3arr,[n1,n2,n3])
+print 'the new dimension of q3arr is', shape(q3arr)
+
+print 'now a print of the array content', q3arr[14,:,12]
+
+z_pos=25 # I fix a certain position along z
+
+
+
+for i in range(1,n1):
+	for k in range(1,n3):
+	#my_mean[:]= q3arr[i,:,z_pos]+my_mean[:]
+	# what is above corresponds to the mean radial profile along z for a certain snapshot of the velocity field
+		my_mean[:]= q3arr[i,:,k]+my_mean[:]
+# now this corresponds to the mean z velocity for a certain snapshot
+#my_mean=my_mean/float(n1)
+my_mean=my_mean/float(n1*n3)
+
+print "my_mean is", my_mean # this is the mean v_z(r) averaged on the cross section for a single velocity  profile for a certain specific z coordinate and a single snapshot of the velocity.
+print 'the dimension of my_mean is', shape(my_mean)
+
+
+data = pylab.load("radstr.out")
+
+
+
+pylab.plot(data[:,0],my_mean)
+pylab.xlabel('Radial Coordinate')
+pylab.ylabel('Physical Coordinate')
+pylab.title('Radial Grid')
+pylab.grid(True)
+pylab.savefig('Average_v_z(r)_on_snapshot')
+
+
+
+print 'So far so good'
\ No newline at end of file