Index: /trunk/LMDZ.COMMON/libf/evolution/adsorption_mod.F90
===================================================================
--- /trunk/LMDZ.COMMON/libf/evolution/adsorption_mod.F90	(revision 2938)
+++ /trunk/LMDZ.COMMON/libf/evolution/adsorption_mod.F90	(revision 2939)
@@ -230,5 +230,4 @@
    m_h2o_completesoil(:,:,:) = dm_h2o_regolith_slope(:,:,:)
 
-
  RETURN
 #endif
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 2938)
+++ /trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90	(revision 2939)
@@ -80,11 +80,10 @@
     print *, "Tendencies on ice increasing=", pos_tend
     print *, "This can be due to the absence of water ice in the PCM run!!"
-    print *,cell_area,qsurf(:,:),ngrid,STOPPING
-      call criterion_waterice_stop(cell_area,1.,qsurf(:,:)*0.,STOPPING,ngrid,qsurf(:,:)*0.)
-      do i=1,ngrid
-         do islope=1,nslope
-          new_tendencies(i,islope)=0
-         enddo
-      enddo
+    call criterion_waterice_stop(cell_area,1.,qsurf(:,:)*0.,STOPPING,ngrid,qsurf(:,:)*0.)
+    do i=1,ngrid
+       do islope=1,nslope
+         new_tendencies(i,islope)=0
+       enddo
+    enddo
   endif
 
Index: /trunk/LMDZ.COMMON/libf/evolution/pem.F90
===================================================================
--- /trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 2938)
+++ /trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 2939)
@@ -749,4 +749,5 @@
 ! II.a.1. Compute updated global pressure
      print *, "Recomputing the new pressure..."
+
      do i=1,ngrid
        do islope=1,nslope
@@ -754,5 +755,5 @@
       enddo
      enddo
-       print *, 'Global average pressure old time step',global_ave_press_old
+     print *, 'Global average pressure old time step',global_ave_press_old
 
      call WRITEDIAGFI(ngrid,'ps_ave','Global average pressure','Pa',0,global_ave_press_new)
Index: /trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90
===================================================================
--- /trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90	(revision 2938)
+++ /trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90	(revision 2939)
@@ -90,5 +90,4 @@
 
 write(*,*)'Is start PEM?',startpem_file
-startpem_file = .true.
 
 !1. Run
@@ -275,4 +274,5 @@
 
 !4. CO2 & H2O Adsorption
+
  if(adsorption_pem) then
   DO islope=1,nslope
@@ -281,15 +281,17 @@
     if((.not.found)) then
        m_co2_regolith_phys(:,:,:) = 0.
+       exit
     endif
-    exit
+    
   ENDDO
 
   DO islope=1,nslope
    write(num,fmt='(i2.2)') islope
-   call get_field("mh2o_reg_ads_slope"//num,m_co2_regolith_phys(:,:,islope),found2)
+   call get_field("mh2o_reg_ads_slope"//num,m_h2o_regolith_phys(:,:,islope),found2)
     if((.not.found2)) then
        m_h2o_regolith_phys(:,:,:) = 0.
+      exit
     endif
-    exit
+    
   ENDDO
 
Index: /trunk/LMDZ.COMMON/libf/evolution/update_soil.F90
===================================================================
--- /trunk/LMDZ.COMMON/libf/evolution/update_soil.F90	(revision 2938)
+++ /trunk/LMDZ.COMMON/libf/evolution/update_soil.F90	(revision 2939)
@@ -44,5 +44,4 @@
    regolith_inertia(:,islope) = inertiedat_PEM(:,1)
    do ig = 1,ngrid
-      write(*,*) 'ig,islope',ig,islope,tendencies_waterice(ig,islope),waterice(ig,islope)
       if((tendencies_waterice(ig,islope).lt.-1e-5).and.(waterice(ig,islope).eq.0)) then
               regolith_inertia(ig,islope) = inertie_averaged