trunk/dyn3d/bernoui.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/bernoui.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
!
      SUBROUTINE bernoui (ngrid,nlay,pphi,pecin,pbern)
      use dimens_m
      use paramet_m
      use logic
      use filtreg_m, only: filtreg
      IMPLICIT NONE

c=======================================================================
c
c   Auteur:   P. Le Van
c   -------
c
c   Objet:
c   ------
c     calcul de la fonction de Bernouilli aux niveaux s  .....
c     phi  et  ecin  sont des arguments d'entree pour le s-pg .......
c          bern       est un  argument de sortie pour le s-pg  ......
c
c    fonction de Bernouilli = bern = filtre de( geopotentiel + 
c                              energ.cinet.)
c
c=======================================================================
c
c-----------------------------------------------------------------------
c   Decalrations:
c   -------------
c
c
c   Arguments:
c   ----------
c
      INTEGER nlay,ngrid
      REAL pphi(ngrid*nlay),pecin(ngrid*nlay),pbern(ngrid*nlay)
c
c   Local:
c   ------
c
      INTEGER   ijl
c
c-----------------------------------------------------------------------
c   calcul de Bernouilli:
c   ---------------------
c
      DO 4 ijl = 1,ngrid*nlay
         pbern( ijl ) =  pphi( ijl ) + pecin( ijl )
   4  CONTINUE
c
c-----------------------------------------------------------------------
c   filtre:
c   -------
c
      CALL filtreg( pbern, jjp1, llm, 2,1, .true., 1 )
c
c-----------------------------------------------------------------------
      RETURN
      END
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/bernoui.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3	!
4	SUBROUTINE bernoui (ngrid,nlay,pphi,pecin,pbern)
5	use dimens_m
6	use paramet_m
7	use logic
8	use filtreg_m, only: filtreg
9	IMPLICIT NONE
10
11	c=======================================================================
12	c
13	c Auteur: P. Le Van
14	c -------
15	c
16	c Objet:
17	c ------
18	c calcul de la fonction de Bernouilli aux niveaux s .....
19	c phi et ecin sont des arguments d'entree pour le s-pg .......
20	c bern est un argument de sortie pour le s-pg ......
21	c
22	c fonction de Bernouilli = bern = filtre de( geopotentiel +
23	c energ.cinet.)
24	c
25	c=======================================================================
26	c
27	c-----------------------------------------------------------------------
28	c Decalrations:
29	c -------------
30	c
31	c
32	c Arguments:
33	c ----------
34	c
35	INTEGER nlay,ngrid
36	REAL pphi(ngridnlay),pecin(ngridnlay),pbern(ngrid*nlay)
37	c
38	c Local:
39	c ------
40	c
41	INTEGER ijl
42	c
43	c-----------------------------------------------------------------------
44	c calcul de Bernouilli:
45	c ---------------------
46	c
47	DO 4 ijl = 1,ngrid*nlay
48	pbern( ijl ) = pphi( ijl ) + pecin( ijl )
49	4 CONTINUE
50	c
51	c-----------------------------------------------------------------------
52	c filtre:
53	c -------
54	c
55	CALL filtreg( pbern, jjp1, llm, 2,1, .true., 1 )
56	c
57	c-----------------------------------------------------------------------
58	RETURN
59	END