trunk/dyn3d/dudv1.f

module dudv1_m

  IMPLICIT NONE

contains

  SUBROUTINE dudv1(vorpot, pbaru, pbarv, du, dv)

    ! From LMDZ4/libf/dyn3d/dudv1.F, version 1.1.1.1, 2004/05/19 12:53:06 

    ! Author: P. Le Van

    ! Objet: calcul du terme de rotation. Ce terme est ajouté à
    ! d(ucov)/dt et à d(vcov)/dt.

    USE dimens_m, ONLY: iim, llm
    USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1

    REAL, intent(in):: vorpot(ip1jm, llm)
    REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)

    real du(ip1jmp1, llm), dv(ip1jm, llm)
    ! du et dv sont des arguments de sortie pour le s-pg 

    ! Local:
    INTEGER l, ij

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

    DO l = 1, llm
       DO ij = iip2, ip1jm - 1
          du(ij, l) = 0.125 * (vorpot(ij - iip1, l) + vorpot(ij, l)) &
               * (pbarv(ij - iip1, l) + pbarv(ij - iim, l) + pbarv(ij, l) &
               + pbarv(ij + 1, l))
       END DO

       DO ij = 1, ip1jm - 1
          dv(ij + 1, l) = - 0.125 * (vorpot(ij, l) + vorpot(ij + 1, l)) &
               * (pbaru(ij, l) + pbaru(ij + 1, l) + pbaru(ij + iip1, l) &
               + pbaru(ij + iip2, l))
       END DO

       ! correction pour dv(1, j, l) 
       ! dv(1, j, l) = dv(iip1, j, l) 
       DO ij = 1, ip1jm, iip1
          dv(ij, l) = dv(ij + iim, l)
       END DO
    END DO

  END SUBROUTINE dudv1

end module dudv1_m
1	guez	88	module dudv1_m
2	guez	3
3	guez	81	IMPLICIT NONE
4
5	guez	88	contains
6	guez	81
7	guez	88	SUBROUTINE dudv1(vorpot, pbaru, pbarv, du, dv)
8	guez	81
9	guez	88	! From LMDZ4/libf/dyn3d/dudv1.F, version 1.1.1.1, 2004/05/19 12:53:06
10	guez	81
11	guez	88	! Author: P. Le Van
12	guez	81
13	guez	88	! Objet: calcul du terme de rotation. Ce terme est ajouté à
14			! d(ucov)/dt et à d(vcov)/dt.
15	guez	81
16	guez	88	USE dimens_m, ONLY: iim, llm
17			USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1
18	guez	81
19	guez	88	REAL, intent(in):: vorpot(ip1jm, llm)
20			REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
21	guez	81
22	guez	88	real du(ip1jmp1, llm), dv(ip1jm, llm)
23			! du et dv sont des arguments de sortie pour le s-pg
24	guez	81
25	guez	88	! Local:
26			INTEGER l, ij
27	guez	81
28	guez	88	!----------------------------------------------------------------------
29	guez	81
30	guez	88	DO l = 1, llm
31			DO ij = iip2, ip1jm - 1
32			du(ij, l) = 0.125 * (vorpot(ij - iip1, l) + vorpot(ij, l)) &
33			* (pbarv(ij - iip1, l) + pbarv(ij - iim, l) + pbarv(ij, l) &
34			+ pbarv(ij + 1, l))
35			END DO
36	guez	81
37	guez	88	DO ij = 1, ip1jm - 1
38			dv(ij + 1, l) = - 0.125 * (vorpot(ij, l) + vorpot(ij + 1, l)) &
39			* (pbaru(ij, l) + pbaru(ij + 1, l) + pbaru(ij + iip1, l) &
40			+ pbaru(ij + iip2, l))
41			END DO
42
43			! correction pour dv(1, j, l)
44			! dv(1, j, l) = dv(iip1, j, l)
45			DO ij = 1, ip1jm, iip1
46			dv(ij, l) = dv(ij + iim, l)
47			END DO
48	guez	81	END DO
49
50	guez	88	END SUBROUTINE dudv1
51
52			end module dudv1_m