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	guez	23	module caldyn0_m
2	guez	3
3	guez	23	IMPLICIT NONE
4	guez	4
5	guez	23	contains
6	guez	3
7	guez	23	SUBROUTINE caldyn0(ucov, vcov, teta, ps, masse, pk, phis, phi, w, &
8			pbaru, pbarv)
9	guez	3
10	guez	23	! 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	guez	3
15	guez	23	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	guez	3
20	guez	23	! 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	guez	3
30	guez	23	! Local:
31	guez	3
32	guez	23	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	guez	3
39	guez	37	INTEGER :: ij, l
40	guez	3
41	guez	23	!-----------------------------------------------------------------------
42	guez	3
43	guez	23	PRINT *, 'Call sequence information: caldyn0'
44	guez	3
45	guez	23	! Calcul des tendances dynamiques:
46	guez	3
47	guez	23	CALL covcont(llm, ucov, vcov, ucont, vcont)
48	guez	37	forall (l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps
49	guez	23	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	guez	3
56	guez	23	DO ij = 1, ip1jmp1
57			dp(ij) = convm(ij, 1)/airesurg(ij)
58			END DO
59	guez	3
60	guez	23	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