Not yet reviewed by Pith; the record is open.

Context: Traditional one-dimensional (1D) hydrostatic model atmospheres introduce systematic modelling errors into spectroscopic analyses of FGK-type stars. Aims: We present an updated version of the STAGGER-grid of 3D model atmospheres, and explore the accuracy of post-processing methods in preparation for spectral synthesis. Methods: New and old models were (re)computed following an updated workflow, including an updated opacity binning technique. Spectroscopic tests were performed in 3D LTE for a grid of 216 fictitious Fe I lines, spanning a wide range in oscillator strength, excitation potential and central wavelength, and eight model atmospheres that cover the stellar atmospheric parameter range (Teff, log g, [Fe/H]) of FGK-type stars. Using this grid, the impact of vertical and horizontal resolution, and temporal sampling of model atmospheres on spectroscopic diagnostics was tested. Results: We find that downsampling the horizontal mesh from its original size of 240 x 240 grid cells to 80 x 80 cells, i.e. sampling every third grid cell, introduces minimal errors on the equivalent width and normalized line flux across the line and stellar parameter space. Regarding temporal sampling, we find that sampling ten statistically independent snapshots is sufficient to accurately model the shape of spectral line profiles. For equivalent widths, a subsample consisting of only two snapshots is sufficient, introducing an abundance error of less than 0.015 dex. Conclusions: We have computed 32 new model atmospheres and recomputed 116 old model atmospheres present in the original grid. The public release of the STAGGER-grid contains 243 models, excluding models with [Fe/H] = -4.00, and the processed snapshots can be used to improve the accuracy of spectroscopic analyses.