/[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.f90 revision 81 by guez, Wed Mar 5 14:38:41 2014 UTC
# Line 1  Line 1 
1  !  
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convmas.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convmas.F,v 1.1.1.1 2004/05/19
3  !  ! 12:53:07 lmdzadmin Exp $
4        SUBROUTINE convmas (pbaru, pbarv, convm )  
5  c  SUBROUTINE convmas(pbaru, pbarv, convm)
6        use dimens_m  
7        use paramet_m    USE dimens_m
8        use comvert    USE paramet_m
9        use logic    USE disvert_m
10        IMPLICIT NONE    USE conf_gcm_m
11      USE filtreg_m, ONLY: filtreg
12  c=======================================================================    IMPLICIT NONE
13  c  
14  c   Auteurs:  P. Le Van , F. Hourdin  .    ! =======================================================================
15  c   -------  
16  c    ! Auteurs:  P. Le Van , F. Hourdin  .
17  c   Objet:    ! -------
18  c   ------  
19  c    ! Objet:
20  c   ********************************************************************    ! ------
21  c   .... calcul de la convergence du flux de masse aux niveaux p ...  
22  c   ********************************************************************    ! ********************************************************************
23  c    ! .... calcul de la convergence du flux de masse aux niveaux p ...
24  c    ! ********************************************************************
25  c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....  
26  c      .....  convm      est  un argument de sortie pour le s-pg  ....  
27  c    ! pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....
28  c    le calcul se fait de haut en bas,    ! .....  convm      est  un argument de sortie pour le s-pg  ....
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 .    ! le calcul se fait de haut en bas,
31  c    ! la convergence de masse au niveau p(llm+1) est egale a 0. et
32  c    ! n'est pas stockee dans le tableau convm .
33  c=======================================================================  
34  c  
35  c   Declarations:    ! =======================================================================
36  c   -------------  
37      ! Declarations:
38      ! -------------
39        REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm(  ip1jmp1,llm )  
40        INTEGER   l,ij  
41      REAL, INTENT (IN) :: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
42      REAL, INTENT (OUT) :: convm(ip1jmp1, llm)
43  c-----------------------------------------------------------------------    INTEGER l, ij
44  c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......  
45    
46        CALL  convflu( pbaru, pbarv, llm, convm )    ! -----------------------------------------------------------------------
47      ! ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
48  c-----------------------------------------------------------------------  
49  c   filtrage:    CALL convflu(pbaru, pbarv, llm, convm)
50  c   ---------  
51      ! -----------------------------------------------------------------------
52         CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )    ! filtrage:
53      ! ---------
54  c    integration de la convergence de masse de haut  en bas ......  
55      CALL filtreg(convm, jjp1, llm, 2, 2, .TRUE.)
56        DO      l      = llmm1, 1, -1  
57          DO    ij     = 1, ip1jmp1    ! integration de la convergence de masse de haut  en bas ......
58           convm(ij,l) = convm(ij,l) + convm(ij,l+1)  
59          ENDDO    DO l = llmm1, 1, -1
60        ENDDO      DO ij = 1, ip1jmp1
61  c        convm(ij, l) = convm(ij, l) + convm(ij, l+1)
62        RETURN      END DO
63        END    END DO
64    
65      RETURN
66    END SUBROUTINE convmas
Legend:
Removed from v.3  
changed lines
  Added in v.81
  ViewVC Help
Powered by ViewVC 1.1.21