trunk/dyn3d/caldyn0.f

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

    !   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, l

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

    PRINT *, 'Call sequence information: caldyn0'

    !   Calcul des tendances dynamiques:

    CALL covcont(llm, ucov, vcov, ucont, vcont)
    forall (l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps
    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
20	! Arguments:
21	REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)
22	REAL :: teta(ip1jmp1, llm)
23	REAL, INTENT (IN) :: ps(ip1jmp1)
24	REAL, INTENT (IN) :: phis(ip1jmp1)
25	REAL, INTENT (IN) :: pk(iip1, jjp1, llm)
26	REAL :: vcont(ip1jm, llm), ucont(ip1jmp1, llm)
27	REAL :: phi(ip1jmp1, llm), masse(ip1jmp1, llm)
28	REAL :: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
29
30	! Local:
31
32	REAL :: p(ip1jmp1, llmp1)
33	REAL :: massebx(ip1jmp1, llm), masseby(ip1jm, llm), psexbarxy(ip1jm)
34	REAL :: vorpot(ip1jm, llm)
35	REAL :: w(ip1jmp1, llm), ecin(ip1jmp1, llm), convm(ip1jmp1, llm)
36	REAL :: bern(ip1jmp1, llm)
37	REAL :: massebxy(ip1jm, llm), dp(ip1jmp1)
38
39	INTEGER :: ij, l
40
41	!-----------------------------------------------------------------------
42
43	PRINT *, 'Call sequence information: caldyn0'
44
45	! Calcul des tendances dynamiques:
46
47	CALL covcont(llm, ucov, vcov, ucont, vcont)
48	forall (l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps
49	CALL psextbar(ps, psexbarxy)
50	CALL massdair(p, masse)
51	CALL massbar(masse, massebx, masseby)
52	CALL massbarxy(masse, massebxy)
53	CALL flumass(massebx, masseby, vcont, ucont, pbaru, pbarv)
54	CALL convmas(pbaru, pbarv, convm)
55
56	DO ij = 1, ip1jmp1
57	dp(ij) = convm(ij, 1)/airesurg(ij)
58	END DO
59
60	CALL vitvert(convm, w)
61	CALL tourpot(vcov, ucov, massebxy, vorpot)
62	CALL enercin(vcov, ucov, vcont, ucont, ecin)
63	CALL bernoui(ip1jmp1, llm, phi, ecin, bern)
64	CALL sortvarc0(ucov, teta, ps, masse, pk, phis, vorpot, phi, bern, dp)
65
66	END SUBROUTINE caldyn0
67
68	end module caldyn0_m