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
    USE paramet_m

    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	43	USE dimens_m
15			USE paramet_m
16	guez	3
17	guez	43	INTEGER, INTENT (IN):: ngrid
18			REAL, INTENT (IN):: teta(ngrid, llm), pks(ngrid)
19			REAL, INTENT (IN) :: phis(ngrid)
20			REAL, INTENT (IN) :: pk(ngrid, llm)
21			REAL, INTENT (out):: phi(ngrid, llm)
22	guez	3
23	guez	43	! Local:
24			INTEGER l
25	guez	3
26	guez	43	! -----------------------------------------------------------------------
27	guez	3
28	guez	43	! Calcul de phi au niveau 1 près du sol :
29			phi(:, 1) = phis + teta(:, 1) * (pks - pk(:, 1))
30
31			! Calcul de phi aux niveaux supérieurs :
32			DO l = 2, llm
33			phi(:, l) = phi(:, l-1) + 0.5 * (teta(:, l) + teta(:, l-1)) &
34			* (pk(:, l-1) - pk(:, l))
35			END DO
36
37			END SUBROUTINE geopot
38
39			end module geopot_m