Changeset 4289 for LMDZ6/trunk

Timestamp:

Oct 2, 2022, 10:32:19 PM (2 years ago)

Author:

evignon

Message:

ajoute de la dissipation de tke par les arbres pour l'option yamada4_num=1 (schema numerique du MAR)
plus un petit peu de nettoyage

File:

• LMDZ6/trunk/libf/phylmd/yamada4.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/yamada4.F90

@@ -444 +444 @@
          tkeprov= tkeprov*                           &
            &  (tkeprov+dt*(shear(ig,k)+max(0.,buoy(ig,k))))/ &
-           &  (tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k)))
+           &  (tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k)+drgpro(ig,k)*tkeprov))
          q2(ig,k)=tkeprov*ydeux
         ENDDO
@@ -497 +497 @@
          DO ig=1,ngrid
          tkeprov=q2(ig,k)/ydeux
-
-!             if(ifl_pbltree .eq. 0) then
-!         disseff=dissip(ig,k)-min(0.,buoy(ig,k))
-!         tkeexp=exp(-dt*disseff/tkeprov)
-!         tkeprov= shear(ig,k)*tkeprov/disseff*(1.-tkeexp)+tkeprov*tkeexp
-!           else
 !FC on ajoute la dissipation due aux arbres
          disseff=dissip(ig,k)-min(0.,buoy(ig,k)) + drgpro(ig,k)*tkeprov
@@ -510 +504 @@
          tkeprov= (shear(ig,k)+ &
           & drgpro(ig,k)*(winds(ig,k))**3)*tkeprov/disseff*(1.-tkeexp)+tkeprov*tkeexp
-!               endif
-
          q2(ig,k)=tkeprov*ydeux
-
          ! En cas stable, on traite la flotabilite comme la
          ! dissipation, en supposant que buoy/q2^3 est constant.
@@ -522 +513 @@
        DO k = 2, klev - 1
          DO ig=1,ngrid
-         !tkeprov=q2(ig,k)/ydeux
-         !tkeprov=tkeprov+max(buoy(ig,k)+shear(ig,k),0.)*dt
-         !disseff=dissip(ig,k)-min(0.,buoy(ig,k)+shear(ig,k)) 
-         !tkeexp=exp(-dt*disseff/tkeprov)
-         !tkeprov= tkeprov*tkeexp
-         !q2(ig,k)=tkeprov*ydeux
          ! En cas stable, on traite la flotabilite comme la
          ! dissipation, en supposant que dissipeff/TKE est constant.
          ! Puis on prend la solution exacte
-         ! 
          ! With drag and dissipation from high vegetation (EV & FC, 05/10/2020)
          winds(ig,k)=sqrt(u(ig,k)**2+v(ig,k)**2)