Index: LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90	(revision 3974)
+++ LMDZ6/trunk/libf/phylmdiso/ocean_forced_mod.F90	(revision 3975)
@@ -262,4 +262,6 @@
 !
     USE dimphy
+    USE geometry_mod, ONLY: longitude,latitude
+
     USE calcul_fluxs_mod
     USE surface_data,     ONLY : calice, calsno
@@ -379,5 +381,7 @@
     IF (soil_model) THEN 
 ! update tsoil and calculate soilcap and soilflux
-       CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, tsoil,soilcap, soilflux)
+       CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, qsol, &
+        & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil,soilcap, soilflux)
+
        cal(1:knon) = RCPD / soilcap(1:knon)
        radsol(1:knon) = radsol(1:knon)  + soilflux(1:knon)
Index: LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90	(revision 3974)
+++ LMDZ6/trunk/libf/phylmdiso/surf_land_bucket_mod.F90	(revision 3975)
@@ -30,5 +30,5 @@
     USE cpl_mod
     USE dimphy
-    USE geometry_mod, ONLY: latitude 
+    USE geometry_mod, ONLY: longitude,latitude 
     USE mod_grid_phy_lmdz
     USE mod_phys_lmdz_para
@@ -176,5 +176,7 @@
     IF (soil_model) THEN 
 !       write(*,*) 'surf_land_bucket 185'
-       CALL soil(dtime, is_ter, knon, snow, tsurf, tsoil, soilcap, soilflux)
+       CALL soil(dtime, is_ter, knon, snow, tsurf, qsol,  &
+      & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
+
 !       write(*,*) 'surf_land_bucket 187'
        DO i=1, knon
Index: LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90	(revision 3974)
+++ LMDZ6/trunk/libf/phylmdiso/surf_landice_mod.F90	(revision 3975)
@@ -44,4 +44,6 @@
 #endif
 #endif
+    USE geometry_mod,     ONLY : longitude,latitude
+
 !FC
     USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -53,4 +55,5 @@
 
     USE indice_sol_mod
+
 
 !    INCLUDE "indicesol.h"
@@ -337,5 +340,7 @@
     ! use soil model and recalculate properly cal
     IF (soil_model) THEN 
-       CALL soil(dtime, is_lic, knon, snow, tsurf, tsoil, soilcap, soilflux)
+       CALL soil(dtime, is_lic, knon, snow, tsurf, qsol, &
+    & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
+
        cal(1:knon) = RCPD / soilcap(1:knon)
        radsol(1:knon)  = radsol(1:knon) + soilflux(1:knon)