velocity.

pydemo/euler/euler.ufl

@@ -18,6 +18,10 @@ class Compressible2DEuler:
     def maxLambda(U,n):
         rho, u,v, p = Compressible2DEuler.toPrim(U)
         return abs(u*n[0]+v*n[1]) + sqrt(gamma*p/rho)
+    def velocity(U):
+        return as_vector( [U[1],U[2]] )
+#    def physical(U)
+#        return conditional( (U[0]>1e-8) )
 
 Model = Compressible2DEuler