From 275bd7d2cd192e05e3084188fd8d730016255318 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Tue, 22 Aug 2023 17:24:04 +0200
Subject: [PATCH] Quentin 22/08/2023: bugfix SSO_DIR reading from PREP_PGD +
 PREP_IDEAL_CASE : prep_ideal calls init_pgd_surf_atm twice. In v9.0, XSSO_DIR
 now allocated at read_pgd_isban. At the second call of
 init_pgd_surf_atm>init_surf_atm, XSSO_DIR is already associated on nature
 points and then is not re-allocated at the full size of the domain as
 expected

---
 src/SURFEX/read_sson.F90 | 1 +
 1 file changed, 1 insertion(+)

diff --git a/src/SURFEX/read_sson.F90 b/src/SURFEX/read_sson.F90
index 60d652ba9..cffa00c81 100644
--- a/src/SURFEX/read_sson.F90
+++ b/src/SURFEX/read_sson.F90
@@ -127,6 +127,7 @@ YRECFM='SSO_ANIS'
  CALL READ_SURF(HPROGRAM,YRECFM,USS%XSSO_ANIS(:),IRESP)
 WHERE (PSEA(:) == 1.) USS%XSSO_ANIS(:) = XUNDEF
 !
+IF (ASSOCIATED(USS%XSSO_DIR)) DEALLOCATE(USS%XSSO_DIR) ! For MesoNH PREP_IDEAL_CASE
 IF(.NOT.ASSOCIATED(USS%XSSO_DIR)) ALLOCATE(USS%XSSO_DIR(KSIZE_FULL))
 YRECFM='SSO_DIR'
  CALL READ_SURF(HPROGRAM,YRECFM,USS%XSSO_DIR(:),IRESP)
-- 
GitLab