Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/callphysiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/callphysiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/callphysiq_mod.F	(revision 2277)
@@ -13,12 +13,7 @@
 
 SUBROUTINE call_physiq(planet_type, klon,llm,nqtot,                       &
-                       debut_split,lafin_split,                           &
-                       jD_cur,jH_cur_split,zdt_split,                     &
-                       zplev_omp,zplay_omp,                               &
-                       zpk_omp,zphi_omp,zphis_omp,                        &
-                       presnivs_omp,                                      &
-                       zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp,      &
-                       flxwfi_omp,                                        &
-                       zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp)
+                       debut_split,lafin_split)
+
+  USE variables_mod
 
   IMPLICIT NONE
@@ -30,25 +25,4 @@
   LOGICAL,INTENT(IN) :: debut_split ! .true. if very first call to physics
   LOGICAL,INTENT(IN) :: lafin_split ! .true. if last call to physics
-  REAL,INTENT(IN) :: JD_cur ! Julian day
-  REAL,INTENT(IN) :: JH_cur_split ! Julian hour (fraction of day)
-  REAL,INTENT(IN) :: zdt_split ! time step over which the physics are evaluated
-  REAL,INTENT(IN) :: zplev_omp(klon,llm+1) ! interlayer pressure (Pa)
-  REAL,INTENT(IN) :: zplay_omp(klon,llm) ! mid-layer pressure (Pa)
-  REAL,INTENT(INOUT) :: zpk_omp(klon,llm)
-  REAL,INTENT(IN) :: zphi_omp(klon,llm) ! geopotential at midlayer
-  REAL,INTENT(INOUT) :: zphis_omp(klon) ! surface geopotential
-  REAL,INTENT(INOUT) :: presnivs_omp(llm) ! approximate pressure of atm. layers
-  REAL,INTENT(IN) :: zufi_omp(klon,llm) ! zonal wind (m/s)
-  REAL,INTENT(IN) :: zvfi_omp(klon,llm) ! meridional wind (m/s)
-  REAL,INTENT(INOUT) :: zrfi_omp(klon,llm) ! relative wind vorticity, in s-1
-  REAL,INTENT(IN) :: ztfi_omp(klon,llm) ! temperature (K)
-  REAL,INTENT(IN) :: zqfi_omp(klon,llm,nqtot) ! tracers (*/kg of air)
-  REAL,INTENT(IN) :: flxwfi_omp(klon,llm) ! Vertical mass flux on lower mesh interfaces (kg/s) 
-  ! tendencies (in */s) from the physics:
-  REAL,INTENT(OUT) :: zdufi_omp(klon,llm) ! tendency on zonal winds
-  REAL,INTENT(OUT) :: zdvfi_omp(klon,llm) ! tendency on meridional winds
-  REAL,INTENT(OUT) :: zdtfi_omp(klon,llm) ! tendency on temperature
-  REAL,INTENT(OUT) :: zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
-  REAL,INTENT(OUT) :: zdpsrf_omp(klon) ! tendency on surface pressure
   ! data file
   INTEGER :: i,k
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/iniphysiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/iniphysiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/iniphysiq_mod.F	(revision 2277)
@@ -3,6 +3,6 @@
 CONTAINS
 
-subroutine iniphysiq(ngrid,nlayer,nq,phour_ini,piphysiq,&
-                     punjours, pdayref,ptimestep, &
+subroutine iniphysiq(ngrid,nlayer,nq,piphysiq,&
+                     punjours, pdayref, &
                      prad,pg,pr,pcpp,iflag_phys)
 
@@ -15,5 +15,4 @@
 real,intent(in) :: punjours ! length (in s) of a standard day [daysec]
 integer,intent(in) :: pdayref ! reference day of for the simulation [day_ini]
-real,intent(in) :: ptimestep !physics time step (s) [dtphys]
 integer,intent(in) :: iflag_phys ! type of physics to be called
 
@@ -21,5 +20,4 @@
 integer,intent(in) :: nlayer ! number of atmospheric layers
 integer,intent(in) :: nq ! number of tracers
-real,intent(in) :: phour_ini   ! start time (fraction of day) of the run 0=<phour_ini<1
 real,intent(in) :: piphysiq   ! call physics every piphysiq dynamical timesteps
 
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_inputs_physiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_inputs_physiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_inputs_physiq_mod.F	(revision 2277)
@@ -19,9 +19,9 @@
             elaps,&
             lct_input,lon_input,ls_input,&
-            ptime,pday,MY)
+            MY)
 
   INTEGER, INTENT(IN) :: JULDAY, JULYR
   REAL, INTENT(IN) :: GMT,elaps,lon_input,ls_input,lct_input
-  REAL,INTENT(OUT) :: pday,ptime,MY
+  REAL,INTENT(OUT) :: MY
 
 END SUBROUTINE update_inputs_physiq_time
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_outputs_physiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_outputs_physiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/update_outputs_physiq_mod.F	(revision 2277)
@@ -87,4 +87,5 @@
      HR_SW,HR_LW,RDUST,VMR_ICE,RICE,CLOUDFRAC,HR_DYN,DT,DTRAD,RAIN,SNOW,&
      LSCEZ,H2OICE_REFF
+   REAL, DIMENSION(kms:kme) :: phr_sw,phr_lw,phr_dyn
    INTEGER :: i,j,subs
 
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/variables_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/variables_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_prescribed_void_lmd_new/variables_mod.F	(revision 2277)
@@ -2,4 +2,28 @@
 
 IMPLICIT NONE
+
+REAL*8 :: JD_cur !pday !JD_cur ! Julian day
+REAL*8 :: JH_cur_split !ptime !JH_cur_split ! Julian hour (fraction of day)
+REAL*8 :: zdt_split !ptimestep !zdt_split ! time step over which the physics are evaluated
+REAL*8 :: phour_ini ! start time (fraction of day) of the run 0=<phour_ini<1
+REAL*8 :: counter
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zplev_omp !pplev !zplev_omp(klon,llm+1) ! interlayer pressure (Pa)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zplay_omp !pplay !zplay_omp(klon,llm) ! mid-layer pressure (Pa)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zpk_omp!(klon,llm)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zphi_omp !pphi !zphi_omp(klon,llm) ! geopotential at midlayer
+REAL*8,DIMENSION(:),ALLOCATABLE :: zphis_omp!(klon) ! surface geopotential
+REAL*8,DIMENSION(:),ALLOCATABLE :: presnivs_omp!(llm) ! approximate pressure of atm. layers
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zufi_omp !pu !zufi_omp(klon,llm) ! zonal wind (m/s)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zvfi_omp !pv !zvfi_omp(klon,llm) ! meridional wind (m/s)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zrfi_omp!(klon,llm) ! relative wind vorticity, in s-1
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: ztfi_omp !pt !ztfi_omp(klon,llm) ! temperature (K)
+REAL*8,DIMENSION(:,:,:),ALLOCATABLE :: zqfi_omp !pq !zqfi_omp(klon,llm,nqtot) ! tracers (*/kg of air)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: flxwfi_omp !flxw !flxwfi_omp(klon,llm) ! Vertical mass flux on lower mesh interfaces (kg/s)  ! tendencies (in */s) from the physics:
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zdufi_omp !pdu !zdufi_omp(klon,llm) ! tendency on zonal winds
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zdvfi_omp !pdv !zdvfi_omp(klon,llm) ! tendency on meridional winds
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zdtfi_omp !pdt !zdtfi_omp(klon,llm) ! tendency on temperature
+REAL*8,DIMENSION(:,:,:),ALLOCATABLE :: zdqfi_omp !pdq !zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
+REAL*8,DIMENSION(:),ALLOCATABLE :: zdpsrf_omp !pdpsrf !zdpsrf_omp(klon) !tendency on surface pressure
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: hrlw,hrsw
 
 CONTAINS
@@ -8,13 +32,53 @@
 
   IMPLICIT NONE
- 
+
   INTEGER,INTENT(IN) :: klon ! (local) number of atmospheric columns
   INTEGER,INTENT(IN) :: llm  ! number of atmospheric layers
   INTEGER,INTENT(IN) :: nqtot ! number of tracers
- 
+
+  IF (.NOT.ALLOCATED(zplev_omp)) ALLOCATE(zplev_omp(klon,llm+1))
+  IF (.NOT.ALLOCATED(zplay_omp)) ALLOCATE(zplay_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zpk_omp)) ALLOCATE(zpk_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zphi_omp)) ALLOCATE(zphi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zphis_omp)) ALLOCATE(zphis_omp(klon))
+  IF (.NOT.ALLOCATED(presnivs_omp)) ALLOCATE(presnivs_omp(llm))
+  IF (.NOT.ALLOCATED(zufi_omp)) ALLOCATE(zufi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zvfi_omp)) ALLOCATE(zvfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zrfi_omp)) ALLOCATE(zrfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(ztfi_omp)) ALLOCATE(ztfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zqfi_omp)) ALLOCATE(zqfi_omp(klon,llm,nqtot))
+  IF (.NOT.ALLOCATED(flxwfi_omp)) ALLOCATE(flxwfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdufi_omp)) ALLOCATE(zdufi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdvfi_omp)) ALLOCATE(zdvfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdtfi_omp)) ALLOCATE(zdtfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdqfi_omp)) ALLOCATE(zdqfi_omp(klon,llm,nqtot))
+  IF (.NOT.ALLOCATED(zdpsrf_omp)) ALLOCATE(zdpsrf_omp(klon))
+  IF (.NOT.ALLOCATED(hrlw)) ALLOCATE(hrlw(klon,llm))
+  IF (.NOT.ALLOCATED(hrsw)) ALLOCATE(hrsw(klon,llm))
+
 END SUBROUTINE allocate_interface
 
 SUBROUTINE deallocate_interface
 
+  DEALLOCATE(zplev_omp)
+  DEALLOCATE(zplay_omp)
+  DEALLOCATE(zpk_omp)
+  DEALLOCATE(zphi_omp)
+  DEALLOCATE(zphis_omp)
+  DEALLOCATE(presnivs_omp)
+  DEALLOCATE(zufi_omp)
+  DEALLOCATE(zvfi_omp)
+  DEALLOCATE(zrfi_omp)
+  DEALLOCATE(ztfi_omp)
+  DEALLOCATE(zqfi_omp)
+  DEALLOCATE(flxwfi_omp)
+  DEALLOCATE(zdufi_omp)
+  DEALLOCATE(zdvfi_omp)
+  DEALLOCATE(zdtfi_omp)
+  DEALLOCATE(zdqfi_omp)
+  DEALLOCATE(zdpsrf_omp)
+  DEALLOCATE(hrlw)
+  DEALLOCATE(hrsw)
+
 END SUBROUTINE deallocate_interface
 
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/callphysiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/callphysiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/callphysiq_mod.F	(revision 2277)
@@ -18,5 +18,4 @@
   USE physiq_mod, ONLY: physiq
   USE module_model_constants, only : p0,rcp,cp
-  use moyzon_mod, only: tmoy
   IMPLICIT NONE
 
@@ -81,6 +80,4 @@
               zqfi_omp,       &
               flxwfi_omp,     &
-              zplevmoy,       &
-              ztmoy,          &
               zdufi_omp,      &
               zdvfi_omp,      &
@@ -88,5 +85,4 @@
               zdqfi_omp,      &
               zdpsrf_omp)
-   deallocate(tmoy)
 END SUBROUTINE call_physiq
 
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/iniphysiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/iniphysiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/iniphysiq_mod.F	(revision 2277)
@@ -62,5 +62,5 @@
       call init_dimphy(ngrid,nlayer)
 
-      Call phys_state_var_init      
+      Call phys_state_var_init(nq)
       ! allocate "surfdat_h" arrays
       !call ini_surfdat_h(ngrid,nq)
@@ -79,5 +79,5 @@
 
       ! Initialize some physical constants
-      call suphec
+      call suphec(pcpp)
 
       cpp=pcpp
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F	(revision 2277)
@@ -60,4 +60,5 @@
     JH_cur_split = JH_cur_split / 2808.
     print*,'JH_cur_split',JH_cur_split
+    JD_cur = INT((3700*(lct_input - lon_input / 15.) +elaps)/88800)
     JD_cur = MODULO(int(JD_cur),669)
     MY = 2024
@@ -70,5 +71,5 @@
 SUBROUTINE update_inputs_physiq_tracers(nq,MARS_MODE)
 
-  use infotrac, only: nqtot,tname
+  use infotrac_phy, only: nqtot,tname
   use logic_mod, only: iflag_trac
   implicit none
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/callphysiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/callphysiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/callphysiq_mod.F	(revision 2277)
@@ -13,17 +13,11 @@
 
 SUBROUTINE call_physiq(planet_type, klon,llm,nqtot,                       &
-                       debut_split,lafin_split,                           &
-                       jD_cur,jH_cur_split,zdt_split,                     &
-                       zplev_omp,zplay_omp,                               &
-                       zpk_omp,zphi_omp,zphis_omp,                        &
-                       presnivs_omp,                                      &
-                       zufi_omp,zvfi_omp,zrfi_omp,ztfi_omp,zqfi_omp,      &
-                       flxwfi_omp,                                        &
-                       zdufi_omp,zdvfi_omp,zdtfi_omp,zdqfi_omp,zdpsrf_omp)
+                       debut_split,lafin_split)
+
+  USE variables_mod
 
   IMPLICIT NONE
 
   character(len=10),INTENT(IN) :: planet_type ! planet type ('earth','mars',...)
-
   INTEGER,INTENT(IN) :: klon ! (local) number of atmospheric columns
   INTEGER,INTENT(IN) :: llm  ! number of atmospheric layers
@@ -31,25 +25,4 @@
   LOGICAL,INTENT(IN) :: debut_split ! .true. if very first call to physics
   LOGICAL,INTENT(IN) :: lafin_split ! .true. if last call to physics
-  REAL,INTENT(IN) :: JD_cur ! Julian day
-  REAL,INTENT(IN) :: JH_cur_split ! Julian hour (fraction of day)
-  REAL,INTENT(IN) :: zdt_split ! time step over which the physics are evaluated
-  REAL,INTENT(IN) :: zplev_omp(klon,llm+1) ! interlayer pressure (Pa)
-  REAL,INTENT(IN) :: zplay_omp(klon,llm) ! mid-layer pressure (Pa)
-  REAL,INTENT(INOUT) :: zpk_omp(klon,llm)
-  REAL,INTENT(IN) :: zphi_omp(klon,llm) ! geopotential at midlayer
-  REAL,INTENT(INOUT) :: zphis_omp(klon) ! surface geopotential
-  REAL,INTENT(INOUT) :: presnivs_omp(llm) ! approximate pressure of atm. layers
-  REAL,INTENT(IN) :: zufi_omp(klon,llm) ! zonal wind (m/s)
-  REAL,INTENT(IN) :: zvfi_omp(klon,llm) ! meridional wind (m/s)
-  REAL,INTENT(INOUT) :: zrfi_omp(klon,llm) ! relative wind vorticity, in s-1
-  REAL,INTENT(IN) :: ztfi_omp(klon,llm) ! temperature (K)
-  REAL,INTENT(IN) :: zqfi_omp(klon,llm,nqtot) ! tracers (*/kg of air)
-  REAL,INTENT(IN) :: flxwfi_omp(klon,llm) ! Vertical mass flux on lower mesh interfaces (kg/s) 
-  ! tendencies (in */s) from the physics:
-  REAL,INTENT(OUT) :: zdufi_omp(klon,llm) ! tendency on zonal winds
-  REAL,INTENT(OUT) :: zdvfi_omp(klon,llm) ! tendency on meridional winds
-  REAL,INTENT(OUT) :: zdtfi_omp(klon,llm) ! tendency on temperature
-  REAL,INTENT(OUT) :: zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
-  REAL,INTENT(OUT) :: zdpsrf_omp(klon) ! tendency on surface pressure
 
 END SUBROUTINE call_physiq
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/iniphysiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/iniphysiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/iniphysiq_mod.F	(revision 2277)
@@ -3,6 +3,6 @@
 CONTAINS
 
-subroutine iniphysiq(ngrid,nlayer,nq,phour_ini,piphysiq,&
-                     punjours, pdayref,ptimestep, &
+subroutine iniphysiq(ngrid,nlayer,nq,piphysiq,&
+                     punjours, pdayref, &
                      prad,pg,pr,pcpp,iflag_phys)
 
@@ -15,5 +15,4 @@
 real,intent(in) :: punjours ! length (in s) of a standard day [daysec]
 integer,intent(in) :: pdayref ! reference day of for the simulation [day_ini]
-real,intent(in) :: ptimestep !physics time step (s) [dtphys]
 integer,intent(in) :: iflag_phys ! type of physics to be called
 
@@ -21,5 +20,4 @@
 integer,intent(in) :: nlayer ! number of atmospheric layers
 integer,intent(in) :: nq ! number of tracers
-real,intent(in) :: phour_ini   ! start time (fraction of day) of the run 0=<phour_ini<1
 real,intent(in) :: piphysiq   ! call physics every piphysiq dynamical timesteps
 
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/update_inputs_physiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/update_inputs_physiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/update_inputs_physiq_mod.F	(revision 2277)
@@ -19,9 +19,9 @@
             elaps,&
             lct_input,lon_input,ls_input,&
-            ptime,pday,MY)
+            MY)
 
   INTEGER, INTENT(IN) :: JULDAY, JULYR
   REAL, INTENT(IN) :: GMT,elaps,lon_input,ls_input,lct_input
-  REAL,INTENT(OUT) :: pday,ptime,MY
+  REAL,INTENT(OUT) :: MY
 
 END SUBROUTINE update_inputs_physiq_time
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/update_outputs_physiq_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/update_outputs_physiq_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/update_outputs_physiq_mod.F	(revision 2277)
@@ -71,6 +71,9 @@
             SWDOWNZ,TAU_DUST,QSURFDUST,&
             MTOT,ICETOT,TAU_ICE,&
-            HR_SW,HR_LW,HR_DYN,DT,DTRAD,DTVDF,DTAJS,&
-            RDUST,VMR_ICE,RICE)
+            HR_SW,HR_LW,HR_DYN,DT,DTRAD,&
+            RDUST,VMR_ICE,RICE,&
+            CLOUDFRAC,TOTCLOUDFRAC,&
+            RAIN,SNOW,REEVAP,SURFRAIN,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,&
+            FLUXSURF_LW,FLXGRD,LSCEZ,H2OICE_REFF,LATENT_HF)
 
    INTEGER, INTENT(IN) :: ims,ime,jms,jme,kms,kme
@@ -78,11 +81,13 @@
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: &
      SWDOWNZ,TAU_DUST,QSURFDUST,&
-     MTOT,ICETOT,TAU_ICE
+     MTOT,ICETOT,TAU_ICE,&
+     TOTCLOUDFRAC,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,&
+     FLUXSURF_LW,FLXGRD,LATENT_HF,REEVAP,SURFRAIN
    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: &
-     HR_SW,HR_LW,HR_DYN,RDUST,VMR_ICE,RICE,DT,DTRAD,DTVDF,DTAJS
+     HR_SW,HR_LW,RDUST,VMR_ICE,RICE,CLOUDFRAC,HR_DYN,DT,DTRAD,RAIN,SNOW,&
+     LSCEZ,H2OICE_REFF
    INTEGER :: i,j,subs
 
 END SUBROUTINE update_outputs_physiq_diag
-
 
 END MODULE update_outputs_physiq_mod
Index: /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/variables_mod.F
===================================================================
--- /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/variables_mod.F	(revision 2276)
+++ /trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_void_lmd_new/variables_mod.F	(revision 2277)
@@ -2,4 +2,28 @@
 
 IMPLICIT NONE
+
+REAL*8 :: JD_cur !pday !JD_cur ! Julian day
+REAL*8 :: JH_cur_split !ptime !JH_cur_split ! Julian hour (fraction of day)
+REAL*8 :: zdt_split !ptimestep !zdt_split ! time step over which the physics are evaluated
+REAL*8 :: phour_ini ! start time (fraction of day) of the run 0=<phour_ini<1
+REAL*8 :: counter
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zplev_omp !pplev !zplev_omp(klon,llm+1) ! interlayer pressure (Pa)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zplay_omp !pplay !zplay_omp(klon,llm) ! mid-layer pressure (Pa)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zpk_omp!(klon,llm)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zphi_omp !pphi !zphi_omp(klon,llm) ! geopotential at midlayer
+REAL*8,DIMENSION(:),ALLOCATABLE :: zphis_omp!(klon) ! surface geopotential
+REAL*8,DIMENSION(:),ALLOCATABLE :: presnivs_omp!(llm) ! approximate pressure of atm. layers
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zufi_omp !pu !zufi_omp(klon,llm) ! zonal wind (m/s)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zvfi_omp !pv !zvfi_omp(klon,llm) ! meridional wind (m/s)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zrfi_omp!(klon,llm) ! relative wind vorticity, in s-1
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: ztfi_omp !pt !ztfi_omp(klon,llm) ! temperature (K)
+REAL*8,DIMENSION(:,:,:),ALLOCATABLE :: zqfi_omp !pq !zqfi_omp(klon,llm,nqtot) ! tracers (*/kg of air)
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: flxwfi_omp !flxw !flxwfi_omp(klon,llm) ! Vertical mass flux on lower mesh interfaces (kg/s)  ! tendencies (in */s) from the physics:
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zdufi_omp !pdu !zdufi_omp(klon,llm) ! tendency on zonal winds
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zdvfi_omp !pdv !zdvfi_omp(klon,llm) ! tendency on meridional winds
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: zdtfi_omp !pdt !zdtfi_omp(klon,llm) ! tendency on temperature
+REAL*8,DIMENSION(:,:,:),ALLOCATABLE :: zdqfi_omp !pdq !zdqfi_omp(klon,llm,nqtot) ! tendency on tracers
+REAL*8,DIMENSION(:),ALLOCATABLE :: zdpsrf_omp !pdpsrf !zdpsrf_omp(klon) !tendency on surface pressure
+REAL*8,DIMENSION(:,:),ALLOCATABLE :: hrlw,hrsw
 
 CONTAINS
@@ -8,13 +32,53 @@
 
   IMPLICIT NONE
- 
+
   INTEGER,INTENT(IN) :: klon ! (local) number of atmospheric columns
   INTEGER,INTENT(IN) :: llm  ! number of atmospheric layers
   INTEGER,INTENT(IN) :: nqtot ! number of tracers
- 
+
+  IF (.NOT.ALLOCATED(zplev_omp)) ALLOCATE(zplev_omp(klon,llm+1))
+  IF (.NOT.ALLOCATED(zplay_omp)) ALLOCATE(zplay_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zpk_omp)) ALLOCATE(zpk_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zphi_omp)) ALLOCATE(zphi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zphis_omp)) ALLOCATE(zphis_omp(klon))
+  IF (.NOT.ALLOCATED(presnivs_omp)) ALLOCATE(presnivs_omp(llm))
+  IF (.NOT.ALLOCATED(zufi_omp)) ALLOCATE(zufi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zvfi_omp)) ALLOCATE(zvfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zrfi_omp)) ALLOCATE(zrfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(ztfi_omp)) ALLOCATE(ztfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zqfi_omp)) ALLOCATE(zqfi_omp(klon,llm,nqtot))
+  IF (.NOT.ALLOCATED(flxwfi_omp)) ALLOCATE(flxwfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdufi_omp)) ALLOCATE(zdufi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdvfi_omp)) ALLOCATE(zdvfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdtfi_omp)) ALLOCATE(zdtfi_omp(klon,llm))
+  IF (.NOT.ALLOCATED(zdqfi_omp)) ALLOCATE(zdqfi_omp(klon,llm,nqtot))
+  IF (.NOT.ALLOCATED(zdpsrf_omp)) ALLOCATE(zdpsrf_omp(klon))
+  IF (.NOT.ALLOCATED(hrlw)) ALLOCATE(hrlw(klon,llm))
+  IF (.NOT.ALLOCATED(hrsw)) ALLOCATE(hrsw(klon,llm))
+
 END SUBROUTINE allocate_interface
 
 SUBROUTINE deallocate_interface
 
+  DEALLOCATE(zplev_omp)
+  DEALLOCATE(zplay_omp)
+  DEALLOCATE(zpk_omp)
+  DEALLOCATE(zphi_omp)
+  DEALLOCATE(zphis_omp)
+  DEALLOCATE(presnivs_omp)
+  DEALLOCATE(zufi_omp)
+  DEALLOCATE(zvfi_omp)
+  DEALLOCATE(zrfi_omp)
+  DEALLOCATE(ztfi_omp)
+  DEALLOCATE(zqfi_omp)
+  DEALLOCATE(flxwfi_omp)
+  DEALLOCATE(zdufi_omp)
+  DEALLOCATE(zdvfi_omp)
+  DEALLOCATE(zdtfi_omp)
+  DEALLOCATE(zdqfi_omp)
+  DEALLOCATE(zdpsrf_omp)
+  DEALLOCATE(hrlw)
+  DEALLOCATE(hrsw)
+
 END SUBROUTINE deallocate_interface