Changeset 1346 for trunk/LMDZ.GENERIC

Timestamp:

Sep 10, 2014, 2:28:04 PM (10 years ago)

Author:

emillour

Message:

Mars and Generic GCM:

• Bug correction on the computation of vertical velocity.

EM

File:

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

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

@@ -833 +833 @@
 
      ! Compute vertical velocity (m/s) from vertical mass flux
-     ! w = F / (rho*area) and rho = r*T/P
+     ! w = F / (rho*area) and rho = P/(r*T)
      ! but first linearly interpolate mass flux to mid-layers
      do l=1,nlayer-1
@@ -840 +840 @@
      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))
+       pw(1:ngrid,l)=(pw(1:ngrid,l)*r*pt(1:ngrid,l)) /  &
+                     (pplay(1:ngrid,l)*area(1:ngrid))
      enddo