Changeset 1346

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

Location:

trunk

Files:

• LMDZ.GENERIC/libf/phystd/physiq.F90 (modified) (2 diffs)
• LMDZ.MARS/libf/phymars/physiq.F (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
 

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