trunk/dyn3d/interpost.f


! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/interpost.F,v 1.1.1.1 2004/05/19
! 12:53:06 lmdzadmin Exp $

SUBROUTINE interpost(q, qppm)

  USE dimens_m
  USE paramet_m
  USE comconst
  USE disvert_m
  USE comgeom
  IMPLICIT NONE


  ! Arguments
  REAL q(iip1, jjp1, llm)
  REAL qppm(iim, jjp1, llm)
  ! Local
  INTEGER l, i, j

  ! RE-INVERSION DES NIVEAUX
  ! le programme ppm3d travaille avec une 3ème coordonnée inversée par
  ! rapport
  ! de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
  ! On passe donc des niveaux de Lin à ceux du LMDZ

  DO l = 1, llm
    DO j = 1, jjp1
      DO i = 1, iim
        q(i, j, l) = qppm(i, j, llm-l+1)
      END DO
    END DO
  END DO

  ! BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D

  DO l = 1, llm
    DO j = 1, jjp1
      q(iip1, j, l) = q(1, j, l)
    END DO
  END DO


  RETURN

END SUBROUTINE interpost
1	guez	3
2	guez	81	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/interpost.F,v 1.1.1.1 2004/05/19
3			! 12:53:06 lmdzadmin Exp $
4	guez	3
5	guez	81	SUBROUTINE interpost(q, qppm)
6	guez	3
7	guez	81	USE dimens_m
8			USE paramet_m
9			USE comconst
10			USE disvert_m
11			USE comgeom
12			IMPLICIT NONE
13	guez	3
14
15
16	guez	81	! Arguments
17			REAL q(iip1, jjp1, llm)
18			REAL qppm(iim, jjp1, llm)
19			! Local
20			INTEGER l, i, j
21
22			! RE-INVERSION DES NIVEAUX
23			! le programme ppm3d travaille avec une 3ème coordonnée inversée par
24			! rapport
25			! de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
26			! On passe donc des niveaux de Lin à ceux du LMDZ
27
28			DO l = 1, llm
29			DO j = 1, jjp1
30			DO i = 1, iim
31			q(i, j, l) = qppm(i, j, llm-l+1)
32			END DO
33			END DO
34			END DO
35
36			! BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
37
38			DO l = 1, llm
39			DO j = 1, jjp1
40			q(iip1, j, l) = q(1, j, l)
41			END DO
42			END DO
43
44
45			RETURN
46
47			END SUBROUTINE interpost