/[lmdze]/trunk/phylmd/nflxtr.f
ViewVC logotype

Annotation of /trunk/phylmd/nflxtr.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 10 - (hide annotations)
Fri Apr 18 14:45:53 2008 UTC (16 years, 1 month ago) by guez
Original Path: trunk/libf/phylmd/nflxtr.f
File size: 5362 byte(s) 
Added NetCDF directory "/home/guez/include" in "g95.mk" and
"nag_tools.mk".

Added some "intent" attributes in "PVtheta", "advtrac", "caladvtrac",
"calfis", "diagedyn", "dissip", "vlspltqs", "aeropt", "ajsec",
"calltherm", "clmain", "cltrac", "cltracrn", "concvl", "conema3",
"conflx", "fisrtilp", "newmicro", "nuage", "diagcld1", "diagcld2",
"drag_noro", "lift_noro", "SUGWD", "physiq", "phytrac", "radlwsw", "thermcell".

Removed the case "ierr == 0" in "abort_gcm"; moved call to "histclo"
and messages for end of run from "abort_gcm" to "gcm"; replaced call
to "abort_gcm" in "leapfrog" by exit from outer loop.

In "calfis": removed argument "pp" and variable "unskap"; changed
"pksurcp" from scalar to rank 2; use "pressure_var"; rewrote
computation of "zplev", "zplay", "ztfi", "pcvgt" using "dyn_phy";
added computation of "pls".

Removed unused variable in "dynredem0".

In "exner_hyb": changed "dellta" from scalar to rank 1; replaced call
to "ssum" by call to "sum"; removed variables "xpn" and "xps";
replaced some loops by array expressions.

In "leapfrog": use "pressure_var"; deleted variables "p", "longcles".

Converted common blocks "YOECUMF", "YOEGWD" to modules.

Removed argument "pplay" in "cvltr", "diagetpq", "nflxtr".

Created module "raddimlw" from include file "raddimlw.h".

Corrected call to "new_unit" in "test_disvert".
1 guez 3 !
2     ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/nflxtr.F,v 1.1.1.1 2004/05/19 12:53:08 lmdzadmin Exp $
3     !
4     SUBROUTINE nflxtr(pdtime,pmfu,pmfd,pen_u,pde_u,pen_d,pde_d,
5 guez 10 . paprs,x,dx)
6 guez 3 use dimens_m
7     use dimphy
8     use YOMCST
9 guez 10 use yoecumf
10 guez 3 IMPLICIT NONE
11     c=====================================================================
12     c Objet : Melange convectif de traceurs a partir des flux de masse
13     c Date : 13/12/1996 -- 13/01/97
14     c Auteur: O. Boucher (LOA) sur inspiration de Z. X. Li (LMD),
15     c Brinkop et Sausen (1996) et Boucher et al. (1996).
16     c ATTENTION : meme si cette routine se veut la plus generale possible,
17     c elle a herite de certaines notations et conventions du
18     c schema de Tiedtke (1993).
19     c --En particulier, les couches sont numerotees de haut en bas !!!
20     c Ceci est valable pour les flux
21 guez 10 c mais pas pour les entrees x, paprs !!!!
22 guez 3 c --pmfu est positif, pmfd est negatif
23     c --Tous les flux d'entrainements et de detrainements sont positifs
24     c contrairement au schema de Tiedtke d'ou les changements de signe!!!!
25     c=====================================================================
26     c
27     c
28 guez 7 REAL, intent(in):: pdtime
29 guez 3 c--les flux sont definis au 1/2 niveaux
30     c--pmfu(klev+1) et pmfd(klev+1) sont implicitement nuls
31     REAL pmfu(klon,klev) ! flux de masse dans le panache montant
32     REAL pmfd(klon,klev) ! flux de masse dans le panache descendant
33     REAL pen_u(klon,klev) ! flux entraine dans le panache montant
34     REAL pde_u(klon,klev) ! flux detraine dans le panache montant
35     REAL pen_d(klon,klev) ! flux entraine dans le panache descendant
36     REAL pde_d(klon,klev) ! flux detraine dans le panache descendant
37    
38     REAL, intent(in):: paprs(klon,klev+1) ! pression aux 1/2 couches (bas en haut)
39     REAL, intent(in):: x(klon,klev) ! q de traceur (bas en haut)
40     REAL dx(klon,klev) ! tendance de traceur (bas en haut)
41     c
42     c--flux convectifs mais en variables locales
43     REAL zmfu(klon,klev+1)
44     REAL zmfd(klon,klev+1)
45     REAL zen_u(klon,klev)
46     REAL zde_u(klon,klev)
47     REAL zen_d(klon,klev)
48     REAL zde_d(klon,klev)
49     real zmfe
50     c
51     c--variables locales
52     c--les flux de x sont definis aux 1/2 niveaux
53     c--xu et xd sont definis aux niveaux complets
54     REAL xu(klon,klev) ! q de traceurs dans le panache montant
55     REAL xd(klon,klev) ! q de traceurs dans le panache descendant
56     REAL zmfux(klon,klev+1) ! flux de x dans le panache montant
57     REAL zmfdx(klon,klev+1) ! flux de x dans le panache descendant
58     REAL zmfex(klon,klev+1) ! flux de x dans l'environnement
59     INTEGER i, k
60     REAL zmfmin
61     PARAMETER (zmfmin=1.E-10)
62    
63     c =========================================
64     c
65     c
66     c Extension des flux UP et DN sur klev+1 niveaux
67     c =========================================
68     do k=1,klev
69     do i=1,klon
70     zmfu(i,k)=pmfu(i,k)
71     zmfd(i,k)=pmfd(i,k)
72     enddo
73     enddo
74     do i=1,klon
75     zmfu(i,klev+1)=0.
76     zmfd(i,klev+1)=0.
77     enddo
78    
79     c--modif pour diagnostiquer les detrainements
80     c =========================================
81     c on privilegie l'ajustement de l'entrainement dans l'ascendance.
82    
83     do k=1, klev
84     do i=1, klon
85     zen_d(i,k)=pen_d(i,k)
86     zde_u(i,k)=pde_u(i,k)
87     zde_d(i,k) =-zmfd(i,k+1)+zmfd(i,k)+zen_d(i,k)
88     zen_u(i,k) = zmfu(i,k+1)-zmfu(i,k)+zde_u(i,k)
89     enddo
90     enddo
91     c
92     c--calcul des flux dans le panache montant
93     c =========================================
94     c
95     c Dans la premiere couche, on prend q comme valeur de qu
96     c
97     do i=1, klon
98     zmfux(i,1)=0.0
99     enddo
100     c
101     c Autres couches
102     do k=1,klev
103     do i=1, klon
104     if ((zmfu(i,k+1)+zde_u(i,k)).lt.zmfmin) THEN
105     xu(i,k)=x(i,k)
106     else
107     xu(i,k)=(zmfux(i,k)+zen_u(i,k)*x(i,k))
108     s /(zmfu(i,k+1)+zde_u(i,k))
109     endif
110     zmfux(i,k+1)=zmfu(i,k+1)*xu(i,k)
111     enddo
112     enddo
113     c
114     c--calcul des flux dans le panache descendant
115     c =========================================
116     c
117     do i=1, klon
118     zmfdx(i,klev+1)=0.0
119     enddo
120     c
121     do k=klev,1,-1
122     do i=1, klon
123     if ((zde_d(i,k)-zmfd(i,k)).lt.zmfmin) THEN
124     xd(i,k)=x(i,k)
125     else
126     xd(i,k)=(zmfdx(i,k+1)-zen_d(i,k)*x(i,k)) /
127     . (zmfd(i,k)-zde_d(i,k))
128     endif
129     zmfdx(i,k)=zmfd(i,k)*xd(i,k)
130     enddo
131     enddo
132     c
133     c--introduction du flux de retour dans l'environnement
134     c =========================================
135     c
136     do k=2, klev
137     do i=1, klon
138     zmfe=-zmfu(i,k)-zmfd(i,k)
139     if (zmfe.le.0.) then
140     zmfex(i,k)= zmfe*x(i,k)
141     else
142     zmfex(i,k)= zmfe*x(i,k-1)
143     endif
144     enddo
145     enddo
146    
147     do i=1, klon
148     zmfex(i,1)=0.
149     zmfex(i,klev+1)=0.
150     enddo
151     c
152     c--calcul final des tendances
153     c
154     do k=1, klev
155     do i=1, klon
156     dx(i,k)=RG/(paprs(i,k)-paprs(i,k+1))*pdtime*
157     . ( zmfux(i,k) - zmfux(i,k+1) +
158     . zmfdx(i,k) - zmfdx(i,k+1) +
159     . zmfex(i,k) - zmfex(i,k+1) )
160     enddo
161     enddo
162     c
163     return
164     end
  ViewVC Help
Powered by ViewVC 1.1.21