trunk/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
      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			IMPLICIT NONE
11
12			c=======================================================================
13			c
14			c Auteurs: P. Le Van , F. Hourdin .
15			c -------
16			c
17			c Objet:
18			c ------
19			c
20			c ********************************************************************
21			c .... calcul de la convergence du flux de masse aux niveaux p ...
22			c ********************************************************************
23			c
24			c
25			c pbaru et pbarv sont des arguments d'entree pour le s-pg ....
26			c ..... convm est un argument de sortie pour le s-pg ....
27			c
28			c le calcul se fait de haut en bas,
29			c la convergence de masse au niveau p(llm+1) est egale a 0. et
30			c n'est pas stockee dans le tableau convm .
31			c
32			c
33			c=======================================================================
34			c
35			c Declarations:
36			c -------------
37
38
39			REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm( ip1jmp1,llm )
40			INTEGER l,ij
41
42
43			c-----------------------------------------------------------------------
44			c .... calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
45
46			CALL convflu( pbaru, pbarv, llm, convm )
47
48			c-----------------------------------------------------------------------
49			c filtrage:
50			c ---------
51
52			CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )
53
54			c integration de la convergence de masse de haut en bas ......
55
56			DO l = llmm1, 1, -1
57			DO ij = 1, ip1jmp1
58			convm(ij,l) = convm(ij,l) + convm(ij,l+1)
59			ENDDO
60			ENDDO
61			c
62			RETURN
63			END