Changeset 5468 for LMDZ6

Timestamp:

Jan 9, 2025, 4:42:44 PM (28 hours ago)

Author:

fhourdin

Message:

New subotpion iflag_pbl_split=20 (or 21).

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1835 +1835 @@
     
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+   CALL iophys_ini(pdtphys,nsurfwind) ! replay automatic include  ! replay automatic include
        CALL wake_ini(rg,rd,rv,prt_level)
        CALL yamada_ini(klon,lunout,prt_level)
@@ -3718 +3719 @@
              !  poches, la tendance moyenne associ\'ee doit etre
              !  multipliee par la fraction surfacique qu'ils couvrent.
+             IF (mod(iflag_pbl_split/10,10) == 1) THEN
+                ! On tient compte du splitting pour modifier les profils deltatq/T des poches
+                DO k=1,klev
+                   DO i=1,klon
+                      d_deltat_the(i,k) = - d_t_ajs(i,k)
+                      d_deltaq_the(i,k) = - d_q_ajs(i,k)
+                   ENDDO
+                ENDDO
+             ELSE
+                d_deltat_the(:,:) = 0.
+                d_deltaq_the(:,:) = 0.
+             ENDIF
+
              DO k=1,klev
                 DO i=1,klon
-                   !
-                   d_deltat_the(i,k) = - d_t_ajs(i,k)
-                   d_deltaq_the(i,k) = - d_q_ajs(i,k)
-                   !
                    d_u_ajs(i,k) = d_u_ajs(i,k)*(1.-wake_s(i))
                    d_v_ajs(i,k) = d_v_ajs(i,k)*(1.-wake_s(i))
                    d_t_ajs(i,k) = d_t_ajs(i,k)*(1.-wake_s(i))
                    d_q_ajs(i,k) = d_q_ajs(i,k)*(1.-wake_s(i))
-                   !
                 ENDDO
              ENDDO