source: LMDZ.3.3/branches/rel-LF/libf/dyn3d/traceur.F @ 279

Last change on this file since 279 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 2.6 KB
Line 
1      SUBROUTINE traceur (iq,iadv,q,teta,pk,w, pbaru, pbarv,dq )
2      IMPLICIT NONE
3
4c=======================================================================
5c
6c   Auteur:  F.Forget / A.Armengaud
7c   -------
8c
9c   ....  Modif  P.Le Van  ,  29/04/97  et 01/06/98    ....
10c
11c   ********************************************************************
12c   Transport des traceurs (dont l'humidite q(1)) par differents shemas
13c   ********************************************************************
14c     iq,q,pbaru et pbarv sont des arguments d'entree  pour le s-pg ....
15c     dq               sont des arguments de sortie pour le s-pg ....
16c
17c     iadv = 1 :  Schema de transport type "humidite specifique  LMD"
18c                 (shema amont (horiz),moy arithmetique (vertic))
19c
20c     iadv = 2  :  Shema amont
21c
22c
23c     iadv = 0  : Indice montrant que q n'est pas vraiment un
24c                traceur, mais une variable pronostique attachee
25c                (par exemple : une pente du shemas de pente)
26c=======================================================================
27
28
29#include "dimensions.h"
30#include "paramet.h"
31#include "comconst.h"
32
33c   Arguments:
34c   ----------
35      INTEGER iq,iadv(iq)
36      REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm)
37      REAL q( ip1jmp1,llm,iq), dq( ip1jmp1,llm,iq)
38      REAL w(ip1jmp1,llm),teta(ip1jmp1,llm),pk(ip1jmp1,llm)
39
40c   Local:
41c   ------
42
43      LOGICAL     qsat
44      PARAMETER ( qsat = .TRUE. )
45
46c
47cc ******     Si  iadv(iq) =  2   et   qsat = .TRUE.  ,   
48cc        ,  on choisit  le shema  amont avec le q sature   ********
49c
50
51      EXTERNAL adv_h2o, amont, amont_qsat
52
53
54       IF( iadv(iq).GT.2 )  THEN
55         PRINT *,' Erreur dans le choix de iadv  pour la routine
56     *  traceur .  Corrigez . '
57         STOP
58       ENDIF
59
60c   ----------------------------------------------------------------
61c   Schema de transport type "humidite specifique  LMD"
62c   (shema amont (horiz),moy arithmetique (vertic))
63c   ----------------------------------------------------------------
64         IF (iadv(iq).EQ.1) THEN
65              CALL adv_h2o (iq,iq,q,w, pbaru, pbarv,dq)       
66c   ----------------------------------------------------------------
67c   Shema amont (dans les 3 dimensions)
68c   ----------------------------------------------------------------
69         ELSE IF (iadv(iq).EQ.2) THEN
70           IF( qsat )   THEN
71              CALL amont_qsat ( iq,iq,q,teta,pk,w, pbaru, pbarv,dq )
72           ELSE
73              CALL amont      ( iq,iq,q,        w, pbaru, pbarv,dq )
74           ENDIF
75c   ----------------------------------------------------------------
76
77         ENDIF
78
79
80      RETURN
81      END
Note: See TracBrowser for help on using the repository browser.