From a39cd6363a3623232cb4e5ddd2efd987d326a93a Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Wed, 30 Jan 2019 15:48:53 -0200 Subject: [PATCH] USER_PROBLEM: Limit blim parameter. Signed-off-by: Grzegorz Kowal --- sources/user_problem.F90 | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index 5687c6a..50054b7 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -99,9 +99,10 @@ module user_problem ! include external procedures and variables ! - use constants , only : pi, pi2 - use parameters , only : get_parameter - use random , only : randomu, randomn + use constants , only : pi, pi2 + use coordinates, only : ng, ady + use parameters , only : get_parameter + use random , only : randomu, randomn ! local variables are not implicit by default ! @@ -187,6 +188,10 @@ module user_problem ! kvec = pi2 * kper +! lower limit for blim +! + blim = max(blim, ng * ady(1)) + ! prepare the random perturbation of velocity ! if (pert == 1) then