trunk/dyn3d/massdair.f

module massdair_m

  IMPLICIT NONE

contains

  SUBROUTINE massdair(p, masse)

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

    ! Calcule la masse d'air dans chaque maille.
    ! Authors: P. Le Van, F. Hourdin

    USE comgeom, ONLY: airesurg
    USE dimens_m, ONLY: iim, llm
    USE paramet_m, ONLY: iip1, ip1jmp1, llmp1

    REAL, intent(in):: p(ip1jmp1, llmp1) ! aux interfaces des llm couches
    real, intent(out):: masse(ip1jmp1, llm)

    ! Variables locales
    INTEGER l, ij

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

    DO l = 1, llm
       masse(:, l) = airesurg * (p(:, l) - p(:, l + 1))

       DO ij = 1, ip1jmp1, iip1
          masse(ij + iim, l) = masse(ij, l)
       ENDDO
    end DO

  END SUBROUTINE massdair

end module massdair_m
1	guez	67	module massdair_m
2	guez	3
3	guez	67	IMPLICIT NONE
4	guez	3
5	guez	67	contains
6	guez	3
7	guez	67	SUBROUTINE massdair(p, masse)
8	guez	3
9	guez	67	! From LMDZ4/libf/dyn3d/massdair.F, version 1.1.1.1 2004/05/19 12:53:07
10
11			! Calcule la masse d'air dans chaque maille.
12	guez	79	! Authors: P. Le Van, F. Hourdin
13	guez	67
14			USE comgeom, ONLY: airesurg
15			USE dimens_m, ONLY: iim, llm
16			USE paramet_m, ONLY: iip1, ip1jmp1, llmp1
17
18			REAL, intent(in):: p(ip1jmp1, llmp1) ! aux interfaces des llm couches
19			real, intent(out):: masse(ip1jmp1, llm)
20
21			! Variables locales
22			INTEGER l, ij
23
24			!----------------------------------------------------------
25
26	guez	78	DO l = 1, llm
27			masse(:, l) = airesurg * (p(:, l) - p(:, l + 1))
28	guez	67
29			DO ij = 1, ip1jmp1, iip1
30			masse(ij + iim, l) = masse(ij, l)
31			ENDDO
32			end DO
33
34			END SUBROUTINE massdair
35
36			end module massdair_m