ROMS mixed layer depth output

[Honolulu_Wang]

Hi All,
May I ask the possible and simple way for obtaining the ROMS mixed layer depth under GLS mixing scheme.
The ROMS only has surface boundary depth for application under LMD scheme, but that scheme mixes almost whole water column. Any suggestion?

Thanks