Jean W. 18/07/2024: SURFEX fix to avoid too violent spin-up in Post-treatment...

Jean W. 18/07/2024: SURFEX fix to avoid too violent spin-up in Post-treatment on fluxes to avoid extra-large positive latent and negative sensible fluxes

src/SURFEX/isba_hveg.F90

@@ -399,6 +399,11 @@ WHERE (PLE_HVEG>PRN_HVEG .AND. PEK%XTV<PTA-1.)
   ZDELTA_H = PRHOA * XCPD * ZVEG * (ZDT ) / PEK%XRESA / PEXNS
   ZDELTA_G = ZVEG * ZCHEATV * ZDT / PTSTEP
   ZDELTA_LE= - ZDELTA_H - ZDELTA_G
+  !
+  ! J. Wurtz & V. Masson Check if ZDELTA_LE > 0 to avoid too violent spin-up
+  ! Residual is increased of remaining energy
+  ZDELTA_LE= MAX(ZDELTA_LE, -PLE_HVEG)
+  ZDELTA_G= - ZDELTA_H - ZDELTA_LE
 END WHERE
 !
 PH_HVEG (:) = PH_HVEG (:) + ZDELTA_H (:)