Index: LMDZ6/trunk/libf/phylmdiso/surf_land_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/surf_land_mod.F90	(revision 4527)
+++ LMDZ6/trunk/libf/phylmdiso/surf_land_mod.F90	(revision 4528)
@@ -171,7 +171,7 @@
 
 IF (ok_bs) THEN
-precip_totsnow=precip_snow+precip_bs
+precip_totsnow(:)=precip_snow(:)+precip_bs(:)
 ELSE
-precip_totsnow=precip_snow
+precip_totsnow(:)=precip_snow(:)
 ENDIF
 !****************************************************************************************
@@ -240,5 +240,5 @@
             precip_rain, precip_totsnow, lwdown_m, swnet, swdown, &
             pref_tmp, q2m, t2m, &
-            evap, fluxsens, fluxlat,fluxbs, &              
+            evap, fluxsens, fluxlat,  &              
             tsol_rad, tsurf_new, alb1_new, alb2_new, &
             emis_new, z0m, z0h, qsurf, &
Index: LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90	(revision 4527)
+++ LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90	(revision 4528)
@@ -490,9 +490,9 @@
 !****************************************************************************************
     IF (ok_bs) THEN
-      precip_totsnow=precip_snow+precip_bs
-      evap_totsnow=evap-fluxbs ! flux bs is positive towards the surface (snow erosion)
+      precip_totsnow(:)=precip_snow(:)+precip_bs(:)
+      evap_totsnow(:)=evap(:)-fluxbs(:) ! flux bs is positive towards the surface (snow erosion)
     ELSE
-      precip_totsnow=precip_snow
-      evap_totsnow=evap
+      precip_totsnow(:)=precip_snow(:)
+      evap_totsnow(:)=evap(:)
     ENDIF
 
Index: LMDZ6/trunk/libf/phylmdiso/surf_ocean_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/surf_ocean_mod.F90	(revision 4527)
+++ LMDZ6/trunk/libf/phylmdiso/surf_ocean_mod.F90	(revision 4528)
@@ -206,7 +206,7 @@
 
     IF (ok_bs) THEN
-       precip_totsnow=precip_snow+precip_bs
+       precip_totsnow(:)=precip_snow(:)+precip_bs(:)
     ELSE
-       precip_totsnow=precip_snow
+       precip_totsnow(:)=precip_snow(:)
     ENDIF