source: LMDZ5/branches/testing/libf/dyn3dpar/convmas.F @ 1663

Last change on this file since 1663 was 1000, checked in by Laurent Fairhead, 16 years ago
  • Modifs sur le parallelisme: masquage dans la physique
  • Inclusion strato
  • mise en coherence etat0
  • le mode offline fonctionne maintenant en parallele,
  • les fichiers de la dynamiques sont correctement sortis et peuvent etre reconstruit avec rebuild
  • la version parallele de la dynamique peut s'executer sans MPI (sur 1 proc)
  • L'OPENMP fonctionne maintenant sans la parallelisation MPI.

YM
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.6 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE convmas (pbaru, pbarv, convm )
5c
6      IMPLICIT NONE
7
8c=======================================================================
9c
10c   Auteurs:  P. Le Van , F. Hourdin  .
11c   -------
12c
13c   Objet:
14c   ------
15c
16c   ********************************************************************
17c   .... calcul de la convergence du flux de masse aux niveaux p ...
18c   ********************************************************************
19c
20c
21c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....
22c      .....  convm      est  un argument de sortie pour le s-pg  ....
23c
24c    le calcul se fait de haut en bas,
25c    la convergence de masse au niveau p(llm+1) est egale a 0. et
26c    n'est pas stockee dans le tableau convm .
27c
28c
29c=======================================================================
30c
31c   Declarations:
32c   -------------
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
43c-----------------------------------------------------------------------
44c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
45
46      CALL  convflu( pbaru, pbarv, llm, convm )
47
48c-----------------------------------------------------------------------
49c   filtrage:
50c   ---------
51
52       CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )
53
54c    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
61c
62      RETURN
63      END
Note: See TracBrowser for help on using the repository browser.