Hello, the YSU PBL scheme has the diffusivity coefficient Km scaling with (1 - z/h_pbl) where h_pbl is the diagnosed boundary layer height. What happens when z > h_pbl, say for shallow boundary layers? How is turbulent mixing modeled above the boundary layer?
Thanks!