parameters with/without morphogen

fixed_boundaries.json

+{
+    "morphogen" : true,
+
+    "resolution" : 64,
+    "system_size_by_2PI" : 1 ,
+
+    "elasticity" : 0.0,
+    "viscosity" : 1.0,
+    "friction" : 1.0,
+    "lamda": 1.75,
+
+    "diffusion_rho" : 1.0,
+    "turnover_rho" : 0.0,
+    "average_rho" : 1.0,
+    "saturation_rho" : 1.0,
+
+    "diffusion_c" : 1.0,
+    "turnover_c" : 0.0,
+    "average_c" : 1.0,
+
+    "noise_level": 0.1,
+
+    "timestep" : 0.05,
+    "maxtime" : 200.0
+
+}

no_c_fixed_bdry_fft.py

@@ -10,13 +10,13 @@ import progressbar
 Lx = 2*np.pi
 Lx = 2*np.pi
 Nx = 200
 Nx = 200
 T = 100
 T = 100
-dt = 0.1
+dt = 0.01
 Nt = int(T/dt)
 Nt = int(T/dt)
 times = np.linspace(0,T,Nt)
 times = np.linspace(0,T,Nt)
 
 
 K = 1
 K = 1
 eta = 1
 eta = 1
-lamda = 1
+lamda = 5.0
 krho = 1
 krho = 1
 rho0 = 1
 rho0 = 1
 Drho = 1
 Drho = 1

with_c_lsa.py

@@ -19,13 +19,13 @@ ax.set_xlabel(r'$q$')
 ax.set_ylabel(r'$Re[\, \lambda(q) \, ]_{\rm max}$')
 ax.set_ylabel(r'$Re[\, \lambda(q) \, ]_{\rm max}$')
 ax.axhline(y=0, color='black')
 ax.axhline(y=0, color='black')
 
 
-lamda0 = 6.0
-K0 = 4.0
+lamda0 = 2.6
+K0 = 1.0
 krho0 = 1.0
 krho0 = 1.0
-eta0 = 4.0
+eta0 = 1.0
 rho00 = 1.0
 rho00 = 1.0
 Dc0 = 1.0
 Dc0 = 1.0
-Drho0 = 0.0
+Drho0 = 1.0
 kc0 = 1.0
 kc0 = 1.0
 c00 = 1.0
 c00 = 1.0
 gamma0 = 1.0
 gamma0 = 1.0