Sources/dyn3d/convmas.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convmas.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
!
      SUBROUTINE convmas (pbaru, pbarv, convm )
c
      use dimens_m
      use paramet_m
      use comvert
      use logic
      use filtreg_m, only: filtreg
      IMPLICIT NONE

c=======================================================================
c
c   Auteurs:  P. Le Van , F. Hourdin  .
c   -------
c
c   Objet:
c   ------
c
c   ********************************************************************
c   .... calcul de la convergence du flux de masse aux niveaux p ...
c   ********************************************************************
c
c
c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....
c      .....  convm      est  un argument de sortie pour le s-pg  ....
c
c    le calcul se fait de haut en bas, 
c    la convergence de masse au niveau p(llm+1) est egale a 0. et
c    n'est pas stockee dans le tableau convm .
c
c
c=======================================================================
c
c   Declarations:
c   -------------


      REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm(  ip1jmp1,llm )
      INTEGER   l,ij


c-----------------------------------------------------------------------
c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......

      CALL  convflu( pbaru, pbarv, llm, convm )

c-----------------------------------------------------------------------
c   filtrage:
c   ---------

       CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )

c    integration de la convergence de masse de haut  en bas ......

      DO      l      = llmm1, 1, -1
        DO    ij     = 1, ip1jmp1
         convm(ij,l) = convm(ij,l) + convm(ij,l+1)
        ENDDO
      ENDDO
c
      RETURN
      END
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convmas.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
3			!
4			SUBROUTINE convmas (pbaru, pbarv, convm )
5			c
6			use dimens_m
7			use paramet_m
8			use comvert
9			use logic
10	guez	27	use filtreg_m, only: filtreg
11	guez	3	IMPLICIT NONE
12
13			c=======================================================================
14			c
15			c Auteurs: P. Le Van , F. Hourdin .
16			c -------
17			c
18			c Objet:
19			c ------
20			c
21			c ********************************************************************
22			c .... calcul de la convergence du flux de masse aux niveaux p ...
23			c ********************************************************************
24			c
25			c
26			c pbaru et pbarv sont des arguments d'entree pour le s-pg ....
27			c ..... convm est un argument de sortie pour le s-pg ....
28			c
29			c le calcul se fait de haut en bas,
30			c la convergence de masse au niveau p(llm+1) est egale a 0. et
31			c n'est pas stockee dans le tableau convm .
32			c
33			c
34			c=======================================================================
35			c
36			c Declarations:
37			c -------------
38
39
40			REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm( ip1jmp1,llm )
41			INTEGER l,ij
42
43
44			c-----------------------------------------------------------------------
45			c .... calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
46
47			CALL convflu( pbaru, pbarv, llm, convm )
48
49			c-----------------------------------------------------------------------
50			c filtrage:
51			c ---------
52
53			CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )
54
55			c integration de la convergence de masse de haut en bas ......
56
57			DO l = llmm1, 1, -1
58			DO ij = 1, ip1jmp1
59			convm(ij,l) = convm(ij,l) + convm(ij,l+1)
60			ENDDO
61			ENDDO
62			c
63			RETURN
64			END