Index: trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90	(revision 1875)
+++ trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90	(revision 1877)
@@ -236,4 +236,5 @@
       real zh(ngrid,nlayer)               ! Potential temperature (K).
       real pw(ngrid,nlayer)               ! Vertical velocity (m/s). (NOTE : >0 WHEN DOWNWARDS !!)
+      real omega(ngrid,nlayer)            ! omega velocity (Pa/s, >0 when downward)
 
       integer l,ig,ierr,iq,nw,isoil
@@ -757,4 +758,12 @@
                        (pplay(1:ngrid,l)*cell_area(1:ngrid))
       enddo
+      ! omega in Pa/s
+      do l=1,nlayer-1
+         omega(1:ngrid,l)=0.5*(flxw(1:ngrid,l)+flxw(1:ngrid,l+1))
+      enddo
+      omega(1:ngrid,nlayer)=0.5*flxw(1:ngrid,nlayer) ! since flxw(nlayer+1)=0
+      do l=1,nlayer
+         omega(1:ngrid,l)=g*omega(1:ngrid,l)/cell_area(1:ngrid)
+      enddo
 
       ! ----------------------------------------------------------------
@@ -2172,4 +2181,8 @@
       CALL send_xios_field("u",zu)
       CALL send_xios_field("v",zv)
+      CALL send_xios_field("omega",omega)
+
+      CALL send_xios_field("ISR",fluxtop_dn)
+      CALL send_xios_field("OLR",fluxtop_lw)
 
       if (lastcall.and.is_omp_master) then