libf/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 comvert
      use comgeom
       implicit none


c Arguments   
      real   q(iip1,jjp1,llm)
      real   qppm(iim,jjp1,llm)
c Local
      integer l,i,j
  
c RE-INVERSION DES NIVEAUX
c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport
c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
c 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)
             enddo
          enddo
         enddo
            
c BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D

         do l=1,llm
           do j=1,jjp1
            q(iip1,j,l)=q(1,j,l)
           enddo
         enddo
  
      
       return

       end
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/interpost.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3	!
4	subroutine interpost(q,qppm)
5
6	use dimens_m
7	use paramet_m
8	use comconst
9	use comvert
10	use comgeom
11	implicit none
12
13
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