Changeset 8274 for trunk/anuga_work/development/2010-projects/anuga_1d/sqpipe/test_sqpipe_domain.py

Timestamp:

Dec 8, 2011, 11:38:18 AM (12 years ago)

Author:

paul

Message:

Fixed loss of mass in sqpipe and rearranged code a little

File:

• trunk/anuga_work/development/2010-projects/anuga_1d/sqpipe/test_sqpipe_domain.py (modified) (1 diff)

trunk/anuga_work/development/2010-projects/anuga_1d/sqpipe/test_sqpipe_domain.py

@@ -5 +5 @@
 
 domain = anuga_1d.sqpipe.parabolic_canal.get_domain(dom)
-finaltime = 10000.0
+finaltime = 1000.0
 yieldstep = 10.0
+#m0 = (domain.get_mass()).sum()
+#for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
+#       m = (domain.get_mass()).sum()
+#       if abs(m - m0) > 0.5:
+#               print m0 - m
+#               
+#       m0 = m
+        
+#outfile="out"
+#anuga_1d.sqpipe.parabolic_canal.write_domain(domain, outfile)
+
 anuga_1d.sqpipe.parabolic_canal.animate_domain(domain, yieldstep, finaltime)
+print "finished"
 anuga_1d.sqpipe.parabolic_canal.plot_domain(domain)