source: LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar/convmas2_p.F @ 1410

Last change on this file since 1410 was 763, checked in by Laurent Fairhead, 17 years ago

Merge entre la version V3_conv et le HEAD
Nouvelles routines parallélisées
YM, JG, LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.4 KB
Line 
1      SUBROUTINE convmas2_p ( 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 )
38      REAL :: convm(  ip1jmp1,llm )
39      INTEGER   l,ij
40      INTEGER ijb,ije,jjb,jje
41 
42c$OMP MASTER
43c    integration de la convergence de masse de haut  en bas ......
44       ijb=ij_begin
45       ije=ij_end+iip1
46       if (pole_sud) ije=ij_end
47           
48      DO      l      = llmm1, 1, -1
49        DO    ij     = ijb, ije
50         convm(ij,l) = convm(ij,l) + convm(ij,l+1)
51        ENDDO
52      ENDDO
53c
54c$OMP END MASTER
55      RETURN
56      END
Note: See TracBrowser for help on using the repository browser.