From b6cba63b72580ed5904c002f2efdac102b9d3c4c Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Wed, 11 Nov 2020 11:19:30 +0100
Subject: [PATCH] Real parameters can now be calculated also for square
mobility.
---
@Ternary_model/calc_real_params.m | 11 ++++++++++-
1 file changed, 10 insertions(+), 1 deletion(-)
diff --git a/@Ternary_model/calc_real_params.m b/@Ternary_model/calc_real_params.m
index be348b2..16aebb3 100644
--- a/@Ternary_model/calc_real_params.m
+++ b/@Ternary_model/calc_real_params.m
@@ -3,11 +3,20 @@ function calc_real_params(T)
% partitioning etc..
T.real_params.phi_in = Ternary_model.phi_tot(-1000, T.a, T.b, T.e, T.u0, 0);
T.real_params.phi_out = Ternary_model.phi_tot(1000, T.a, T.b, T.e, T.u0, 0);
+T.real_params.partitioning = T.real_params.phi_in/T.real_params.phi_out;
+
+if strcmp(T.mode, 'Constituent') || strcmp(T.mode, 'Client')
T.real_params.Gamma_in = Ternary_model.gamma0(-1000, T.a, T.b, T.e_g0,...
T.u_g0, 0);
T.real_params.Gamma_out = Ternary_model.gamma0(1000, T.a, T.b, T.e_g0,...
T.u_g0, 0);
T.real_params.D_in = (1-T.real_params.phi_in)*T.real_params.Gamma_in;
T.real_params.D_out = (1-T.real_params.phi_out)*T.real_params.Gamma_out;
-T.real_params.partitioning = T.real_params.phi_in/T.real_params.phi_out;
+
+elseif strcmp(T.mode, 'Const_mob')
+ T.real_params.D_in = (1-T.real_params.phi_in+...
+ T.e_g0*T.real_params.phi_in^2);
+ T.real_params.D_out = (1-T.real_params.phi_out+...
+ T.e_g0*T.real_params.phi_out^2);
+end
end
\ No newline at end of file
--
GitLab