Index: /issm/workshop/2012/Talks/06_Mesh_generation/Code/circles.m
===================================================================
--- /issm/workshop/2012/Talks/06_Mesh_generation/Code/circles.m	(revision 14122)
+++ /issm/workshop/2012/Talks/06_Mesh_generation/Code/circles.m	(revision 14122)
@@ -0,0 +1,7 @@
+function vel=circles(x,y),
+
+	u=4*x-2; v=4*y-2;
+
+	vel=tanh(30*(u.^2+v.^2-0.25)) ...
+		+tanh(30*((u-0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u-0.75).^2+(v+0.75).^2-0.25)) ...
+		+tanh(30*((u+0.75).^2+(v-0.75).^2-0.25)) +tanh(30*((u+0.75).^2+(v+0.75).^2-0.25)) ;