Context Navigation

cv_driver.F @ 717

Last change on this file since 717 was 619, checked in by lmdzadmin, 20 years ago
Rajout convection Kerry Emanuel pour traceurs- MAF+JYG
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 23.7 KB

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

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

Download in other formats:

Original Format