Context Navigation

cv_driver.F @ 2052

Last change on this file since 2052 was 546, checked in by lmdzadmin, 20 years ago
Suite des adaptations au portage sur SGI et VPP pour Prism (Caubel & Demory) LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 22.7 KB

Line
1	SUBROUTINE cv_driver(len,nd,ndp1,ntra,iflag_con,
2	& t1,q1,qs1,u1,v1,tra1,
3	& p1,ph1,iflag1,ft1,fq1,fu1,fv1,ftra1,
4	& precip1,
5	& cbmf1,sig1,w01,
6	& delt,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1)
7	C
8	implicit none
9	C
10	C.............................START PROLOGUE............................
11	C
12	C PARAMETERS:
13	C Name Type Usage Description
14	C ---------- ---------- ------- ----------------------------
15	C
16	C len Integer Input first (i) dimension
17	C nd Integer Input vertical (k) dimension
18	C ndp1 Integer Input nd + 1
19	C ntra Integer Input number of tracors
20	C iflag_con Integer Input version of convect (3/4)
21	C t1 Real Input temperature
22	C q1 Real Input specific hum
23	C qs1 Real Input sat specific hum
24	C u1 Real Input u-wind
25	C v1 Real Input v-wind
26	C tra1 Real Input tracors
27	C p1 Real Input full level pressure
28	C ph1 Real Input half level pressure
29	C iflag1 Integer Output flag for Emanuel conditions
30	C ft1 Real Output temp tend
31	C fq1 Real Output spec hum tend
32	C fu1 Real Output u-wind tend
33	C fv1 Real Output v-wind tend
34	C ftra1 Real Output tracor tend
35	C precip1 Real Output precipitation
36	C cbmf1 Real Output cloud base mass flux
37	C sig1 Real In/Out section adiabatic updraft
38	C w01 Real In/Out vertical velocity within adiab updraft
39	C delt Real Input time step
40	C Ma1 Real Output mass flux adiabatic updraft
41	C upwd1 Real Output total upward mass flux (adiab+mixed)
42	C dnwd1 Real Output saturated downward mass flux (mixed)
43	C dnwd01 Real Output unsaturated downward mass flux
44	C qcondc1 Real Output in-cld mixing ratio of condensed water
45	C wd1 Real Output downdraft velocity scale for sfc fluxes
46	C cape1 Real Output CAPE
47	C
48	C S. Bony, Mar 2002:
49	C * Several modules corresponding to different physical processes
50	C * Several versions of convect may be used:
51	C - iflag_con=3: version lmd (previously named convect3)
52	C - iflag_con=4: version 4.3b (vect. version, previously convect1/2)
53	C + tard: - iflag_con=5: version lmd with ice (previously named convectg)
54	C S. Bony, Oct 2002:
55	C * Vectorization of convect3 (ie version lmd)
56	C
57	C..............................END PROLOGUE.............................
58	c
59	c
60	#include "dimensions.h"
61	#include "dimphy.h"
62
63	integer len
64	integer nd
65	integer ndp1
66	integer noff
67	integer iflag_con
68	integer ntra
69	real t1(len,nd)
70	real q1(len,nd)
71	real qs1(len,nd)
72	real u1(len,nd)
73	real v1(len,nd)
74	real p1(len,nd)
75	real ph1(len,ndp1)
76	integer iflag1(len)
77	real ft1(len,nd)
78	real fq1(len,nd)
79	real fu1(len,nd)
80	real fv1(len,nd)
81	real precip1(len)
82	real cbmf1(len)
83	real Ma1(len,nd)
84	real upwd1(len,nd)
85	real dnwd1(len,nd)
86	real dnwd01(len,nd)
87
88	real qcondc1(len,nd) ! cld
89	real wd1(len) ! gust
90	real cape1(len)
91
92	real tra1(len,nd,ntra)
93	real ftra1(len,nd,ntra)
94
95	real delt
96
97	!-------------------------------------------------------------------
98	! --- ARGUMENTS
99	!-------------------------------------------------------------------
100	! --- On input:
101	!
102	! t: Array of absolute temperature (K) of dimension ND, with first
103	! index corresponding to lowest model level. Note that this array
104	! will be altered by the subroutine if dry convective adjustment
105	! occurs and if IPBL is not equal to 0.
106	!
107	! q: Array of specific humidity (gm/gm) of dimension ND, with first
108	! index corresponding to lowest model level. Must be defined
109	! at same grid levels as T. Note that this array will be altered
110	! if dry convective adjustment occurs and if IPBL is not equal to 0.
111	!
112	! qs: Array of saturation specific humidity of dimension ND, with first
113	! index corresponding to lowest model level. Must be defined
114	! at same grid levels as T. Note that this array will be altered
115	! if dry convective adjustment occurs and if IPBL is not equal to 0.
116	!
117	! u: Array of zonal wind velocity (m/s) of dimension ND, witth first
118	! index corresponding with the lowest model level. Defined at
119	! same levels as T. Note that this array will be altered if
120	! dry convective adjustment occurs and if IPBL is not equal to 0.
121	!
122	! v: Same as u but for meridional velocity.
123	!
124	! tra: Array of passive tracer mixing ratio, of dimensions (ND,NTRA),
125	! where NTRA is the number of different tracers. If no
126	! convective tracer transport is needed, define a dummy
127	! input array of dimension (ND,1). Tracers are defined at
128	! same vertical levels as T. Note that this array will be altered
129	! if dry convective adjustment occurs and if IPBL is not equal to 0.
130	!
131	! p: Array of pressure (mb) of dimension ND, with first
132	! index corresponding to lowest model level. Must be defined
133	! at same grid levels as T.
134	!
135	! ph: Array of pressure (mb) of dimension ND+1, with first index
136	! corresponding to lowest level. These pressures are defined at
137	! levels intermediate between those of P, T, Q and QS. The first
138	! value of PH should be greater than (i.e. at a lower level than)
139	! the first value of the array P.
140	!
141	! nl: The maximum number of levels to which convection can penetrate, plus 1.
142	! NL MUST be less than or equal to ND-1.
143	!
144	! delt: The model time step (sec) between calls to CONVECT
145	!
146	!----------------------------------------------------------------------------
147	! --- On Output:
148	!
149	! iflag: An output integer whose value denotes the following:
150	! VALUE INTERPRETATION
151	! ----- --------------
152	! 0 Moist convection occurs.
153	! 1 Moist convection occurs, but a CFL condition
154	! on the subsidence warming is violated. This
155	! does not cause the scheme to terminate.
156	! 2 Moist convection, but no precip because ep(inb) lt 0.0001
157	! 3 No moist convection because new cbmf is 0 and old cbmf is 0.
158	! 4 No moist convection; atmosphere is not
159	! unstable
160	! 6 No moist convection because ihmin le minorig.
161	! 7 No moist convection because unreasonable
162	! parcel level temperature or specific humidity.
163	! 8 No moist convection: lifted condensation
164	! level is above the 200 mb level.
165	! 9 No moist convection: cloud base is higher
166	! then the level NL-1.
167	!
168	! ft: Array of temperature tendency (K/s) of dimension ND, defined at same
169	! grid levels as T, Q, QS and P.
170	!
171	! fq: Array of specific humidity tendencies ((gm/gm)/s) of dimension ND,
172	! defined at same grid levels as T, Q, QS and P.
173	!
174	! fu: Array of forcing of zonal velocity (m/s^2) of dimension ND,
175	! defined at same grid levels as T.
176	!
177	! fv: Same as FU, but for forcing of meridional velocity.
178	!
179	! ftra: Array of forcing of tracer content, in tracer mixing ratio per
180	! second, defined at same levels as T. Dimensioned (ND,NTRA).
181	!
182	! precip: Scalar convective precipitation rate (mm/day).
183	!
184	! wd: A convective downdraft velocity scale. For use in surface
185	! flux parameterizations. See convect.ps file for details.
186	!
187	! tprime: A convective downdraft temperature perturbation scale (K).
188	! For use in surface flux parameterizations. See convect.ps
189	! file for details.
190	!
191	! qprime: A convective downdraft specific humidity
192	! perturbation scale (gm/gm).
193	! For use in surface flux parameterizations. See convect.ps
194	! file for details.
195	!
196	! cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST
197	! BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT
198	! ITS NEXT CALL. That is, the value of CBMF must be "remembered"
199	! by the calling program between calls to CONVECT.
200	!
201	! det: Array of detrainment mass flux of dimension ND.
202	!
203	!-------------------------------------------------------------------
204	c
205	c Local arrays
206	c
207
208	integer i,k,n,il,j
209	integer icbmax
210	integer nk1(klon)
211	integer icb1(klon)
212	integer icbs1(klon)
213
214	real plcl1(klon)
215	real tnk1(klon)
216	real qnk1(klon)
217	real gznk1(klon)
218	real pnk1(klon)
219	real qsnk1(klon)
220	real pbase1(klon)
221	real buoybase1(klon)
222
223	real lv1(klon,klev)
224	real cpn1(klon,klev)
225	real tv1(klon,klev)
226	real gz1(klon,klev)
227	real hm1(klon,klev)
228	real h1(klon,klev)
229	real tp1(klon,klev)
230	real tvp1(klon,klev)
231	real clw1(klon,klev)
232	real sig1(klon,klev)
233	real w01(klon,klev)
234	real th1(klon,klev)
235	c
236	integer ncum
237	c
238	c (local) compressed fields:
239	c
240	integer nloc
241	parameter (nloc=klon) ! pour l'instant
242
243	integer idcum(nloc)
244	integer iflag(nloc),nk(nloc),icb(nloc)
245	integer nent(nloc,klev)
246	integer icbs(nloc)
247	integer inb(nloc), inbis(nloc)
248
249	real cbmf(nloc),plcl(nloc),tnk(nloc),qnk(nloc),gznk(nloc)
250	real t(nloc,klev),q(nloc,klev),qs(nloc,klev)
251	real u(nloc,klev),v(nloc,klev)
252	real gz(nloc,klev),h(nloc,klev),lv(nloc,klev),cpn(nloc,klev)
253	real p(nloc,klev),ph(nloc,klev+1),tv(nloc,klev),tp(nloc,klev)
254	real clw(nloc,klev)
255	real dph(nloc,klev)
256	real pbase(nloc), buoybase(nloc), th(nloc,klev)
257	real tvp(nloc,klev)
258	real sig(nloc,klev), w0(nloc,klev)
259	real hp(nloc,klev), ep(nloc,klev), sigp(nloc,klev)
260	real frac(nloc), buoy(nloc,klev)
261	real cape(nloc)
262	real m(nloc,klev), ment(nloc,klev,klev), qent(nloc,klev,klev)
263	real uent(nloc,klev,klev), vent(nloc,klev,klev)
264	real ments(nloc,klev,klev), qents(nloc,klev,klev)
265	real sij(nloc,klev,klev), elij(nloc,klev,klev)
266	real mp(nloc,klev), qp(nloc,klev), up(nloc,klev), vp(nloc,klev)
267	real wt(nloc,klev), water(nloc,klev), evap(nloc,klev)
268	real b(nloc,klev), ft(nloc,klev), fq(nloc,klev)
269	real fu(nloc,klev), fv(nloc,klev)
270	real upwd(nloc,klev), dnwd(nloc,klev), dnwd0(nloc,klev)
271	real Ma(nloc,klev), mike(nloc,klev), tls(nloc,klev)
272	real tps(nloc,klev), qprime(nloc), tprime(nloc)
273	real precip(nloc)
274	real tra(nloc,klev,ntra), trap(nloc,klev,ntra)
275	real ftra(nloc,klev,ntra), traent(nloc,klev,klev,ntra)
276	real qcondc(nloc,klev) ! cld
277	real wd(nloc) ! gust
278
279	!-------------------------------------------------------------------
280	! --- SET CONSTANTS AND PARAMETERS
281	!-------------------------------------------------------------------
282
283	c -- set simulation flags:
284	c (common cvflag)
285
286	CALL cv_flag
287
288	c -- set thermodynamical constants:
289	c (common cvthermo)
290
291	CALL cv_thermo(iflag_con)
292
293	c -- set convect parameters
294	c
295	c includes microphysical parameters and parameters that
296	c control the rate of approach to quasi-equilibrium)
297	c (common cvparam)
298
299	if (iflag_con.eq.3) then
300	CALL cv3_param(nd,delt)
301	endif
302
303	if (iflag_con.eq.4) then
304	CALL cv_param(nd)
305	endif
306
307	!---------------------------------------------------------------------
308	! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS
309	!---------------------------------------------------------------------
310
311	!LF
312	qs(:,:) = 0.
313	qp(:,:) = 0.
314	ep = 0.
315	clw = 0.
316	qcondc = 0.
317	nent = 0
318
319	do 20 k=1,nd
320	do 10 i=1,len
321	ft1(i,k)=0.0
322	fq1(i,k)=0.0
323	fu1(i,k)=0.0
324	fv1(i,k)=0.0
325	tvp1(i,k)=0.0
326	tp1(i,k)=0.0
327	clw1(i,k)=0.0
328	gz1(i,k) = 0.
329
330	Ma1(i,k)=0.0
331	upwd1(i,k)=0.0
332	dnwd1(i,k)=0.0
333	dnwd01(i,k)=0.0
334	qcondc1(i,k)=0.0
335	10 continue
336	20 continue
337
338	do 30 j=1,ntra
339	do 31 k=1,nd
340	do 32 i=1,len
341	ftra1(i,k,j)=0.0
342	32 continue
343	31 continue
344	30 continue
345
346	do 60 i=1,len
347	precip1(i)=0.0
348	iflag1(i)=0
349	wd1(i)=0.0
350	cape1(i)=0.0
351	60 continue
352
353	if (iflag_con.eq.3) then
354	do il=1,len
355	sig1(il,nd)=sig1(il,nd)+1.
356	sig1(il,nd)=amin1(sig1(il,nd),12.1)
357	enddo
358	endif
359
360	!--------------------------------------------------------------------
361	! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
362	!--------------------------------------------------------------------
363
364	if (iflag_con.eq.3) then
365	CALL cv3_prelim(len,nd,ndp1,t1,q1,p1,ph1 ! nd->na
366	o ,lv1,cpn1,tv1,gz1,h1,hm1,th1)
367	endif
368
369	if (iflag_con.eq.4) then
370	CALL cv_prelim(len,nd,ndp1,t1,q1,p1,ph1
371	o ,lv1,cpn1,tv1,gz1,h1,hm1)
372	endif
373
374	!--------------------------------------------------------------------
375	! --- CONVECTIVE FEED
376	!--------------------------------------------------------------------
377
378	if (iflag_con.eq.3) then
379	CALL cv3_feed(len,nd,t1,q1,qs1,p1,ph1,hm1,gz1 ! nd->na
380	o ,nk1,icb1,icbmax,iflag1,tnk1,qnk1,gznk1,plcl1)
381	endif
382
383	if (iflag_con.eq.4) then
384	CALL cv_feed(len,nd,t1,q1,qs1,p1,hm1,gz1
385	o ,nk1,icb1,icbmax,iflag1,tnk1,qnk1,gznk1,plcl1)
386	endif
387
388	!--------------------------------------------------------------------
389	! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part
390	! (up through ICB for convect4, up through ICB+1 for convect3)
391	! Calculates the lifted parcel virtual temperature at nk, the
392	! actual temperature, and the adiabatic liquid water content.
393	!--------------------------------------------------------------------
394
395	if (iflag_con.eq.3) then
396	CALL cv3_undilute1(len,nd,t1,q1,qs1,gz1,plcl1,p1,nk1,icb1 ! nd->na
397	o ,tp1,tvp1,clw1,icbs1)
398	endif
399
400	if (iflag_con.eq.4) then
401	CALL cv_undilute1(len,nd,t1,q1,qs1,gz1,p1,nk1,icb1,icbmax
402	: ,tp1,tvp1,clw1)
403	endif
404
405	!-------------------------------------------------------------------
406	! --- TRIGGERING
407	!-------------------------------------------------------------------
408
409	if (iflag_con.eq.3) then
410	CALL cv3_trigger(len,nd,icb1,plcl1,p1,th1,tv1,tvp1 ! nd->na
411	o ,pbase1,buoybase1,iflag1,sig1,w01)
412	endif
413
414	if (iflag_con.eq.4) then
415	CALL cv_trigger(len,nd,icb1,cbmf1,tv1,tvp1,iflag1)
416	endif
417
418	!=====================================================================
419	! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY
420	!=====================================================================
421
422	ncum=0
423	do 400 i=1,len
424	if(iflag1(i).eq.0)then
425	ncum=ncum+1
426	idcum(ncum)=i
427	endif
428	400 continue
429
430	c print*,'klon, ncum = ',len,ncum
431
432	IF (ncum.gt.0) THEN
433
434	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
435	! --- COMPRESS THE FIELDS
436	! (-> vectorization over convective gridpoints)
437	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
438
439	if (iflag_con.eq.3) then
440	CALL cv3_compress( len,nloc,ncum,nd,ntra
441	: ,iflag1,nk1,icb1,icbs1
442	: ,plcl1,tnk1,qnk1,gznk1,pbase1,buoybase1
443	: ,t1,q1,qs1,u1,v1,gz1,th1
444	: ,tra1
445	: ,h1,lv1,cpn1,p1,ph1,tv1,tp1,tvp1,clw1
446	: ,sig1,w01
447	o ,iflag,nk,icb,icbs
448	o ,plcl,tnk,qnk,gznk,pbase,buoybase
449	o ,t,q,qs,u,v,gz,th
450	o ,tra
451	o ,h,lv,cpn,p,ph,tv,tp,tvp,clw
452	o ,sig,w0 )
453	endif
454
455	if (iflag_con.eq.4) then
456	CALL cv_compress( len,nloc,ncum,nd
457	: ,iflag1,nk1,icb1
458	: ,cbmf1,plcl1,tnk1,qnk1,gznk1
459	: ,t1,q1,qs1,u1,v1,gz1
460	: ,h1,lv1,cpn1,p1,ph1,tv1,tp1,tvp1,clw1
461	o ,iflag,nk,icb
462	o ,cbmf,plcl,tnk,qnk,gznk
463	o ,t,q,qs,u,v,gz,h,lv,cpn,p,ph,tv,tp,tvp,clw
464	o ,dph )
465	endif
466
467	!-------------------------------------------------------------------
468	! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
469	! --- FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
470	! --- &
471	! --- COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
472	! --- FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
473	! --- &
474	! --- FIND THE LEVEL OF NEUTRAL BUOYANCY
475	!-------------------------------------------------------------------
476
477	if (iflag_con.eq.3) then
478	CALL cv3_undilute2(nloc,ncum,nd,icb,icbs,nk !na->nd
479	: ,tnk,qnk,gznk,t,q,qs,gz
480	: ,p,h,tv,lv,pbase,buoybase,plcl
481	o ,inb,tp,tvp,clw,hp,ep,sigp,buoy)
482	endif
483
484	if (iflag_con.eq.4) then
485	CALL cv_undilute2(nloc,ncum,nd,icb,nk
486	: ,tnk,qnk,gznk,t,q,qs,gz
487	: ,p,dph,h,tv,lv
488	o ,inb,inbis,tp,tvp,clw,hp,ep,sigp,frac)
489	endif
490
491	!-------------------------------------------------------------------
492	! --- CLOSURE
493	!-------------------------------------------------------------------
494
495	if (iflag_con.eq.3) then
496	CALL cv3_closure(nloc,ncum,nd,icb,inb ! na->nd
497	: ,pbase,p,ph,tv,buoy
498	o ,sig,w0,cape,m)
499	endif
500
501	if (iflag_con.eq.4) then
502	CALL cv_closure(nloc,ncum,nd,nk,icb
503	: ,tv,tvp,p,ph,dph,plcl,cpn
504	o ,iflag,cbmf)
505	endif
506
507	!-------------------------------------------------------------------
508	! --- MIXING
509	!-------------------------------------------------------------------
510
511	if (iflag_con.eq.3) then
512	CALL cv3_mixing(nloc,ncum,nd,nd,ntra,icb,nk,inb ! na->nd
513	: ,ph,t,q,qs,u,v,tra,h,lv,qnk
514	: ,hp,tv,tvp,ep,clw,m,sig
515	o ,ment,qent,uent,vent,sij,elij,ments,qents,traent)
516	endif
517
518	if (iflag_con.eq.4) then
519	CALL cv_mixing(nloc,ncum,nd,icb,nk,inb,inbis
520	: ,ph,t,q,qs,u,v,h,lv,qnk
521	: ,hp,tv,tvp,ep,clw,cbmf
522	o ,m,ment,qent,uent,vent,nent,sij,elij)
523	endif
524
525	!-------------------------------------------------------------------
526	! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
527	!-------------------------------------------------------------------
528
529	if (iflag_con.eq.3) then
530	CALL cv3_unsat(nloc,ncum,nd,nd,ntra,icb,inb ! na->nd
531	: ,t,q,qs,gz,u,v,tra,p,ph
532	: ,th,tv,lv,cpn,ep,sigp,clw
533	: ,m,ment,elij,delt,plcl
534	o ,mp,qp,up,vp,trap,wt,water,evap,b)
535	endif
536
537	if (iflag_con.eq.4) then
538	CALL cv_unsat(nloc,ncum,nd,inb,t,q,qs,gz,u,v,p,ph
539	: ,h,lv,ep,sigp,clw,m,ment,elij
540	o ,iflag,mp,qp,up,vp,wt,water,evap)
541	endif
542
543	!-------------------------------------------------------------------
544	! --- YIELD
545	! (tendencies, precipitation, variables of interface with other
546	! processes, etc)
547	!-------------------------------------------------------------------
548
549	if (iflag_con.eq.3) then
550	CALL cv3_yield(nloc,ncum,nd,nd,ntra ! na->nd
551	: ,icb,inb,delt
552	: ,t,q,u,v,tra,gz,p,ph,h,hp,lv,cpn,th
553	: ,ep,clw,m,tp,mp,qp,up,vp,trap
554	: ,wt,water,evap,b
555	: ,ment,qent,uent,vent,nent,elij,traent,sig
556	: ,tv,tvp
557	o ,iflag,precip,ft,fq,fu,fv,ftra
558	o ,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd)
559	endif
560
561	if (iflag_con.eq.4) then
562	CALL cv_yield(nloc,ncum,nd,nk,icb,inb,delt
563	: ,t,q,u,v,gz,p,ph,h,hp,lv,cpn
564	: ,ep,clw,frac,m,mp,qp,up,vp
565	: ,wt,water,evap
566	: ,ment,qent,uent,vent,nent,elij
567	: ,tv,tvp
568	o ,iflag,wd,qprime,tprime
569	o ,precip,cbmf,ft,fq,fu,fv,Ma,qcondc)
570	endif
571
572	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
573	! --- UNCOMPRESS THE FIELDS
574	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
575
576
577	if (iflag_con.eq.3) then
578	CALL cv3_uncompress(nloc,len,ncum,nd,ntra,idcum
579	: ,iflag
580	: ,precip,sig,w0
581	: ,ft,fq,fu,fv,ftra
582	: ,Ma,upwd,dnwd,dnwd0,qcondc,wd,cape
583	o ,iflag1
584	o ,precip1,sig1,w01
585	o ,ft1,fq1,fu1,fv1,ftra1
586	o ,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1 )
587	endif
588
589	if (iflag_con.eq.4) then
590	CALL cv_uncompress(nloc,len,ncum,nd,idcum
591	: ,iflag
592	: ,precip,cbmf
593	: ,ft,fq,fu,fv
594	: ,Ma,qcondc
595	o ,iflag1
596	o ,precip1,cbmf1
597	o ,ft1,fq1,fu1,fv1
598	o ,Ma1,qcondc1 )
599	endif
600
601	ENDIF ! ncum>0
602
603	9999 continue
604
605	return
606	end
607
608	!==================================================================
609	SUBROUTINE cv_flag
610	implicit none
611
612	#include "cvflag.h"
613
614	c -- si .TRUE., on rend la gravite plus explicite et eventuellement
615	c differente de 10.0 dans convect3:
616	cvflag_grav = .FALSE.
617
618	return
619	end
620
621	!==================================================================
622	SUBROUTINE cv_thermo(iflag_con)
623	implicit none
624
625	c-------------------------------------------------------------
626	c Set thermodynamical constants for convectL
627	c-------------------------------------------------------------
628
629	#include "YOMCST.h"
630	#include "cvthermo.h"
631
632	integer iflag_con
633
634
635	c original set from convect:
636	if (iflag_con.eq.4) then
637	cpd=1005.7
638	cpv=1870.0
639	cl=4190.0
640	rrv=461.5
641	rrd=287.04
642	lv0=2.501E6
643	g=9.8
644	t0=273.15
645	grav=g
646	endif
647
648	c constants consistent with LMDZ:
649	if (iflag_con.eq.3) then
650	cpd = RCPD
651	cpv = RCPV
652	cl = RCW
653	rrv = RV
654	rrd = RD
655	lv0 = RLVTT
656	g = RG ! not used in convect3
657	c ori t0 = RTT
658	t0 = 273.15 ! convect3 (RTT=273.16)
659	grav= 10. ! implicitely or explicitely used in convect3
660	endif
661
662	rowl=1000.0 !(a quelle variable de YOMCST cela correspond-il?)
663
664	clmcpv=cl-cpv
665	clmcpd=cl-cpd
666	cpdmcp=cpd-cpv
667	cpvmcpd=cpv-cpd
668	cpvmcl=cl-cpv ! for convect3
669	eps=rrd/rrv
670	epsi=1.0/eps
671	epsim1=epsi-1.0
672	c ginv=1.0/g
673	ginv=1.0/grav
674	hrd=0.5*rrd
675
676	return
677	end
678

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

Download in other formats:

Original Format