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

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/dyn3d/calfis.F (modified) (3 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)