My Project
 All Classes Files Functions Variables Macros
convmas.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE convmas (pbaru, pbarv, convm )
5 c
6  IMPLICIT NONE
7 
8 c=======================================================================
9 c
10 c Auteurs: P. Le Van , F. Hourdin .
11 c -------
12 c
13 c Objet:
14 c ------
15 c
16 c ********************************************************************
17 c .... calcul de la convergence du flux de masse aux niveaux p ...
18 c ********************************************************************
19 c
20 c
21 c pbaru et pbarv sont des arguments d'entree pour le s-pg ....
22 c ..... convm est un argument de sortie pour le s-pg ....
23 c
24 c le calcul se fait de haut en bas,
25 c la convergence de masse au niveau p(llm+1) est egale a 0. et
26 c n'est pas stockee dans le tableau convm .
27 c
28 c
29 c=======================================================================
30 c
31 c Declarations:
32 c -------------
33 
34 #include "dimensions.h"
35 #include "paramet.h"
36 #include "comvert.h"
37 #include "logic.h"
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