This work represents a first step in the definition of a framework aimed at finding, by means of efficient global optimization based on metamodels, an optimal configuration of physical parameters for the ICON (ICOsahedral Nonhydrostatic) Limited Area Mode at high resolution (about 1.1 km) over Southern Italy, to be used for operational runs. The objective of the optimization is to reduce the distance between observed meteorological variables and modeled data. This distance is measured by an opportunely designed objective function. This work represents a preparatory step, since the input parameters considered are only a reduced number with respect to the huge amount of parameters potentially involved. First, domain size sensitivity was performed to choose the optimal domain. Then, the optimization was conducted by means of an Efficient Global Optimization algorithm relying on a Gaussian-based metamodel. The four parameters considered control the heat transfer in the turbulent layer, the laminar resistance and the snow vertical velocity. They were optimized over a week in November 2018, a period characterized by extreme events in the region considered. The results demonstrated the effectiveness of the proposed approach, reducing the distance from observed data, and the method can be considered promising from the perspective taking into account a larger set of physical parameters, and validation over a wider time-window.