Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys	(revision 2742)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys	(revision 2743)
@@ -108,9 +108,9 @@
 state  real  GSNOW  ikj  misc  1  -  -  "SNOW"     "SNOW kg/kg.s-1"           ""
 state  real  ALBEQ  ij misc  1  -  -  "ALBEQ"    "ALBEDO EQUIVALENT"   ""
-state  real  FLUXTOP_DN ij misc  1  -  - "FLUXTOP_DN"   "Incoming SW radiation at the TOA" "W.m-2"
-state  real  FLUXABS_SW ij misc  1  -  - "FLUXABS_SW"  "Absorbed SW flux" "W.m-2"
-state  real  FLUXTOP_LW ij misc  1  -  - "FLUXTOP_LW"  "Outgoing LW flux TOA" "W.m-2"
-state  real  FLUXSURF_SW ij misc  1  -  - "FLUXSURF_SW" "SW Surface flux" "W.m-2"
-state  real  FLUXSURF_LW ij misc  1  -  -  "FLUXSURF_LW" "LW Surface flux" "W.m-2"
+state  real  FLUXTOP_DN ij misc  1  -  h "FLUXTOP_DN"   "Incoming SW radiation at the TOA" "W.m-2"
+state  real  FLUXABS_SW ij misc  1  -  h "FLUXABS_SW"  "Absorbed SW flux" "W.m-2"
+state  real  FLUXTOP_LW ij misc  1  -  h "FLUXTOP_LW"  "Outgoing LW flux TOA" "W.m-2"
+state  real  FLUXSURF_SW ij misc  1  -  h "FLUXSURF_SW" "SW Surface flux" "W.m-2"
+state  real  FLUXSURF_LW ij misc  1  -  h  "FLUXSURF_LW" "LW Surface flux" "W.m-2"
 state  real  FLXGRD ij misc  1  -  h "FLXGRD" "heat flux from ground"  "W.m-2"
 state  real  LSCEZ ikj misc  1  -  - "LSCEZ" "heat from largescale" "W.m-2"
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F	(revision 2742)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F	(revision 2743)
@@ -176,5 +176,5 @@
    END SELECT
 
-   IF (planet == "mars") THEN
+   IF (planet == "mars" .or. planet == "titan") THEN
      stretch_grid = .false.
      !! FOR LES, set stretch to false
@@ -259,5 +259,5 @@
      !!MARS
      !!MARS
-     IF (planet .ne. "mars") THEN
+     IF (planet .ne. "mars" .and. planet .ne. "titan") THEN
        open(unit=12,file='levels',form='formatted',status='old')
        rewind(12)
@@ -276,5 +276,5 @@
 
    !! SPECIFIC FOR LES PBL MARS
-   IF (planet == "mars") THEN
+   IF (planet == "mars" .or. planet == "titan") THEN
      !!!MARS
      grid%znw(1)=1.000
@@ -286,4 +286,13 @@
      ENDDO
    ENDIF
+
+   open(unit=12,file='levels',form='formatted',status='old')
+   rewind(12)
+   DO k=1, kde
+     write(12,*) grid%znw(k)
+     write(6,*) 'for update_inputs_physiq_mod (e.g. Titan, generic)',k,grid%znw(k)
+   ENDDO
+   close(12)
+
 
 
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/update_inputs_physiq_mod.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/update_inputs_physiq_mod.F	(revision 2742)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/update_inputs_physiq_mod.F	(revision 2743)
@@ -187,4 +187,6 @@
    use comgeomfi_h, only: totarea, totarea_planet
    use planete_mod, only: ini_planete_mod
+   ! in LMD-WRF interface
+   USE variables_mod, only: ptopwrf
 
    INTEGER, INTENT(IN) :: ims,ime,jms,jme
@@ -198,5 +200,4 @@
    REAL, DIMENSION(nlayer+1) :: znw
    REAL*8, DIMENSION(nlayer+1) :: apdyn,bpdyn
-   REAL :: ptop
    REAL*8 :: PP0
    DO j = jps,jpe
@@ -280,6 +281,5 @@
    close(12)
    
-   ptop=0.5
-   apdyn=ptop*(1-znw)
+   apdyn=ptopwrf*(1-znw)
    bpdyn=znw
    PP0=p0
Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/variables_mod.F
===================================================================
--- trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/variables_mod.F	(revision 2742)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/variables_mod.F	(revision 2743)
@@ -7,4 +7,5 @@
 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 :: ptopwrf ! pressure at top of model in WRF [should be equal to grid%p_top]
 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)
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 2742)
+++ trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new	(revision 2743)
@@ -264,4 +264,5 @@
 wappel_phys = RADT
 zdt_split = dt*wappel_phys            ! physical timestep (s)
+ptopwrf = ptop                        ! top pressure in WRF coordinates (in variables_mod)
 ngrid=(ipe-ips+1)*(jpe-jps+1)         ! size of the horizontal grid
 nlayer = kpe-kps+1                    ! number of vertical layers: nlayermx