Changeset 1312 for trunk/LMDZ.GENERIC/libf

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

Location:

trunk/LMDZ.GENERIC/libf

Files:

• dyn3d/calfis.F (modified) (3 diffs)
• phystd/physiq.F90 (modified) (5 diffs)

trunk/LMDZ.GENERIC/libf/dyn3d/calfis.F

@@ -43 +43 @@
 c       pmasse          masse d'air dans chaque maille
 c       pts             surface temperature  (K)
-c       pw              flux vertical (kg m-2)
+c       pw              flux vertical (kg/s)
 c
 c    Output :
@@ -308 +308 @@
 c On interpole "lineairement" la temperature entre les couches(FF,10/95)
 
-      DO ig=1,ngridmx
-         zvervel(ig,1)=0.
-      END DO
-      DO l=2,llm
-        zvervel(1,l)=(pw(1,1,l)/apoln)
-     &  * r *0.5*(ztfi(1,l)+ztfi(1,l-1)) /zplev(1,l)              
-        ig0=2
-       DO j=2,jjm
-           DO i = 1, iim
-              zvervel(ig0,l) = pw(i,j,l) * unsaire(i,j)
-     &        * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
-              ig0 = ig0 + 1
-           ENDDO
-       ENDDO
-        zvervel(ig0,l)=(pw(1,jjp1,l)/apols)
-     &  * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
-      ENDDO
+!      DO ig=1,ngridmx
+!         zvervel(ig,1)=0.
+!      END DO
+!      DO l=2,llm
+!        zvervel(1,l)=(pw(1,1,l)/apoln)
+!     &  * r *0.5*(ztfi(1,l)+ztfi(1,l-1)) /zplev(1,l)              
+!        ig0=2
+!       DO j=2,jjm
+!           DO i = 1, iim
+!              zvervel(ig0,l) = pw(i,j,l) * unsaire(i,j)
+!     &        * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
+!              ig0 = ig0 + 1
+!           ENDDO
+!       ENDDO
+!        zvervel(ig0,l)=(pw(1,jjp1,l)/apols)
+!     &  * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
+!      ENDDO
 
 c    .........  Reindexation : calcul de zvervel au MILIEU des couches
-       DO l=1,llm-1
-              DO ig=1,ngridmx
-                     zvervel(ig,l) = 0.5*(zvervel(ig,l)+zvervel(ig,l+1))
-          END DO 
-       END DO 
+!       DO l=1,llm-1
+!             DO ig=1,ngridmx
+!                    zvervel(ig,l) = 0.5*(zvervel(ig,l)+zvervel(ig,l+1))
+!          END DO 
+!       END DO 
 c      (dans la couche llm, on garde la valeur à la limite inférieure llm)
 
@@ -430 +430 @@
      ,     zplev,zplay,zphi,
      ,     zufi, zvfi,ztfi, zqfi,  
-     ,     zvervel,
+!     ,     zvervel,
+     ,     pw,
 C - sorties
      s     zdufi, zdvfi, zdtfi, zdqfi,zdpsrf,tracer)

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
 
 !-----------------------------------------------------------------------