Changeset 3267 for LMDZ6/trunk/libf

Timestamp:

Mar 15, 2018, 10:32:39 AM (7 years ago)

Author:

jyg

Message:

Refinement of commit 3253

File:

• LMDZ6/trunk/libf/phylmd/cvltr_noscav.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/cvltr_noscav.F90

@@ -39 +39 @@
   INTEGER                         :: i,k,j 
   REAL                            :: pdtimeRG
+  REAL                            :: smallest_mp
   real conserv
   real smfd
@@ -47 +48 @@
 ! calcul des tendances liees au downdraft
 ! =========================================
+!
+  smallest_mp = tiny(mp(1,1))
 !cdir collapse
   qfeed(:,it) = 0.
@@ -88 +91 @@
 !!              zmd(i,j,k)=zmd(i,j+1,k)*min(1.,mp(i,j)/mp(i,j+1))
 !!           ENDif
-           zmd(i,j,k)=zmd(i,j+1,k)*mp(i,j)/max(mp(i,j),max(mp(i,j+1),1.e-20))
+           zmd(i,j,k)=zmd(i,j+1,k)*mp(i,j)/max(mp(i,j),mp(i,j+1),smallest_mp)
         END DO
      END DO