trunk/dyn3d/geopot.f

module geopot_m

  IMPLICIT NONE

contains

  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'intégration se fait de bas en haut.

    USE dimens_m, ONLY: llm

    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

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

    ! Calcul de phi au niveau 1 près du sol :
    phi(:, 1) = phis + teta(:, 1) * (pks - pk(:, 1))

    ! Calcul de phi aux niveaux supérieurs :
    DO l = 2, llm
       phi(:, l) = phi(:, l-1) + 0.5 * (teta(:, l) + teta(:, l-1)) &
            * (pk(:, l-1) - pk(:, l))
    END DO

  END SUBROUTINE geopot

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