subroutine mfpbltq	(	integer	im,
		integer	ix,
		integer	km,
		integer	kmpbl,
		integer	ntcw,
		integer	ntrac1,
		real(kind=kind_phys)	delt
	)

GFS mfpblt General Algorithm

• Compute thermal excess
• Compute entrainment rate
• Compute buoyancy for updraft air parcel
• Compute updraft velocity square(wu2, eqn 13 in Han et al.(2019) [73])
• 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
• Compute updraft fraction as a function of mean entrainment rate (Grell and Freitas (2014) [70]
• Compute scale-aware function based on Arakawa and Wu (2013) [9]
• Final scale-aware updraft mass flux
• Compute updraft property using updated entranment rate

References physcons::con_cp, physcons::con_eps, physcons::con_epsm1, physcons::con_fvirt, physcons::con_g, physcons::con_hvap, physcons::con_rv, and funcphys::fpvs().

Referenced by satmedmfvdifq_run().

Here is the call graph for this function: