/[lmdze]/trunk/Sources/dyn3d/convmas.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/convmas.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/dyn3d/convmas.f revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC trunk/dyn3d/convmas.f revision 104 by guez, Thu Sep 4 10:05:52 2014 UTC
# Line 1  Line 1 
1  !  module convmas_m
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convmas.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $  
3  !    IMPLICIT NONE
4        SUBROUTINE convmas (pbaru, pbarv, convm )  
5  c  contains
6        use dimens_m  
7        use paramet_m    SUBROUTINE convmas(pbaru, pbarv, convm)
8        use comvert  
9        use logic      ! From LMDZ4/libf/dyn3d/convmas.F, version 1.1.1.1, 2004/05/19 12:53:07
10        IMPLICIT NONE  
11        USE dimens_m, ONLY: iim, jjm, llm
12  c=======================================================================      USE paramet_m, ONLY: ip1jm, ip1jmp1, jjp1, llmm1
13  c      USE filtreg_m, ONLY: filtreg
14  c   Auteurs:  P. Le Van , F. Hourdin  .  
15  c   -------      ! Authors: P. Le Van, F. Hourdin
16  c      ! Objet: calcul de la convergence du flux de masse aux niveaux p
17  c   Objet:  
18  c   ------      ! Le calcul se fait de haut en bas, la convergence de masse au
19  c      ! niveau p(llm+1) est égale à 0 et n'est pas stockée dans le
20  c   ********************************************************************      ! tableau convm.
21  c   .... calcul de la convergence du flux de masse aux niveaux p ...  
22  c   ********************************************************************      REAL, INTENT(IN):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
23  c      REAL, INTENT(OUT):: convm(iim + 1, jjm + 1, llm)
24  c  
25  c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....      ! Local:
26  c      .....  convm      est  un argument de sortie pour le s-pg  ....      INTEGER l
27  c  
28  c    le calcul se fait de haut en bas,      !-----------------------------------------------------------------------
29  c    la convergence de masse au niveau p(llm+1) est egale a 0. et  
30  c    n'est pas stockee dans le tableau convm .      ! Calcul de - (d(pbaru)/dx + d(pbarv)/dy) :
31  c      CALL convflu(pbaru, pbarv, llm, convm)
32  c  
33  c=======================================================================      ! Filtrage :
34  c      CALL filtreg(convm, jjp1, llm, 2, 2, .TRUE.)
35  c   Declarations:  
36  c   -------------      ! Intégration de la convergence de masse de haut en bas :
37        DO l = llmm1, 1, -1
38           convm(:, :, l) = convm(:, :, l) + convm(:, :, l+1)
39        REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm(  ip1jmp1,llm )      END DO
40        INTEGER   l,ij  
41      END SUBROUTINE convmas
42    
43  c-----------------------------------------------------------------------  end module convmas_m
 c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......  
   
       CALL  convflu( pbaru, pbarv, llm, convm )  
   
 c-----------------------------------------------------------------------  
 c   filtrage:  
 c   ---------  
   
        CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )  
   
 c    integration de la convergence de masse de haut  en bas ......  
   
       DO      l      = llmm1, 1, -1  
         DO    ij     = 1, ip1jmp1  
          convm(ij,l) = convm(ij,l) + convm(ij,l+1)  
         ENDDO  
       ENDDO  
 c  
       RETURN  
       END  
Legend:
Removed from v.3  
changed lines
  Added in v.104
  ViewVC Help
Powered by ViewVC 1.1.21