Context Navigation

cv_driver.F @ 520

Last change on this file since 520 was 418, checked in by (none), 22 years ago
This commit was manufactured by cvs2svn to create branch 'rel-LF'.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 22.6 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	do 20 k=1,nd
312	do 10 i=1,len
313	ft1(i,k)=0.0
314	fq1(i,k)=0.0
315	fu1(i,k)=0.0
316	fv1(i,k)=0.0
317	tvp1(i,k)=0.0
318	tp1(i,k)=0.0
319	clw1(i,k)=0.0
320	gz1(i,k) = 0.
321
322	Ma1(i,k)=0.0
323	upwd1(i,k)=0.0
324	dnwd1(i,k)=0.0
325	dnwd01(i,k)=0.0
326	qcondc1(i,k)=0.0
327	10 continue
328	20 continue
329
330	do 30 j=1,ntra
331	do 31 k=1,nd
332	do 32 i=1,len
333	ftra1(i,k,j)=0.0
334	32 continue
335	31 continue
336	30 continue
337
338	do 60 i=1,len
339	precip1(i)=0.0
340	iflag1(i)=0
341	wd1(i)=0.0
342	cape1(i)=0.0
343	60 continue
344
345	if (iflag_con.eq.3) then
346	do il=1,len
347	sig1(il,nd)=sig1(il,nd)+1.
348	sig1(il,nd)=amin1(sig1(il,nd),12.1)
349	enddo
350	endif
351
352	!--------------------------------------------------------------------
353	! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
354	!--------------------------------------------------------------------
355
356	if (iflag_con.eq.3) then
357	CALL cv3_prelim(len,nd,ndp1,t1,q1,p1,ph1 ! nd->na
358	o ,lv1,cpn1,tv1,gz1,h1,hm1,th1)
359	endif
360
361	if (iflag_con.eq.4) then
362	CALL cv_prelim(len,nd,ndp1,t1,q1,p1,ph1
363	o ,lv1,cpn1,tv1,gz1,h1,hm1)
364	endif
365
366	!--------------------------------------------------------------------
367	! --- CONVECTIVE FEED
368	!--------------------------------------------------------------------
369
370	if (iflag_con.eq.3) then
371	CALL cv3_feed(len,nd,t1,q1,qs1,p1,ph1,hm1,gz1 ! nd->na
372	o ,nk1,icb1,icbmax,iflag1,tnk1,qnk1,gznk1,plcl1)
373	endif
374
375	if (iflag_con.eq.4) then
376	CALL cv_feed(len,nd,t1,q1,qs1,p1,hm1,gz1
377	o ,nk1,icb1,icbmax,iflag1,tnk1,qnk1,gznk1,plcl1)
378	endif
379
380	!--------------------------------------------------------------------
381	! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part
382	! (up through ICB for convect4, up through ICB+1 for convect3)
383	! Calculates the lifted parcel virtual temperature at nk, the
384	! actual temperature, and the adiabatic liquid water content.
385	!--------------------------------------------------------------------
386
387	if (iflag_con.eq.3) then
388	CALL cv3_undilute1(len,nd,t1,q1,qs1,gz1,plcl1,p1,nk1,icb1 ! nd->na
389	o ,tp1,tvp1,clw1,icbs1)
390	endif
391
392	if (iflag_con.eq.4) then
393	CALL cv_undilute1(len,nd,t1,q1,qs1,gz1,p1,nk1,icb1,icbmax
394	: ,tp1,tvp1,clw1)
395	endif
396
397	!-------------------------------------------------------------------
398	! --- TRIGGERING
399	!-------------------------------------------------------------------
400
401	if (iflag_con.eq.3) then
402	CALL cv3_trigger(len,nd,icb1,plcl1,p1,th1,tv1,tvp1 ! nd->na
403	o ,pbase1,buoybase1,iflag1,sig1,w01)
404	endif
405
406	if (iflag_con.eq.4) then
407	CALL cv_trigger(len,nd,icb1,cbmf1,tv1,tvp1,iflag1)
408	endif
409
410	!=====================================================================
411	! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY
412	!=====================================================================
413
414	ncum=0
415	do 400 i=1,len
416	if(iflag1(i).eq.0)then
417	ncum=ncum+1
418	idcum(ncum)=i
419	endif
420	400 continue
421
422	c print*,'klon, ncum = ',len,ncum
423
424	IF (ncum.gt.0) THEN
425
426	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
427	! --- COMPRESS THE FIELDS
428	! (-> vectorization over convective gridpoints)
429	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
430
431	if (iflag_con.eq.3) then
432	CALL cv3_compress( len,nloc,ncum,nd,ntra
433	: ,iflag1,nk1,icb1,icbs1
434	: ,plcl1,tnk1,qnk1,gznk1,pbase1,buoybase1
435	: ,t1,q1,qs1,u1,v1,gz1,th1
436	: ,tra1
437	: ,h1,lv1,cpn1,p1,ph1,tv1,tp1,tvp1,clw1
438	: ,sig1,w01
439	o ,iflag,nk,icb,icbs
440	o ,plcl,tnk,qnk,gznk,pbase,buoybase
441	o ,t,q,qs,u,v,gz,th
442	o ,tra
443	o ,h,lv,cpn,p,ph,tv,tp,tvp,clw
444	o ,sig,w0 )
445	endif
446
447	if (iflag_con.eq.4) then
448	CALL cv_compress( len,nloc,ncum,nd
449	: ,iflag1,nk1,icb1
450	: ,cbmf1,plcl1,tnk1,qnk1,gznk1
451	: ,t1,q1,qs1,u1,v1,gz1
452	: ,h1,lv1,cpn1,p1,ph1,tv1,tp1,tvp1,clw1
453	o ,iflag,nk,icb
454	o ,cbmf,plcl,tnk,qnk,gznk
455	o ,t,q,qs,u,v,gz,h,lv,cpn,p,ph,tv,tp,tvp,clw
456	o ,dph )
457	endif
458
459	!-------------------------------------------------------------------
460	! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
461	! --- FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
462	! --- &
463	! --- COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
464	! --- FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
465	! --- &
466	! --- FIND THE LEVEL OF NEUTRAL BUOYANCY
467	!-------------------------------------------------------------------
468
469	if (iflag_con.eq.3) then
470	CALL cv3_undilute2(nloc,ncum,nd,icb,icbs,nk !na->nd
471	: ,tnk,qnk,gznk,t,q,qs,gz
472	: ,p,h,tv,lv,pbase,buoybase,plcl
473	o ,inb,tp,tvp,clw,hp,ep,sigp,buoy)
474	endif
475
476	if (iflag_con.eq.4) then
477	CALL cv_undilute2(nloc,ncum,nd,icb,nk
478	: ,tnk,qnk,gznk,t,q,qs,gz
479	: ,p,dph,h,tv,lv
480	o ,inb,inbis,tp,tvp,clw,hp,ep,sigp,frac)
481	endif
482
483	!-------------------------------------------------------------------
484	! --- CLOSURE
485	!-------------------------------------------------------------------
486
487	if (iflag_con.eq.3) then
488	CALL cv3_closure(nloc,ncum,nd,icb,inb ! na->nd
489	: ,pbase,p,ph,tv,buoy
490	o ,sig,w0,cape,m)
491	endif
492
493	if (iflag_con.eq.4) then
494	CALL cv_closure(nloc,ncum,nd,nk,icb
495	: ,tv,tvp,p,ph,dph,plcl,cpn
496	o ,iflag,cbmf)
497	endif
498
499	!-------------------------------------------------------------------
500	! --- MIXING
501	!-------------------------------------------------------------------
502
503	if (iflag_con.eq.3) then
504	CALL cv3_mixing(nloc,ncum,nd,nd,ntra,icb,nk,inb ! na->nd
505	: ,ph,t,q,qs,u,v,tra,h,lv,qnk
506	: ,hp,tv,tvp,ep,clw,m,sig
507	o ,ment,qent,uent,vent,sij,elij,ments,qents,traent)
508	endif
509
510	if (iflag_con.eq.4) then
511	CALL cv_mixing(nloc,ncum,nd,icb,nk,inb,inbis
512	: ,ph,t,q,qs,u,v,h,lv,qnk
513	: ,hp,tv,tvp,ep,clw,cbmf
514	o ,m,ment,qent,uent,vent,nent,sij,elij)
515	endif
516
517	!-------------------------------------------------------------------
518	! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
519	!-------------------------------------------------------------------
520
521	if (iflag_con.eq.3) then
522	CALL cv3_unsat(nloc,ncum,nd,nd,ntra,icb,inb ! na->nd
523	: ,t,q,qs,gz,u,v,tra,p,ph
524	: ,th,tv,lv,cpn,ep,sigp,clw
525	: ,m,ment,elij,delt,plcl
526	o ,mp,qp,up,vp,trap,wt,water,evap,b)
527	endif
528
529	if (iflag_con.eq.4) then
530	CALL cv_unsat(nloc,ncum,nd,inb,t,q,qs,gz,u,v,p,ph
531	: ,h,lv,ep,sigp,clw,m,ment,elij
532	o ,iflag,mp,qp,up,vp,wt,water,evap)
533	endif
534
535	!-------------------------------------------------------------------
536	! --- YIELD
537	! (tendencies, precipitation, variables of interface with other
538	! processes, etc)
539	!-------------------------------------------------------------------
540
541	if (iflag_con.eq.3) then
542	CALL cv3_yield(nloc,ncum,nd,nd,ntra ! na->nd
543	: ,icb,inb,delt
544	: ,t,q,u,v,tra,gz,p,ph,h,hp,lv,cpn,th
545	: ,ep,clw,m,tp,mp,qp,up,vp,trap
546	: ,wt,water,evap,b
547	: ,ment,qent,uent,vent,nent,elij,traent,sig
548	: ,tv,tvp
549	o ,iflag,precip,ft,fq,fu,fv,ftra
550	o ,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd)
551	endif
552
553	if (iflag_con.eq.4) then
554	CALL cv_yield(nloc,ncum,nd,nk,icb,inb,delt
555	: ,t,q,u,v,gz,p,ph,h,hp,lv,cpn
556	: ,ep,clw,frac,m,mp,qp,up,vp
557	: ,wt,water,evap
558	: ,ment,qent,uent,vent,nent,elij
559	: ,tv,tvp
560	o ,iflag,wd,qprime,tprime
561	o ,precip,cbmf,ft,fq,fu,fv,Ma,qcondc)
562	endif
563
564	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
565	! --- UNCOMPRESS THE FIELDS
566	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
567
568
569	if (iflag_con.eq.3) then
570	CALL cv3_uncompress(nloc,len,ncum,nd,ntra,idcum
571	: ,iflag
572	: ,precip,sig,w0
573	: ,ft,fq,fu,fv,ftra
574	: ,Ma,upwd,dnwd,dnwd0,qcondc,wd,cape
575	o ,iflag1
576	o ,precip1,sig1,w01
577	o ,ft1,fq1,fu1,fv1,ftra1
578	o ,Ma1,upwd1,dnwd1,dnwd01,qcondc1,wd1,cape1 )
579	endif
580
581	if (iflag_con.eq.4) then
582	CALL cv_uncompress(nloc,len,ncum,nd,idcum
583	: ,iflag
584	: ,precip,cbmf
585	: ,ft,fq,fu,fv
586	: ,Ma,qcondc
587	o ,iflag1
588	o ,precip1,cbmf1
589	o ,ft1,fq1,fu1,fv1
590	o ,Ma1,qcondc1 )
591	endif
592
593	ENDIF ! ncum>0
594
595	9999 continue
596
597	return
598	end
599
600	!==================================================================
601	SUBROUTINE cv_flag
602	implicit none
603
604	#include "cvflag.h"
605
606	c -- si .TRUE., on rend la gravite plus explicite et eventuellement
607	c differente de 10.0 dans convect3:
608	cvflag_grav = .FALSE.
609
610	return
611	end
612
613	!==================================================================
614	SUBROUTINE cv_thermo(iflag_con)
615	implicit none
616
617	c-------------------------------------------------------------
618	c Set thermodynamical constants for convectL
619	c-------------------------------------------------------------
620
621	#include "YOMCST.h"
622	#include "cvthermo.h"
623
624	integer iflag_con
625
626
627	c original set from convect:
628	if (iflag_con.eq.4) then
629	cpd=1005.7
630	cpv=1870.0
631	cl=4190.0
632	rrv=461.5
633	rrd=287.04
634	lv0=2.501E6
635	g=9.8
636	t0=273.15
637	grav=g
638	endif
639
640	c constants consistent with LMDZ:
641	if (iflag_con.eq.3) then
642	cpd = RCPD
643	cpv = RCPV
644	cl = RCW
645	rrv = RV
646	rrd = RD
647	lv0 = RLVTT
648	g = RG ! not used in convect3
649	c ori t0 = RTT
650	t0 = 273.15 ! convect3 (RTT=273.16)
651	grav= 10. ! implicitely or explicitely used in convect3
652	endif
653
654	rowl=1000.0 !(a quelle variable de YOMCST cela correspond-il?)
655
656	clmcpv=cl-cpv
657	clmcpd=cl-cpd
658	cpdmcp=cpd-cpv
659	cpvmcpd=cpv-cpd
660	cpvmcl=cl-cpv ! for convect3
661	eps=rrd/rrv
662	epsi=1.0/eps
663	epsim1=epsi-1.0
664	c ginv=1.0/g
665	ginv=1.0/grav
666	hrd=0.5*rrd
667
668	return
669	end
670

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

Download in other formats:

Original Format