interpost.F

Go to the documentation of this file.

!
! $Header$
!
        subroutine interpost(q,qppm)

       implicit none


#include "dimensions.h"
#include "paramet.h"
#include "comconst.h"
#include "comvert.h"
#include "comgeom2.h"

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