source: LMDZ4/trunk/libf/dyn3dpar/convmas_p.F @ 706

Last change on this file since 706 was 630, checked in by Laurent Fairhead, 20 years ago

Import d'une version parallele de la dynamique YM
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.9 KB
RevLine 
[630]1      SUBROUTINE convmas_p (pbaru, pbarv, convm )
2c
3      USE parallel
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#include "comvert.h"
35#include "logic.h"
36
37      REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm(  ip1jmp1,llm )
38      INTEGER   l,ij
39
40      EXTERNAL   filtreg_p
41      EXTERNAL   convflu_p
42     
43      INTEGER ijb,ije,jjb,jje
44     
45     
46
47c-----------------------------------------------------------------------
48c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
49
50      CALL  convflu_p( pbaru, pbarv, llm, convm )
51
52c-----------------------------------------------------------------------
53c   filtrage:
54c   ---------
55       
56       jjb=jj_begin
57       jje=jj_end+1
58       if (pole_sud) jje=jj_end
59 
60       CALL filtreg_p( convm, jjb, jje, jjp1, llm, 2, 2, .true., 1 )
61
62c    integration de la convergence de masse de haut  en bas ......
63       ijb=ij_begin
64       ije=ij_end+iip1
65       if (pole_sud) ije=ij_end
66      DO      l      = llmm1, 1, -1
67        DO    ij     = ijb, ije
68         convm(ij,l) = convm(ij,l) + convm(ij,l+1)
69        ENDDO
70      ENDDO
71c
72      RETURN
73      END
Note: See TracBrowser for help on using the repository browser.