Index: trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90	(revision 2945)
+++ trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90	(revision 2950)
@@ -42,5 +42,5 @@
   pos_tend=0.
   neg_tend=0.
-
+if (ngrid.NE.1) then ! to make sure we are not in 1D 
   do i=1,ngrid
      do islope=1,nslope
@@ -54,6 +54,5 @@
      enddo
   enddo
-
-  if(neg_tend.GT.pos_tend .and. pos_tend.GT.0) then
+   if(neg_tend.GT.pos_tend .and. pos_tend.GT.0) then
      do i=1,ngrid
        do islope=1,nslope
@@ -65,5 +64,5 @@
        enddo
      enddo
-  elseif(neg_tend.LT.pos_tend .and. neg_tend.GT.0) then
+   elseif(neg_tend.LT.pos_tend .and. neg_tend.GT.0) then
      do i=1,ngrid
        do islope=1,nslope
@@ -75,5 +74,5 @@
        enddo
      enddo
-  elseif(pos_tend.EQ.0 .OR. neg_tend.EQ.0) then 
+   elseif(pos_tend.EQ.0 .OR. neg_tend.EQ.0) then 
     print *, "Reason of stopping : There is either no water ice sublimating or no water ice increasing !!"
     print *, "Tendencies on ice sublimating=", neg_tend
@@ -86,6 +85,6 @@
        enddo
     enddo
+   endif
   endif
-
   negative_part = 0.
 
@@ -115,6 +114,7 @@
     enddo
   enddo
-
-
+else ! ngrid==1;
+        qsurf(i,islope)=qsurf(i,islope)+new_tendencies(i,islope)*dt_pem
+endif
 
 END SUBROUTINE evol_h2o_ice_s_slope