Index: /LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 5947)
+++ /LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 5948)
@@ -325,5 +325,4 @@
   END SUBROUTINE pbl_surface_init_iso
 #endif
-
 !
 !****************************************************************************************
@@ -1495,5 +1494,4 @@
 
 
-
   SUBROUTINE pbl_surface_uncompress_pre( &
        itap,          &
@@ -1542,4 +1540,6 @@
 #endif      
      &   )
+!$gpum horizontal klon
+
 !****************************************************************************************
 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
@@ -1624,5 +1624,4 @@
     use hbtm_mod, only: hbtm
     USE indice_sol_mod
-    USE time_phylmdz_mod,   ONLY : 
     USE mod_grid_phy_lmdz,  ONLY :  grid1dto2d_glo
     USE print_control_mod,  ONLY : prt_level
@@ -1643,6 +1642,4 @@
     USE ioipsl_getin_p_mod, ONLY : getin_p
     use phys_state_var_mod, only:  frac_tersrf, albedo_tersrf !AM 
-    use phys_output_var_mod, only: 
-    use lmdz_blowing_snow_ini, only : 
     USE dimsoil_mod_h, ONLY: nsoilmx
     USE surf_param_mod, ONLY: eff_surf_param  !AM
@@ -1871,5 +1868,5 @@
     REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: trmb3_w
 
-    REAL, PARAMETER                    :: facteur=2./sqrt(3.14)
+    REAL, PARAMETER                    :: facteur = 2. / 1.772  ! ( == 2. / SQRT(3.14))
     REAL, PARAMETER                    :: inertia=2000.
     REAL, DIMENSION(klon,nbsrf), INTENT(OUT)        :: pblh         ! height of the planetary boundary layer
@@ -2258,14 +2255,14 @@
 ! >> PC
    IF (carbon_cycle_cpl .AND. carbon_cycle_tr .AND. nbcf_out.GT.0 ) THEN
-       r_co2_ppm(:) = co2_send(:)
+       r_co2_ppm(1:klon) = co2_send(1:klon)
        DO n=1,nbcf_out
-           IF (cfname_out(n) == "atmco2" ) fields_out(:,n) = co2_send(:)
+           IF (cfname_out(n) == "atmco2" ) fields_out(1:klon,n) = co2_send(1:klon)
        ENDDO
    ENDIF
 
    IF ( .NOT. carbon_cycle_tr .AND. nbcf_out.GT.0 ) THEN
-       r_co2_ppm(:) = co2_ppm     ! Constant field
+       r_co2_ppm(1:klon) = co2_ppm     ! Constant field
        DO n=1,nbcf_out
-           IF (cfname_out(n) == "atmco2" ) fields_out(:,n) = co2_ppm
+           IF (cfname_out(n) == "atmco2" ) fields_out(1:klon,n) = co2_ppm
        ENDDO
    ENDIF
@@ -2881,5 +2878,5 @@
     REAL, DIMENSION(knon)              :: ydTs_ins, ydqs_ins
 !
-    REAL, PARAMETER                    :: facteur=2./sqrt(3.14)
+    REAL, PARAMETER                    :: facteur = 2. / 1.772  ! ( == 2. / SQRT(3.14))
     REAL, PARAMETER                    :: inertia=2000.
     REAL, DIMENSION(knon)              :: ydtsurf_th
@@ -5431,5 +5428,5 @@
     REAL, DIMENSION(klon, nbsrf), INTENT(IN)       :: pblh_w
     REAL, DIMENSION(klon, nbsrf), INTENT(IN)       :: plcl_w
-    REAL, PARAMETER                    :: facteur=2./sqrt(3.14)
+    REAL, PARAMETER                    :: facteur = 2. / 1.772  ! ( == 2. / SQRT(3.14))
     REAL, PARAMETER                    :: inertia=2000.
 
@@ -5613,5 +5610,6 @@
     zxtsol(:) = zxtsol(:)+1.5*meansqT(:)/zxtsol(:)
    ENDIF   ! iflag_order2_sollw == 1
-       
+
+!$gpum nocall       
        CALL checksum("n2mout", n2mout) 
        CALL checksum("n2mout_x", n2mout_x)