Two different mobilities and space-dependent reaction rate in

Hammerhead/Ribozyme system.

flory_hugg_pde.m

@@ -27,12 +27,14 @@ elseif strcmp(mode, 'Mobility_ratio')
 elseif strcmp(mode, 'Hammer_Ribo')
     % Use code from client (flux) and constituent (mobility)
     c = [1; 1];
-    g0 = Ternary_model.gamma0(x, a+t*v, b, e_g0, u_g0, v*t, e, u0);
-    chi_phi = -4;
-    chi_phi1 = -10;
+    g01 = Ternary_model.gamma0(x, a+t*v, b, e_g0, u_g0, v*t, e, u0);
+    g02 = Ternary_model.gamma0(x, a+t*v, b, e_g0, 30*u_g0, v*t, e, u0);
+    chi_phi = -1;
+    chi_phi1 = -2;
     k_H = 0.2;
-    f = [1;1].*g0.*(dudx+[chi_phi; chi_phi1].*u*gra_a);
-    F = k_H*pt*u(1);
+    f = [g01;g02].*(dudx+[chi_phi; chi_phi1].*u*gra_a);
+    reaction_rate = Ternary_model.phi_tot(x, a, b, 0.499, u0, v*t);
+    F = k_H*reaction_rate*u(1);
     s = [-F; F];
 end
 end
\ No newline at end of file