From 06bfd1684b859df720e4890683beef152fa94615 Mon Sep 17 00:00:00 2001
From: Grzegorz Kowal <grzegorz@amuncode.org>
Date: Fri, 26 Nov 2021 11:39:39 -0300
Subject: [PATCH] USER_PROBLEM: Fix the exponential decay in the absoption
 layer.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
---
 sources/user_problem.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/sources/user_problem.F90 b/sources/user_problem.F90
index 6868d74..62a6185 100644
--- a/sources/user_problem.F90
+++ b/sources/user_problem.F90
@@ -1008,7 +1008,7 @@ module user_problem
       dh(3) = adz(pdata%meta%level)
 #endif /* NDIMS == 3 */
 
-      fa = exp(- dt / tabs)
+      fa = 1.0d+00 - exp(- dt / tabs)
       fb = 5.0d-01 * adif * minval(dh(1:NDIMS))**2
 
       call laplace(dh, pdata%q(ivy,:,:,:), q2)