Jean W. 18/07/2024: SURFEX add init of old pressure to actual pressure....

Jean W. 18/07/2024: SURFEX add init of old pressure to actual pressure. Previous pressure is used in BEM but is default is wrong leading to possible crash at first timestep.

src/SURFEX/coupling_tebn.F90

@@ -1061,6 +1061,14 @@ DO JP = 1,TOP%NTEB_PATCH
      CALL ABOR1_SFX("COUPLING_TEBN: Negative humidity in canyon")
   ENDIF
   !
+    !CHECK if old pressure is initialized
+  
+  IF (ANY(NB%AL(JP)%XPSOLD(:)==XUNDEF)) THEN
+  NB%AL(JP)%XPSOLD(:)=PPS(:)
+  
+  END IF
+  
+  
   CALL TOWN_ENERGY_BALANCE(DTCO, G, TOP, SPAOP, NT%AL(JP), BOP, NB%AL(JP), TPN, TIR, TD%NDMT%AL(JP), GDM, GRM, &
        HM, SB, CT, JP, HPROGRAM, CIMPLICIT_WIND, PTSUN, ZT_CAN, ZQ_CAN, ZU_CANYON, ZT_LOWCAN,  &
        ZQ_LOWCAN, ZU_LOWCAN, ZZ_LOWCAN, ZTA_HVEG, ZQA_HVEG,                                    &