Commit: 08a73f562b8ef9ea912481f8cc3ed6658967ca6a - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revision: 08a73f562b8ef9ea912481f8cc3ed6658967ca6a (tree)
• Zeit: 2007-03-07 19:57:33
• Autor: iselllo
• Commiter: iselllo

Log Message

I modified the code Python-codes/postprocess.py by adding some extra
functionalities to read & postprocess the output of the DNS simulations.

Ändern Zusammenfassung

• modified: Python-codes/postprocess.py

Diff

diff -r 896c4a6f209d -r 08a73f562b8e Python-codes/postprocess.py
--- a/Python-codes/postprocess.py	Mon Mar 05 17:55:03 2007 +0000
+++ b/Python-codes/postprocess.py	Wed Mar 07 10:57:33 2007 +0000
@@ -2,13 +2,13 @@
 from scipy import *
 import pylab  # used to read the .csv file
 
-from pylab import *
+#from pylab import *
 
 # now I want to try reading some .csv file
 data = pylab.load("radjet.out")
 #print "data now is", data
 
-hold(False) # so that the plots are NOT superimposed
+pylab.hold(False) # so that the plots are NOT superimposed
 
 skip_this=0
 if (skip_this !=1):
@@ -25,15 +25,14 @@
 
 
 x=linspace(0.,1.,len(data[:,0]))
-print 'the length of data[:,0] is', len(data[:,0])
-print 'the length of x is', len(x)
+
 
 
 rext=5.
 strr=3.5 # stretching parameter
 istr_0=tanh(x*strr)/tanh(strr)*rext
 
-print 'istr_0 is', istr_0
+#print 'istr_0 is', istr_0
 
 print 'the length of istr_0 is', len(istr_0)
 
@@ -47,13 +46,7 @@
 pylab.savefig('radial_physical_grid_istrr=0')
 #pylab.show() # uncomment them to plot a distribution
 
-#s = file( './field0160.dat','rb' ).read( )
 
-#newarr = fromstring(s ,Float)
-
-#print 'newarr is', newarr
-
-#print 'the dimensions of newarr are', shape(newarr)
 
 data2 = pylab.load("rg2.out")
 
@@ -65,7 +58,7 @@
 
 pylab.plot(seq2,data2[:,0],seq2,data2[:,1],seq2,data2[:,2],seq2,data2[:,3],linewidth=2.0)
 #hold(True) # now I want to overlap the plots
-legend(('rc','rm','g2rc','g2rm'))
+pylab.legend(('rc','rm','g2rc','g2rm'))
 pylab.xlabel('grid point number')
 pylab.ylabel('physical coordinate')
 #pylab.title('Radial Grid ISTR=0')
@@ -89,7 +82,7 @@
 
 pylab.plot(seq2,data2[:,0],seq2,data2[:,1],seq2,data2[:,2],linewidth=2.0)
 #hold(True) # now I want to overlap the plots
-legend(('x2','yd','rc'))
+pylab.legend(('x2','yd','rc'))
 pylab.xlabel('grid point number')
 pylab.ylabel('physical coordinate')
 #pylab.title('Radial Grid ISTR=0')
@@ -98,3 +91,103 @@
 #pylab.plot(seq2,data2[:,2])
 #pylab.plot(seq2,data2[:,3])
 pylab.savefig('radstrout_plot')
+
+
+# now I try reading the binary file
+
+#s = file( './field0160.dat','rb' ).read( )
+
+#newarr = fromstring(s ,Float)
+
+#print 'newarr is', newarr
+
+#print 'the dimensions of newarr are', shape(newarr)
+
+#lenarr=len(newarr)
+#print max(newarr)
+
+pylab.hold(False)
+
+v_prof = pylab.load("v_prof")
+
+pylab.plot(data2[:,2],v_prof)
+pylab.xlabel('radial distance from axis')
+pylab.ylabel('Velocity profile')
+pylab.grid(True)
+pylab.savefig('velocity_profile')
+
+
+
+my_mean=zeros(len(data2[:,0]))
+
+
+pylab.hold(True) # s e t t h e ho ld s t a t e t o be on
+for i in range (1,20 ) : # s t a r t a t 1 , end a t 100
+	fname = 'vpro%03d'%i # %03d pads t h e i n t e g e r s wi th z e r o s
+	new_v_prof = pylab.load(fname)
+	pylab.plot(data2[:,2],new_v_prof,linewidth=1.5)
+	pylab.xlabel('radial distance from axis')
+	pylab.ylabel('Velocity profile')
+	pylab.grid(True)
+new_v_prof = pylab.load('vpro001')
+pylab.plot(data2[:,2],new_v_prof,linewidth=1.5)
+pylab.savefig('velocity_profile_different_times')
+#pylab.show()
+
+pylab.hold(False)
+
+ini_read=15
+end_read=20
+
+#hold(True) # s e t t h e ho ld s t a t e t o be on
+for i in range (ini_read,(end_read+1) ) : # s t a r t a t 1 , end a t 100
+	#print 'i is', i
+	fname = 'vpro%03d'%i # %03d pads t h e i n t e g e r s wi th z e r o s
+	new_v_prof = pylab.load(fname)
+	my_mean=my_mean+new_v_prof
+	pylab.plot(data2[:,2],new_v_prof,linewidth=1.5)
+	pylab.hold(True)
+	pylab.xlabel('radial distance from axis')
+	pylab.ylabel('Velocity profile')
+	pylab.grid(True)
+
+my_mean=my_mean/(end_read-ini_read+1)
+#print 'mymat[] is', mymat[:,1]
+
+ini_v_prof = pylab.load('vpro001')
+pylab.plot(data2[:,2],ini_v_prof,linewidth=1.5)
+
+
+
+
+
+
+pylab.legend(('t=150','t=160','t=170','t=180','t=190','t=200','initial profile'))
+pylab.savefig('velocity_profile_late_times')
+#pylab.show()
+
+pylab.hold(False)
+
+
+pylab.plot(data2[:,2],my_mean,data2[:,2],ini_v_prof,linewidth=1.5,)
+pylab.legend(('average profile at late t','initial profile'))
+pylab.savefig('late-average-vs-initial')
+
+# now I also try plotting the initial profile I would get using an empirical correlation
+
+def v_appr(r):
+	return U*(1-r/R)**(1./6.)
+
+U=my_mean[0]
+R=1.
+v_emp=map(v_appr,data2[:,2])
+pylab.hold(False)
+
+pylab.plot(data2[:,2],my_mean,data2[:,2],v_emp,linewidth=1.5,)
+pylab.legend(('average profile at late t','empirical_profile'))
+pylab.grid(True)
+pylab.xlabel('radial distance from axis')
+pylab.ylabel('Velocity profile')
+pylab.savefig('late-average-vs-empirical')
+
+print 'So far so good' 
\ No newline at end of file