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_2d
    USE dimens_m, ONLY: iim, jjm, llm

    REAL, intent(in):: p(iim + 1, jjm + 1, llm + 1)
    ! aux interfaces des llm couches

    real, intent(out):: masse(iim + 1, jjm + 1, llm)

    ! Local:
    INTEGER l

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

    forall (l = 1: llm) masse(:iim, :, l) = airesurg_2d(:iim, :) &
         * (p(:iim, :, l) - p(:iim, :, l + 1))
    masse(iim + 1, :, :) = masse(1, :, :)

  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_2d
15	USE dimens_m, ONLY: iim, jjm, llm
16
17	REAL, intent(in):: p(iim + 1, jjm + 1, llm + 1)
18	! aux interfaces des llm couches
19
20	real, intent(out):: masse(iim + 1, jjm + 1, llm)
21
22	! Local:
23	INTEGER l
24
25	!----------------------------------------------------------
26
27	forall (l = 1: llm) masse(:iim, :, l) = airesurg_2d(:iim, :) &
28	* (p(:iim, :, l) - p(:iim, :, l + 1))
29	masse(iim + 1, :, :) = masse(1, :, :)
30
31	END SUBROUTINE massdair
32
33	end module massdair_m