source: LMDZ5/trunk/libf/dyn3dmem/bernoui_loc.F @ 5362

Last change on this file since 5362 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
File size: 1.8 KB
Line 
1      SUBROUTINE bernoui_loc (ngrid,nlay,pphi,pecin,pbern)
2      USE parallel_lmdz
3      USE mod_filtreg_p
4      IMPLICIT NONE
5
6c=======================================================================
7c
8c   Auteur:   P. Le Van
9c   -------
10c
11c   Objet:
12c   ------
13c     calcul de la fonction de Bernouilli aux niveaux s  .....
14c     phi  et  ecin  sont des arguments d'entree pour le s-pg .......
15c          bern       est un  argument de sortie pour le s-pg  ......
16c
17c    fonction de Bernouilli = bern = filtre de( geopotentiel +
18c                              energ.cinet.)
19c
20c=======================================================================
21c
22c-----------------------------------------------------------------------
23c   Decalrations:
24c   -------------
25c
26      include "dimensions.h"
27      include "paramet.h"
28c
29c   Arguments:
30c   ----------
31c
32      INTEGER nlay,ngrid
33      REAL pphi(ijb_u:ije_u,nlay),pecin(ijb_u:ije_u,nlay)
34      REAL pbern(ijb_u:ije_u,nlay)
35c
36c   Local:
37c   ------
38c
39      INTEGER   ij,l,ijb,ije,jjb,jje
40c
41c-----------------------------------------------------------------------
42c   calcul de Bernouilli:
43c   ---------------------
44c
45      ijb=ij_begin
46      ije=ij_end+iip1
47      if (pole_sud) ije=ij_end
48
49      jjb=jj_begin
50      jje=jj_end+1
51      if (pole_sud) jje=jj_end
52
53c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)               
54      DO l=1,llm
55   
56        DO 4 ij = ijb,ije
57          pbern( ij,l ) =  pphi( ij,l ) + pecin( ij,l )
58   4    CONTINUE
59       
60       ENDDO
61c$OMP END DO NOWAIT
62c
63c-----------------------------------------------------------------------
64c   filtre:
65c   -------
66c
67
68       
69        CALL filtreg_p( pbern,jjb_u,jje_u,jjb,jje, jjp1, llm,
70     &                  2,1, .true., 1 )
71c
72c-----------------------------------------------------------------------
73     
74     
75      RETURN
76      END
Note: See TracBrowser for help on using the repository browser.