There is growing evidence, from experiments and numerical simulations, that a key feature of sufficiently disordered superconductors is the spatial inhomogeneity of the order parameter. However not much is known analytically about the details of its spatial distribution or the associated global critical temperature that signals the breaking of long-range order. Here we address this problem for disordered systems around an Anderson transition characterized by multifractal one-body eigenstates. In the limit of weak multifractality and for weakly coupled superconductors we compute the superconducting order parameter analytically, including its energy dependence and statistical distribution in space. The spatial distribution of the order parameter is found to be always lognormal. The global critical temperature, computed by percolation techniques and neglecting phase fluctuations, is enhanced with respect to the clean limit only for very weakly coupled superconductors. Some enhancement still persists even in the presence of moderate phase fluctuations crudely modelled by increasing the percolation threshold. Our results are also consistent with experiments, where enhancement of the critical temperature is observed in Al thin films, a very weakly coupled metallic superconductor, but not in more strongly coupled materials.PACS numbers: 74.78. Na, 75.10.Pq For many years the role of disorder in superconductivity was believed to be well understood. According to the so called Anderson theorem [1], also stated independently by Gor'kov and Abrikosov [2], the critical temperature of a conventional weakly-coupled superconductor is not affected by weak non-magnetic impurity scattering. These results are based on the assumption that the local density of states in the material is unaffected by weak disorder [3,4]. However with the development of the Bogoliubovde Gennes theory of superconductivity [5] it became clear that the order parameter becomes increasingly inhomogeneous with increasing disorder.Experimentally it is well established [6][7][8][9][10][11][12][13][14], especially for conventional superconducting thin films, that the critical temperature decreases monotonically as disorder increases. Analytic results [15,16], obtained using mesoscopic techniques, confirmed that the interplay between weak disorder and Coulomb interactions could explain this suppression of the critical temperature. For stronger disorder around the superconductor insulator transition there is recent numerical [17,18] evidence that, even in the absence of Coulomb interactions, phase fluctuations are enhanced [19] and the superconducting order parameter becomes highly inhomogeneous [20,21]. Close to the Berezinski-Kosterlitz-Thouless transition phase correlation only persist along a ramified network, reminiscent of a percolation transition [22]. This is consistent with experimental observations of a universal scaling of the order parameter amplitude distribution function [23], emergent granularity [24,25] and reports of glassy features [26], with...