From 48b26323a7541c024804e85001b625214395ea11 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Wed, 11 May 2022 17:57:54 +0200
Subject: [PATCH] Quentin : 11/05/2022: bugfix GRIB2 ARPEGE since Epygram 1.4.8
read the Land-sea mask wrongly (always =1) because a field with the same
parameterNumber and level was read
---
src/SURFEX/mode_read_grib.F90 | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/SURFEX/mode_read_grib.F90 b/src/SURFEX/mode_read_grib.F90
index f92149d19..50ba15ed1 100644
--- a/src/SURFEX/mode_read_grib.F90
+++ b/src/SURFEX/mode_read_grib.F90
@@ -446,6 +446,7 @@ INTEGER :: ILEV ! level
INTEGER :: INUM_ZS,ISIZE,ICOUNT,JLOOP,IPARAM,IGRIB,IPARAM2
CHARACTER(LEN=24) :: YLTYPELU
REAL(KIND=JPRB) :: ZHOOK_HANDLE
+CHARACTER(LEN=50) :: CNAME ! name of the parameter (for ARPEGE GRIB2 converted with Epygram 1.4.8)
!-------------------------------------------------------------------
IF (LHOOK) CALL DR_HOOK('MODE_READ_GRIB:READ_GRIB_LAND_MASK',0,ZHOOK_HANDLE)
WRITE (KLUOUT,'(A)') 'MODE_READ_GRIB:READ_GRIB_LAND_MASK: | Reading land mask from ',HINMODEL
@@ -458,7 +459,8 @@ SELECT CASE (HINMODEL)
CASE ('ARPEGE','ALADIN','MOCAGE')
IF(HINMODEL=='ARPEGE' .AND. NGRIB_VERSION==2) THEN
ILTYPE=1
- CALL READ_GRIB(HGRIB,HINMODEL,KLUOUT,0,IRET,PMASK,KLTYPE=ILTYPE)
+ CNAME='Land-sea mask'
+ CALL READ_GRIB(HGRIB,HINMODEL,KLUOUT,0,IRET,PMASK,KLTYPE=ILTYPE,HNAME=CNAME)
ELSE
CALL READ_GRIB(HGRIB,HINMODEL,KLUOUT,81,IRET,PMASK)
ENDIF
--
GitLab