Context Navigation

concvl.F90 @ 2850

Last change on this file since 2850 was 2824, checked in by Laurent Fairhead, 8 years ago
Introduction of convoccur variable and modification of the wbeff variable in the convection so that wbeff/convoccur is the mean value of wbeff when convection occurs Jyg
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 18.7 KB

Line
1	SUBROUTINE concvl(iflag_clos, &
2	dtime, paprs, pplay, k_upper_cv, &
3	t, q, t_wake, q_wake, s_wake, u, v, tra, ntra, &
4	Ale, Alp, sig1, w01, &
5	d_t, d_q, d_u, d_v, d_tra, &
6	rain, snow, kbas, ktop, sigd, &
7	cbmf, plcl, plfc, wbeff, convoccur, &
8	upwd, dnwd, dnwdbis, &
9	Ma, mip, Vprecip, &
10	cape, cin, tvp, Tconv, iflag, &
11	pbase, bbase, dtvpdt1, dtvpdq1, dplcldt, dplcldr, &
12	qcondc, wd, pmflxr, pmflxs, &
13	!RomP >>>
14	!! . da,phi,mp,dd_t,dd_q,lalim_conv,wght_th)
15	da, phi, mp, phi2, d1a, dam, sij, clw, elij, & ! RomP
16	dd_t, dd_q, lalim_conv, wght_th, & ! RomP
17	evap, ep, epmlmMm, eplaMm, & ! RomP
18	wdtrainA, wdtrainM, wght, qtc, sigt, &
19	tau_cld_cv, coefw_cld_cv, & ! RomP+RL, AJ
20	!RomP <<<
21	epmax_diag) ! epmax_cape
22	! **************************************************************
23	! *
24	! CONCVL *
25	! *
26	! *
27	! written by : Sandrine Bony-Lena, 17/05/2003, 11.16.04 *
28	! modified by : *
29	! **************************************************************
30
31
32	USE dimphy
33	USE infotrac_phy, ONLY: nbtr
34	USE phys_local_var_mod, ONLY: omega
35	USE print_control_mod, ONLY: prt_level, lunout
36	IMPLICIT NONE
37	! ======================================================================
38	! Auteur(s): S. Bony-Lena (LMD/CNRS) date: ???
39	! Objet: schema de convection de Emanuel (1991) interface
40	! ======================================================================
41	! Arguments:
42	! dtime--input-R-pas d'integration (s)
43	! s-------input-R-la vAleur "s" pour chaque couche
44	! sigs----input-R-la vAleur "sigma" de chaque couche
45	! sig-----input-R-la vAleur de "sigma" pour chaque niveau
46	! psolpa--input-R-la pression au sol (en Pa)
47	! pskapa--input-R-exponentiel kappa de psolpa
48	! h-------input-R-enthAlpie potentielle (CpT/P*kappa)
49	! q-------input-R-vapeur d'eau (en kg/kg)
50
51	! work*: input et output: deux variables de travail,
52	! on peut les mettre a 0 au debut
53	! ALE--------input-R-energie disponible pour soulevement
54	! ALP--------input-R-puissance disponible pour soulevement
55
56	! d_h--------output-R-increment de l'enthAlpie potentielle (h)
57	! d_q--------output-R-increment de la vapeur d'eau
58	! rain-------output-R-la pluie (mm/s)
59	! snow-------output-R-la neige (mm/s)
60	! upwd-------output-R-saturated updraft mass flux (kg/m**2/s)
61	! dnwd-------output-R-saturated downdraft mass flux (kg/m**2/s)
62	! dnwd0------output-R-unsaturated downdraft mass flux (kg/m**2/s)
63	! Ma---------output-R-adiabatic ascent mass flux (kg/m2/s)
64	! mip--------output-R-mass flux shed by adiabatic ascent (kg/m2/s)
65	! Vprecip----output-R-vertical profile of total precipitation (kg/m2/s)
66	! Tconv------output-R-environment temperature seen by convective scheme (K)
67	! Cape-------output-R-CAPE (J/kg)
68	! Cin -------output-R-CIN (J/kg)
69	! Tvp--------output-R-Temperature virtuelle d'une parcelle soulevee
70	! adiabatiquement a partir du niveau 1 (K)
71	! deltapb----output-R-distance entre LCL et base de la colonne (<0 ; Pa)
72	! Ice_flag---input-L-TRUE->prise en compte de la thermodynamique de la glace
73	! dd_t-------output-R-increment de la temperature du aux descentes precipitantes
74	! dd_q-------output-R-increment de la vapeur d'eau du aux desc precip
75	! lalim_conv-
76	! wght_th----
77	! evap-------output-R
78	! ep---------output-R
79	! epmlmMm----output-R
80	! eplaMm-----output-R
81	! wdtrainA---output-R
82	! wdtrainM---output-R
83	! wght-------output-R
84	! ======================================================================
85
86
87	include "clesphys.h"
88
89	INTEGER iflag_clos
90
91	REAL dtime, paprs(klon, klev+1), pplay(klon, klev)
92	INTEGER k_upper_cv
93	REAL t(klon, klev), q(klon, klev), u(klon, klev), v(klon, klev)
94	REAL t_wake(klon, klev), q_wake(klon, klev)
95	REAL s_wake(klon)
96	REAL tra(klon, klev, nbtr)
97	INTEGER ntra
98	REAL sig1(klon, klev), w01(klon, klev), ptop2(klon)
99	REAL pmflxr(klon, klev+1), pmflxs(klon, klev+1)
100	REAL Ale(klon), Alp(klon)
101
102	REAL d_t(klon, klev), d_q(klon, klev), d_u(klon, klev), d_v(klon, klev)
103	REAL dd_t(klon, klev), dd_q(klon, klev)
104	REAL d_tra(klon, klev, nbtr)
105	REAL rain(klon), snow(klon)
106
107	INTEGER kbas(klon), ktop(klon)
108	REAL em_ph(klon, klev+1), em_p(klon, klev)
109	REAL upwd(klon, klev), dnwd(klon, klev), dnwdbis(klon, klev)
110
111	!! REAL Ma(klon,klev), mip(klon,klev),Vprecip(klon,klev) !jyg
112	REAL Ma(klon, klev), mip(klon, klev), Vprecip(klon, klev+1) !jyg
113	REAL Vprecipi(klon, klev+1) !jyg
114	REAL wght(klon, klev) !RL
115
116	REAL da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
117	! RomP >>>
118	REAL phi2(klon, klev, klev)
119	REAL d1a(klon, klev), dam(klon, klev)
120	REAL sij(klon, klev, klev), clw(klon, klev), elij(klon, klev, klev)
121	REAL wdtrainA(klon, klev), wdtrainM(klon, klev)
122	REAL evap(klon, klev), ep(klon, klev)
123	REAL epmlmMm(klon, klev, klev), eplaMm(klon, klev)
124	! RomP <<<
125	REAL cape(klon), cin(klon), tvp(klon, klev)
126	REAL Tconv(klon, klev)
127
128	!CR:test: on passe lentr et alim_star des thermiques
129	INTEGER lalim_conv(klon)
130	REAL wght_th(klon, klev)
131	REAL em_sig1feed ! sigma at lower bound of feeding layer
132	REAL em_sig2feed ! sigma at upper bound of feeding layer
133	REAL em_wght(klev) ! weight density determining the feeding mixture
134	!on enleve le save
135	! SAVE em_sig1feed,em_sig2feed,em_wght
136
137	INTEGER iflag(klon)
138	REAL rflag(klon)
139	REAL pbase(klon), bbase(klon)
140	REAL dtvpdt1(klon, klev), dtvpdq1(klon, klev)
141	REAL dplcldt(klon), dplcldr(klon)
142	REAL qcondc(klon, klev)
143	REAL qtc(klon, klev)
144	REAL sigt(klon, klev)
145	REAL wd(klon)
146	REAL plim1(klon), plim2(klon), asupmax(klon, klev)
147	REAL supmax0(klon), asupmaxmin(klon)
148
149	REAL sigd(klon)
150	REAL zx_t, zdelta, zx_qs, zcor
151	REAL tau_cld_cv, coefw_cld_cv
152	REAL epmax_diag(klon) ! epmax_cape
153
154	! INTEGER iflag_mix
155	! SAVE iflag_mix
156	INTEGER noff, minorig
157	INTEGER i, k, itra
158	REAL qs(klon, klev), qs_wake(klon, klev)
159	REAL cbmf(klon), plcl(klon), plfc(klon), wbeff(klon)
160	REAL convoccur(klon)
161	!LF SAVE cbmf
162	!IM/JYG REAL, SAVE, ALLOCATABLE :: cbmf(:)
163	!!!$OMP THREADPRIVATE(cbmf)!
164	REAL cbmflast(klon)
165	INTEGER ifrst
166	SAVE ifrst
167	DATA ifrst/0/
168	!$OMP THREADPRIVATE(ifrst)
169
170
171	! Variables supplementaires liees au bilan d'energie
172	! Real paire(klon)
173	!LF Real ql(klon,klev)
174	! Save paire
175	!LF Save ql
176	!LF Real t1(klon,klev),q1(klon,klev)
177	!LF Save t1,q1
178	! Data paire /1./
179	REAL, SAVE, ALLOCATABLE :: ql(:, :), q1(:, :), t1(:, :)
180	!$OMP THREADPRIVATE(ql, q1, t1)
181
182	! Variables liees au bilan d'energie et d'enthAlpi
183	REAL ztsol(klon)
184	REAL h_vcol_tot, h_dair_tot, h_qw_tot, h_ql_tot, &
185	h_qs_tot, qw_tot, ql_tot, qs_tot, ec_tot
186	SAVE h_vcol_tot, h_dair_tot, h_qw_tot, h_ql_tot, &
187	h_qs_tot, qw_tot, ql_tot, qs_tot, ec_tot
188	!$OMP THREADPRIVATE(h_vcol_tot, h_dair_tot, h_qw_tot, h_ql_tot)
189	!$OMP THREADPRIVATE(h_qs_tot, qw_tot, ql_tot, qs_tot , ec_tot)
190	REAL d_h_vcol, d_h_dair, d_qt, d_qw, d_ql, d_qs, d_ec
191	REAL d_h_vcol_phy
192	REAL fs_bound, fq_bound
193	SAVE d_h_vcol_phy
194	!$OMP THREADPRIVATE(d_h_vcol_phy)
195	REAL zero_v(klon)
196	CHARACTER *15 ztit
197	INTEGER ip_ebil ! PRINT level for energy conserv. diag.
198	SAVE ip_ebil
199	DATA ip_ebil/2/
200	!$OMP THREADPRIVATE(ip_ebil)
201	INTEGER if_ebil ! level for energy conserv. dignostics
202	SAVE if_ebil
203	DATA if_ebil/2/
204	!$OMP THREADPRIVATE(if_ebil)
205	!+jld ec_conser
206	REAL d_t_ec(klon, klev) ! tendance du a la conersion Ec -> E thermique
207	REAL zrcpd
208	!-jld ec_conser
209	!LF
210	INTEGER nloc
211	LOGICAL, SAVE :: first = .TRUE.
212	!$OMP THREADPRIVATE(first)
213	INTEGER, SAVE :: itap, igout
214	!$OMP THREADPRIVATE(itap, igout)
215
216
217	include "YOMCST.h"
218	include "YOMCST2.h"
219	include "YOETHF.h"
220	include "FCTTRE.h"
221	!jyg<
222	include "conema3.h"
223	!>jyg
224
225	IF (first) THEN
226	! Allocate some variables LF 04/2008
227
228	!IM/JYG allocate(cbmf(klon))
229	ALLOCATE (ql(klon,klev))
230	ALLOCATE (t1(klon,klev))
231	ALLOCATE (q1(klon,klev))
232	!
233	convoccur(:) = 0.
234	!
235	itap = 0
236	igout = klon/2 + 1/klon
237	END IF
238	! Incrementer le compteur de la physique
239	itap = itap + 1
240
241	! Copy T into Tconv
242	DO k = 1, klev
243	DO i = 1, klon
244	Tconv(i, k) = t(i, k)
245	END DO
246	END DO
247
248	IF (if_ebil>=1) THEN
249	DO i = 1, klon
250	ztsol(i) = t(i, 1)
251	zero_v(i) = 0.
252	DO k = 1, klev
253	ql(i, k) = 0.
254	END DO
255	END DO
256	END IF
257
258	! ym
259	snow(:) = 0
260
261	! IF (ifrst .EQ. 0) THEN
262	! ifrst = 1
263	IF (first) THEN
264	first = .FALSE.
265
266	! ===========================================================================
267	! READ IN PARAMETERS FOR THE CLOSURE AND THE MIXING DISTRIBUTION
268	! ===========================================================================
269
270	IF (iflag_con==3) THEN
271	! CALL cv3_inicp()
272	CALL cv3_inip()
273	END IF
274
275	! ===========================================================================
276	! READ IN PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS
277	! ===========================================================================
278
279	! c$$$ open (56,file='supcrit.data')
280	! c$$$ read (56,*) Supcrit1, Supcrit2
281	! c$$$ close (56)
282
283	IF (prt_level>=10) WRITE (lunout, *) 'supcrit1, supcrit2', supcrit1, supcrit2
284
285	! ===========================================================================
286	! Initialisation pour les bilans d'eau et d'energie
287	! ===========================================================================
288	IF (if_ebil>=1) d_h_vcol_phy = 0.
289
290	DO i = 1, klon
291	cbmf(i) = 0.
292	!! plcl(i) = 0.
293	sigd(i) = 0.
294	END DO
295	END IF !(ifrst .EQ. 0)
296
297	! Initialisation a chaque pas de temps
298	plfc(:) = 0.
299	wbeff(:) = 100.
300	plcl(:) = 0.
301
302	DO k = 1, klev + 1
303	DO i = 1, klon
304	em_ph(i, k) = paprs(i, k)/100.0
305	pmflxr(i, k) = 0.
306	pmflxs(i, k) = 0.
307	END DO
308	END DO
309
310	DO k = 1, klev
311	DO i = 1, klon
312	em_p(i, k) = pplay(i, k)/100.0
313	END DO
314	END DO
315
316
317	! Feeding layer
318
319	em_sig1feed = 1.
320	!jyg<
321	! em_sig2feed = 0.97
322	em_sig2feed = cvl_sig2feed
323	!>jyg
324	! em_sig2feed = 0.8
325	! Relative Weight densities
326	DO k = 1, klev
327	em_wght(k) = 1.
328	END DO
329	!CRtest: couche alim des tehrmiques ponderee par a*
330	! DO i = 1, klon
331	! do k=1,lalim_conv(i)
332	! em_wght(k)=wght_th(i,k)
333	! print*,'em_wght=',em_wght(k),wght_th(i,k)
334	! end do
335	! END DO
336
337	IF (iflag_con==4) THEN
338	DO k = 1, klev
339	DO i = 1, klon
340	zx_t = t(i, k)
341	zdelta = max(0., sign(1.,rtt-zx_t))
342	zx_qs = min(0.5, r2es*foeew(zx_t,zdelta)/em_p(i,k)/100.0)
343	zcor = 1./(1.-retv*zx_qs)
344	qs(i, k) = zx_qs*zcor
345	END DO
346	DO i = 1, klon
347	zx_t = t_wake(i, k)
348	zdelta = max(0., sign(1.,rtt-zx_t))
349	zx_qs = min(0.5, r2es*foeew(zx_t,zdelta)/em_p(i,k)/100.0)
350	zcor = 1./(1.-retv*zx_qs)
351	qs_wake(i, k) = zx_qs*zcor
352	END DO
353	END DO
354	ELSE ! iflag_con=3 (modif de puristes qui fait la diffce pour la convergence numerique)
355	DO k = 1, klev
356	DO i = 1, klon
357	zx_t = t(i, k)
358	zdelta = max(0., sign(1.,rtt-zx_t))
359	zx_qs = r2es*foeew(zx_t, zdelta)/em_p(i, k)/100.0
360	zx_qs = min(0.5, zx_qs)
361	zcor = 1./(1.-retv*zx_qs)
362	zx_qs = zx_qs*zcor
363	qs(i, k) = zx_qs
364	END DO
365	DO i = 1, klon
366	zx_t = t_wake(i, k)
367	zdelta = max(0., sign(1.,rtt-zx_t))
368	zx_qs = r2es*foeew(zx_t, zdelta)/em_p(i, k)/100.0
369	zx_qs = min(0.5, zx_qs)
370	zcor = 1./(1.-retv*zx_qs)
371	zx_qs = zx_qs*zcor
372	qs_wake(i, k) = zx_qs
373	END DO
374	END DO
375	END IF ! iflag_con
376
377	! ------------------------------------------------------------------
378
379	! Main driver for convection:
380	! iflag_con=3 -> nvlle version de KE (JYG)
381	! iflag_con = 30 -> equivAlent to convect3
382	! iflag_con = 4 -> equivAlent to convect1/2
383
384
385	IF (iflag_con==30) THEN
386
387	! print *, '-> cv_driver' !jyg
388	CALL cv_driver(klon, klev, klevp1, ntra, iflag_con, &
389	t, q, qs, u, v, tra, &
390	em_p, em_ph, iflag, &
391	d_t, d_q, d_u, d_v, d_tra, rain, &
392	Vprecip, cbmf, sig1, w01, & !jyg
393	kbas, ktop, &
394	dtime, Ma, upwd, dnwd, dnwdbis, qcondc, wd, cape, &
395	da, phi, mp, phi2, d1a, dam, sij, clw, elij, & !RomP
396	evap, ep, epmlmMm, eplaMm, & !RomP
397	wdtrainA, wdtrainM, & !RomP
398	epmax_diag) ! epmax_cape
399	! print *, 'cv_driver ->' !jyg
400
401	DO i = 1, klon
402	cbmf(i) = Ma(i, kbas(i))
403	END DO
404
405	!RL
406	wght(:, :) = 0.
407	DO i = 1, klon
408	wght(i, 1) = 1.
409	END DO
410	!RL
411
412	ELSE
413
414	!LF necessary for gathered fields
415	nloc = klon
416	CALL cva_driver(klon, klev, klev+1, ntra, nloc, k_upper_cv, &
417	iflag_con, iflag_mix, iflag_ice_thermo, &
418	iflag_clos, ok_conserv_q, dtime, cvl_comp_threshold, &
419	t, q, qs, t_wake, q_wake, qs_wake, s_wake, u, v, tra, &
420	em_p, em_ph, &
421	Ale, Alp, omega, &
422	em_sig1feed, em_sig2feed, em_wght, &
423	iflag, d_t, d_q, d_u, d_v, d_tra, rain, kbas, ktop, &
424	cbmf, plcl, plfc, wbeff, sig1, w01, ptop2, sigd, &
425	Ma, mip, Vprecip, Vprecipi, upwd, dnwd, dnwdbis, qcondc, wd, &
426	cape, cin, tvp, &
427	dd_t, dd_q, plim1, plim2, asupmax, supmax0, &
428	asupmaxmin, lalim_conv, &
429	!AC!+!RomP+jyg
430	!! da,phi,mp,phi2,d1a,dam,sij,clw,elij, & ! RomP
431	!! evap,ep,epmlmMm,eplaMm, ! RomP
432	da, phi, mp, phi2, d1a, dam, sij, wght, & ! RomP+RL
433	clw, elij, evap, ep, epmlmMm, eplaMm, & ! RomP+RL
434	wdtrainA, wdtrainM, qtc, sigt, &
435	tau_cld_cv, coefw_cld_cv, & ! RomP,AJ
436	!AC!+!RomP+jyg
437	epmax_diag) ! epmax_cape
438	END IF
439	! ------------------------------------------------------------------
440	IF (prt_level>=10) WRITE (lunout, *) ' cva_driver -> cbmf,plcl,plfc,wbeff ', &
441	cbmf(1), plcl(1), plfc(1), wbeff(1)
442
443	DO i = 1, klon
444	rain(i) = rain(i)/86400.
445	rflag(i) = iflag(i)
446	END DO
447
448	DO k = 1, klev
449	DO i = 1, klon
450	d_t(i, k) = dtime*d_t(i, k)
451	d_q(i, k) = dtime*d_q(i, k)
452	d_u(i, k) = dtime*d_u(i, k)
453	d_v(i, k) = dtime*d_v(i, k)
454	END DO
455	END DO
456
457	IF (iflag_con==3) THEN
458	DO i = 1,klon
459	IF (wbeff(i) > 100. .OR. wbeff(i) == 0 .OR. iflag(i) > 3) THEN
460	wbeff(i) = 0.
461	convoccur(i) = 0.
462	ELSE
463	convoccur(i) = 1.
464	ENDIF
465	ENDDO
466	ENDIF
467
468	IF (iflag_con==30) THEN
469	DO itra = 1, ntra
470	DO k = 1, klev
471	DO i = 1, klon
472	!RL! d_tra(i,k,itra) =dtime*d_tra(i,k,itra)
473	d_tra(i, k, itra) = 0.
474	END DO
475	END DO
476	END DO
477	END IF
478
479	!!AC!
480	IF (iflag_con==3) THEN
481	DO itra = 1, ntra
482	DO k = 1, klev
483	DO i = 1, klon
484	!RL! d_tra(i,k,itra) =dtime*d_tra(i,k,itra)
485	d_tra(i, k, itra) = 0.
486	END DO
487	END DO
488	END DO
489	END IF
490	!!AC!
491
492	DO k = 1, klev
493	DO i = 1, klon
494	t1(i, k) = t(i, k) + d_t(i, k)
495	q1(i, k) = q(i, k) + d_q(i, k)
496	END DO
497	END DO
498	! !jyg
499	IF (iflag_con == 30 .OR. iflag_ice_thermo ==0) THEN
500	! --Separation neige/pluie (pour diagnostics) !jyg
501	DO k = 1, klev !jyg
502	DO i = 1, klon !jyg
503	IF (t1(i,k)<rtt) THEN !jyg
504	pmflxs(i, k) = Vprecip(i, k) !jyg
505	ELSE !jyg
506	pmflxr(i, k) = Vprecip(i, k) !jyg
507	END IF !jyg
508	END DO !jyg
509	END DO !jyg
510	ELSE
511	DO k = 1, klev !jyg
512	DO i = 1, klon !jyg
513	pmflxs(i, k) = Vprecipi(i, k) !jyg
514	pmflxr(i, k) = Vprecip(i, k)-Vprecipi(i, k) !jyg
515	END DO !jyg
516	END DO !jyg
517	ENDIF
518
519	! c IF (if_ebil.ge.2) THEN
520	! c ztit='after convect'
521	! c CALL diagetpq(paire,ztit,ip_ebil,2,2,dtime
522	! c e , t1,q1,ql,qs,u,v,paprs,pplay
523	! c s , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)
524	! c call diagphy(paire,ztit,ip_ebil
525	! c e , zero_v, zero_v, zero_v, zero_v, zero_v
526	! c e , zero_v, rain, zero_v, ztsol
527	! c e , d_h_vcol, d_qt, d_ec
528	! c s , fs_bound, fq_bound )
529	! c END IF
530
531
532	! les traceurs ne sont pas mis dans cette version de convect4:
533	IF (iflag_con==4) THEN
534	DO itra = 1, ntra
535	DO k = 1, klev
536	DO i = 1, klon
537	d_tra(i, k, itra) = 0.
538	END DO
539	END DO
540	END DO
541	END IF
542	! print*, 'concvl->: dd_t,dd_q ',dd_t(1,1),dd_q(1,1)
543
544	DO k = 1, klev
545	DO i = 1, klon
546	dtvpdt1(i, k) = 0.
547	dtvpdq1(i, k) = 0.
548	END DO
549	END DO
550	DO i = 1, klon
551	dplcldt(i) = 0.
552	dplcldr(i) = 0.
553	END DO
554
555	IF (prt_level>=20) THEN
556	DO k = 1, klev
557	! print*,'physiq apres_add_con i k it d_u d_v d_t d_q qdl0',igout, &
558	! k,itap,d_u_con(igout,k) ,d_v_con(igout,k), d_t_con(igout,k), &
559	! d_q_con(igout,k),dql0(igout,k)
560	! print*,'phys apres_add_con itap Ma cin ALE ALP wak t q undi t q', &
561	! itap,Ma(igout,k),cin(igout),ALE(igout), ALP(igout), &
562	! t_wake(igout,k),q_wake(igout,k),t_undi(igout,k),q_undi(igout,k)
563	! print*,'phy apres_add_con itap CON rain snow EMA wk1 wk2 Vpp mip', &
564	! itap,rain_con(igout),snow_con(igout),ema_work1(igout,k), &
565	! ema_work2(igout,k),Vprecip(igout,k), mip(igout,k)
566	! print*,'phy apres_add_con itap upwd dnwd dnwd0 cape tvp Tconv ', &
567	! itap,upwd(igout,k),dnwd(igout,k),dnwd0(igout,k),cape(igout), &
568	! tvp(igout,k),Tconv(igout,k)
569	! print*,'phy apres_add_con itap dtvpdt dtvdq dplcl dplcldr qcondc', &
570	! itap,dtvpdt1(igout,k),dtvpdq1(igout,k),dplcldt(igout), &
571	! dplcldr(igout),qcondc(igout,k)
572	! print*,'phy apres_add_con itap wd pmflxr Kpmflxr Kp1 Kpmflxs Kp1', &
573	! itap,wd(igout),pmflxr(igout,k),pmflxr(igout,k+1),pmflxs(igout,k), &
574	! pmflxs(igout,k+1)
575	! print*,'phy apres_add_con itap da phi mp ftd fqd lalim wgth', &
576	! itap,da(igout,k),phi(igout,k,k),mp(igout,k),ftd(igout,k), &
577	! fqd(igout,k),lalim_conv(igout),wght_th(igout,k)
578	END DO
579	END IF !(prt_level.EQ.20) THEN
580
581	RETURN
582	END SUBROUTINE concvl
583

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format