Changeset 2952 for LMDZ5/trunk/libf/phylmd/yamada4.F90

Timestamp:

Jul 14, 2017, 10:23:40 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Parametrization of drag by copses
Need version 4465 of ORCHIDEE at least

6. Cheruy

Location:

LMDZ5/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/yamada4.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/yamada4.F90

@@ -2 +2 @@
 
 SUBROUTINE yamada4(ni, nsrf, ngrid, dt, g, rconst, plev, temp, zlev, zlay, u, v, teta, &
-    cd, tke, km, kn, kq, ustar, iflag_pbl)
+    cd, tke, km, kn, kq, ustar, iflag_pbl, drgpro)
 
   USE dimphy
@@ -101 +101 @@
   INTEGER ni(klon)
 
+!FC
+  REAL drgpro(klon,klev)
+  REAL winds(klon,klev)
 
   ! Local
@@ -479 +482 @@
          DO ig=1,ngrid
          tkeprov=q2(ig,k)/ydeux
-         disseff=dissip(ig,k)-min(0.,buoy(ig,k))
+
+!             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
          tkeexp=exp(-dt*disseff/tkeprov)
-         tkeprov= shear(ig,k)*tkeprov/disseff*(1.-tkeexp)+tkeprov*tkeexp
+! on prend en compte la tke cree par les arbres
+         winds(ig,k)=sqrt(u(ig,k)**2+v(ig,k)**2)
+         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.