libf/dyn3d/caldyn0.f90

module caldyn0_m

  IMPLICIT NONE

contains

  SUBROUTINE caldyn0(ucov, vcov, teta, ps, masse, pk, phis, phi, w, &
       pbaru, pbarv)

    ! From dyn3d/caldyn0.F, v 1.1.1.1 2004/05/19 12:53:07
    ! Auteur :  P. Le Van
    ! Objet: calcul des tendances dynamiques
    ! Modif 04/93 F.Forget

    USE dimens_m, ONLY : llm
    USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1
    USE comvert, ONLY : ap, bp
    USE comgeom, ONLY : airesurg
    USE pression_m, ONLY : pression

    !   Arguments:
    REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)
    REAL :: teta(ip1jmp1, llm)
    REAL, INTENT (IN) :: ps(ip1jmp1)
    REAL, INTENT (IN) :: phis(ip1jmp1)
    REAL, INTENT (IN) :: pk(iip1, jjp1, llm)
    REAL :: vcont(ip1jm, llm), ucont(ip1jmp1, llm)
    REAL :: phi(ip1jmp1, llm), masse(ip1jmp1, llm)
    REAL :: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)

    !   Local:

    REAL :: p(ip1jmp1, llmp1)
    REAL :: massebx(ip1jmp1, llm), masseby(ip1jm, llm), psexbarxy(ip1jm)
    REAL :: vorpot(ip1jm, llm)
    REAL :: w(ip1jmp1, llm), ecin(ip1jmp1, llm), convm(ip1jmp1, llm)
    REAL :: bern(ip1jmp1, llm)
    REAL :: massebxy(ip1jm, llm), dp(ip1jmp1)

    INTEGER :: ij

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

    PRINT *, 'Call sequence information: caldyn0'

    !   Calcul des tendances dynamiques:

    CALL covcont(llm, ucov, vcov, ucont, vcont)
    CALL pression(ip1jmp1, ap, bp, ps, p)
    CALL psextbar(ps, psexbarxy)
    CALL massdair(p, masse)
    CALL massbar(masse, massebx, masseby)
    CALL massbarxy(masse, massebxy)
    CALL flumass(massebx, masseby, vcont, ucont, pbaru, pbarv)
    CALL convmas(pbaru, pbarv, convm)

    DO ij = 1, ip1jmp1
       dp(ij) = convm(ij, 1)/airesurg(ij)
    END DO

    CALL vitvert(convm, w)
    CALL tourpot(vcov, ucov, massebxy, vorpot)
    CALL enercin(vcov, ucov, vcont, ucont, ecin)
    CALL bernoui(ip1jmp1, llm, phi, ecin, bern)
    CALL sortvarc0(ucov, teta, ps, masse, pk, phis, vorpot, phi, bern, dp)

  END SUBROUTINE caldyn0

end module caldyn0_m
1	module caldyn0_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE caldyn0(ucov, vcov, teta, ps, masse, pk, phis, phi, w, &
8	pbaru, pbarv)
9
10	! From dyn3d/caldyn0.F, v 1.1.1.1 2004/05/19 12:53:07
11	! Auteur : P. Le Van
12	! Objet: calcul des tendances dynamiques
13	! Modif 04/93 F.Forget
14
15	USE dimens_m, ONLY : llm
16	USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1
17	USE comvert, ONLY : ap, bp
18	USE comgeom, ONLY : airesurg
19	USE pression_m, ONLY : pression
20
21	! Arguments:
22	REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)
23	REAL :: teta(ip1jmp1, llm)
24	REAL, INTENT (IN) :: ps(ip1jmp1)
25	REAL, INTENT (IN) :: phis(ip1jmp1)
26	REAL, INTENT (IN) :: pk(iip1, jjp1, llm)
27	REAL :: vcont(ip1jm, llm), ucont(ip1jmp1, llm)
28	REAL :: phi(ip1jmp1, llm), masse(ip1jmp1, llm)
29	REAL :: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
30
31	! Local:
32
33	REAL :: p(ip1jmp1, llmp1)
34	REAL :: massebx(ip1jmp1, llm), masseby(ip1jm, llm), psexbarxy(ip1jm)
35	REAL :: vorpot(ip1jm, llm)
36	REAL :: w(ip1jmp1, llm), ecin(ip1jmp1, llm), convm(ip1jmp1, llm)
37	REAL :: bern(ip1jmp1, llm)
38	REAL :: massebxy(ip1jm, llm), dp(ip1jmp1)
39
40	INTEGER :: ij
41
42	!-----------------------------------------------------------------------
43
44	PRINT *, 'Call sequence information: caldyn0'
45
46	! Calcul des tendances dynamiques:
47
48	CALL covcont(llm, ucov, vcov, ucont, vcont)
49	CALL pression(ip1jmp1, ap, bp, ps, p)
50	CALL psextbar(ps, psexbarxy)
51	CALL massdair(p, masse)
52	CALL massbar(masse, massebx, masseby)
53	CALL massbarxy(masse, massebxy)
54	CALL flumass(massebx, masseby, vcont, ucont, pbaru, pbarv)
55	CALL convmas(pbaru, pbarv, convm)
56
57	DO ij = 1, ip1jmp1
58	dp(ij) = convm(ij, 1)/airesurg(ij)
59	END DO
60
61	CALL vitvert(convm, w)
62	CALL tourpot(vcov, ucov, massebxy, vorpot)
63	CALL enercin(vcov, ucov, vcont, ucont, ecin)
64	CALL bernoui(ip1jmp1, llm, phi, ecin, bern)
65	CALL sortvarc0(ucov, teta, ps, masse, pk, phis, vorpot, phi, bern, dp)
66
67	END SUBROUTINE caldyn0
68
69	end module caldyn0_m