diff --git a/@Ternary_model/create_mesh.m b/@Ternary_model/create_mesh.m
index 35884d83ce52d336a2efc9d8ed0bed5cfed22c08..6646f2f05eb1c93f22ad3978f2693d6c7ad671a8 100644
--- a/@Ternary_model/create_mesh.m
+++ b/@Ternary_model/create_mesh.m
@@ -33,5 +33,12 @@ elseif strcmp(T.mode, 'Client') && T.v ~= 0
 %     T.x = linspace(0, T.system_size, 24000);
     T.x = [linspace(0, -T.a-1, 1000), linspace(-T.a-1+0.0001, -T.a+1, 4800),...
            linspace(-T.a+1+0.0001, T.system_size, 1000)];
+    % Fine mesh close to x0 if ic = 'Gauss'
+    [~, ind] = min(abs(T.x-T.x0));
+    [~, ind_plus] = min(abs(T.x-(T.x0+0.6)));
+    [~, ind_minus] = min(abs(T.x-(T.x0-0.6)));
+    T.x = [T.x(1:ind_minus),...
+           linspace(T.x(ind_minus)+0.0001, T.x(ind_plus)-0.0001, 1000),...
+           T.x(ind_plus:end)];
 end
 end