Commit: 5e56b77d4446fae9d45bf176488783696dd4a75b - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revision: 5e56b77d4446fae9d45bf176488783696dd4a75b (tree)
• Zeit: 2007-07-04 21:38:16
• Autor: iselllo
• Commiter: iselllo

Log Message

I added a simple python code temperature_comparison.py to read and plot
a couple of quantities
of
interest (as output by several simulations) to see the effect of the
temperature

Ändern Zusammenfassung

• add: Python-codes/temperature_comparison.py

Diff

diff -r 8926a1f71b2c -r 5e56b77d4446 Python-codes/temperature_comparison.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Python-codes/temperature_comparison.py	Wed Jul 04 12:38:16 2007 +0000
@@ -0,0 +1,39 @@
+#! /usr/bin/env python
+from scipy import *
+import pylab  # used to read the .csv file
+
+mu_350=pylab.load("log_mu_evolution")
+mu_400=pylab.load("../D_f-19-T400/log_mu_evolution")
+mu_450=pylab.load("../D_f-19-T450/log_mu_evolution")
+
+print 'mu_350 is', mu_350
+
+pylab.plot(mu_350[:,0],mu_350[:,1],mu_400[:,0],mu_400[:,1],mu_450[:,0],mu_450[:,1])
+pylab.xlabel('time')
+pylab.ylabel('geometric mean')
+pylab.legend(('T=350','T=400',"T=450"))
+pylab.title('Effect of T on mean')
+pylab.grid(True)
+pylab.savefig('log_mu_vs_T')
+
+pylab.hold(False)
+
+
+mu_350=pylab.load("tot_population_evolution")
+mu_400=pylab.load("../D_f-19-T400/tot_population_evolution")
+mu_450=pylab.load("../D_f-19-T450/tot_population_evolution")
+
+print 'now mu_350 is', mu_350
+
+pylab.plot(mu_350[:,0],mu_350[:,1],mu_400[:,0],mu_400[:,1],mu_450[:,0],mu_450[:,1])
+pylab.xlabel('time')
+pylab.ylabel('total particle concentration')
+pylab.legend(('T=350','T=400',"T=450"))
+pylab.title('Effect of T on mean')
+pylab.grid(True)
+pylab.savefig('N_tot_vs_T')
+
+pylab.hold(False)
+
+
+print 'So far so good'