Index: /issm/trunk-jpl/examples/Mesh/circles.py
===================================================================
--- /issm/trunk-jpl/examples/Mesh/circles.py	(revision 26562)
+++ /issm/trunk-jpl/examples/Mesh/circles.py	(revision 26562)
@@ -0,0 +1,11 @@
+import numpy as np
+
+
+def circles(x, y):
+
+    u = 4 * x - 2
+    v = 4 * y - 2
+
+    vel = np.tanh(30 * (u**2 + v**2 - 0.25)) + np.tanh(30 * ((u - 0.75)**2 + (v - 0.75)**2 - 0.25)) + np.tanh(30 * ((u - 0.75)**2 + (v + 0.75)**2 - 0.25)) + np.tanh(30 * ((u + 0.75)**2 + (v - 0.75)**2 - 0.25)) + np.tanh(30 * ((u + 0.75)**2 + (v + 0.75)**2 - 0.25))
+
+    return vel
Index: /issm/trunk-jpl/examples/Mesh/shock.py
===================================================================
--- /issm/trunk-jpl/examples/Mesh/shock.py	(revision 26562)
+++ /issm/trunk-jpl/examples/Mesh/shock.py	(revision 26562)
@@ -0,0 +1,7 @@
+import numpy as np
+
+
+def shock(x, y):
+
+    vel = np.exp(-(np.sqrt((x + 0.1)**2 + (y + 0.1)**2) - 0.75)**2 * 1e6) + ((x + 0.1)**2 + (y + 0.1)**2) / 2
+    return vel