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 conf_gcm_m
  USE filtreg_m, ONLY: filtreg
  IMPLICIT NONE

  ! =======================================================================

  ! Auteur:   P. Le Van
  ! -------

  ! Objet:
  ! ------
  ! calcul de la fonction de Bernouilli aux niveaux s  .....
  ! phi  et  ecin  sont des arguments d'entree pour le s-pg .......
  ! bern       est un  argument de sortie pour le s-pg  ......

  ! fonction de Bernouilli = bern = filtre de( geopotentiel +
  ! energ.cinet.)

  ! =======================================================================

  ! -----------------------------------------------------------------------
  ! Decalrations:
  ! -------------


  ! Arguments:
  ! ----------

  INTEGER nlay, ngrid
  REAL, INTENT (IN) :: pphi(ngrid*nlay), pecin(ngrid*nlay)
  REAL pbern(ngrid*nlay)

  ! Local:
  ! ------

  INTEGER ijl

  ! -----------------------------------------------------------------------
  ! calcul de Bernouilli:
  ! ---------------------

  DO ijl = 1, ngrid*nlay
    pbern(ijl) = pphi(ijl) + pecin(ijl)
  END DO

  ! -----------------------------------------------------------------------
  ! filtre:
  ! -------

  CALL filtreg(pbern, jjp1, llm, 2, 1, .TRUE.)

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