mfpblt_mod Module Reference

This module contains the subroutine that calculates mass flux and updraft parcel properties for thermals driven by surface heating for use in the TKE-EDMF PBL scheme. More...

Functions/Subroutines
subroutine	mfpblt (im, ix, km, kmpbl, ntcw, ntrac1, delt, cnvflg, zl, zm, q1, t1, u1, v1, plyr, pix, thlx, thvx, gdx, hpbl, kpbl, vpert, buo, xmf, tcko, qcko, ucko, vcko, xlamue)
	This subroutine computes mass flux and updraft parcel properties for thermals driven by surface heating.

Detailed Description

This module contains the subroutine that calculates mass flux and updraft parcel properties for thermals driven by surface heating for use in the TKE-EDMF PBL scheme.

Function/Subroutine Documentation

mfpblt()

subroutine mfpblt_mod::mfpblt	(	integer	im,
		integer	ix,
		integer	km,
		integer	kmpbl,
		integer	ntcw,
		integer	ntrac1,
		real(kind=kind_phys)	delt,
		logical, dimension(im)	cnvflg,
		real(kind=kind_phys), dimension(im,km)	zl,
		real(kind=kind_phys), dimension(im,km)	zm,
		real(kind=kind_phys), dimension(ix,km,ntrac1)	q1,
		real(kind=kind_phys), dimension(ix,km)	t1,
		real(kind=kind_phys), dimension(ix,km)	u1,
		real(kind=kind_phys), dimension(ix,km)	v1,
		real(kind=kind_phys), dimension(im,km)	plyr,
		real(kind=kind_phys), dimension(im,km)	pix,
		real(kind=kind_phys), dimension(im,km)	thlx,
		real(kind=kind_phys), dimension(im,km)	thvx,
		real(kind=kind_phys), dimension(im)	gdx,
		real(kind=kind_phys), dimension(im)	hpbl,
		integer, dimension(im)	kpbl,
		real(kind=kind_phys), dimension(im)	vpert,
		real(kind=kind_phys), dimension(im,km)	buo,
		real(kind=kind_phys), dimension(im,km)	xmf,
		real(kind=kind_phys), dimension(im,km)	tcko,
		real(kind=kind_phys), dimension(im,km,ntrac1)	qcko,
		real(kind=kind_phys), dimension(im,km)	ucko,
		real(kind=kind_phys), dimension(im,km)	vcko,
		real(kind=kind_phys), dimension(im,km-1)	xlamue
	)

This subroutine computes mass flux and updraft parcel properties for thermals driven by surface heating.

GFS mfpblt General Algorithm

• Compute thermal excess

• Compute entrainment rate (eqn 30 of Han et al. (2019) [77])

• Compute buoyancy for updraft air parcel

• Compute updraft velocity square(wu2, eqn 13 in Han et al.(2019) [77])

• Update pbl height as the height where updraft velocity vanishes

• Update entrainment rate

• Compute entrainment rate averaged over the whole pbl

• Updraft mass flux as a function of updraft velocity profile (eqn 12 of Han et al.(2019) [77])

• Compute updraft fraction as a function of mean entrainment rate (Grell and Freitas (2014) [73]

• Compute scale-aware function based on Arakawa and Wu (2013) [7]

• Final scale-aware updraft mass flux

• Compute updraft property using updated entranment rate

Definition at line 16 of file mfpblt.f.