From 6680ed594f6939cb1b410b1dac9458f539878c1f Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Mon, 9 Aug 2021 16:55:35 +0200
Subject: [PATCH] Starting work on dip in mobility. Not finished.
---
@Ternary_model/gamma0.m | 12 +++++++++++-
1 file changed, 11 insertions(+), 1 deletion(-)
diff --git a/@Ternary_model/gamma0.m b/@Ternary_model/gamma0.m
index c5bb0ff..a1d0456 100644
--- a/@Ternary_model/gamma0.m
+++ b/@Ternary_model/gamma0.m
@@ -3,7 +3,17 @@ function g0 = gamma0(x, a, b, e_g0, u_g0, vt, e, u0)
% Based on a tanh profile.
if isnumeric(e_g0)
g0 = e_g0*(tanh((x+a+vt)/b)+1)/2+u_g0;
-else
+elseif strcmp(e_g0, 'const')
g0 = u_g0./(1-Ternary_model.phi_tot(x, a, b, e, u0, vt));
+elseif ~isempty(str2double(e_g0))
+ eg = str2double(e_g0);
+ al = 2.7;
+ p1 = eg*(tanh((x+a+vt)/b)+1)/2+u_g0;
+ p2 = eg/2/b*(1-tanh((x+a+vt)/b).^2)/al;
+ g0 = p1-p2;
+% figure; hold on;
+% plot(x, p1); plot(x, p2); plot(x, g0);
+% pause();
+% close all;
end
end
\ No newline at end of file
--
GitLab