Hi Andres,
The key parameter to set in run_Aem3d for BBL mixing is ibblmixing. Set this to one.
Without this bottom energy is generated and added to the bottom cell of each column but is only used in mixing of the bottom two cells. WIth ibblmixing switched on there is an additional upsweep of TKE mixing.
Chris