Mosaicing of multitemporal Sentinel-2 images

My aim is to combine several Sentinel-2 Level2A images of the same location/orbit number but at different dates of a year.

I would like to have different options to create that mosaicing (e.g., if multiple data per pixel is available than create mean/median/std etc., if no data is available than keep Nodata, of only one dataset has data than keep them).

Is there any tool in ENVI which can be used for such mosaicing processes? Alternativly I could also work in ERDAS, SNAP or ESRI ArcMap/ArcPro.