  We introduce a new class of inhomogeneous cosmological models as solutions to
the Einstein-Maxwell equations in electrovacuum. The new models can be
considered to be nonlinear perturbations, through an electromagnetic field, of
the previously studied `smooth Gowdy-symmetric generalised Taub-NUT solutions'
in vacuum. Utilising methods from soliton theory, we analyse the effects of the
Maxwell field on global properties of the solutions. In particular, we show
existence of regular Cauchy horizons, and we investigate special singular cases
in which curvature singularities form.
