Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/inigeomphy_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/inigeomphy_mod.F90	(revision 1621)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/inigeomphy_mod.F90	(revision 1682)
@@ -73,5 +73,6 @@
   REAL,ALLOCATABLE,SAVE :: boundslonfi(:,:)
   REAL,ALLOCATABLE,SAVE :: boundslatfi(:,:)
-!$OMP THREADPRIVATE (latfi,lonfi,cufi,cvfi,airefi,boundslonfi,boundslatfi)
+  INTEGER,ALLOCATABLE,SAVE :: ind_cell_glo_fi(:)
+!$OMP THREADPRIVATE (latfi,lonfi,cufi,cvfi,airefi,boundslonfi,boundslatfi,ind_cell_glo_fi)
 
   ! Initialize Physics distibution and parameters and interface with dynamics
@@ -201,5 +202,6 @@
   ALLOCATE(boundslonfi(klon_omp,4))
   ALLOCATE(boundslatfi(klon_omp,4))
-!  CALL initcomgeomphy
+  ALLOCATE(ind_cell_glo_fi(klon_omp))
+
 
   offset = klon_mpi_begin - 1
@@ -211,8 +213,9 @@
   boundslonfi(1:klon_omp,:) = boundslonfi_glo(offset+klon_omp_begin:offset+klon_omp_end,:)
   boundslatfi(1:klon_omp,:) = boundslatfi_glo(offset+klon_omp_begin:offset+klon_omp_end,:)
+  ind_cell_glo_fi(1:klon_omp)=(/ (i,i=offset+klon_omp_begin,offset+klon_omp_end) /)
 
   ! copy over local grid longitudes and latitudes
   CALL init_geometry(klon_omp,lonfi,latfi,boundslonfi,boundslatfi, &
-                     airefi,cufi,cvfi)
+                     airefi,ind_cell_glo_fi,cufi,cvfi)
 
   ! copy over preff , ap(), bp(), etc 
Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/mod_interface_dyn_phys.F90
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/mod_interface_dyn_phys.F90	(revision 1621)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/mod_interface_dyn_phys.F90	(revision 1682)
@@ -22,5 +22,5 @@
     
     k=1
-    IF (is_north_pole) THEN
+    IF (is_north_pole_dyn) THEN
       index_i(k)=1
       index_j(k)=1
@@ -42,5 +42,5 @@
     ENDDO
     
-    IF (is_south_pole) THEN
+    IF (is_south_pole_dyn) THEN
       index_i(k)=1
       index_j(k)=jj_end
Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/iniphysiq_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/iniphysiq_mod.F90	(revision 1621)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/iniphysiq_mod.F90	(revision 1682)
@@ -66,6 +66,6 @@
 ! and do some initializations 
 
-! Initialize dimphy module
-call init_dimphy(klon_omp,nlayer)
+! Initialize dimphy module => Now done in physics_distribution_mod
+!call init_dimphy(klon_omp,nlayer)
 
 call phys_state_var_init(klon_omp,nlayer,nqtot,tname, &