The polar atmospheric surface layer is often stably stratified, which strongly influences turbulent transport processes between the atmosphere and sea ice/ocean. Transport is usually parametrized applying Monin Obukhov Similarity Theory (MOST) which delivers transfer coefficients as a function of stability parameters (see below). In a series of papers (Gryanik and Lüpkes, 2018; Gryanik et al., 2020,2021; Gryanik and Lüpkes, 2022) it has been shown that differences between existing parametrizations are large, especially for strong stability. One reason is that they are based on different data sets, for which the origin of differences is still unclear. In this situation Gryanik et al. (2021) as well as Gryanik and Lüpkes (2022) proposed a numerically efficient method, which can be used for most of the existing data sets and their specific stability dependences. A package of parametrization resulted that is suitable for its application in weather prediction and climate models. Especially, calculation of fluxes over sea ice were improved. Combined with latest parametrizations of surface roughness it has a large impact on large scale fields as shown recently by Schneider et al. (2021) who applied some members of the package.