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 comgeom, ONLY : airesurg
    USE dimens_m, ONLY : llm
    USE disvert_m, ONLY : ap, bp
    use massdair_m, only: massdair
    USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1

    !   Arguments:
    REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)
    REAL, INTENT(IN):: 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)
    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 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 comgeom, ONLY : airesurg
16	USE dimens_m, ONLY : llm
17	USE disvert_m, ONLY : ap, bp
18	use massdair_m, only: massdair
19	USE paramet_m, ONLY : iip1, ip1jm, ip1jmp1, jjp1, llmp1
20
21	! Arguments:
22	REAL, INTENT (IN) :: vcov(ip1jm, llm), ucov(ip1jmp1, llm)
23	REAL, INTENT(IN):: 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)
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, l
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	forall (l = 1: llm + 1) p(:, l) = ap(l) + bp(l) * ps
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