!

! $Id: caldyn.F 2600 2016-07-23 05:45:38Z emillour $

!

     $ (itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,

     $  phi,conser,du,dv,dteta,dp,w,pbaru,pbarv,time )

      USE comvert_mod, ONLY: ap, bp

      IMPLICIT NONE

!=======================================================================

!

!  Auteur :  P. Le Van

!

!   Objet:

!   ------

!

!   Calcul des tendances dynamiques.

!

! Modif 04/93 F.Forget

!=======================================================================

!-----------------------------------------------------------------------

!   0. Declarations:

!   ----------------

      include "dimensions.h"

      include "paramet.h"

      include "comgeom.h"

!   Arguments:

!   ----------

      LOGICAL,INTENT(IN) :: conser ! triggers printing some diagnostics

      INTEGER,INTENT(IN) :: itau ! time step index

      REAL,INTENT(IN) :: vcov(ip1jm,llm) ! covariant meridional wind

      REAL,INTENT(IN) :: ucov(ip1jmp1,llm) ! covariant zonal wind

      REAL,INTENT(IN) :: teta(ip1jmp1,llm) ! potential temperature

      REAL,INTENT(IN) :: ps(ip1jmp1) ! surface pressure

      REAL,INTENT(IN) :: phis(ip1jmp1) ! geopotential at the surface

      REAL,INTENT(IN) :: pk(ip1jmp1,llm) ! Exner at mid-layer

      REAL,INTENT(IN) :: pkf(ip1jmp1,llm) ! filtered Exner

      REAL,INTENT(IN) :: phi(ip1jmp1,llm) ! geopotential

      REAL,INTENT(OUT) :: masse(ip1jmp1,llm) ! air mass

      REAL,INTENT(OUT) :: dv(ip1jm,llm) ! tendency on vcov

      REAL,INTENT(OUT) :: du(ip1jmp1,llm) ! tendency on ucov

      REAL,INTENT(OUT) :: dteta(ip1jmp1,llm) ! tenddency on teta

      REAL,INTENT(OUT) :: dp(ip1jmp1) ! tendency on ps

      REAL,INTENT(OUT) :: w(ip1jmp1,llm) ! vertical velocity

      REAL,INTENT(OUT) :: pbaru(ip1jmp1,llm) ! mass flux in the zonal direction

      REAL,INTENT(OUT) :: pbarv(ip1jm,llm) ! mass flux in the meridional direction

      REAL,INTENT(IN) :: time ! current time

!   Local:

!   ------

      REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)

      REAL ang(ip1jmp1,llm),p(ip1jmp1,llmp1)

      REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm),psexbarxy(ip1jm)

      REAL vorpot(ip1jm,llm)

      REAL ecin(ip1jmp1,llm),convm(ip1jmp1,llm)

      REAL bern(ip1jmp1,llm)

      REAL massebxy(ip1jm,llm)

      INTEGER   ij,l

!-----------------------------------------------------------------------

!   Compute dynamical tendencies:

!--------------------------------

      ! compute contravariant winds ucont() and vcont

      ! compute pressure p()

      ! compute psexbarxy() XY-area weighted-averaged surface pressure (what for?)

      ! compute mass in each atmospheric mesh: masse()

      ! compute X and Y-averages of mass, massebx() and masseby()

      ! compute XY-average of mass, massebxy()

      ! compute mass fluxes pbaru() and pbarv()

      ! compute dteta() , horizontal converging flux of theta

      ! compute convm(), horizontal converging flux of mass

      ! compute pressure variation due to mass convergence

      ENDDO

      ! compute vertical velocity w()

      ! compute potential vorticity vorpot()

      ! compute rotation induced du() and dv()

      ! compute kinetic energy ecin()

      ! compute Bernouilli function bern()

      ! compute and add du() and dv() contributions from Bernouilli and pressure

         ENDDO

      ENDDO

      ! compute vertical advection contributions to du(), dv() and dteta()

!  WARNING probleme de peridocite de dv sur les PC/linux. Pb d'arrondi

!          probablement. Observe sur le code compile avec pgf90 3.0-1

!         PRINT *,'!!!ATTENTION!!! probleme de periodicite sur vcov',

!    ,   ' dans caldyn'

!         PRINT *,' l,  ij = ', l, ij, ij+iim,dv(ij+iim,l),dv(ij,l)

           ENDIF

         ENDDO

      ENDDO

!-----------------------------------------------------------------------

!   Output some control variables:

!---------------------------------

        CALL sortvarc

      ENDIF