gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

USER_PROBLEM: Fix the exponential decay in the absoption layer.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2021-11-26 11:39:39 -03:00
parent 8484f52154
commit 06bfd1684b
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
sources/user_problem.F90
View File

@ -1008,7 +1008,7 @@ module user_problem
dh(3) = adz(pdata%meta%level) dh(3) = adz(pdata%meta%level)
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
fa = exp(- dt / tabs) fa = 1.0d+00 - exp(- dt / tabs)
fb = 5.0d-01 * adif * minval(dh(1:NDIMS))**2 fb = 5.0d-01 * adif * minval(dh(1:NDIMS))**2
call laplace(dh, pdata%q(ivy,:,:,:), q2) call laplace(dh, pdata%q(ivy,:,:,:), q2)