Changeset 1346 for trunk/LMDZ.MARS/libf

Timestamp:

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

Author:

emillour

Message:

Mars and Generic GCM:

• Bug correction on the computation of vertical velocity.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/physiq.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -580 +580 @@
 
       ! 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
@@ -587 +587 @@
       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))
-       ! NB: here we use r and nor rnew since this diagnostic comes
+       pw(1:ngrid,l)=(pw(1:ngrid,l)*r*pt(1:ngrid,l)) / 
+     &               (pplay(1:ngrid,l)*area(1:ngrid))
+       ! NB: here we use r and not rnew since this diagnostic comes
        ! from the dynamics
       enddo