From 57af99f98c9ace8a1835871bb032cef05f8d2ca9 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 6 May 2022 12:20:30 +0200
Subject: [PATCH] M.Taufour: 06/05/2022: correct position of test on small
value of BHMIE IN argument
---
src/MNH/mode_fscatter.f90 | 14 +++++++-------
1 file changed, 7 insertions(+), 7 deletions(-)
diff --git a/src/MNH/mode_fscatter.f90 b/src/MNH/mode_fscatter.f90
index 6c5bdbb67..3e66e3a43 100644
--- a/src/MNH/mode_fscatter.f90
+++ b/src/MNH/mode_fscatter.f90
@@ -153,6 +153,12 @@ CONTAINS
!
! Modification (C.Lac) 04/2014 : exclude very small values of x
!
+ IF (X <= 1.E-07) THEN
+ QEXT = 0.
+ QSCA = 0.
+ QBACK = 0.
+ RETURN
+ ELSE
! -----------------------------------------------------------
! del is the inner sphere convergence criterion
! -----------------------------------------------------------
@@ -169,13 +175,7 @@ CONTAINS
qext = 0.0
xback = (0.0,0.0)
n=1
- IF (x <= 1.E-07) THEN
- QEXT = 0.
- QSCA = 0.
- QBACK = 0.
- RETURN
- ELSE
- do while(n<=nstop)
+ do while(n<=nstop)
! DO n=1,nstop
DX = 1.0/(n/x-dx) - n/x
DY = 1.0/(n/y-dy) - n/y
--
GitLab