Changeset 6063

Timestamp:

Feb 2, 2026, 12:15:47 PM (2 days ago)

Author:

evignon

Message:

correction nsrf en nbsrf

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_call_gwd.f90 (modified) (1 diff)
• lmdz_gwd_tendtotke.f90 (modified) (4 diffs)
• physiq_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/lmdz_call_gwd.f90

@@ -488 +488 @@
          forall (k=1:klev) exner(:, k) = (pplay(:, k)/paprs(:, 1))**rkappa
 
-   CALL tend_to_tke(phys_tstep, klon, klev, nbsrf, RG, RCPD, paprs, exner, t_seri, u_seri, v_seri, dtadd, duadd, dvadd, pctsrf, tke)
+         CALL tend_to_tke(phys_tstep, klon, klev, nbsrf, paprs, exner, t_seri, u_seri, v_seri, dtadd, duadd, dvadd, pctsrf, tke)
 
          ! Prevent pbl_tke_w from becoming negative

LMDZ6/trunk/libf/phylmd/lmdz_gwd_tendtotke.f90

@@ -39 +39 @@
   CONTAINS
 
- SUBROUTINE tend_to_tke(dt,klon,klev,nbsrf,rg,rcpd,plev,exner,temp,windu,windv,dt_a,du_a,dv_a,pctsrf,tke)
+ SUBROUTINE tend_to_tke(dt,klon,klev,nbsrf,plev,exner,temp,windu,windv,dt_a,du_a,dv_a,pctsrf,tke)
 
+USE lmdz_gwd_ini, ONLY: RG, RCPD
 
 IMPLICIT NONE
@@ -54 +55 @@
   INTEGER, INTENT(IN)  :: klon,klev ! horizontal and vertical dimensions
   INTEGER, INTENT(IN)  :: nbsrf ! number of subsurfaces
-  REAL, INTENT(IN) :: RG ! gravity acceleration [m/s2]
-  REAL, INTENT(IN) :: RCPD ! dry air cp [J/K/kg]
   REAL, DIMENSION(klon,klev+1), INTENT(IN) :: plev  ! inter-layer pressure [Pa]
   REAL, DIMENSION(klon,klev), INTENT(IN) :: temp      ! temperature [K], grid-cell average or for a one subsurface
@@ -132 +131 @@
 
  DO isrf=1,nbsrf
-    DO k=1,klev
+   DO k=1,klev
        DO i=1,klon
           IF (pctsrf(i,isrf)>0.) THEN
@@ -142 +141 @@
  ENDDO
 
+! The tendency on mean TKE (nbsrf+1) must be added as well
+! or mesh-averaged TKE should be updated here
+
+
 
  END SUBROUTINE tend_to_tke

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -4724 +4724 @@
     !
 
-    CALL call_gwd(klon,klev,nsrf,is_ter,is_lic,abortphy,flag_inhib_tend,itap, JD_cur, JD_ref, JH_cur, &
+    CALL call_gwd(klon,klev,nbsrf,is_ter,is_lic,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, &