Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/Makefile
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/Makefile	(revision 1578)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/Makefile	(revision 1579)
@@ -11,4 +11,5 @@
         update_inputs_physiq_mod.o \
         update_outputs_physiq_mod.o \
+        iniphysiq_mod.o \
 	module_lmd_driver.o 
  
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new	(revision 1578)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new	(revision 1579)
@@ -62,7 +62,6 @@
    USE update_outputs_physiq_mod !! to get outputs from physiq
    USE comm_wrf !! to get fields to be written from physiq
+   USE iniphysiq_mod !! to get iniphysiq subroutine
    USE physiq_mod, only : physiq
-   USE phys_state_var_init_mod, only : phys_state_var_init
-   !!! ajouter init_time (voir iniphysiq. ne plus le faire dans la physiq)
 #endif
 
@@ -454,13 +453,11 @@
             lct_input,lon_input,ls_input,&
             ptime,pday,MY)
-  !!
-  CALL phys_state_var_init(ngrid,nlayer,nq,&
-                           floor(pday),ptime,&
-                           wdaysec,ptimestep,1/reradius,g,r_d,cp)
   !! Fill planetary parameters in modules
   !! Values defined in the module_model_constants.F WRF routine
   CALL update_inputs_physiq_constants
-  !! Read callphys.def !!
-  call conf_phys(ngrid,nlayer,nq)
+  !! Initialize physics
+  CALL iniphysiq(ngrid,nlayer,nq,ptime,wappel_phys,&
+                     wdaysec,floor(pday),ptimestep, &
+                     1./reradius,g,r_d,cp,1)
 ENDIF allocation_firstcall
 #endif
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/update_inputs_physiq_mod.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/update_inputs_physiq_mod.F	(revision 1578)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/update_inputs_physiq_mod.F	(revision 1579)
@@ -227,6 +227,13 @@
             XLAT,XLONG)
 
+   ! in WRF (share)
    USE module_model_constants, only: DEGRAD
+   ! in LMD (phymars)
    use comgeomfi_h, only: ini_fillgeom
+   ! in LMD (phy_common)
+   USE mod_grid_phy_lmdz, ONLY: init_grid_phy_lmdz
+   USE geometry_mod, ONLY: latitude,latitude_deg,&
+                           longitude,longitude_deg,&
+                           cell_area
 
    INTEGER, INTENT(IN) :: ims,ime,jms,jme
@@ -273,4 +280,25 @@
    !! FILL GEOMETRICAL ARRAYS !!
    call ini_fillgeom(ngrid,plat,plon,parea)
+
+   !!! ----------------------------------------------------------
+   !!! --- initializing geometry in phy_common
+   !!! --- (this is quite planet-independent)
+   !!! ----------------------------------------------------------
+   ! initialize mod_grid_phy_lmdz
+   CALL init_grid_phy_lmdz(1,1,ipe-ips+1,jpe-jps+1,nlayer)
+   ! fill in geometry_mod variables
+   ! ... copy over local grid longitudes and latitudes
+   ! ... partly what is done in init_geometry
+   IF(.not.ALLOCATED(longitude)) ALLOCATE(longitude(ngrid))
+   IF(.not.ALLOCATED(longitude_deg)) ALLOCATE(longitude_deg(ngrid))
+   IF(.not.ALLOCATED(latitude)) ALLOCATE(latitude(ngrid))
+   IF(.not.ALLOCATED(latitude_deg)) ALLOCATE(latitude_deg(ngrid))
+   IF(.not.ALLOCATED(cell_area)) ALLOCATE(cell_area(ngrid))
+   longitude(:) = plon(:)
+   latitude(:) = plat(:)
+   longitude_deg(:) = plon(:)/DEGRAD
+   latitude_deg(:) = plat(:)/DEGRAD
+   cell_area(:) = parea(:)
+   !!! ----------------------------------------------------------
 
 END SUBROUTINE update_inputs_physiq_geom