Index: trunk/LMDZ.COMMON/libf/evolution/xios_data.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/xios_data.F90	(revision 4136)
+++ trunk/LMDZ.COMMON/libf/evolution/xios_data.F90	(revision 4138)
@@ -194,6 +194,7 @@
             end do
         end do
-        ! Keep deep layers consistent with PCM bottom layer when nsoil > nsoil_PCM.
-        h2o_soildensity_avg(:,nsoil_PCM + 1:nsoil,islope) = h2o_soildensity_avg(:,nsoil_PCM,islope)
+        do isoil = nsoil_PCM + 1,nsoil
+            h2o_soildensity_avg(:,isoil,islope) = h2o_soildensity_avg(:,nsoil_PCM,islope) ! Explicit initialization because dimension with size nsoil > nsoil_PCM
+        end do
     end do
     h2o_soildensity_avg(:,:,:) = h2o_soildensity_avg(:,:,:)/real(nday,dp)