Height distribution of ionospheric plasma parameters in the F2-region is closely related to height distribution of the main thermospheric parameters. Therefore, they can be extracted from ionospheric observations solving an inverse problem of aeronomy. A self-consistent approach to the Ne(h) modeling at the F2-region heights has been applied to solve the problem. Using routine incoherent scatter radar observations (N e(h), Te(h), Ti(h), Vi(h) profiles) the method yields a self-consistent set of main aeronomic parameters responsible for the F2-region formation. The list of derived parameters includes: neutral temperature profile Tn(h) depending on the exospheric temperature Tex, the temperature at 120 km T120 and the shape-parameter S, which determine the temperature profile, concentration of neutral species [O], [O2], [N2], vertical plasma drift W, which may be converted to the meridional thermospheric wind Vnx, total solar EUV flux and ion composition (O + , O2 + , NO + , N2 + , N + ) as a result of Ne(h) fitting. Therefore, the method gives a complete description of the upper atmosphere condition in the vicinity of incoherent scatter facility for the periods of observation. Analysis of all available EIS-CAT (CP-1, CP-2) observations has shown wide deviations from MSIS-86 model predictions for geomagnetically disturbed conditions while the retrieved parameters are close to the model ones for quiet periods. The approach turns out to be very useful for physical analyses of the F2-layer disturbance mechanisms giving a complete picture of the phenomenon in question. Limitations and problems related to method application are discussed. Under existing conditions when thermospheric observations are not conducted currently the proposed method may be considered a real tool for thermosphere investigation and monitoring at least for the periods of ISR observations.