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	module massdair_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE massdair(p, masse)
8
9	! 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	! Authors: P. Le Van, F. Hourdin
13
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	DO l = 1, llm
27	masse(:, l) = airesurg * (p(:, l) - p(:, l + 1))
28
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