Ignore:
Timestamp:
Jun 20, 2001, 3:53:15 PM (23 years ago)
Author:
lmdzadmin
Message:

Merge par rapport a la branche principale
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/tracvl.F

    r2 r232  
    11      SUBROUTINE tracvl(numvanle,iapp_tracvl,nq,pbaru,pbarv ,
    2      *                            p,  masse , q, iapptrac    )
     2     *                    p, masse , q, iapptrac, iadv1, teta, pk  )
    33c
    44c     Auteur :  F. Hourdin
     
    66c
    77ccc   ..   Modif. P. Le Van  ( 20/12/97 )  ...
     8c                 F. Codron     (10/99)
     9
    810c
    911      IMPLICIT NONE
     
    1517#include "comgeom.h"
    1618
    17       INTEGER nq,iapp_tracvl
     19c     .... Arguments  ....
     20c
     21      INTEGER numvanle, nq, iapp_tracvl, iapptrac, iadv1
    1822
    1923      REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
    2024      REAL q(ip1jmp1,llm,nq),masse(ip1jmp1,llm)
    21       REAL p( ip1jmp1,llmp1 )
     25      REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
     26      REAL pk(ip1jmp1,llm)
    2227
     28c     ....  var. locales  .....
     29c
    2330      REAL pbaruc(ip1jmp1,llm),pbarvc(ip1jm,llm)
    2431      REAL massem(ip1jmp1,llm),zdp(ip1jmp1)
     
    2633      REAL pbarug(ip1jmp1,llm),pbarvg(ip1jm,llm),wg(ip1jmp1,llm)
    2734
    28       INTEGER iapptrac
    2935
    30       INTEGER iadvtr, numvanle
     36      INTEGER iadvtr, numvan
    3137      INTEGER ij,l,iq
    3238      REAL zdpmin, zdpmax
    3339      EXTERNAL  minmax
    34  
    35       SAVE iadvtr, massem,pbaruc,pbarvc
     40      SAVE iadvtr, massem, pbaruc, pbarvc, numvan
    3641      DATA iadvtr/0/
     42
     43      numvan = numvanle
    3744
    3845      IF(iadvtr.EQ.0) THEN
     
    99106
    100107c   Advection proprement dite.
    101          DO iq = numvanle, nq
     108c
     109c   test sur iadv1 pour le schema de vapeur d'eau
     110c
     111         IF (numvanle.EQ.1.AND.iadv1.EQ.4) THEN
     112           CALL vlspltqs( q(1,1,1), 2., massem, wg ,
     113     *                 pbarug,pbarvg,dtvr,p,pk,teta )
     114           numvan = 2
     115         ENDIF
     116
     117         DO iq = numvan, nq
    102118          CALL vlsplt( q(1,1,iq), 2. ,massem,wg,pbarug,pbarvg,dtvr )
    103119         ENDDO
Note: See TracChangeset for help on using the changeset viewer.