Changeset 3784 for LMDZ6/trunk/libf/phylmd/yamada4.F90

Timestamp:

Nov 9, 2020, 3:09:22 PM (4 years ago)

Author:

evignon

Message:

correction v3783 pour convergence et compilation ancienne physique, Etienne aide par Ehouarn

File:

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

LMDZ6/trunk/libf/phylmd/yamada4.F90

@@ -189 +189 @@
        viscom=1.46E-5
        viscoh=2.06E-5
+       !lmixmin=1.0E-3
        lmixmin=0.
        yamada4_num=5
@@ -432 +433 @@
       shear(ig,k)=km(ig, k)*m2(ig, k)
       buoy(ig,k)=km(ig, k)*m2(ig, k)*(-1.*rif(ig,k))
-      dissip(ig,k)=min(max(((sqrt(q2(ig,k)))**3)/(b1*l(ig,k)),1.E-12),1.E4)
+!      dissip(ig,k)=min(max(((sqrt(q2(ig,k)))**3)/(b1*l(ig,k)),1.E-12),1.E4)
+      dissip(ig,k)=((sqrt(q2(ig,k)))**3)/(b1*l(ig,k)) 
      ENDDO
     ENDDO
@@ -537 +539 @@
          tkeexp=exp(-dt*disseff/tkeprov)
          tkeprov= tkeprov*tkeexp
-         q2(ig,k)=tkeprov*ydeux        
+         q2(ig,k)=tkeprov*ydeux
          
         ENDDO
@@ -556 +558 @@
 !     q2(1:ngrid, k) = q2(1:ngrid, k) + dt*km(1:ngrid, k)*m2(1:ngrid, k)*(1.-rif(1:ngrid,k))
       q2(1:ngrid, k) = min(max(q2(1:ngrid,k),1.E-10), 1.E4)
-      q2(1:ngrid, k) = 1./(1./sqrt(q2(1:ngrid,k))+dt/(yun*l(1:ngrid,k)*b1))
+       q2(1:ngrid, k) = 1./(1./sqrt(q2(1:ngrid,k))+dt/(yun*l(1:ngrid,k)*b1))
 !     q2(1:ngrid, k) = 1./(1./sqrt(q2(1:ngrid,k))+dt/(2*l(1:ngrid,k)*b1))
       q2(1:ngrid, k) = q2(1:ngrid, k)*q2(1:ngrid, k)
@@ -753 +755 @@
 ! Diagnostics
  tke_dissip(1:ngrid,:,nsrf)=0.
- DO k=2,klev
-    DO ig=1,ngrid
-       jg=ni(ig)
-       tke_dissip(jg,k,nsrf)=dissip(ig,k)
-    ENDDO
- ENDDO
+! DO k=2,klev
+!    DO ig=1,ngrid
+!       jg=ni(ig)
+!       tke_dissip(jg,k,nsrf)=dissip(ig,k)
+!    ENDDO
+! ENDDO
  
 !=============================================================================
@@ -1105 +1107 @@
    l2(1:ngrid,:)=0.0
    l_mixmin(1:ngrid,:,nsrf)=0.
-   l_mix(1:ngrid,:,nsrf)=1.E-5
+   l_mix(1:ngrid,:,nsrf)=0.
 
    IF (nsrf .EQ. 1) THEN
@@ -1174 +1176 @@
     DO ig=1,ngrid
        lmix(ig,k)=MAX(MAX(l1(ig,k), l2(ig,k)),lmixmin)
-       lmix(ig,k)=MAX(lmix(ig,k),1.E-5)
    ENDDO
  ENDDO
@@ -1180 +1181 @@
 ! Diagnostics
 
- DO k=2,klev+1
+ DO k=2,klev
     DO ig=1,ngrid
        jg=ni(ig)