Changeset 1312 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Jul 11, 2014, 11:48:09 AM (11 years ago)

Author:

emillour

Message:

Generic and Mars GCM:
LMDZ.COMMON dynamics sends mass flux to physics and not vertical velocity.
Add the computation of vertical velocity from input mass flux in the physics,
and also modify "old" LMDZ.GENERIC and LMDZ.MARS dynamics to be consistent.
EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq.F90 (modified) (5 diffs)

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

@@ -5 +5 @@
                   pplev,pplay,pphi,        &
                   pu,pv,pt,pq,             &
-                  pw,                      &
+                  flxw,                    &
                   pdu,pdv,pdt,pdq,pdpsrf,tracerdyn)
  
@@ -98 +98 @@
 !    ptdyn(ngrid,nlayer)     / corresponding variables
 !    pqdyn(ngrid,nlayer,nq) /
-!    pw(ngrid,?)           vertical velocity
+!    flxw(ngrid,nlayer)      vertical mass flux (kg/s) at layer lower boundary
 !
 !   output
@@ -161 +161 @@
       real,intent(in) :: pt(ngrid,nlayer) ! temperature (K)
       real,intent(in) :: pq(ngrid,nlayer,nq) ! tracers (.../kg_of_air)
-      real,intent(in) :: pw(ngrid,nlayer)    ! vertical velocity (m/s)
-
+      real,intent(in) :: flxw(ngrid,nlayer) ! vertical mass flux (ks/s)
+                                            ! at lower boundary of layer
 
 
@@ -210 +210 @@
       real aerosol(ngrid,nlayer,naerkind) 
       real zh(ngrid,nlayer)      ! potential temperature (K)
+      real pw(ngrid,nlayer) ! vertical velocity (m/s) (>0 when downwards)
 
       character*80 fichier 
@@ -803 +804 @@
       enddo
 
+     ! Compute vertical velocity (m/s) from vertical mass flux
+     ! w = F / (rho*area) and rho = r*T/P
+     ! but first linearly interpolate mass flux to mid-layers
+     do l=1,nlayer-1
+       pw(1:ngrid,l)=0.5*(flxw(1:ngrid,l)+flxw(1:ngrid,l+1))
+     enddo
+     pw(1:ngrid,nlayer)=0.5*flxw(1:ngrid,nlayer) ! since flxw(nlayer+1)=0
+     do l=1,nlayer
+       pw(1:ngrid,l)=(pw(1:ngrid,l)*pplay(1:ngrid,l)) /  &
+                     (r*pt(1:ngrid,l)*area(1:ngrid))
+     enddo
 
 !-----------------------------------------------------------------------