trunk/dyn3d/geopot.f

SUBROUTINE geopot(ngrid,teta,pk,pks,phis,phi)

  ! From libf/dyn3d/geopot.F, version 1.1.1.1 2004/05/19
  ! Author: P. Le Van
  ! Objet : calcul du géopotentiel aux milieux des couches
  ! L'integration se fait de bas en haut.

  USE dimens_m
  USE paramet_m

  IMPLICIT NONE

  ! Arguments:
  INTEGER, INTENT (IN):: ngrid
  REAL, INTENT (IN):: teta(ngrid,llm), pks(ngrid)
  REAL, INTENT (IN) :: phis(ngrid)
  REAL, INTENT (IN) :: pk(ngrid,llm)
  REAL, INTENT (out)::  phi(ngrid,llm)

  ! Local:
  INTEGER l, ij

  ! -----------------------------------------------------------------------

  ! calcul de phi au niveau 1 pres du sol
  DO  ij = 1, ngrid
     phi(ij,1) = phis(ij) + teta(ij,1)*(pks(ij)-pk(ij,1))
  end DO

  ! calcul de phi aux niveaux superieurs
  DO l = 2, llm
     DO ij = 1, ngrid
        phi(ij,l) = phi(ij,l-1) + 0.5 * (teta(ij,l) + teta(ij,l-1)) &
             * (pk(ij,l-1) - pk(ij,l))
     END DO
  END DO

END SUBROUTINE geopot
1	guez	3	SUBROUTINE geopot(ngrid,teta,pk,pks,phis,phi)
2
3	guez	32	! From libf/dyn3d/geopot.F, version 1.1.1.1 2004/05/19
4			! Author: P. Le Van
5			! Objet : calcul du géopotentiel aux milieux des couches
6			! L'integration se fait de bas en haut.
7	guez	3
8			USE dimens_m
9			USE paramet_m
10
11			IMPLICIT NONE
12
13			! Arguments:
14			INTEGER, INTENT (IN):: ngrid
15			REAL, INTENT (IN):: teta(ngrid,llm), pks(ngrid)
16			REAL, INTENT (IN) :: phis(ngrid)
17			REAL, INTENT (IN) :: pk(ngrid,llm)
18			REAL, INTENT (out):: phi(ngrid,llm)
19
20			! Local:
21			INTEGER l, ij
22
23			! -----------------------------------------------------------------------
24
25	guez	32	! calcul de phi au niveau 1 pres du sol
26	guez	3	DO ij = 1, ngrid
27			phi(ij,1) = phis(ij) + teta(ij,1)*(pks(ij)-pk(ij,1))
28			end DO
29
30	guez	32	! calcul de phi aux niveaux superieurs
31	guez	3	DO l = 2, llm
32			DO ij = 1, ngrid
33			phi(ij,l) = phi(ij,l-1) + 0.5 * (teta(ij,l) + teta(ij,l-1)) &
34			* (pk(ij,l-1) - pk(ij,l))
35			END DO
36			END DO
37
38			END SUBROUTINE geopot