libf/dyn3d/geopot.f90

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	guez	43	module geopot_m
2	guez	3
3	guez	43	IMPLICIT NONE
4	guez	3
5	guez	43	contains
6	guez	3
7	guez	43	SUBROUTINE geopot(ngrid, teta, pk, pks, phis, phi)
8	guez	3
9	guez	43	! 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	guez	3
14	guez	44	USE dimens_m, ONLY: llm
15	guez	3
16	guez	43	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	guez	3
22	guez	43	! Local:
23			INTEGER l
24	guez	3
25	guez	43	! -----------------------------------------------------------------------
26	guez	3
27	guez	43	! 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