Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2332)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2333)
@@ -697,5 +697,5 @@
         write(*,*) "physiq firstcall: call initialize_xios_output"
         call initialize_xios_output(pday,ptime,ptimestep,daysec,
-     &                              presnivs,pseudoalt)
+     &                              presnivs,pseudoalt,mlayer)
 #endif
       ENDIF        !  (end of "if firstcall")
@@ -3745,8 +3745,21 @@
       CALL send_xios_field("ps",ps)
       CALL send_xios_field("area",cell_area)
+
+      CALL send_xios_field("ISR",fluxtop_sw_tot)
+      CALL send_xios_field("OLR",fluxtop_lw)
+
+      CALL send_xios_field("tsurf",tsurf)
+      CALL send_xios_field("inertiedat",inertiedat)
+      CALL send_xios_field("tsoil",tsoil)
+      CALL send_xios_field("co2ice",co2ice)
       
       CALL send_xios_field("temperature",zt)
       CALL send_xios_field("u",zu)
       CALL send_xios_field("v",zv)
+      !CREATE IF CO2CYCLE
+      !if (tracer.and.(igcm_co2.ne.0)) then
+      !   CALL send_xios_field("co2",zq(1,1,igcm_co2))
+      !endif
+
 
       if (lastcall.and.is_omp_master) then
Index: trunk/LMDZ.MARS/libf/phymars/xios_output_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/xios_output_mod.F90	(revision 2332)
+++ trunk/LMDZ.MARS/libf/phymars/xios_output_mod.F90	(revision 2333)
@@ -18,5 +18,5 @@
 
   SUBROUTINE initialize_xios_output(day,timeofday,dtphys,daysec,&
-                                    presnivs,pseudoalt)
+                                    presnivs,pseudoalt,mlayer)
   USE mod_phys_lmdz_para, ONLY: jj_nb, jj_begin, jj_end, ii_begin, ii_end, &
                                 mpi_size, mpi_rank, klon_mpi, &
@@ -39,4 +39,5 @@
   REAL,INTENT(IN) :: presnivs(:) ! vertical grid approximate pressure (Pa)
   REAL,INTENT(IN) :: pseudoalt(:) ! vertical grid approximate altitude (km)
+  REAL,INTENT(IN) :: mlayer(:) ! soil layer depth at intermediate level (m)
   
   
@@ -56,4 +57,7 @@
     CALL xios_set_axis_attr("altitude", n_glo=size(pseudoalt), value=pseudoalt,&
                             unit="km",positive="up")
+    IF (prt_level>=10) WRITE(lunout,*) "initialize_xios_output: call xios_set_axis_attr for soil"
+    CALL xios_set_axis_attr("soil_layers", n_glo=size(mlayer), value=mlayer,&
+                            unit="m",positive="down")
     
     ! 2. Declare horizontal domain