USER_PROBLEM: Limit blim parameter.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2019-01-30 15:48:53 -02:00
parent 6087e1b833
commit a39cd6363a

View File

@ -100,6 +100,7 @@ module user_problem
! include external procedures and variables ! include external procedures and variables
! !
use constants , only : pi, pi2 use constants , only : pi, pi2
use coordinates, only : ng, ady
use parameters , only : get_parameter use parameters , only : get_parameter
use random , only : randomu, randomn use random , only : randomu, randomn
@ -187,6 +188,10 @@ module user_problem
! !
kvec = pi2 * kper kvec = pi2 * kper
! lower limit for blim
!
blim = max(blim, ng * ady(1))
! prepare the random perturbation of velocity ! prepare the random perturbation of velocity
! !
if (pert == 1) then if (pert == 1) then