diff --git a/src/SURFEX/coupling_tebn.F90 b/src/SURFEX/coupling_tebn.F90 index 3626e37642f404c87e3d89aed12dd4ada8bb2a5c..73615be56125e7b522def1edae61efd13cc68406 100644 --- a/src/SURFEX/coupling_tebn.F90 +++ b/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, &