Changeset 6068 for LMDZ6/trunk/libf/phylmd/lmdz_call_gwd.f90

Timestamp:

Feb 5, 2026, 4:11:56 PM (2 weeks ago)

Author:

evignon

Message:

update of mesh-averaged TKE after tend_to_tke. This only impact the writting of restartphy and ener_conserv

File:

• LMDZ6/trunk/libf/phylmd/lmdz_call_gwd.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_call_gwd.f90

@@ -20 +20 @@
 !===================================================================
 
-   SUBROUTINE call_gwd(klon, klev, nbsrf, is_ter, is_lic, abortphy, flag_inhib_tend, itap, JD_cur, JD_ref, JH_cur, &
+   SUBROUTINE call_gwd(klon, klev, nbsrf, is_ter, is_lic, is_ave, abortphy, flag_inhib_tend, itap, JD_cur, JD_ref, JH_cur, &
                        pctsrf, is_sequential, phys_tstep, cell_area, longitude_deg, latitude_deg, pphis, &
                        zstd, zpic, zmea, zval, zsig, zgam, zthe, pplay, paprs, presnivs, rain_fall, snow_fall, &
@@ -61 +61 @@
       REAL, INTENT(IN)     :: JD_ref ! start day of the run
       REAL, INTENT(IN)     :: JH_cur ! time of the day in seconds
-      INTEGER, INTENT(IN)  :: is_ter, is_lic ! indices for land and landice subsurfaces
+      INTEGER, INTENT(IN)  :: is_ter, is_lic, is_ave ! indices for land and landice subsurfaces and mesh-averaged
       LOGICAL, INTENT(IN)  :: is_sequential ! sequential or parallel model
       REAL, INTENT(IN)  :: phys_tstep ! time step [s]
@@ -488 +488 @@
          forall (k=1:klev) exner(:, k) = (pplay(:, k)/paprs(:, 1))**rkappa
 
-         CALL tend_to_tke(phys_tstep, klon, klev, nbsrf, paprs, exner, t_seri, u_seri, v_seri, dtadd, duadd, dvadd, pctsrf, tke)
+         CALL tend_to_tke(phys_tstep, klon, klev, nbsrf, is_ave, paprs, exner, t_seri, u_seri, v_seri, dtadd, duadd, dvadd, pctsrf, tke)
 
          ! Prevent pbl_tke_w from becoming negative