diff --git a/src/MNH/mode_fscatter.f90 b/src/MNH/mode_fscatter.f90 index 6c5bdbb670eba3e73137a8e01323136edfa1b4f9..3e66e3a4394b708b49acc9dbe6c7c4d65a37ef1c 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