Meteorological centres compare global and large-scale regional model output versus radiosonde data more or less continuously for specific models and also make inter-model comparisons. They normally provide comparisons of meteorological variables at standard pressure levels, and similar comparisons for finer-scale models may be found as output from field tests and other experiments. Additional means to evaluate vertical profiles derived from model output more fully should help provide a more complete evaluation of model output compared with observations. Here, two versions of a program to produce vertical profiles of meteorological variables in terms of user input heights or pressure levels are developed, and integrated mean value profiles of the variables for layers defined by those heights or pressures are generated. Output profiles were entered into spreadsheets for calculation of the derived variables of density and vector wind speed. The level and mean layer values from the Weather Research and Forecasting model (version 3.6.1) were compared with co-incident World Meteorological Organization radiosonde observation data and the differences were computed for several variables. The methods developed here may be used for a variety of model-generated and observed vertical profiles with only minimal changes, primarily to the input function.