We evaluate validity of NLFFF extrapolation performed with Optimization class (OPTI) codes. While explaining inevitable for OPTI partial non-solenoidality caused by the gas pressure notable role in pressure balance at photospheric heights and by mathematical aspects related to optimization and BVP (boundary value problem), we justify elimination of the non-solenoidal component (postprocessing) from the OPTI result obtained. In essence, postprocessing converts the entire non-solenoidal part into a solenoidal force part, which possibly reflects factual deviation of magnetic field from its force-free pproximation on the photosphere and in the solar corona. Two forms of postprocessing have been analyzed in this paper.Postprocessing I eliminates the non-solenoidal component without changing transverse field at the measurement level, and Postprocessing II leaves the field normal component unchanged. Extrapolation, postprocessing, and then comparison of metric and energy characteristics are performed over AR 11158 active region for a small fragment of its evolution containing the February X-class flare.Our version of OPTI code showed that free energy decreased by $\sim 10^{32}$ erg within 1 hour, which corresponds to theoretical estimations of the flare-caused magnetic energy loss. This result differs significantly from the one in Sun et al. (2012).Therefore, we also comment on some features of our OPTI code implementation, which may cause significant differences between our results and those obtained using the Wiegelmann (2004) version of OPTI code in study by Sun et al. (2012).

In this paper, we came to conclusion that there is a significant systematic error in the SDO/HMI vector magnetic data, which reveals itself in a significant deviation of the lines of the knot magnetic fields from the radial direction. The value of this deviation demonstrates a clear dependence on the distance to the disk center. This paper suggests a method for correction of the vector magnetograms that eliminates the detected systematic error.

A method of construction of solution for acoustic-gravity waves (AGW) above a wave source, taking dissipation throughout the atmosphere into account (Dissipative Solution above Source, DSAS), is proposed. The method is to combine three solutions for three parts of the atmosphere: an analytical solution for the upper isothermal part and numerical solutions for the real non-isothermal dissipative atmosphere in the middle part and for the real non-isothermal small dissipation atmosphere in the lower one. In this paper the method has been carried out for the atmosphere with thermal conductivity but without viscosity. The heights of strong dissipation and the total absorption index in the regions of weak and average dissipation are found. For internal gravity waves the results of test calculations for an isothermal atmosphere and calculations for a real non-isothermal atmosphere are shown in graphical form. An algorithm and appropriate code to calculate DSAS, taking dissipation due to finite thermal conductivity into account throughout the atmosphere, are developed. The results of test DSAS calculations for an everywhere isothermal atmosphere are given. The calculation results for DSAS for the real non-isothermal atmosphere are also presented. A method for constructing of the 2x2 Green's matrix fully taking dissipation into account and allowing to find disturbance from some source of AGW in the atmosphere is proposed.

Trapped atmosphere waves, such as IGW waveguide modes and Lamb modes, are described using dissipative solution above source (DSAS) (Dmitrienko and Rudenko, 2016). Accordingly this description, the modes are disturbances penetrating without limit in the upper atmosphere and dissipating their energy throughout the atmosphere; leakage from a trapping region to the upper atmosphere is taken in consideration. The DSAS results are compared to those based on both accurate and WKB approximated dissipationless equations. It is shown that the spatial and frequency characteristics of modes in the upper atmosphere calculated by any of the methods are close to each other and are in good agreement with the observed characteristics of traveling ionospheric disturbances.