source: LMDZ5/branches/LF-private/libf/dyn3dmem/convmas2_loc.F @ 5440

Last change on this file since 5440 was 1823, checked in by Ehouarn Millour, 11 years ago

Remplacement de parallel.F90 (en conflit avec orchidée) par parallel_lmdz.F90.
UG
.........................................
Renaming parallel.F90 (conflicting with orchidée) into parallel_lmdz.F90.
UG

File size: 1.4 KB
Line 
1      SUBROUTINE convmas2_loc ( 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#include "comvert.h"
35#include "logic.h"
36
37      REAL pbaru( ijb_u:ije_u,llm ),pbarv( ijb_v:ije_v,llm )
38      REAL :: convm(  ijb_u:ije_u,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.