readDiag is a Python package that provides a class to read and plot the Gridpoint Statistical Interpolation (GSI) diagnostics files. It can be used to retrieve and investigate important information from the data assimilation process