Within the framework of free interaction theory numerical methods are used to investigate the occurrence of supersonic zones with shocks in the outer inviscid region for flow past roughness in the lower viscous sublayer, with and without the formation of local separation zones.