periodic jet initialization

examples/bubble/periodic_jet_levelset.py

@@ -24,7 +24,7 @@ def init_phi(data, coords, L):
     Di = np.empty_like(phi)
     X,Y    = coords 
     Ly, Lx = L
-    R = Lx/16*(2.5 - 0.75*np.sin(2*np.pi*Y))
+    R = Lx*(0.15 - 0.05*np.sin(2*np.pi*Y))
     Di[...] = (X-0.5*Lx)**2 - R**2
     mask = (np.abs(Di)<np.abs(phi))
     phi[mask] = Di[mask]