source: LMDZ5/trunk/libf/dyn3dpar/convmas_p.F @ 4735

Last change on this file since 4735 was 2603, checked in by Ehouarn Millour, 8 years ago

Cleanup in the dynamics: turn logic.h into module logic_mod.F90
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.8 KB
Line 
1      SUBROUTINE convmas_p (pbaru, pbarv, convm )
2c
3      USE parallel_lmdz
4      IMPLICIT NONE
5
6c=======================================================================
7c
8c   Auteurs:  P. Le Van , F. Hourdin  .
9c   -------
10c
11c   Objet:
12c   ------
13c
14c   ********************************************************************
15c   .... calcul de la convergence du flux de masse aux niveaux p ...
16c   ********************************************************************
17c
18c
19c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....
20c      .....  convm      est  un argument de sortie pour le s-pg  ....
21c
22c    le calcul se fait de haut en bas,
23c    la convergence de masse au niveau p(llm+1) est egale a 0. et
24c    n'est pas stockee dans le tableau convm .
25c
26c
27c=======================================================================
28c
29c   Declarations:
30c   -------------
31
32#include "dimensions.h"
33#include "paramet.h"
34
35      REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm )
36      REAL, target :: convm(  ip1jmp1,llm )
37      INTEGER   l,ij
38
39      INTEGER ijb,ije,jjb,jje
40 
41     
42c-----------------------------------------------------------------------
43c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
44
45      CALL  convflu_p( pbaru, pbarv, llm, convm )
46
47c-----------------------------------------------------------------------
48c   filtrage:
49c   ---------
50       
51       jjb=jj_begin
52       jje=jj_end+1
53       if (pole_sud) jje=jj_end
54 
55       CALL filtreg_p( convm, jjb, jje, jjp1, llm, 2, 2, .true., 1 )
56
57c    integration de la convergence de masse de haut  en bas ......
58       ijb=ij_begin
59       ije=ij_end+iip1
60       if (pole_sud) ije=ij_end
61           
62      DO      l      = llmm1, 1, -1
63        DO    ij     = ijb, ije
64         convm(ij,l) = convm(ij,l) + convm(ij,l+1)
65        ENDDO
66      ENDDO
67c
68      RETURN
69      END
Note: See TracBrowser for help on using the repository browser.