When parsing a CASTEP formatted grid with python, you may be tempted to use numpy.loadtxt – in particular, as it allows for fast access to the numerical data afterwards. However, you need to know the number of grid points in each dimension in order to work out the proper arguments […]