Context Navigation

newmicro.F90 @ 5407

Last change on this file since 5407 was 2641, checked in by Laurent Fairhead, 8 years ago
Merged trunk changes r2593:2640 into testing branch
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 20.5 KB

Rev	Line
[1279]	1	! $Id: newmicro.F90 2641 2016-09-29 21:26:46Z abarral $
[1523]	2
	3
[1712]	4
[1992]	5	SUBROUTINE newmicro(ok_cdnc, bl95_b0, bl95_b1, paprs, pplay, t, pqlwp, pclc, &
	6	pcltau, pclemi, pch, pcl, pcm, pct, pctlwp, xflwp, xfiwp, xflwc, xfiwc, &
	7	mass_solu_aero, mass_solu_aero_pi, pcldtaupi, re, fl, reliq, reice, &
	8	reliq_pi, reice_pi)
[524]	9
[1992]	10	USE dimphy
	11	USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
	12	reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra
	13	USE phys_state_var_mod, ONLY: rnebcon, clwcon
[2160]	14	USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
[1992]	15	IMPLICIT NONE
	16	! ======================================================================
	17	! Auteur(s): Z.X. Li (LMD/CNRS) date: 19930910
	18	! O. Boucher (LMD/CNRS) mise a jour en 201212
[2641]	19	! I. Musat (LMD/CNRS) : prise en compte de la meme hypothese de recouvrement
	20	! pour les nuages que pour le rayonnement rrtm via
	21	! le parametre novlp de radopt.h : 20160721
[1992]	22	! Objet: Calculer epaisseur optique et emmissivite des nuages
	23	! ======================================================================
	24	! Arguments:
	25	! ok_cdnc-input-L-flag pour calculer les rayons a partir des aerosols
[1146]	26
[1992]	27	! t-------input-R-temperature
	28	! pqlwp---input-R-eau liquide nuageuse dans l'atmosphere dans la partie
	29	! nuageuse (kg/kg)
	30	! pclc----input-R-couverture nuageuse pour le rayonnement (0 a 1)
	31	! mass_solu_aero-----input-R-total mass concentration for all soluble
	32	! aerosols[ug/m^3]
	33	! mass_solu_aero_pi--input-R-ditto, pre-industrial value
[1279]	34
[1992]	35	! bl95_b0-input-R-a PARAMETER, may be varied for tests (s-sea, l-land)
	36	! bl95_b1-input-R-a PARAMETER, may be varied for tests ( -"- )
	37
	38	! re------output-R-Cloud droplet effective radius multiplied by fl [um]
	39	! fl------output-R-Denominator to re, introduced to avoid problems in
	40	! the averaging of the output. fl is the fraction of liquid
	41	! water clouds within a grid cell
	42
	43	! pcltau--output-R-epaisseur optique des nuages
	44	! pclemi--output-R-emissivite des nuages (0 a 1)
	45	! pcldtaupi-output-R-pre-industrial value of cloud optical thickness,
	46
	47	! pcl-output-R-2D low-level cloud cover
	48	! pcm-output-R-2D mid-level cloud cover
	49	! pch-output-R-2D high-level cloud cover
	50	! pct-output-R-2D total cloud cover
	51	! ======================================================================
	52
	53	include "YOMCST.h"
	54	include "nuage.h"
	55	include "radepsi.h"
	56	include "radopt.h"
	57
[2641]	58	! choix de l'hypothese de recouvrement nuageuse via radopt.h (IM, 19.07.2016)
	59	! !novlp=1: max-random
	60	! !novlp=2: maximum
	61	! !novlp=3: random
	62	! LOGICAL random, maximum_random, maximum
	63	! PARAMETER (random=.FALSE., maximum_random=.TRUE., maximum=.FALSE.)
[1992]	64
	65	LOGICAL, SAVE :: first = .TRUE.
	66	!$OMP THREADPRIVATE(FIRST)
	67	INTEGER flag_max
	68
	69	! threshold PARAMETERs
	70	REAL thres_tau, thres_neb
	71	PARAMETER (thres_tau=0.3, thres_neb=0.001)
	72
	73	REAL phase3d(klon, klev)
	74	REAL tcc(klon), ftmp(klon), lcc_integrat(klon), height(klon)
	75
	76	REAL paprs(klon, klev+1)
	77	REAL pplay(klon, klev)
	78	REAL t(klon, klev)
	79	REAL pclc(klon, klev)
	80	REAL pqlwp(klon, klev)
	81	REAL pcltau(klon, klev)
	82	REAL pclemi(klon, klev)
	83	REAL pcldtaupi(klon, klev)
	84
	85	REAL pct(klon)
	86	REAL pcl(klon)
	87	REAL pcm(klon)
	88	REAL pch(klon)
	89	REAL pctlwp(klon)
	90
	91	LOGICAL lo
	92
	93	! !Abderr modif JL mail du 19.01.2011 18:31
	94	! REAL cetahb, cetamb
	95	! PARAMETER (cetahb = 0.45, cetamb = 0.80)
	96	! Remplacer
	97	! cetahb*paprs(i,1) par prmhc
	98	! cetamb*paprs(i,1) par prlmc
	99	REAL prmhc ! Pressure between medium and high level cloud in Pa
	100	REAL prlmc ! Pressure between low and medium level cloud in Pa
	101	PARAMETER (prmhc=440.100., prlmc=680.100.)
	102
	103	INTEGER i, k
	104	REAL xflwp(klon), xfiwp(klon)
	105	REAL xflwc(klon, klev), xfiwc(klon, klev)
	106
	107	REAL radius
	108
	109	REAL coef_froi, coef_chau
	110	PARAMETER (coef_chau=0.13, coef_froi=0.09)
	111
	112	REAL seuil_neb
	113	PARAMETER (seuil_neb=0.001)
	114
[2056]	115	! JBM (3/14) nexpo is replaced by exposant_glace
	116	! INTEGER nexpo ! exponentiel pour glace/eau
	117	! PARAMETER (nexpo=6)
	118	! PARAMETER (nexpo=1)
	119	! if iflag_t_glace=0, the old values are used:
	120	REAL, PARAMETER :: t_glace_min_old = 258.
	121	REAL, PARAMETER :: t_glace_max_old = 273.13
[1992]	122
	123	REAL rel, tc, rei
	124	REAL k_ice0, k_ice, df
	125	PARAMETER (k_ice0=0.005) ! units=m2/g
	126	PARAMETER (df=1.66) ! diffusivity factor
	127
	128	! jq for the aerosol indirect effect
	129	! jq introduced by Johannes Quaas (quaas@lmd.jussieu.fr), 27/11/2003
	130	! jq
	131	REAL mass_solu_aero(klon, klev) ! total mass concentration for all soluble aerosols [ug m-3]
	132	REAL mass_solu_aero_pi(klon, klev) ! - " - (pre-industrial value)
	133	REAL cdnc(klon, klev) ! cloud droplet number concentration [m-3]
	134	REAL re(klon, klev) ! cloud droplet effective radius [um]
	135	REAL cdnc_pi(klon, klev) ! cloud droplet number concentration [m-3] (pi value)
	136	REAL re_pi(klon, klev) ! cloud droplet effective radius [um] (pi value)
	137
	138	REAL fl(klon, klev) ! xliq * rneb (denominator to re; fraction of liquid water clouds
	139	! within the grid cell)
	140
	141	LOGICAL ok_cdnc
	142	REAL bl95_b0, bl95_b1 ! Parameter in B&L 95-Formula
	143
	144	! jq-end
	145	! IM cf. CR:parametres supplementaires
	146	REAL zclear(klon)
	147	REAL zcloud(klon)
	148	REAL zcloudh(klon)
	149	REAL zcloudm(klon)
	150	REAL zcloudl(klon)
	151	REAL rhodz(klon, klev) !--rho*dz pour la couche
	152	REAL zrho(klon, klev) !--rho pour la couche
	153	REAL dh(klon, klev) !--dz pour la couche
	154	REAL zfice(klon, klev)
	155	REAL rad_chaud(klon, klev) !--rayon pour les nuages chauds
	156	REAL rad_chaud_pi(klon, klev) !--rayon pour les nuages chauds pre-industriels
	157	REAL zflwp_var, zfiwp_var
	158	REAL d_rei_dt
	159
	160	! Abderrahmane oct 2009
	161	REAL reliq(klon, klev), reice(klon, klev)
	162	REAL reliq_pi(klon, klev), reice_pi(klon, klev)
	163
	164	! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	165	! FH : 2011/05/24
	166
	167	! rei = ( rei_max - rei_min ) * T(°C) / 81.4 + rei_max
	168	! to be used for a temperature in celcius T(°C) < 0
	169	! rei=rei_min for T(°C) < -81.4
	170
	171	! Calcul de la pente de la relation entre rayon effective des cristaux
	172	! et la température.
	173	! Pour retrouver les résultats numériques de la version d'origine,
	174	! on impose 0.71 quand on est proche de 0.71
	175
	176	d_rei_dt = (rei_max-rei_min)/81.4
	177	IF (abs(d_rei_dt-0.71)<1.E-4) d_rei_dt = 0.71
	178	! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	179
	180	! Calculer l'epaisseur optique et l'emmissivite des nuages
	181	! IM inversion des DO
	182
	183	xflwp = 0.D0
	184	xfiwp = 0.D0
	185	xflwc = 0.D0
	186	xfiwc = 0.D0
	187
	188	reliq = 0.
	189	reice = 0.
	190	reliq_pi = 0.
	191	reice_pi = 0.
	192
[2056]	193	IF (iflag_t_glace.EQ.0) THEN
	194	DO k = 1, klev
	195	DO i = 1, klon
	196	! -layer calculation
	197	rhodz(i, k) = (paprs(i,k)-paprs(i,k+1))/rg ! kg/m2
	198	zrho(i, k) = pplay(i, k)/t(i, k)/rd ! kg/m3
	199	dh(i, k) = rhodz(i, k)/zrho(i, k) ! m
	200	! -Fraction of ice in cloud using a linear transition
	201	zfice(i, k) = 1.0 - (t(i,k)-t_glace_min_old)/(t_glace_max_old-t_glace_min_old)
	202	zfice(i, k) = min(max(zfice(i,k),0.0), 1.0)
	203	! -IM Total Liquid/Ice water content
	204	xflwc(i, k) = (1.-zfice(i,k))*pqlwp(i, k)
	205	xfiwc(i, k) = zfice(i, k)*pqlwp(i, k)
	206	END DO
[1992]	207	END DO
[2056]	208	ELSE ! of IF (iflag_t_glace.EQ.0)
	209	DO k = 1, klev
[2160]	210	CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k))
	211
	212
	213	! JBM: icefrac_lsc is now contained icefrac_lsc_mod
	214	! zfice(i, k) = icefrac_lsc(t(i,k), t_glace_min, &
	215	! t_glace_max, exposant_glace)
[2056]	216	DO i = 1, klon
	217	! -layer calculation
	218	rhodz(i, k) = (paprs(i,k)-paprs(i,k+1))/rg ! kg/m2
	219	zrho(i, k) = pplay(i, k)/t(i, k)/rd ! kg/m3
	220	dh(i, k) = rhodz(i, k)/zrho(i, k) ! m
	221	! -IM Total Liquid/Ice water content
	222	xflwc(i, k) = (1.-zfice(i,k))*pqlwp(i, k)
	223	xfiwc(i, k) = zfice(i, k)*pqlwp(i, k)
	224	END DO
	225	END DO
	226	ENDIF
[1992]	227
	228	IF (ok_cdnc) THEN
	229
	230	! --we compute cloud properties as a function of the aerosol load
	231
	232	DO k = 1, klev
	233	DO i = 1, klon
	234
	235	! Formula "D" of Boucher and Lohmann, Tellus, 1995
	236	! Cloud droplet number concentration (CDNC) is restricted
	237	! to be within [20, 1000 cm^3]
	238
	239	! --present-day case
	240	cdnc(i, k) = 10.*(bl95_b0+bl95_b1log(max(mass_solu_aero(i,k), &
	241	1.E-4))/log(10.))*1.E6 !-m-3
	242	cdnc(i, k) = min(1000.E6, max(20.E6,cdnc(i,k)))
	243
	244	! --pre-industrial case
	245	cdnc_pi(i, k) = 10.*(bl95_b0+bl95_b1log(max(mass_solu_aero_pi(i,k), &
	246	1.E-4))/log(10.))*1.E6 !-m-3
	247	cdnc_pi(i, k) = min(1000.E6, max(20.E6,cdnc_pi(i,k)))
	248
	249	! --present-day case
	250	rad_chaud(i, k) = 1.1((pqlwp(i,k)pplay(i, &
	251	k)/(rdt(i,k)))/(4./3rpi1000.cdnc(i,k)))**(1./3.)
	252	rad_chaud(i, k) = max(rad_chaud(i,k)*1.E6, 5.)
	253
	254	! --pre-industrial case
	255	rad_chaud_pi(i, k) = 1.1((pqlwp(i,k)pplay(i, &
	256	k)/(rdt(i,k)))/(4./3.rpi1000.cdnc_pi(i,k)))**(1./3.)
	257	rad_chaud_pi(i, k) = max(rad_chaud_pi(i,k)*1.E6, 5.)
	258
	259	! --pre-industrial case
	260	! --liquid/ice cloud water paths:
	261	IF (pclc(i,k)<=seuil_neb) THEN
	262
	263	pcldtaupi(i, k) = 0.0
	264
	265	ELSE
	266
	267	zflwp_var = 1000.(1.-zfice(i,k))pqlwp(i, k)/pclc(i, k)* &
	268	rhodz(i, k)
	269	zfiwp_var = 1000.zfice(i, k)pqlwp(i, k)/pclc(i, k)*rhodz(i, k)
	270	tc = t(i, k) - 273.15
	271	rei = d_rei_dt*tc + rei_max
	272	IF (tc<=-81.4) rei = rei_min
	273
	274	! -- cloud optical thickness :
	275	! [for liquid clouds, traditional formula,
	276	! for ice clouds, Ebert & Curry (1992)]
	277
	278	IF (zfiwp_var==0. .OR. rei<=0.) rei = 1.
	279	pcldtaupi(i, k) = 3.0/2.0*zflwp_var/rad_chaud_pi(i, k) + &
	280	zfiwp_var*(3.448E-03+2.431/rei)
	281
	282	END IF
	283
	284	END DO
	285	END DO
	286
	287	ELSE !--not ok_cdnc
	288
	289	! -prescribed cloud droplet radius
	290
	291	DO k = 1, min(3, klev)
	292	DO i = 1, klon
	293	rad_chaud(i, k) = rad_chau2
	294	rad_chaud_pi(i, k) = rad_chau2
	295	END DO
	296	END DO
	297	DO k = min(3, klev) + 1, klev
	298	DO i = 1, klon
	299	rad_chaud(i, k) = rad_chau1
	300	rad_chaud_pi(i, k) = rad_chau1
	301	END DO
	302	END DO
	303
	304	END IF !--ok_cdnc
	305
	306	! --computation of cloud optical depth and emissivity
	307	! --in the general case
	308
	309	DO k = 1, klev
	310	DO i = 1, klon
	311
	312	IF (pclc(i,k)<=seuil_neb) THEN
	313
	314	! effective cloud droplet radius (microns) for liquid water clouds:
	315	! For output diagnostics cloud droplet effective radius [um]
	316	! we multiply here with f * xl (fraction of liquid water
	317	! clouds in the grid cell) to avoid problems in the averaging of the
	318	! output.
	319	! In the output of IOIPSL, derive the REAL cloud droplet
	320	! effective radius as re/fl
	321
	322	fl(i, k) = seuil_neb*(1.-zfice(i,k))
	323	re(i, k) = rad_chaud(i, k)*fl(i, k)
	324	rel = 0.
	325	rei = 0.
	326	pclc(i, k) = 0.0
	327	pcltau(i, k) = 0.0
	328	pclemi(i, k) = 0.0
	329
	330	ELSE
	331
	332	! -- liquid/ice cloud water paths:
	333
	334	zflwp_var = 1000.(1.-zfice(i,k))pqlwp(i, k)/pclc(i, k)*rhodz(i, k)
	335	zfiwp_var = 1000.zfice(i, k)pqlwp(i, k)/pclc(i, k)*rhodz(i, k)
	336
	337	! effective cloud droplet radius (microns) for liquid water clouds:
	338	! For output diagnostics cloud droplet effective radius [um]
	339	! we multiply here with f * xl (fraction of liquid water
	340	! clouds in the grid cell) to avoid problems in the averaging of the
	341	! output.
	342	! In the output of IOIPSL, derive the REAL cloud droplet
	343	! effective radius as re/fl
	344
	345	fl(i, k) = pclc(i, k)*(1.-zfice(i,k))
	346	re(i, k) = rad_chaud(i, k)*fl(i, k)
	347
	348	rel = rad_chaud(i, k)
	349
	350	! for ice clouds: as a function of the ambiant temperature
	351	! [formula used by Iacobellis and Somerville (2000), with an
	352	! asymptotical value of 3.5 microns at T<-81.4 C added to be
	353	! consistent with observations of Heymsfield et al. 1986]:
	354	! 2011/05/24 : rei_min = 3.5 becomes a free PARAMETER as well as
	355	! rei_max=61.29
	356
	357	tc = t(i, k) - 273.15
	358	rei = d_rei_dt*tc + rei_max
	359	IF (tc<=-81.4) rei = rei_min
	360
	361	! -- cloud optical thickness :
	362	! [for liquid clouds, traditional formula,
	363	! for ice clouds, Ebert & Curry (1992)]
	364
	365	IF (zflwp_var==0.) rel = 1.
	366	IF (zfiwp_var==0. .OR. rei<=0.) rei = 1.
	367	pcltau(i, k) = 3.0/2.0(zflwp_var/rel) + zfiwp_var(3.448E-03+2.431/ &
	368	rei)
	369
	370	! -- cloud infrared emissivity:
	371	! [the broadband infrared absorption coefficient is PARAMETERized
	372	! as a function of the effective cld droplet radius]
	373	! Ebert and Curry (1992) formula as used by Kiehl & Zender (1995):
	374
	375	k_ice = k_ice0 + 1.0/rei
	376
	377	pclemi(i, k) = 1.0 - exp(-coef_chauzflwp_var-dfk_ice*zfiwp_var)
	378
	379	END IF
	380
	381	reice(i, k) = rei
	382
	383	xflwp(i) = xflwp(i) + xflwc(i, k)*rhodz(i, k)
	384	xfiwp(i) = xfiwp(i) + xfiwc(i, k)*rhodz(i, k)
	385
	386	END DO
	387	END DO
	388
	389	! --if cloud droplet radius is fixed, then pcldtaupi=pcltau
	390
	391	IF (.NOT. ok_cdnc) THEN
	392	DO k = 1, klev
	393	DO i = 1, klon
	394	pcldtaupi(i, k) = pcltau(i, k)
	395	reice_pi(i, k) = reice(i, k)
	396	END DO
	397	END DO
	398	END IF
	399
	400	DO k = 1, klev
	401	DO i = 1, klon
	402	reliq(i, k) = rad_chaud(i, k)
	403	reliq_pi(i, k) = rad_chaud_pi(i, k)
	404	reice_pi(i, k) = reice(i, k)
	405	END DO
	406	END DO
	407
	408	! COMPUTE CLOUD LIQUID PATH AND TOTAL CLOUDINESS
	409	! IM cf. CR:test: calcul prenant ou non en compte le recouvrement
	410	! initialisations
	411
	412	DO i = 1, klon
	413	zclear(i) = 1.
	414	zcloud(i) = 0.
	415	zcloudh(i) = 0.
	416	zcloudm(i) = 0.
	417	zcloudl(i) = 0.
	418	pch(i) = 1.0
	419	pcm(i) = 1.0
	420	pcl(i) = 1.0
	421	pctlwp(i) = 0.0
	422	END DO
	423
	424	! --calculation of liquid water path
	425
	426	DO k = klev, 1, -1
	427	DO i = 1, klon
	428	pctlwp(i) = pctlwp(i) + pqlwp(i, k)*rhodz(i, k)
	429	END DO
	430	END DO
	431
	432	! --calculation of cloud properties with cloud overlap
	433
	434	IF (novlp==1) THEN
	435	DO k = klev, 1, -1
	436	DO i = 1, klon
	437	zclear(i) = zclear(i)*(1.-max(pclc(i,k),zcloud(i)))/(1.-min(real( &
	438	zcloud(i),kind=8),1.-zepsec))
	439	pct(i) = 1. - zclear(i)
	440	IF (paprs(i,k)<prmhc) THEN
	441	pch(i) = pch(i)*(1.-max(pclc(i,k),zcloudh(i)))/(1.-min(real(zcloudh &
	442	(i),kind=8),1.-zepsec))
	443	zcloudh(i) = pclc(i, k)
	444	ELSE IF (paprs(i,k)>=prmhc .AND. paprs(i,k)<prlmc) THEN
	445	pcm(i) = pcm(i)*(1.-max(pclc(i,k),zcloudm(i)))/(1.-min(real(zcloudm &
	446	(i),kind=8),1.-zepsec))
	447	zcloudm(i) = pclc(i, k)
	448	ELSE IF (paprs(i,k)>=prlmc) THEN
	449	pcl(i) = pcl(i)*(1.-max(pclc(i,k),zcloudl(i)))/(1.-min(real(zcloudl &
	450	(i),kind=8),1.-zepsec))
	451	zcloudl(i) = pclc(i, k)
	452	END IF
	453	zcloud(i) = pclc(i, k)
	454	END DO
	455	END DO
	456	ELSE IF (novlp==2) THEN
	457	DO k = klev, 1, -1
	458	DO i = 1, klon
	459	zcloud(i) = max(pclc(i,k), zcloud(i))
	460	pct(i) = zcloud(i)
	461	IF (paprs(i,k)<prmhc) THEN
	462	pch(i) = min(pclc(i,k), pch(i))
	463	ELSE IF (paprs(i,k)>=prmhc .AND. paprs(i,k)<prlmc) THEN
	464	pcm(i) = min(pclc(i,k), pcm(i))
	465	ELSE IF (paprs(i,k)>=prlmc) THEN
	466	pcl(i) = min(pclc(i,k), pcl(i))
	467	END IF
	468	END DO
	469	END DO
	470	ELSE IF (novlp==3) THEN
	471	DO k = klev, 1, -1
	472	DO i = 1, klon
	473	zclear(i) = zclear(i)*(1.-pclc(i,k))
	474	pct(i) = 1 - zclear(i)
	475	IF (paprs(i,k)<prmhc) THEN
	476	pch(i) = pch(i)*(1.0-pclc(i,k))
	477	ELSE IF (paprs(i,k)>=prmhc .AND. paprs(i,k)<prlmc) THEN
	478	pcm(i) = pcm(i)*(1.0-pclc(i,k))
	479	ELSE IF (paprs(i,k)>=prlmc) THEN
	480	pcl(i) = pcl(i)*(1.0-pclc(i,k))
	481	END IF
	482	END DO
	483	END DO
	484	END IF
	485
	486	DO i = 1, klon
	487	pch(i) = 1. - pch(i)
	488	pcm(i) = 1. - pcm(i)
	489	pcl(i) = 1. - pcl(i)
	490	END DO
	491
	492	! ========================================================
	493	! DIAGNOSTICS CALCULATION FOR CMIP5 PROTOCOL
	494	! ========================================================
	495	! change by Nicolas Yan (LSCE)
	496	! Cloud Droplet Number Concentration (CDNC) : 3D variable
	497	! Fractionnal cover by liquid water cloud (LCC3D) : 3D variable
	498	! Cloud Droplet Number Concentration at top of cloud (CLDNCL) : 2D variable
	499	! Droplet effective radius at top of cloud (REFFCLWTOP) : 2D variable
	500	! Fractionnal cover by liquid water at top of clouds (LCC) : 2D variable
	501
	502	IF (ok_cdnc) THEN
	503
	504	DO k = 1, klev
	505	DO i = 1, klon
	506	phase3d(i, k) = 1 - zfice(i, k)
	507	IF (pclc(i,k)<=seuil_neb) THEN
	508	lcc3d(i, k) = seuil_neb*phase3d(i, k)
	509	ELSE
	510	lcc3d(i, k) = pclc(i, k)*phase3d(i, k)
	511	END IF
	512	scdnc(i, k) = lcc3d(i, k)*cdnc(i, k) ! m-3
	513	END DO
	514	END DO
	515
	516	DO i = 1, klon
	517	lcc(i) = 0.
	518	reffclwtop(i) = 0.
	519	cldncl(i) = 0.
[2641]	520	IF (novlp.EQ.3 .OR. novlp.EQ.1) tcc(i) = 1.
	521	IF (novlp.EQ.2) tcc(i) = 0.
[1992]	522	END DO
	523
	524	DO i = 1, klon
	525	DO k = klev - 1, 1, -1 !From TOA down
	526
	527	! Test, if the cloud optical depth exceeds the necessary
	528	! threshold:
	529
	530	IF (pcltau(i,k)>thres_tau .AND. pclc(i,k)>thres_neb) THEN
	531
[2641]	532	IF (novlp.EQ.2) THEN
[1992]	533	IF (first) THEN
	534	WRITE (, ) 'Hypothese de recouvrement: MAXIMUM'
	535	first = .FALSE.
	536	END IF
	537	flag_max = -1.
	538	ftmp(i) = max(tcc(i), pclc(i,k))
	539	END IF
	540
[2641]	541	IF (novlp.EQ.3) THEN
[1992]	542	IF (first) THEN
	543	WRITE (, ) 'Hypothese de recouvrement: RANDOM'
	544	first = .FALSE.
	545	END IF
	546	flag_max = 1.
	547	ftmp(i) = tcc(i)*(1-pclc(i,k))
	548	END IF
	549
[2641]	550	IF (novlp.EQ.1) THEN
[1992]	551	IF (first) THEN
	552	WRITE (, ) 'Hypothese de recouvrement: MAXIMUM_ &
	553	& &
	554	& RANDOM'
	555	first = .FALSE.
	556	END IF
	557	flag_max = 1.
	558	ftmp(i) = tcc(i)*(1.-max(pclc(i,k),pclc(i,k+1)))/(1.-min(pclc(i, &
	559	k+1),1.-thres_neb))
	560	END IF
	561	! Effective radius of cloud droplet at top of cloud (m)
	562	reffclwtop(i) = reffclwtop(i) + rad_chaud(i, k)1.0E-06phase3d(i, &
	563	k)(tcc(i)-ftmp(i))flag_max
	564	! CDNC at top of cloud (m-3)
	565	cldncl(i) = cldncl(i) + cdnc(i, k)phase3d(i, k)(tcc(i)-ftmp(i))* &
	566	flag_max
	567	! Liquid Cloud Content at top of cloud
	568	lcc(i) = lcc(i) + phase3d(i, k)(tcc(i)-ftmp(i))flag_max
	569	! Total Cloud Content at top of cloud
	570	tcc(i) = ftmp(i)
	571
	572	END IF ! is there a visible, not-too-small cloud?
	573	END DO ! loop over k
	574
[2641]	575	IF (novlp.EQ.3 .OR. novlp.EQ.1) tcc(i) = 1. - tcc(i)
[1992]	576
	577	END DO ! loop over i
	578
	579	! ! Convective and Stratiform Cloud Droplet Effective Radius (REFFCLWC
	580	! REFFCLWS)
	581	DO i = 1, klon
[524]	582	DO k = 1, klev
[1992]	583	! Weight to be used for outputs: eau_liquide*couverture nuageuse
	584	lcc3dcon(i, k) = rnebcon(i, k)phase3d(i, k)clwcon(i, k) ! eau liquide convective
	585	lcc3dstra(i, k) = pclc(i, k)pqlwp(i, k)phase3d(i, k)
	586	lcc3dstra(i, k) = lcc3dstra(i, k) - lcc3dcon(i, k) ! eau liquide stratiforme
	587	lcc3dstra(i, k) = max(lcc3dstra(i,k), 0.0)
	588	! Compute cloud droplet radius as above in meter
	589	radius = 1.1((pqlwp(i,k)pplay(i,k)/(rdt(i,k)))/(4./3rpi1000. &
	590	cdnc(i,k)))**(1./3.)
	591	radius = max(radius, 5.E-6)
	592	! Convective Cloud Droplet Effective Radius (REFFCLWC) : variable 3D
	593	reffclwc(i, k) = radius
	594	reffclwc(i, k) = reffclwc(i, k)*lcc3dcon(i, k)
	595	! Stratiform Cloud Droplet Effective Radius (REFFCLWS) : variable 3D
	596	reffclws(i, k) = radius
	597	reffclws(i, k) = reffclws(i, k)*lcc3dstra(i, k)
	598	END DO !klev
	599	END DO !klon
[524]	600
[1992]	601	! Column Integrated Cloud Droplet Number (CLDNVI) : variable 2D
[524]	602
[1992]	603	DO i = 1, klon
	604	cldnvi(i) = 0.
	605	lcc_integrat(i) = 0.
	606	height(i) = 0.
[1989]	607	DO k = 1, klev
[1992]	608	cldnvi(i) = cldnvi(i) + cdnc(i, k)lcc3d(i, k)dh(i, k)
	609	lcc_integrat(i) = lcc_integrat(i) + lcc3d(i, k)*dh(i, k)
	610	height(i) = height(i) + dh(i, k)
	611	END DO ! klev
	612	lcc_integrat(i) = lcc_integrat(i)/height(i)
	613	IF (lcc_integrat(i)<=1.0E-03) THEN
	614	cldnvi(i) = cldnvi(i)*lcc(i)/seuil_neb
	615	ELSE
	616	cldnvi(i) = cldnvi(i)*lcc(i)/lcc_integrat(i)
	617	END IF
	618	END DO ! klon
[1337]	619
[1992]	620	DO i = 1, klon
	621	DO k = 1, klev
	622	IF (scdnc(i,k)<=0.0) scdnc(i, k) = 0.0
	623	IF (reffclws(i,k)<=0.0) reffclws(i, k) = 0.0
	624	IF (reffclwc(i,k)<=0.0) reffclwc(i, k) = 0.0
	625	IF (lcc3d(i,k)<=0.0) lcc3d(i, k) = 0.0
	626	IF (lcc3dcon(i,k)<=0.0) lcc3dcon(i, k) = 0.0
	627	IF (lcc3dstra(i,k)<=0.0) lcc3dstra(i, k) = 0.0
	628	END DO
	629	IF (reffclwtop(i)<=0.0) reffclwtop(i) = 0.0
	630	IF (cldncl(i)<=0.0) cldncl(i) = 0.0
	631	IF (cldnvi(i)<=0.0) cldnvi(i) = 0.0
	632	IF (lcc(i)<=0.0) lcc(i) = 0.0
	633	END DO
[1337]	634
[1992]	635	END IF !ok_cdnc
[1337]	636
[1992]	637	RETURN
[1337]	638
[1992]	639	END SUBROUTINE newmicro

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

Download in other formats:

Original Format