trunk/dyn3d/interpost.f90


! $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
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/interpost.F,v 1.1.1.1 2004/05/19
3	! 12:53:06 lmdzadmin Exp $
4
5	SUBROUTINE interpost(q, qppm)
6
7	USE dimens_m
8	USE paramet_m
9	USE comconst
10	USE disvert_m
11	USE comgeom
12	IMPLICIT NONE
13
14
15
16	! 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