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 |