From ae3319ce75c7cd53cf1443153b8114067cd70174 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 24 Jul 2020 12:10:39 +0200
Subject: [PATCH] Marie M. 24/07/2020: SURFEX bugfix : correction of IF
condition impacting the interpolation of tiles fractions at PGD step (when
fractions are prescribed in a file and when some points don't have data) from
commit f70f6454 of SURFEX repository 29/04/2019
---
src/SURFEX/pgd_fieldin.F90 | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/SURFEX/pgd_fieldin.F90 b/src/SURFEX/pgd_fieldin.F90
index 683611a1c..990727e15 100644
--- a/src/SURFEX/pgd_fieldin.F90
+++ b/src/SURFEX/pgd_fieldin.F90
@@ -57,7 +57,7 @@ USE MODD_PGDWORK, ONLY : XALL, NSIZE_ALL, CATYPE, NSIZE, XSUMVAL, &
USE MODD_SURF_PAR, ONLY : XUNDEF
USE MODD_PGD_GRID, ONLY : NL
!
-USE MODD_DATA_COVER_PAR, ONLY : NTYPE, LVEG_PRES, NVEGTYPE
+USE MODD_DATA_COVER_PAR, ONLY : NTYPE, LVEG_PRES, NVEGTYPE, NVEGTYPE_OLD
!
USE MODI_GET_LUOUT
USE MODI_TREAT_FIELD
@@ -293,8 +293,8 @@ IF (LEN_TRIM(HFILE)/=0) THEN
!
DO JT=1,SIZE(NSIZE,2)
- IF (.NOT.U%LECOSG.AND.JT>NVEGTYPE) EXIT
-
+ IF (.NOT.U%LECOSG.AND.JT>NVEGTYPE_OLD) EXIT
+
!multitype input file
IF (SIZE(ZFIELD,2)>1) THEN
--
GitLab