Software Open Access
markwieczorek; Matthias Meschede; Ilya Oshchepkov; Elliott Sales de Andrade; mreineck; Armin Corbin; xoviat; Benda Xu; Stefan Schröder; Akihisa Hattori; Aaryaman Vasishta; Andrew Walker; Erik Schnetter; Katrin Leinweber
Implemented the option to use a different backend when performing certain operations requiring spherical harmonic transforms. At present, only 'shtools' (default) and the Distinctly Useful Code Collection ('ducc') are supported.
backendsthat has functions allowing the user to control which backend is used. To set the backend for all subsequent operations, use
nthreads('ducc' only) to all methods of the pyshtools classes that allow multiple backends (such as
backendsweb documentation page that describes the use of the new module.
SHGrid.plot_histogram()for generating area-weighted histograms.
titlesizefor consistency. This affects the methods
Iwith the alias
IO and datasets
ras an alias for backwards compatibility).
dotterm. When a line in the data section starts with an unknown key, a warning is now printed to the screen (which can be turned off by specifying
quiet=True). For files formatted as
icgem2.0, time variable terms were simply ignored if the specified epoch fell outside of the allowed range. Now, the routine will instead raise an error. Finally, the documentation was improved by describing the allowable keyword entries of the header and data section of the file.
encodingfor all routines and methods that read or write text-based spherical harmonic files.
utf-8in order to avoid problems with the XGM2019E dataset that has a character that can not be decoded by the GBK encoding that is the default in some Chinese installations.
SHCoeffs.gradient(). The original behavior was to compute the gradient on a sphere of radius
rwas the degree 0 coefficient of the function. The new behaviour is to compute the gradient on the unit sphere. This radius can be modified by supplying the optional argument
coeffs[1, lmax, lmax]was in error when
lmaxwas odd and when using
DHgrids. Second, when using
SHGrid.expand()with grid type
sampling=2was not passed to the Fortran routine.
makefile installto include example data files, and to place them in the correct directories when installing with homebrew.
SHGridthat could arise if the value input to arccos was greater than 1.
np.complex128to avoid numpy deprectation warning.
threadsafeto numpy signature files.
M. A. Wieczorek, M. Meschede, T. Brugere, A. Corbin, A. Hattori, K. Leinweber, I. Oshchepkov, M. Reinecke, E. Sales de Andrade, E. Schnetter, S. Schröder, A. Vasishta, A. Walker, B. Xu (2021). SHTOOLS: Version 4.9, Zenodo, doi:10.5281/zenodo.592762