source: LMDZ5/branches/LMDZ5V1.0-dev/libf/dyn3dpar/bernoui_p.F @ 1523

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

Merge entre la version V3_conv et le HEAD
YM, JG, LF

  • 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
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"
27#include "logic.h"
28c
29c   Arguments:
30c   ----------
31c
32      INTEGER nlay,ngrid
33      REAL pphi(ngrid,nlay),pecin(ngrid,nlay),pbern(ngrid,nlay)
34c
35c   Local:
36c   ------
37c
38      INTEGER   ij,l,ijb,ije,jjb,jje
39c
40c-----------------------------------------------------------------------
41c   calcul de Bernouilli:
42c   ---------------------
43c
44      ijb=ij_begin
45      ije=ij_end+iip1
46      if (pole_sud) ije=ij_end
47
48      jjb=jj_begin
49      jje=jj_end+1
50      if (pole_sud) jje=jj_end
51
52c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)               
53      DO l=1,llm
54   
55        DO 4 ij = ijb,ije
56          pbern( ij,l ) =  pphi( ij,l ) + pecin( ij,l )
57   4    CONTINUE
58       
59       ENDDO
60c$OMP END DO NOWAIT
61c
62c-----------------------------------------------------------------------
63c   filtre:
64c   -------
65c
66
67       
68        CALL filtreg_p( pbern,jjb,jje, jjp1, llm, 2,1, .true., 1 )
69c
70c-----------------------------------------------------------------------
71     
72     
73      RETURN
74      END
Note: See TracBrowser for help on using the repository browser.