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

USER_PROBLEM: Fix a mistake in the absoption layer profile.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2021-11-05 23:17:24 -03:00
parent 05f46eeb7e
commit 9b7f2f1044
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
sources/user_problem.F90
View File

@ -1111,9 +1111,9 @@ module user_problem
fc(:) = (abs(yc(:)) - yabs) / adec fc(:) = (abs(yc(:)) - yabs) / adec
do j = 1, nn do j = 1, nn
if (fc(j) >= 0.0d+00) then if (fc(j) > 0.0d+00) then
if (fc(j) <= 1.0d+00) then if (fc(j) < 1.0d+00) then
fr = fa * sin(5.0d-01 * pi * yc(j))**2 fr = fa * sin(5.0d-01 * pi * fc(j))**2
else else
fr = fa fr = fa
end if end if