Index: LMDZ6/trunk/libf/phylmd/oasis.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/oasis.F90	(revision 4595)
+++ LMDZ6/trunk/libf/phylmd/oasis.F90	(revision 4596)
@@ -141,4 +141,5 @@
     CHARACTER (len = 80)               :: abort_message 
     LOGICAL, SAVE                      :: cpl_current_omp
+    INTEGER, DIMENSION(klon_mpi)       :: ind_cell_glo_mpi
 
 !*    1. Initializations
@@ -175,4 +176,9 @@
     cpl_current = cpl_current_omp
     WRITE(lunout,*) 'Couple ocean currents, cpl_current = ',cpl_current 
+
+!************************************************************************************
+! Gather global index to be used for oasis decomposition
+!************************************************************************************
+    CALL gather_omp(ind_cell_glo,ind_cell_glo_mpi)
 
 !************************************************************************************
@@ -291,5 +297,5 @@
 
       DO jf=1, klon_mpi_para_nb(mpi_rank)
-        ig_paral(2+jf) = ind_cell_glo(jf)
+        ig_paral(2+jf) = ind_cell_glo_mpi(jf)
       ENDDO
 
Index: LMDZ6/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 4595)
+++ LMDZ6/trunk/libf/phylmd/physiq_mod.F90	(revision 4596)
@@ -1315,6 +1315,6 @@
        irneb= strIdx(tracers(:)%name, addPhase('H2O', 'r'))
        ibs  = strIdx(tracers(:)%name, addPhase('H2O', 'b'))
-       CALL init_etat0_limit_unstruct
-       IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
+!       CALL init_etat0_limit_unstruct
+!       IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
        !CR:nvelles variables convection/poches froides
 
@@ -1339,4 +1339,6 @@
             read_climoz, &
             alp_offset)
+       CALL init_etat0_limit_unstruct
+       IF (.NOT. create_etat0_limit) CALL init_limit_read(days_elapsed)
        CALL phys_state_var_init(read_climoz)
        CALL phys_output_var_init