source: LMDZ.3.3/trunk/libf/dyn3d/convmas.F @ 61

Last change on this file since 61 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.6 KB
Line 
1      SUBROUTINE convmas (pbaru, pbarv, convm )
2c
3      IMPLICIT NONE
4
5c=======================================================================
6c
7c   Auteurs:  P. Le Van , F. Hourdin  .
8c   -------
9c
10c   Objet:
11c   ------
12c
13c   ********************************************************************
14c   .... calcul de la convergence du flux de masse aux niveaux p ...
15c   ********************************************************************
16c
17c
18c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....
19c      .....  convm      est  un argument de sortie pour le s-pg  ....
20c
21c    le calcul se fait de haut en bas,
22c    la convergence de masse au niveau p(llm+1) est egale a 0. et
23c    n'est pas stockee dans le tableau convm .
24c
25c
26c=======================================================================
27c
28c   Declarations:
29c   -------------
30
31#include "dimensions.h"
32#include "paramet.h"
33#include "comvert.h"
34#include "logic.h"
35
36      REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm(  ip1jmp1,llm )
37      INTEGER   l,ij
38
39      EXTERNAL   filtreg
40      EXTERNAL   convflu
41
42c-----------------------------------------------------------------------
43c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
44
45      CALL  convflu( pbaru, pbarv, llm, convm )
46
47c-----------------------------------------------------------------------
48c   filtrage:
49c   ---------
50
51       CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )
52
53c    integration de la convergence de masse de haut  en bas ......
54
55      DO      l      = llmm1, 1, -1
56        DO    ij     = 1, ip1jmp1
57         convm(ij,l) = convm(ij,l) + convm(ij,l+1)
58        ENDDO
59      ENDDO
60c
61      RETURN
62      END
Note: See TracBrowser for help on using the repository browser.