Resolved NMR spectra from samples in inhomogeneous B0 and B1 fields can be obtained with the so-called "ex situ" methodology, employing a train of composite or adiabatic z-rotation RF pulses to periodically refocus the inhomogeneous broadening during the detection of the time-domain signal. Earlier schemes relied on a linear correlation between the inhomogeneous B0 and B1 fields. Here the pulse length, bandwidth, and amplitude of the adiabatic pulses of the hyperbolic secant type are adjusted to improve the refocusing for a setup with non-linear correlation. The field correlation is measured using a two-dimensional nutation experiment augmented with a third dimension with varying RF carrier frequency accounting for off-resonance effects. The pulse optimization is performed with a computer algorithm using the experimentally determined field correlation and a standard adiabatic z-rotation pulse as a starting point for the iterative optimization procedure. The shape of the z-rotation RF pulse is manipulated to provide refocusing for the conditions given by the sample-, magnet-, and RF-coil geometry.