Changeset 3253 for LMDZ6/trunk/libf

Timestamp:

Mar 12, 2018, 6:49:23 PM (7 years ago)

Author:

jyg

Message:

Regularization of a division by 0

File:

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

LMDZ6/trunk/libf/phylmd/cvltr_noscav.F90

@@ -85 +85 @@
      DO j=k-1,1,-1
         DO i=1,klon
-           if(mp(i,j+1).ne.0) then
-              zmd(i,j,k)=zmd(i,j+1,k)*min(1.,mp(i,j)/mp(i,j+1))
-           ENDif
+!!           if(mp(i,j+1).ne.0) then
+!!              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))
         END DO
      END DO