Index: LMDZ4/trunk/libf/phylmd/calltherm.F90
===================================================================
--- LMDZ4/trunk/libf/phylmd/calltherm.F90	(revision 1425)
+++ LMDZ4/trunk/libf/phylmd/calltherm.F90	(revision 1428)
@@ -373,4 +373,12 @@
 
 
+!jyg  Protection contre les temperatures nulles
+          do i=1,klon
+             do k=1,klev
+                if (ztla(i,k) .lt. 1.e-10) fraca(i,k) =0.
+             enddo
+          enddo
+
+
       return
 
Index: LMDZ4/trunk/libf/phylmd/physiq.F
===================================================================
--- LMDZ4/trunk/libf/phylmd/physiq.F	(revision 1425)
+++ LMDZ4/trunk/libf/phylmd/physiq.F	(revision 1428)
@@ -2232,4 +2232,9 @@
 
           DO i = 1, klon
+           itop_con(i) = min(max(itop_con(i),1),klev)
+           ibas_con(i) = min(max(ibas_con(i),1),itop_con(i))
+          ENDDO
+
+          DO i = 1, klon
             ema_pcb(i)  = paprs(i,ibas_con(i))
           ENDDO