diff --git a/prob_laplace.m b/prob_laplace.m
index 1d23534442660391c2aef6a4eed68e5b8ef29483..fc7e071c215d9e0485c778494c7b5fa6e76ef5b0 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -307,12 +307,12 @@ s = '~/Nextcloud/Langevin_vs_MeanField/Data_Figs_FokkPla/';
 %                                   norm_fac*T_mov.sol(end, :)'])
 %% Flux and entropy change for moving boundary
 chi_phi = -4.530864768482371;
+x_interp = (T_mov.x(1:end-1)+T_mov.x(2:end))/2;
+g0 = 0.5;
 s_dot = zeros(1, length(T_mov.t));
 for i = 1:length(T_mov.t)
-    x_interp = (T_mov.x(1:end-1)+T_mov.x(2:end))/2;
     u = norm_fac*T_mov.sol(i, :);
     u_interp = (u(1:end-1)+u(2:end))/2;
-    g0 = 0.5;
     gra_a = Ternary_model.gradient_analytical(x_interp, T_mov.a, T_mov.b,...
                                               T_mov.e, T_mov.v*T_mov.t(i));
     dudx = diff(u)./diff(T_mov.x);