Context Navigation

cva_driver.F @ 5440

Last change on this file since 5440 was 1403, checked in by Laurent Fairhead, 14 years ago

Merged LMDZ4V5.0-dev branch changes r1292:r1399 to trunk.

Validation:
Validation consisted in compiling the HEAD revision of the trunk,
LMDZ4V5.0-dev branch and the merged sources and running different
configurations on local and SX8 machines comparing results.

Local machine: bench configuration, 32x24x11, gfortran

IPSLCM5A configuration (comparison between trunk and merged sources):
- numerical convergence on dynamical fields over 3 days
- start files are equivalent (except for RN and PB fields)
- daily history files equivalent

MH07 configuration, new physics package (comparison between LMDZ4V5.0-dev branch and merged sources):
- numerical convergence on dynamical fields over 3 days
- start files are equivalent (except for RN and PB fields)
- daily history files equivalent

SX8 machine (brodie), 96x95x39 on 4 processors:

IPSLCM5A configuration:
- start files are equivalent (except for RN and PB fields)
- monthly history files equivalent

MH07 configuration:
- start files are equivalent (except for RN and PB fields)
- monthly history files equivalent

Changes to the makegcm and create_make_gcm scripts to take into account
main programs in F90 files

Fusion de la branche LMDZ4V5.0-dev (r1292:r1399) au tronc principal

Validation:
La validation a consisté à compiler la HEAD de le trunk et de la banche
LMDZ4V5.0-dev et les sources fusionnées et de faire tourner le modéle selon
différentes configurations en local et sur SX8 et de comparer les résultats

En local: 32x24x11, config bench/gfortran

pour une config IPSLCM5A (comparaison tronc/fusion):
- convergence numérique sur les champs dynamiques après 3 jours
- restart et restartphy égaux (à part sur RN et Pb)
- fichiers histoire égaux

pour une config nlle physique (MH07) (comparaison LMDZ4v5.0-dev/fusion):
- convergence numérique sur les champs dynamiques après 3 jours
- restart et restartphy égaux
- fichiers histoire équivalents

Sur brodie, 96x95x39 sur 4 proc:

pour une config IPSLCM5A:
- restart et restartphy égaux (à part sur RN et PB)
- pas de différence dans les fichiers histmth.nc

pour une config MH07
- restart et restartphy égaux (à part sur RN et PB)
- pas de différence dans les fichiers histmth.nc

Changement sur makegcm et create_make-gcm pour pouvoir prendre en compte des
programmes principaux en *F90

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 34.4 KB

Rev	Line
[1403]	1	!
	2	! $Id: cva_driver.F 1403 2010-07-01 09:02:53Z evignon $
	3	!
[940]	4	SUBROUTINE cva_driver(len,nd,ndp1,ntra,nloc,
	5	& iflag_con,iflag_mix,
[879]	6	& iflag_clos,delt,
[1146]	7	& t1,q1,qs1,t1_wake,q1_wake,qs1_wake,s1_wake,
[879]	8	& u1,v1,tra1,
	9	& p1,ph1,
	10	& ALE1,ALP1,
	11	& sig1feed1,sig2feed1,wght1,
	12	o iflag1,ft1,fq1,fu1,fv1,ftra1,
	13	& precip1,kbas1,ktop1,cbmf1,
	14	& sig1,w01, !input/output
	15	& ptop21,sigd,
	16	& Ma1,mip1,Vprecip1,upwd1,dnwd1,dnwd01,
	17	& qcondc1,wd1,
	18	& cape1,cin1,tvp1,
	19	& ftd1,fqd1,
	20	& Plim11,Plim21,asupmax1,supmax01,asupmaxmin1
	21	& ,lalim_conv)
	22	***************************************************************
	23	* *
	24	* CV_DRIVER *
	25	* *
	26	* *
	27	* written by : Sandrine Bony-Lena , 17/05/2003, 11.19.41 *
	28	* modified by : *
	29	***************************************************************
	30	***************************************************************
	31	C
[940]	32	USE dimphy
[879]	33	implicit none
	34	C
	35	C.............................START PROLOGUE............................
	36	C
	37	C PARAMETERS:
	38	C Name Type Usage Description
	39	C ---------- ---------- ------- ----------------------------
	40	C
	41	C len Integer Input first (i) dimension
	42	C nd Integer Input vertical (k) dimension
	43	C ndp1 Integer Input nd + 1
	44	C ntra Integer Input number of tracors
	45	C iflag_con Integer Input version of convect (3/4)
	46	C iflag_mix Integer Input version of mixing (0/1/2)
	47	C iflag_clos Integer Input version of closure (0/1)
	48	C delt Real Input time step
	49	C t1 Real Input temperature (sat draught envt)
	50	C q1 Real Input specific hum (sat draught envt)
	51	C qs1 Real Input sat specific hum (sat draught envt)
	52	C t1_wake Real Input temperature (unsat draught envt)
	53	C q1_wake Real Input specific hum(unsat draught envt)
	54	C qs1_wake Real Input sat specific hum(unsat draughts envt)
[1146]	55	C s1_wake Real Input fractionnal area covered by wakes
[879]	56	C u1 Real Input u-wind
	57	C v1 Real Input v-wind
	58	C tra1 Real Input tracors
	59	C p1 Real Input full level pressure
	60	C ph1 Real Input half level pressure
	61	C ALE1 Real Input Available lifting Energy
	62	C ALP1 Real Input Available lifting Power
	63	C sig1feed1 Real Input sigma coord at lower bound of feeding layer
	64	C sig2feed1 Real Input sigma coord at upper bound of feeding layer
	65	C wght1 Real Input weight density determining the feeding mixture
	66	C iflag1 Integer Output flag for Emanuel conditions
	67	C ft1 Real Output temp tend
	68	C fq1 Real Output spec hum tend
	69	C fu1 Real Output u-wind tend
	70	C fv1 Real Output v-wind tend
	71	C ftra1 Real Output tracor tend
	72	C precip1 Real Output precipitation
	73	C kbas1 Integer Output cloud base level
	74	C ktop1 Integer Output cloud top level
	75	C cbmf1 Real Output cloud base mass flux
	76	C sig1 Real In/Out section adiabatic updraft
	77	C w01 Real In/Out vertical velocity within adiab updraft
	78	C ptop21 Real In/Out top of entraining zone
	79	C Ma1 Real Output mass flux adiabatic updraft
	80	C mip1 Real Output mass flux shed by the adiabatic updraft
	81	C Vprecip1 Real Output vertical profile of precipitations
	82	C upwd1 Real Output total upward mass flux (adiab+mixed)
	83	C dnwd1 Real Output saturated downward mass flux (mixed)
	84	C dnwd01 Real Output unsaturated downward mass flux
	85	C qcondc1 Real Output in-cld mixing ratio of condensed water
	86	C wd1 Real Output downdraft velocity scale for sfc fluxes
	87	C cape1 Real Output CAPE
	88	C cin1 Real Output CIN
	89	C tvp1 Real Output adiab lifted parcell virt temp
	90	C ftd1 Real Output precip temp tend
	91	C fqt1 Real Output precip spec hum tend
	92	C Plim11 Real Output
	93	C Plim21 Real Output
	94	C asupmax1 Real Output
	95	C supmax01 Real Output
	96	C asupmaxmin1 Real Output
	97	C S. Bony, Mar 2002:
	98	C * Several modules corresponding to different physical processes
	99	C * Several versions of convect may be used:
	100	C - iflag_con=3: version lmd (previously named convect3)
	101	C - iflag_con=4: version 4.3b (vect. version, previously convect1/2)
	102	C + tard: - iflag_con=5: version lmd with ice (previously named convectg)
	103	C S. Bony, Oct 2002:
	104	C * Vectorization of convect3 (ie version lmd)
	105	C
	106	C..............................END PROLOGUE.............................
	107	c
	108	c
	109	#include "dimensions.h"
[940]	110	ccccc#include "dimphy.h"
[1403]	111	include 'iniprint.h'
	112
[879]	113	c
	114	c Input
	115	integer len
	116	integer nd
	117	integer ndp1
	118	integer ntra
	119	integer iflag_con
	120	integer iflag_mix
	121	integer iflag_clos
	122	real delt
	123	real t1(len,nd)
	124	real q1(len,nd)
	125	real qs1(len,nd)
	126	real t1_wake(len,nd)
	127	real q1_wake(len,nd)
	128	real qs1_wake(len,nd)
[1146]	129	real s1_wake(len)
[879]	130	real u1(len,nd)
	131	real v1(len,nd)
	132	real tra1(len,nd,ntra)
	133	real p1(len,nd)
	134	real ph1(len,ndp1)
	135	real ALE1(len)
	136	real ALP1(len)
	137	real sig1feed1 ! pressure at lower bound of feeding layer
	138	real sig2feed1 ! pressure at upper bound of feeding layer
	139	real wght1(nd) ! weight density determining the feeding mixture
	140	c
	141	c Output
	142	integer iflag1(len)
	143	real ft1(len,nd)
	144	real fq1(len,nd)
	145	real fu1(len,nd)
	146	real fv1(len,nd)
	147	real ftra1(len,nd,ntra)
	148	real precip1(len)
	149	integer kbas1(len)
	150	integer ktop1(len)
	151	real cbmf1(len)
	152	! real cbmflast(len)
	153	real sig1(len,klev) !input/output
	154	real w01(len,klev) !input/output
	155	real ptop21(len)
	156	real Ma1(len,nd)
	157	real mip1(len,nd)
[1403]	158	! real Vprecip1(len,nd)
[1334]	159	real Vprecip1(len,nd+1)
[879]	160	real upwd1(len,nd)
	161	real dnwd1(len,nd)
	162	real dnwd01(len,nd)
	163	real qcondc1(len,nd) ! cld
	164	real wd1(len) ! gust
	165	real cape1(len)
	166	real cin1(len)
	167	real tvp1(len,nd)
	168	c
	169	real ftd1(len,nd)
	170	real fqd1(len,nd)
	171	real Plim11(len)
	172	real Plim21(len)
	173	real asupmax1(len,nd)
	174	real supmax01(len)
	175	real asupmaxmin1(len)
	176	integer lalim_conv(len)
	177	!-------------------------------------------------------------------
	178	! --- ARGUMENTS
	179	!-------------------------------------------------------------------
	180	! --- On input:
	181	!
	182	! t: Array of absolute temperature (K) of dimension ND, with first
	183	! index corresponding to lowest model level. Note that this array
	184	! will be altered by the subroutine if dry convective adjustment
	185	! occurs and if IPBL is not equal to 0.
	186	!
	187	! q: Array of specific humidity (gm/gm) of dimension ND, with first
	188	! index corresponding to lowest model level. Must be defined
	189	! at same grid levels as T. Note that this array will be altered
	190	! if dry convective adjustment occurs and if IPBL is not equal to 0.
	191	!
	192	! qs: Array of saturation specific humidity of dimension ND, with first
	193	! index corresponding to lowest model level. Must be defined
	194	! at same grid levels as T. Note that this array will be altered
	195	! if dry convective adjustment occurs and if IPBL is not equal to 0.
	196	!
	197	! t_wake: Array of absolute temperature (K), seen by unsaturated draughts,
	198	! of dimension ND, with first index corresponding to lowest model level.
	199	!
	200	! q_wake: Array of specific humidity (gm/gm), seen by unsaturated draughts,
	201	! of dimension ND, with first index corresponding to lowest model level.
	202	! Must be defined at same grid levels as T.
	203	!
	204	!qs_wake: Array of saturation specific humidity, seen by unsaturated draughts,
	205	! of dimension ND, with first index corresponding to lowest model level.
	206	! Must be defined at same grid levels as T.
	207	!
[1146]	208	!s_wake: Array of fractionnal area occupied by the wakes.
	209	!
[879]	210	! u: Array of zonal wind velocity (m/s) of dimension ND, witth first
	211	! index corresponding with the lowest model level. Defined at
	212	! same levels as T. Note that this array will be altered if
	213	! dry convective adjustment occurs and if IPBL is not equal to 0.
	214	!
	215	! v: Same as u but for meridional velocity.
	216	!
	217	! tra: Array of passive tracer mixing ratio, of dimensions (ND,NTRA),
	218	! where NTRA is the number of different tracers. If no
	219	! convective tracer transport is needed, define a dummy
	220	! input array of dimension (ND,1). Tracers are defined at
	221	! same vertical levels as T. Note that this array will be altered
	222	! if dry convective adjustment occurs and if IPBL is not equal to 0.
	223	!
	224	! p: Array of pressure (mb) of dimension ND, with first
	225	! index corresponding to lowest model level. Must be defined
	226	! at same grid levels as T.
	227	!
	228	! ph: Array of pressure (mb) of dimension ND+1, with first index
	229	! corresponding to lowest level. These pressures are defined at
	230	! levels intermediate between those of P, T, Q and QS. The first
	231	! value of PH should be greater than (i.e. at a lower level than)
	232	! the first value of the array P.
	233	!
	234	! ALE: Available lifting Energy
	235	!
	236	! ALP: Available lifting Power
	237	!
	238	! nl: The maximum number of levels to which convection can penetrate, plus 1.
	239	! NL MUST be less than or equal to ND-1.
	240	!
	241	! delt: The model time step (sec) between calls to CONVECT
	242	!
	243	!----------------------------------------------------------------------------
	244	! --- On Output:
	245	!
	246	! iflag: An output integer whose value denotes the following:
	247	! VALUE INTERPRETATION
	248	! ----- --------------
	249	! 0 Moist convection occurs.
	250	! 1 Moist convection occurs, but a CFL condition
	251	! on the subsidence warming is violated. This
	252	! does not cause the scheme to terminate.
	253	! 2 Moist convection, but no precip because ep(inb) lt 0.0001
	254	! 3 No moist convection because new cbmf is 0 and old cbmf is 0.
	255	! 4 No moist convection; atmosphere is not
	256	! unstable
	257	! 6 No moist convection because ihmin le minorig.
	258	! 7 No moist convection because unreasonable
	259	! parcel level temperature or specific humidity.
	260	! 8 No moist convection: lifted condensation
	261	! level is above the 200 mb level.
	262	! 9 No moist convection: cloud base is higher
	263	! then the level NL-1.
	264	!
	265	! ft: Array of temperature tendency (K/s) of dimension ND, defined at same
	266	! grid levels as T, Q, QS and P.
	267	!
	268	! fq: Array of specific humidity tendencies ((gm/gm)/s) of dimension ND,
	269	! defined at same grid levels as T, Q, QS and P.
	270	!
	271	! fu: Array of forcing of zonal velocity (m/s^2) of dimension ND,
	272	! defined at same grid levels as T.
	273	!
	274	! fv: Same as FU, but for forcing of meridional velocity.
	275	!
	276	! ftra: Array of forcing of tracer content, in tracer mixing ratio per
	277	! second, defined at same levels as T. Dimensioned (ND,NTRA).
	278	!
	279	! precip: Scalar convective precipitation rate (mm/day).
	280	!
	281	! wd: A convective downdraft velocity scale. For use in surface
	282	! flux parameterizations. See convect.ps file for details.
	283	!
	284	! tprime: A convective downdraft temperature perturbation scale (K).
	285	! For use in surface flux parameterizations. See convect.ps
	286	! file for details.
	287	!
	288	! qprime: A convective downdraft specific humidity
	289	! perturbation scale (gm/gm).
	290	! For use in surface flux parameterizations. See convect.ps
	291	! file for details.
	292	!
	293	! cbmf: The cloud base mass flux ((kg/m**2)/s). THIS SCALAR VALUE MUST
	294	! BE STORED BY THE CALLING PROGRAM AND RETURNED TO CONVECT AT
	295	! ITS NEXT CALL. That is, the value of CBMF must be "remembered"
	296	! by the calling program between calls to CONVECT.
	297	!
	298	! det: Array of detrainment mass flux of dimension ND.
	299	!
	300	! ftd: Array of temperature tendency due to precipitations (K/s) of dimension ND,
	301	! defined at same grid levels as T, Q, QS and P.
	302	!
	303	! fqd: Array of specific humidity tendencies due to precipitations ((gm/gm)/s)
	304	! of dimension ND, defined at same grid levels as T, Q, QS and P.
	305	!
	306	!-------------------------------------------------------------------
	307	c
	308	c Local arrays
	309	c
	310
	311	integer i,k,n,il,j
	312	integer nword1,nword2,nword3,nword4
	313	integer icbmax
	314	integer nk1(klon)
	315	integer icb1(klon)
	316	integer icbs1(klon)
	317
	318	logical ok_inhib ! True => possible inhibition of convection by dryness
[938]	319	logical, save :: debut=.true.
[987]	320	c$OMP THREADPRIVATE(debut)
[879]	321
	322	real plcl1(klon)
	323	real tnk1(klon)
	324	real thnk1(klon)
	325	real qnk1(klon)
	326	real gznk1(klon)
	327	real pnk1(klon)
	328	real qsnk1(klon)
	329	real unk1(klon)
	330	real vnk1(klon)
	331	real cpnk1(klon)
	332	real hnk1(klon)
	333	real pbase1(klon)
	334	real buoybase1(klon)
	335
	336	real lv1(klon,klev) ,lv1_wake(klon,klev)
	337	real cpn1(klon,klev),cpn1_wake(klon,klev)
	338	real tv1(klon,klev) ,tv1_wake(klon,klev)
	339	real gz1(klon,klev) ,gz1_wake(klon,klev)
	340	real hm1(klon,klev) ,hm1_wake(klon,klev)
	341	real h1(klon,klev) ,h1_wake(klon,klev)
	342	real tp1(klon,klev)
	343	real clw1(klon,klev)
	344	real th1(klon,klev) ,th1_wake(klon,klev)
	345	c
	346	real bid(klon,klev) ! dummy array
	347	c
	348	integer ncum
	349	c
	350	integer j1feed(klon)
	351	integer j2feed(klon)
	352	real p1feed1(len) ! pressure at lower bound of feeding layer
	353	real p2feed1(len) ! pressure at upper bound of feeding layer
	354	real wghti1(len,nd) ! weights of the feeding layers
	355	c
	356	c (local) compressed fields:
	357	c
	358	integer nloc
[940]	359	c parameter (nloc=klon) ! pour l'instant
[879]	360
	361	integer idcum(nloc)
	362	integer iflag(nloc),nk(nloc),icb(nloc)
	363	integer nent(nloc,klev)
	364	integer icbs(nloc)
	365	integer inb(nloc), inbis(nloc)
	366
	367	real cbmf(nloc),plcl(nloc)
	368	real t(nloc,klev),q(nloc,klev),qs(nloc,klev)
	369	real t_wake(nloc,klev),q_wake(nloc,klev),qs_wake(nloc,klev)
[1146]	370	real s_wake(nloc)
[879]	371	real u(nloc,klev),v(nloc,klev)
	372	real gz(nloc,klev),h(nloc,klev) ,hm(nloc,klev)
	373	real h_wake(nloc,klev),hm_wake(nloc,klev)
	374	real lv(nloc,klev) ,cpn(nloc,klev)
	375	real lv_wake(nloc,klev),cpn_wake(nloc,klev)
	376	real p(nloc,klev),ph(nloc,klev+1),tv(nloc,klev) ,tp(nloc,klev)
	377	real tv_wake(nloc,klev)
	378	real clw(nloc,klev)
	379	real dph(nloc,klev)
	380	real pbase(nloc), buoybase(nloc), th(nloc,klev)
	381	real th_wake(nloc,klev)
	382	real tvp(nloc,klev)
	383	real sig(nloc,klev), w0(nloc,klev), ptop2(nloc)
	384	real hp(nloc,klev), ep(nloc,klev), sigp(nloc,klev)
	385	real frac(nloc), buoy(nloc,klev)
	386	real cape(nloc)
	387	real cin(nloc)
	388	real m(nloc,klev)
	389	real ment(nloc,klev,klev), sij(nloc,klev,klev)
	390	real qent(nloc,klev,klev)
	391	real hent(nloc,klev,klev)
	392	real uent(nloc,klev,klev), vent(nloc,klev,klev)
	393	real ments(nloc,klev,klev), qents(nloc,klev,klev)
	394	real elij(nloc,klev,klev)
	395	real supmax(nloc,klev)
	396	real ale(nloc),alp(nloc),coef_clos(nloc)
[940]	397	real sigd(nloc)
	398	! real mp(nloc,klev), qp(nloc,klev), up(nloc,klev), vp(nloc,klev)
	399	! real wt(nloc,klev), water(nloc,klev), evap(nloc,klev)
	400	! real b(nloc,klev), sigd(nloc)
	401	! save mp,qp,up,vp,wt,water,evap,b
	402	real, save, allocatable :: mp(:,:),qp(:,:),up(:,:),vp(:,:)
	403	real, save, allocatable :: wt(:,:),water(:,:),evap(:,:), b(:,:)
	404	c$OMP THREADPRIVATE(mp,qp,up,vp,wt,water,evap,b)
[879]	405	real ft(nloc,klev), fq(nloc,klev)
	406	real ftd(nloc,klev), fqd(nloc,klev)
	407	real fu(nloc,klev), fv(nloc,klev)
	408	real upwd(nloc,klev), dnwd(nloc,klev), dnwd0(nloc,klev)
	409	real Ma(nloc,klev), mip(nloc,klev), tls(nloc,klev)
	410	real tps(nloc,klev), qprime(nloc), tprime(nloc)
	411	real precip(nloc)
[1334]	412	! real Vprecip(nloc,klev)
	413	real Vprecip(nloc,klev+1)
[879]	414	real tra(nloc,klev,ntra), trap(nloc,klev,ntra)
	415	real ftra(nloc,klev,ntra), traent(nloc,klev,klev,ntra)
	416	real qcondc(nloc,klev) ! cld
	417	real wd(nloc) ! gust
	418	real Plim1(nloc),Plim2(nloc)
	419	real asupmax(nloc,klev)
	420	real supmax0(nloc)
	421	real asupmaxmin(nloc)
	422	c
	423	real tnk(nloc),qnk(nloc),gznk(nloc)
	424	real wghti(nloc,nd)
	425	real hnk(nloc),unk(nloc),vnk(nloc)
[940]	426	logical, save :: first=.true.
[987]	427	c$OMP THREADPRIVATE(first)
[1403]	428	CHARACTER (LEN=20) :: modname='cva_driver'
	429	CHARACTER (LEN=80) :: abort_message
[879]	430
	431	c
	432	! print *, 't1, t1_wake ',(k,t1(1,k),t1_wake(1,k),k=1,klev)
	433	! print *, 'q1, q1_wake ',(k,q1(1,k),q1_wake(1,k),k=1,klev)
	434
	435	!-------------------------------------------------------------------
	436	! --- SET CONSTANTS AND PARAMETERS
	437	!-------------------------------------------------------------------
	438
[940]	439	if (first) then
	440	allocate(mp(nloc,klev), qp(nloc,klev), up(nloc,klev))
	441	allocate(vp(nloc,klev), wt(nloc,klev), water(nloc,klev))
	442	allocate(evap(nloc,klev), b(nloc,klev))
	443	first=.false.
	444	endif
[879]	445	c -- set simulation flags:
	446	c (common cvflag)
	447
	448	CALL cv_flag
	449
	450	c -- set thermodynamical constants:
	451	c (common cvthermo)
	452
	453	CALL cv_thermo(iflag_con)
	454
	455	c -- set convect parameters
	456	c
	457	c includes microphysical parameters and parameters that
	458	c control the rate of approach to quasi-equilibrium)
	459	c (common cvparam)
	460
	461	if (iflag_con.eq.3) then
	462	CALL cv3_param(nd,delt)
	463
	464	endif
	465
	466	if (iflag_con.eq.4) then
	467	CALL cv_param(nd)
	468	endif
	469
	470	!---------------------------------------------------------------------
	471	! --- INITIALIZE OUTPUT ARRAYS AND PARAMETERS
	472	!---------------------------------------------------------------------
	473	nword1=len
	474	nword2=len*nd
	475	nword3=lenndntra
	476	nword4=lenndnd
	477
	478	! call izilch(iflag1 ,nword1)
	479	! call zilch(iflag1 ,nword1)
	480	do i=1,len
	481	iflag1(i)=0
	482	ktop1(i)=0
	483	kbas1(i)=0
	484	enddo
	485	call zilch(ft1 ,nword2)
	486	call zilch(fq1 ,nword2)
	487	call zilch(fu1 ,nword2)
	488	call zilch(fv1 ,nword2)
	489	call zilch(ftra1 ,nword3)
	490	call zilch(precip1 ,nword1)
	491	! call izilch(kbas1 ,nword1)
	492	! call zilch(kbas1 ,nword1)
	493	! call izilch(ktop1 ,nword1)
	494	! call zilch(ktop1 ,nword1)
	495	call zilch(cbmf1 ,nword1)
	496	call zilch(ptop21 ,nword1)
	497	call zilch(Ma1 ,nword2)
	498	call zilch(mip1 ,nword2)
[1334]	499	! call zilch(Vprecip1,nword2)
	500	Vprecip1=0.
[879]	501	call zilch(upwd1 ,nword2)
	502	call zilch(dnwd1 ,nword2)
	503	call zilch(dnwd01 ,nword2)
	504	call zilch(qcondc1 ,nword2)
	505	!test
	506	! call zilch(qcondc ,nword2)
	507	call zilch(wd1 ,nword1)
	508	call zilch(cape1 ,nword1)
	509	call zilch(cin1 ,nword1)
	510	call zilch(tvp1 ,nword2)
	511	call zilch(ftd1 ,nword2)
	512	call zilch(fqd1 ,nword2)
	513	call zilch(Plim11 ,nword1)
	514	call zilch(Plim21 ,nword1)
	515	call zilch(asupmax1,nword2)
	516	call zilch(supmax01,nword1)
	517	call zilch(asupmaxmin1,nword1)
	518	c
[973]	519	DO il = 1,len
	520	cin1(il) = -100000.
	521	cape1(il) = -1.
	522	ENDDO
	523	c
[879]	524	if (iflag_con.eq.3) then
	525	do il=1,len
	526	sig1(il,nd)=sig1(il,nd)+1.
	527	sig1(il,nd)=amin1(sig1(il,nd),12.1)
	528	enddo
	529	endif
	530
	531	!---------------------------------------------------------------------
	532	! --- INITIALIZE LOCAL ARRAYS AND PARAMETERS
	533	!---------------------------------------------------------------------
	534	!
	535	do il = 1,nloc
	536	coef_clos(il)=1.
	537	enddo
	538
	539	!--------------------------------------------------------------------
	540	! --- CALCULATE ARRAYS OF GEOPOTENTIAL, HEAT CAPACITY & STATIC ENERGY
	541	!--------------------------------------------------------------------
	542
	543	if (iflag_con.eq.3) then
[938]	544
	545	if (debut) THEN
	546	print*,'Emanuel version 3 nouvelle'
	547	endif
[1146]	548	! print*,'t1, q1 ',t1,q1
[879]	549	CALL cv3_prelim(len,nd,ndp1,t1,q1,p1,ph1 ! nd->na
	550	o ,lv1,cpn1,tv1,gz1,h1,hm1,th1)
	551
	552	c
	553	CALL cv3_prelim(len,nd,ndp1,t1_wake,q1_wake,p1,ph1 ! nd->na
	554	o ,lv1_wake,cpn1_wake,tv1_wake,gz1_wake
	555	o ,h1_wake,bid,th1_wake)
	556
	557	endif
	558	c
	559	if (iflag_con.eq.4) then
	560	print*,'Emanuel version 4 '
	561	CALL cv_prelim(len,nd,ndp1,t1,q1,p1,ph1
	562	o ,lv1,cpn1,tv1,gz1,h1,hm1)
	563	endif
	564
	565	!--------------------------------------------------------------------
	566	! --- CONVECTIVE FEED
	567	!--------------------------------------------------------------------
	568	!
	569	! compute feeding layer potential temperature and mixing ratio :
	570	!
	571	! get bounds of feeding layer
	572	!
	573	c test niveaux couche alimentation KE
	574	if(sig1feed1.eq.sig2feed1) then
[1403]	575	write(lunout,*)'impossible de choisir sig1feed=sig2feed'
	576	write(lunout,*)'changer la valeur de sig2feed dans physiq.def'
	577	abort_message = ''
	578	CALL abort_gcm (modname,abort_message,1)
[879]	579	endif
	580	c
	581	do i=1,len
	582	p1feed1(i)=sig1feed1*ph1(i,1)
	583	p2feed1(i)=sig2feed1*ph1(i,1)
	584	ctest maf
	585	c p1feed1(i)=ph1(i,1)
	586	c p2feed1(i)=ph1(i,2)
	587	c p2feed1(i)=ph1(i,3)
	588	ctestCR: on prend la couche alim des thermiques
	589	c p2feed1(i)=ph1(i,lalim_conv(i)+1)
	590	c print*,'lentr=',lentr(i),ph1(i,lentr(i)+1),ph1(i,2)
	591	end do
	592	!
	593	if (iflag_con.eq.3) then
	594	endif
	595	do i=1,len
	596	! print*,'avant cv3_feed plim',p1feed1(i),p2feed1(i)
	597	enddo
	598	if (iflag_con.eq.3) then
	599
	600	c print*, 'IFLAG1 avant cv3_feed'
	601	c print*,'len,nd',len,nd
	602	c write(*,'(64i1)') iflag1(2:klon-1)
	603
	604	CALL cv3_feed(len,nd,t1,q1,u1,v1,p1,ph1,hm1,gz1 ! nd->na
	605	i ,p1feed1,p2feed1,wght1
	606	o ,wghti1,tnk1,thnk1,qnk1,qsnk1,unk1,vnk1
	607	o ,cpnk1,hnk1,nk1,icb1,icbmax,iflag1,gznk1,plcl1)
	608	endif
	609
	610	c print*, 'IFLAG1 apres cv3_feed'
	611	c print*,'len,nd',len,nd
	612	c write(*,'(64i1)') iflag1(2:klon-1)
	613
	614	if (iflag_con.eq.4) then
	615	CALL cv_feed(len,nd,t1,q1,qs1,p1,hm1,gz1
	616	o ,nk1,icb1,icbmax,iflag1,tnk1,qnk1,gznk1,plcl1)
	617	endif
	618	c
	619	! print *, 'cv3_feed-> iflag1, plcl1 ',iflag1(1),plcl1(1)
	620	c
	621	!--------------------------------------------------------------------
	622	! --- UNDILUTE (ADIABATIC) UPDRAFT / 1st part
	623	! (up through ICB for convect4, up through ICB+1 for convect3)
	624	! Calculates the lifted parcel virtual temperature at nk, the
	625	! actual temperature, and the adiabatic liquid water content.
	626	!--------------------------------------------------------------------
	627
	628	if (iflag_con.eq.3) then
	629
	630	CALL cv3_undilute1(len,nd,t1,qs1,gz1,plcl1,p1,icb1,tnk1,qnk1 ! nd->na
	631	o ,gznk1,tp1,tvp1,clw1,icbs1)
	632	endif
	633
	634
	635	if (iflag_con.eq.4) then
	636	CALL cv_undilute1(len,nd,t1,q1,qs1,gz1,p1,nk1,icb1,icbmax
	637	: ,tp1,tvp1,clw1)
	638	endif
	639	c
	640	!-------------------------------------------------------------------
	641	! --- TRIGGERING
	642	!-------------------------------------------------------------------
	643	c
	644	! print *,' avant triggering, iflag_con ',iflag_con
	645	c
	646	if (iflag_con.eq.3) then
	647
	648	CALL cv3_trigger(len,nd,icb1,plcl1,p1,th1,tv1,tvp1,thnk1 ! nd->na
	649	o ,pbase1,buoybase1,iflag1,sig1,w01)
	650
	651
	652	c print*, 'IFLAG1 apres cv3_triger'
	653	c print*,'len,nd',len,nd
	654	c write(*,'(64i1)') iflag1(2:klon-1)
	655
	656	c call dump2d(iim,jjm-1,sig1(2)
	657	endif
	658
	659	if (iflag_con.eq.4) then
	660	CALL cv_trigger(len,nd,icb1,cbmf1,tv1,tvp1,iflag1)
	661	endif
	662	c
	663	c
	664	!=====================================================================
	665	! --- IF THIS POINT IS REACHED, MOIST CONVECTIVE ADJUSTMENT IS NECESSARY
	666	!=====================================================================
	667
	668	ncum=0
	669	do 400 i=1,len
	670	if(iflag1(i).eq.0)then
	671	ncum=ncum+1
	672	idcum(ncum)=i
	673	endif
	674	400 continue
	675	c
	676	! print*,'klon, ncum = ',len,ncum
	677	c
	678	IF (ncum.gt.0) THEN
	679
	680	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	681	! --- COMPRESS THE FIELDS
	682	! (-> vectorization over convective gridpoints)
	683	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	684
	685	if (iflag_con.eq.3) then
[1146]	686	! print*,'ncum tv1 ',ncum,tv1
	687	! print*,'tvp1 ',tvp1
[879]	688	CALL cv3a_compress( len,nloc,ncum,nd,ntra
	689	: ,iflag1,nk1,icb1,icbs1
	690	: ,plcl1,tnk1,qnk1,gznk1,hnk1,unk1,vnk1
	691	: ,wghti1,pbase1,buoybase1
[1146]	692	: ,t1,q1,qs1,t1_wake,q1_wake,qs1_wake,s1_wake
	693	: ,u1,v1,gz1,th1,th1_wake
[879]	694	: ,tra1
	695	: ,h1 ,lv1 ,cpn1 ,p1,ph1,tv1 ,tp1,tvp1,clw1
	696	: ,h1_wake,lv1_wake,cpn1_wake ,tv1_wake
	697	: ,sig1,w01,ptop21
	698	: ,Ale1,Alp1
	699	o ,iflag,nk,icb,icbs
	700	o ,plcl,tnk,qnk,gznk,hnk,unk,vnk
	701	o ,wghti,pbase,buoybase
[1146]	702	o ,t,q,qs,t_wake,q_wake,qs_wake,s_wake
	703	o ,u,v,gz,th,th_wake
[879]	704	o ,tra
	705	o ,h ,lv ,cpn ,p,ph,tv ,tp,tvp,clw
	706	o ,h_wake,lv_wake,cpn_wake ,tv_wake
	707	o ,sig,w0,ptop2
	708	o ,Ale,Alp )
	709
[1146]	710	! print*,'tv ',tv
	711	! print*,'tvp ',tvp
	712
[879]	713	endif
	714
	715	if (iflag_con.eq.4) then
	716	CALL cv_compress( len,nloc,ncum,nd
	717	: ,iflag1,nk1,icb1
	718	: ,cbmf1,plcl1,tnk1,qnk1,gznk1
	719	: ,t1,q1,qs1,u1,v1,gz1
	720	: ,h1,lv1,cpn1,p1,ph1,tv1,tp1,tvp1,clw1
	721	o ,iflag,nk,icb
	722	o ,cbmf,plcl,tnk,qnk,gznk
	723	o ,t,q,qs,u,v,gz,h,lv,cpn,p,ph,tv,tp,tvp,clw
	724	o ,dph )
	725	endif
	726
	727	!-------------------------------------------------------------------
	728	! --- UNDILUTE (ADIABATIC) UPDRAFT / second part :
	729	! --- FIND THE REST OF THE LIFTED PARCEL TEMPERATURES
	730	! --- &
	731	! --- COMPUTE THE PRECIPITATION EFFICIENCIES AND THE
	732	! --- FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
	733	! --- &
	734	! --- FIND THE LEVEL OF NEUTRAL BUOYANCY
	735	!-------------------------------------------------------------------
	736
	737	if (iflag_con.eq.3) then
	738	CALL cv3_undilute2(nloc,ncum,nd,icb,icbs,nk !na->nd
	739	: ,tnk,qnk,gznk,hnk,t,q,qs,gz
	740	: ,p,h,tv,lv,pbase,buoybase,plcl
	741	o ,inb,tp,tvp,clw,hp,ep,sigp,buoy)
	742
	743	endif
	744
	745	if (iflag_con.eq.4) then
	746	CALL cv_undilute2(nloc,ncum,nd,icb,nk
	747	: ,tnk,qnk,gznk,t,q,qs,gz
	748	: ,p,dph,h,tv,lv
	749	o ,inb,inbis,tp,tvp,clw,hp,ep,sigp,frac)
	750	endif
	751	c
	752	!-------------------------------------------------------------------
	753	! --- MIXING(1) (if iflag_mix .ge. 1)
	754	!-------------------------------------------------------------------
	755	IF (iflag_con .eq. 3) THEN
	756	IF (iflag_mix .ge. 1 ) THEN
[1062]	757	CALL zilch(supmax,nloc*klev)
[879]	758	CALL cv3p_mixing(nloc,ncum,nd,nd,ntra,icb,nk,inb ! na->nd
	759	: ,ph,t,q,qs,u,v,tra,h,lv,qnk
	760	: ,unk,vnk,hp,tv,tvp,ep,clw,sig
[991]	761	: ,ment,qent,hent,uent,vent,nent
[879]	762	: ,sij,elij,supmax,ments,qents,traent)
	763	! print*, 'cv3p_mixing-> supmax ', (supmax(1,k), k=1,nd)
	764
	765	ELSE
	766	CALL zilch(supmax,nloc*klev)
	767	ENDIF
	768	ENDIF
	769	!-------------------------------------------------------------------
	770	! --- CLOSURE
	771	!-------------------------------------------------------------------
	772
	773	c
	774	if (iflag_con.eq.3) then
	775	IF (iflag_clos .eq. 0) THEN
	776	CALL cv3_closure(nloc,ncum,nd,icb,inb ! na->nd
	777	: ,pbase,p,ph,tv,buoy
	778	o ,sig,w0,cape,m,iflag)
	779	ENDIF
	780	c
	781	ok_inhib = iflag_mix .EQ. 2
	782	c
	783	IF (iflag_clos .eq. 1) THEN
	784	print *,' pas d appel cv3p_closure'
	785	cc CALL cv3p_closure(nloc,ncum,nd,icb,inb ! na->nd
	786	cc : ,pbase,plcl,p,ph,tv,tvp,buoy
	787	cc : ,supmax
	788	cc o ,sig,w0,ptop2,cape,cin,m)
	789	ENDIF
	790	IF (iflag_clos .eq. 2) THEN
	791	CALL cv3p1_closure(nloc,ncum,nd,icb,inb ! na->nd
	792	: ,pbase,plcl,p,ph,tv,tvp,buoy
	793	: ,supmax,ok_inhib,Ale,Alp
	794	o ,sig,w0,ptop2,cape,cin,m,iflag,coef_clos
	795	: ,Plim1,Plim2,asupmax,supmax0
[1279]	796	: ,asupmaxmin,cbmf)
[879]	797	ENDIF
	798	endif ! iflag_con.eq.3
	799
	800	if (iflag_con.eq.4) then
	801	CALL cv_closure(nloc,ncum,nd,nk,icb
	802	: ,tv,tvp,p,ph,dph,plcl,cpn
	803	o ,iflag,cbmf)
	804	endif
	805	c
	806	! print *,'cv_closure-> cape ',cape(1)
	807	c
	808	!-------------------------------------------------------------------
	809	! --- MIXING(2)
	810	!-------------------------------------------------------------------
	811
	812	if (iflag_con.eq.3) then
	813	IF (iflag_mix.eq.0) THEN
	814	CALL cv3_mixing(nloc,ncum,nd,nd,ntra,icb,nk,inb ! na->nd
	815	: ,ph,t,q,qs,u,v,tra,h,lv,qnk
	816	: ,unk,vnk,hp,tv,tvp,ep,clw,m,sig
[991]	817	o ,ment,qent,uent,vent,nent,sij,elij,ments,qents,traent)
[879]	818	CALL zilch(hent,nlocklevklev)
	819	ELSE
	820	CALL cv3_mixscale(nloc,ncum,nd,ment,m)
[938]	821	if (debut) THEN
	822	print *,' cv3_mixscale-> '
	823	endif !(debut) THEN
[879]	824	ENDIF
	825	endif
	826
	827	if (iflag_con.eq.4) then
	828	CALL cv_mixing(nloc,ncum,nd,icb,nk,inb,inbis
	829	: ,ph,t,q,qs,u,v,h,lv,qnk
	830	: ,hp,tv,tvp,ep,clw,cbmf
	831	o ,m,ment,qent,uent,vent,nent,sij,elij)
	832	endif
	833	c
[938]	834	if (debut) THEN
	835	print *,' cv_mixing ->'
	836	endif !(debut) THEN
[879]	837	c do i = 1,klev
	838	c print*,'cv_mixing-> i,ment ',i,(ment(1,i,j),j=1,klev)
	839	c enddo
	840	c
	841	!-------------------------------------------------------------------
	842	! --- UNSATURATED (PRECIPITATING) DOWNDRAFTS
	843	!-------------------------------------------------------------------
	844	if (iflag_con.eq.3) then
[938]	845	if (debut) THEN
	846	print *,' cva_driver -> cv3_unsat '
	847	endif !(debut) THEN
[879]	848
	849	CALL cv3_unsat(nloc,ncum,nd,nd,ntra,icb,inb,iflag ! na->nd
	850	: ,t_wake,q_wake,qs_wake,gz,u,v,tra,p,ph
	851	: ,th_wake,tv_wake,lv_wake,cpn_wake
	852	: ,ep,sigp,clw
	853	: ,m,ment,elij,delt,plcl,coef_clos
	854	o ,mp,qp,up,vp,trap,wt,water,evap,b,sigd)
	855	endif
	856
	857	if (iflag_con.eq.4) then
	858	CALL cv_unsat(nloc,ncum,nd,inb,t,q,qs,gz,u,v,p,ph
	859	: ,h,lv,ep,sigp,clw,m,ment,elij
	860	o ,iflag,mp,qp,up,vp,wt,water,evap)
	861	endif
	862	c
[938]	863	if (debut) THEN
	864	print *,'cv_unsat-> '
	865	debut=.FALSE.
	866	endif !(debut) THEN
[879]	867	!
	868	c print *,'cv_unsat-> mp ',mp
	869	c print *,'cv_unsat-> water ',water
	870	!-------------------------------------------------------------------
	871	! --- YIELD
	872	! (tendencies, precipitation, variables of interface with other
	873	! processes, etc)
	874	!-------------------------------------------------------------------
	875
	876	if (iflag_con.eq.3) then
	877
	878	CALL cv3_yield(nloc,ncum,nd,nd,ntra ! na->nd
	879	: ,icb,inb,delt
[1146]	880	: ,t,q,t_wake,q_wake,s_wake,u,v,tra
[1044]	881	: ,gz,p,ph,h,hp,lv,cpn,th,th_wake
[879]	882	: ,ep,clw,m,tp,mp,qp,up,vp,trap
	883	: ,wt,water,evap,b,sigd
	884	: ,ment,qent,hent,iflag_mix,uent,vent
	885	: ,nent,elij,traent,sig
	886	: ,tv,tvp,wghti
	887	: ,iflag,precip,Vprecip,ft,fq,fu,fv,ftra
	888	: ,cbmf,upwd,dnwd,dnwd0,ma,mip
	889	: ,tls,tps,qcondc,wd
	890	: ,ftd,fqd)
	891	! print *,' cv3_yield -> fqd(1) = ',fqd(1,1)
	892	endif
	893
	894	if (iflag_con.eq.4) then
	895	CALL cv_yield(nloc,ncum,nd,nk,icb,inb,delt
	896	: ,t,q,t_wake,q_wake,u,v,tra
	897	: ,gz,p,ph,h,hp,lv,cpn,th
	898	: ,ep,clw,frac,m,mp,qp,up,vp
	899	: ,wt,water,evap
	900	: ,ment,qent,uent,vent,nent,elij
	901	: ,tv,tvp
	902	o ,iflag,wd,qprime,tprime
	903	o ,precip,cbmf,ft,fq,fu,fv,Ma,qcondc)
	904	endif
	905
	906	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	907	! --- UNCOMPRESS THE FIELDS
	908	!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
	909
	910
	911	if (iflag_con.eq.3) then
	912	CALL cv3a_uncompress(nloc,len,ncum,nd,ntra,idcum
	913	: ,iflag,icb,inb
[1398]	914	: ,precip,cbmf,sig,w0,ptop2
[879]	915	: ,ft,fq,fu,fv,ftra
	916	: ,Ma,mip,Vprecip,upwd,dnwd,dnwd0
	917	; ,qcondc,wd,cape,cin
	918	: ,tvp
	919	: ,ftd,fqd
	920	: ,Plim1,Plim2,asupmax,supmax0
	921	: ,asupmaxmin
	922	o ,iflag1,kbas1,ktop1
[1398]	923	o ,precip1,cbmf1,sig1,w01,ptop21
[879]	924	o ,ft1,fq1,fu1,fv1,ftra1
	925	o ,Ma1,mip1,Vprecip1,upwd1,dnwd1,dnwd01
	926	o ,qcondc1,wd1,cape1,cin1
	927	o ,tvp1
	928	o ,ftd1,fqd1
	929	o ,Plim11,Plim21,asupmax1,supmax01
	930	o ,asupmaxmin1)
	931	endif
	932
	933	if (iflag_con.eq.4) then
	934	CALL cv_uncompress(nloc,len,ncum,nd,idcum
	935	: ,iflag
	936	: ,precip,cbmf
	937	: ,ft,fq,fu,fv
	938	: ,Ma,qcondc
	939	o ,iflag1
	940	o ,precip1,cbmf1
	941	o ,ft1,fq1,fu1,fv1
	942	o ,Ma1,qcondc1 )
	943	endif
	944
	945	ENDIF ! ncum>0
	946
	947	9999 continue
	948	return
	949	end
	950

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

Download in other formats:

Original Format