Index: LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 1063)
+++ LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 1064)
@@ -1,4 +1,2 @@
-!
-! $Header$
 !
 MODULE pbl_surface_mod
@@ -152,8 +150,8 @@
 
     IF (type_ocean /= 'slab  ' .AND. type_ocean /= 'force ' .AND. type_ocean /= 'couple') THEN
-      WRITE(lunout,*)' *** Warning ***'
-      WRITE(lunout,*)'Option couplage pour l''ocean = ', type_ocean
-      abort_message='option pour l''ocean non valable'
-      CALL abort_gcm(modname,abort_message,1)
+       WRITE(lunout,*)' *** Warning ***'
+       WRITE(lunout,*)'Option couplage pour l''ocean = ', type_ocean
+       abort_message='option pour l''ocean non valable'
+       CALL abort_gcm(modname,abort_message,1)
     ENDIF
 
@@ -255,5 +253,4 @@
     INCLUDE "temps.h"
     INCLUDE "control.h"
-
 ! Input variables
 !****************************************************************************************
@@ -545,4 +542,6 @@
     ytsoil = 999999. 
 
+    rh2m(:)        = 0.
+    qsat2m(:)      = 0.
 !****************************************************************************************
 ! 3) - Calculate pressure thickness of each layer
@@ -642,5 +641,5 @@
 ! 
 !****************************************************************************************
-    
+   
     loop_nbsrf: DO nsrf = 1, nbsrf
 
@@ -988,6 +987,4 @@
        trmb2(:,nsrf)  = 0.        ! inhibition 
        trmb3(:,nsrf)  = 0.        ! Point Omega
-       rh2m(:)        = 0.
-       qsat2m(:)      = 0.
 
 #undef T2m