Changeset 1877 for trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

Timestamp:

Dec 21, 2017, 11:29:56 PM (7 years ago)

Author:

aslmd

Message:

addition of the omega velocity in physiq_mod (awaiting improvements by JL)

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -236 +236 @@
       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 +758 @@
                        (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 +2181 @@
       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