interpre.F

Go to the documentation of this file.

!
! $Id: interpre.F 1403 2010-07-01 09:02:53Z fairhead $
!
       subroutine interpre(q,qppm,w,fluxwppm,masse,
     s            apppm,bpppm,massebx,masseby,pbaru,pbarv,
     s            unatppm,vnatppm,psppm)

      USE control_mod

       implicit none

#include "dimensions.h"
c#include "paramr2.h"
#include "paramet.h"
#include "comconst.h"
#include "comdissip.h"
#include "comvert.h"
#include "comgeom2.h"
#include "logic.h"
#include "temps.h"
#include "ener.h"
#include "description.h"

c---------------------------------------------------
c arguments     
      real   apppm(llm+1),bpppm(llm+1)
      real   q(iip1,jjp1,llm),qppm(iim,jjp1,llm)
c---------------------------------------------------
      real   masse(iip1,jjp1,llm) 
      real   massebx(iip1,jjp1,llm),masseby(iip1,jjm,llm)      
      real   w(iip1,jjp1,llm+1)
      real   fluxwppm(iim,jjp1,llm)
      real   pbaru(iip1,jjp1,llm )
      real   pbarv(iip1,jjm,llm)
      real   unatppm(iim,jjp1,llm)
      real   vnatppm(iim,jjp1,llm)
      real   psppm(iim,jjp1)
c---------------------------------------------------
c local
      real   vnat(iip1,jjp1,llm)
      real   unat(iip1,jjp1,llm)
      real   fluxw(iip1,jjp1,llm)
      real   smass(iip1,jjp1)
c----------------------------------------------------
      integer l,ij,i,j

c       calcul de la pression de surface
c       les coefficients ap et bp sont passés en common 
c       calcul de la pression au sol en mb optimisée pour 
c       la vectorialisation
                   
         do j=1,jjp1
             do i=1,iip1
                smass(i,j)=0.
             enddo
         enddo

         do l=1,llm
             do j=1,jjp1
                 do i=1,iip1
                    smass(i,j)=smass(i,j)+masse(i,j,l)
                 enddo
             enddo
         enddo
      
         do j=1,jjp1
             do i=1,iim
                 psppm(i,j)=smass(i,j)/aire(i,j)*g*0.01
             end do
         end do                        
       
c reconstruction des champs contravariants
c le programme ppm3d travaille avec les composantes
c de vitesse et pas les flux, on doit donc passer de l'un à l'autre
c dans le même temps, on fait le changement d























'orientation du vent en v      do l=1,llm          do j=1,jjm              do i=1,iip1                  vnat(i,j,l)=-pbarv(i,j,l)/masseby(i,j,l)*cv(i,j)                           enddo          enddo          do  i=1,iim          vnat(i,jjp1,l)=0.          enddo          do j=1,jjp1              do i=1,iip1                  unat(i,j,l)=pbaru(i,j,l)/massebx(i,j,l)*cu(i,j)              enddo          enddo      enddo              c CALCUL DU FLUX MASSIQUE VERTICALc Flux en l=1 (sol) nul      fluxw=0.              do l=1,llm           do j=1,jjp1              do i=1,iip1                             fluxw(i,j,l)=w(i,j,l)*g*0.01/aire(i,j)C               print*,i,j,l,'fluxw(i,j,l)=
',fluxw(i,j,l),C     c                      'w(i,j,l)=