Changeset 1766 for LMDZ5/trunk/libf

Timestamp:

Jun 10, 2013, 4:51:15 PM (11 years ago)

Author:

Laurent Fairhead

Message:

Correction pour éviter une division par zéro
JYG

---

Correction to avoid division by zero
JYG

File:

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

LMDZ5/trunk/libf/phylmd/cv3_routines.F

@@ -3602 +3602 @@
             if(k.ge.icb(i).and.k.le.inb(i).and.
      &         j.ge.k.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