| Rev | Line | |
|---|
| [1] | 1 | ! |
|---|
| 2 | ! $Header$ |
|---|
| 3 | ! |
|---|
| 4 | subroutine interpost(q,qppm) |
|---|
| 5 | |
|---|
| 6 | implicit none |
|---|
| 7 | |
|---|
| 8 | |
|---|
| 9 | #include "dimensions.h" |
|---|
| 10 | #include "paramet.h" |
|---|
| 11 | #include "comconst.h" |
|---|
| 12 | #include "comvert.h" |
|---|
| 13 | #include "comgeom2.h" |
|---|
| 14 | |
|---|
| 15 | c Arguments |
|---|
| 16 | real q(iip1,jjp1,llm) |
|---|
| 17 | real qppm(iim,jjp1,llm) |
|---|
| 18 | c Local |
|---|
| 19 | integer l,i,j |
|---|
| 20 | |
|---|
| 21 | c RE-INVERSION DES NIVEAUX |
|---|
| 22 | c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport |
|---|
| 23 | c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface |
|---|
| 24 | c On passe donc des niveaux de Lin à ceux du LMDZ |
|---|
| 25 | |
|---|
| 26 | do l=1,llm |
|---|
| 27 | do j=1,jjp1 |
|---|
| 28 | do i=1,iim |
|---|
| 29 | q(i,j,l)=qppm(i,j,llm-l+1) |
|---|
| 30 | enddo |
|---|
| 31 | enddo |
|---|
| 32 | enddo |
|---|
| 33 | |
|---|
| 34 | c BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D |
|---|
| 35 | |
|---|
| 36 | do l=1,llm |
|---|
| 37 | do j=1,jjp1 |
|---|
| 38 | q(iip1,j,l)=q(1,j,l) |
|---|
| 39 | enddo |
|---|
| 40 | enddo |
|---|
| 41 | |
|---|
| 42 | |
|---|
| 43 | return |
|---|
| 44 | |
|---|
| 45 | end |
|---|
Note: See
TracBrowser
for help on using the repository browser.
