From ee3743dffbb0104dd22bcb40d833e9af09aed917 Mon Sep 17 00:00:00 2001
From: Gaelle DELAUTIER <gaelle.delautier@meteo.fr>
Date: Mon, 18 Sep 2017 10:21:27 +0200
Subject: [PATCH] S.Faroux 09/2017 : correction bug PGD
---
src/SURFEX/interpol_npts.F90 | 27 ++++++++++++++-------------
1 file changed, 14 insertions(+), 13 deletions(-)
diff --git a/src/SURFEX/interpol_npts.F90 b/src/SURFEX/interpol_npts.F90
index d368600a8..bc63cedc9 100644
--- a/src/SURFEX/interpol_npts.F90
+++ b/src/SURFEX/interpol_npts.F90
@@ -164,6 +164,7 @@ DO JL=1,IL
END IF
INPTS = KNPTS
ISCAN = ISCAN_ALL
+ ICOUNT = ISCAN
IINDEX(:) = IINDEX_ALL(:)
!
ELSE
@@ -178,19 +179,19 @@ DO JL=1,IL
END IF
END DO
!
- !IF (ICOUNT>=1) THEN
- IF (ICOUNT>=KNPTS) THEN
- ISCAN = ICOUNT
- !INPTS = MIN(ICOUNT,KNPTS)
- INPTS = KNPTS
- ELSEIF (KNEAR_NBR>=U%NDIM_FULL .AND. ICOUNT>=1) THEN
- ISCAN = ICOUNT
- INPTS = ICOUNT
- ELSE
- KCODE(JL) = -4
- CYCLE
- END IF
- !
+ ENDIF
+ !
+ !IF (ICOUNT>=1) THEN
+ IF (ICOUNT>=KNPTS) THEN
+ ISCAN = ICOUNT
+ !INPTS = MIN(ICOUNT,KNPTS)
+ INPTS = KNPTS
+ ELSEIF (KNEAR_NBR>=U%NDIM_FULL .AND. ICOUNT>=1) THEN
+ ISCAN = ICOUNT
+ INPTS = ICOUNT
+ ELSE
+ KCODE(JL) = -4
+ CYCLE
END IF
!
DO JS=1,ISCAN
--
GitLab