Abstract: Recent progress in materials synthesis enabled fabrication of superconducting atomically flat single-crystalline metallic nanofilms with thicknesses down to a few monolayers. Interest in such nano-thin systems is attracted by the dimensional 3D-2D crossover in their coherent properties which occurs with decreasing the film thickness. The first fundamental aspect of this crossover is dictated by the Mermin-Wagner-Hohenberg theorem and concerns frustration of the long-range order due to superconductive fluctuations and the possibility to track its impact with an unprecedented level of control. The second important aspect is related to the Fabri-Perot modes of the electronic motion strongly bound in the direction perpendicular to the nanofilm. The formation of such modes results in a pronounced multiband structure that changes with the nanofilm thickness and affects both the mean-field behavior and superconductive fluctuations. Though the subject is very rich in physics, it is scarcely investigated to date. The main obstacle is that there are no manageable models to study a complex magnetic response in this case. Full microscopic consideration is rather time consuming, if practicable at all, while the standard Ginzburg-Landau theory is not applicable. In the present work we review the main achievements in the subject to date, and construct and justify an efficient multiband mean-field formalism which allows for numerical and even analytical treatment of nano-thin superconductors in applied magnetic fields.