Changeset 1744

Timestamp:

Apr 9, 2013, 8:31:40 AM (11 years ago)

Author:

Ehouarn Millour

Message:

Fix to avoid dividing by zero.
JYG

File:

• LMDZ5/trunk/libf/phylmd/cv30_routines.F (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/cv30_routines.F

@@ -3106 +3106 @@
             if(k.ge.icb(i).and.k.le.inb(i).and.
      &         j.le.inb(i)) then
-             epm(i,j,k)=1.-(1.-ep(i,j))*clw(i,j)/elij(i,k,j)
+!!jyg             epm(i,j,k)=1.-(1.-ep(i,j))*clw(i,j)/elij(i,k,j)
+             epm(i,j,k)=1.-(1.-ep(i,j))*clw(i,j)/
+     &                     max(elij(i,k,j),1.e-16)
+!!
              epm(i,j,k)=max(epm(i,j,k),0.0)
             endif