trunk/dyn3d/massbar.f

module massbar_m

  IMPLICIT NONE

contains

  SUBROUTINE massbar(masse, massebx, masseby)

    ! From LMDZ4/libf/dyn3d/massbar.F, version 1.1.1.1, 2004/05/19 12:53:05
    ! Authors: P. Le Van, F. Hourdin.

    ! Calcule les moyennes en x et y de la masse d'air dans chaque
    ! maille. Cf. "inigeom.txt" et "massbar.txt".

    USE comgeom, ONLY: alpha1p2, alpha1p4, alpha2p3, alpha3p4
    USE dimens_m, ONLY: iim, llm
    USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1

    REAL, intent(in):: masse(ip1jmp1, llm)
    real, intent(out):: massebx(ip1jmp1, llm), masseby(ip1jm, llm)

    ! Local:
    integer l, ij

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

    DO l = 1, llm
       DO ij = 1, ip1jmp1 - 1
          massebx(ij, l) = masse(ij, l) * alpha1p2(ij) + &
               masse(ij+1, l) * alpha3p4(ij+1)
       ENDDO

       ! correction pour massebx(iip1, j) 
       ! massebx(iip1, j)= massebx(1, j) 
       DO ij = iip1, ip1jmp1, iip1
          massebx(ij, l) = massebx(ij - iim, l)
       ENDDO

       DO ij = 1, ip1jm
          masseby(ij, l) = masse(ij, l) * alpha2p3(ij) + &
               masse(ij+iip1, l) * alpha1p4(ij+iip1)
       ENDDO
    end DO

  END SUBROUTINE massbar

end module massbar_m
1	guez	91	module massbar_m
2	guez	3
3	guez	91	IMPLICIT NONE
4	guez	3
5	guez	91	contains
6	guez	3
7	guez	91	SUBROUTINE massbar(masse, massebx, masseby)
8	guez	3
9	guez	91	! From LMDZ4/libf/dyn3d/massbar.F, version 1.1.1.1, 2004/05/19 12:53:05
10			! Authors: P. Le Van, F. Hourdin.
11	guez	3
12	guez	91	! Calcule les moyennes en x et y de la masse d'air dans chaque
13			! maille. Cf. "inigeom.txt" et "massbar.txt".
14	guez	78
15	guez	164	USE comgeom, ONLY: alpha1p2, alpha1p4, alpha2p3, alpha3p4
16	guez	91	USE dimens_m, ONLY: iim, llm
17			USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1
18	guez	78
19	guez	164	REAL, intent(in):: masse(ip1jmp1, llm)
20			real, intent(out):: massebx(ip1jmp1, llm), masseby(ip1jm, llm)
21	guez	78
22	guez	91	! Local:
23			integer l, ij
24	guez	78
25	guez	91	!--------------------------------------------------------------
26	guez	78
27	guez	91	DO l = 1, llm
28			DO ij = 1, ip1jmp1 - 1
29	guez	164	massebx(ij, l) = masse(ij, l) * alpha1p2(ij) + &
30	guez	91	masse(ij+1, l) * alpha3p4(ij+1)
31			ENDDO
32	guez	78
33	guez	164	! correction pour massebx(iip1, j)
34			! massebx(iip1, j)= massebx(1, j)
35	guez	91	DO ij = iip1, ip1jmp1, iip1
36	guez	164	massebx(ij, l) = massebx(ij - iim, l)
37	guez	91	ENDDO
38	guez	78
39	guez	164	DO ij = 1, ip1jm
40			masseby(ij, l) = masse(ij, l) * alpha2p3(ij) + &
41	guez	91	masse(ij+iip1, l) * alpha1p4(ij+iip1)
42			ENDDO
43			end DO
44
45			END SUBROUTINE massbar
46
47			end module massbar_m