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

USER_PROBLEM: Fix amplitude calculation for one-mode perturbations.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2024-04-21 12:07:08 -03:00
parent 9e45431c89
commit da0cd1b34e
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

@ -674,7 +674,7 @@ module user_problem
! prepare decaying factors
!
fv = 0.5d+00 * pi
fv = 5.0d-01 * pi
do i = 1, nn
xp = fv * min(1.0d+00, max(0.0d+00, abs(xc(i)) - xcut) / xdec)
fx(i) = cos(xp)**2
@ -820,7 +820,7 @@ module user_problem
! calculate the perturbation factor
!
fv = vper * ylen / (pi2 * pi2 * pi * dh(1) * dh(2) * kper)
fv = vper * ylen / (pi2 * pi2 * pi * dh(1) * dh(2) * kper * kper)
! calculate the perturbation profile
!
@ -863,7 +863,7 @@ module user_problem
! calculate the perturbation factor
!
fv = bper * ylen / (pi2 * pi2 * pi * dh(1) * dh(2) * kper)
fv = bper * ylen / (pi2 * pi2 * pi * dh(1) * dh(2) * kper * kper)
! calculate the perturbation profile
!