Changeset 4675

Timestamp:

Sep 5, 2023, 2:29:43 PM (9 months ago)

Author:

evignon

Message:

correction option numerique 1 pour ATKE

File:

• LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90

@@ -240 +240 @@
         DO igrid=1,ngrid
            qq=max(sqrt(2.*tke(igrid,ilay)),1.e-10)
-           delta=1.+4.*dtime/cepsilon/l_exchange(igrid,ilay)/(2.**(3/2)) * &
-                (qq+dtime*l_exchange(igrid,ilay)/sqrt(2.)*Sm(igrid,ilay)*shear2(igrid,ilay) &
-                *(1.-Ri(igrid,ilay)/Prandtl(igrid,ilay)))
-           qq=(-1. + sqrt(delta))/dtime*cepsilon*sqrt(2.)*l_exchange(igrid,ilay)
+           delta=(2.*sqrt(2.)*cepsilon*l_exchange(igrid,ilay)/dtime)**2. &
+                 +4.*(2.*sqrt(2.)*cepsilon*l_exchange(igrid,ilay)/dtime*qq + &
+                 2.*l_exchange(igrid,ilay)*l_exchange(igrid,ilay)*cepsilon*Sm(igrid,ilay) &
+                 *shear2(igrid,ilay) * (1. - Ri(igrid,ilay) / Prandtl(igrid,ilay)))
+           qq=(-2.*sqrt(2.)*cepsilon*l_exchange(igrid,ilay)/dtime + sqrt(delta))/2.
+           qq=max(0.,qq)
            tke(igrid,ilay)=0.5*(qq**2) 
         ENDDO