source: LMDZ5/trunk/libf/dyn3dpar/bernoui_p.F @ 4735

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