this one: https://anaconda.org/eumetsat/fftw3
btw: there are other
fftw3 libs on anaconda that don't have fortran bindings and therefore don't work ...
Comments are welcome. The idea is to start with the shtools functions that already exist, and then to add more complex things (like different ellipsoids, etc.)
Cool! I'm looking forward for this one. I don't have a lot of time this summer, but will try to participate as much as I can. I've written some code for the normal gravity here: https://github.com/ioshchepkov/pygeoid/blob/master/pygeoid/reduction/ellipsoid.py, so it can be implemented in SHTOOLS too, I believe.
It looks like I could probably use a lot of that. One problem with my fortran code is that I only give the option to calculate gravity on a sphere, or a flattened ellipsoid. I'd like to give the option of choosing a different ellipsoid, but that would take a bit more work. Let me get the classes up an running and then I'll see how much time I have left over...
It should be in
/usr/lib/x86_64-linux-gnu; https://packages.ubuntu.com/trusty/amd64/libfftw3-dev/filelist
But numpy should be able to figure this out from the pkg-config file, and that directory should be in the default library search path.
As part of my work on implementing gravity and magnetic field classes, I am also going to redo the planetary constants submodule. The idea is to use the astropy
Constants class that allows you to add units, uncertainties, and much more. If you want me to add anything in specific, let me know here: SHTOOLS/SHTOOLS#156
Here's a notebook describing how it works: https://github.com/MarkWieczorek/SHTOOLS/blob/gravmag/examples/notebooks/Introduction-3.ipynb
I might move the tests to pytest at some point soon Mark
Then I was going to think about make a class SHSlepianCoeffs, for the coefficients of a function expanded using a local Slepian basis.
Here is a new PR that mostly concentrates on make better and publication quality figures. The function
pyshtools.utils.figstyle() sets matplotlib parameters, and also allows to set the maximum width of usable area of a journal page, the relative width of the figure with respect to the maximum size, and the screen resolution in dpi so that matplotlib "inches" really correspond to inches. The plotting routines all have optional parameters for setting tick intervals, font sizes, etc.
I will make a new release in about a week that combines the grav/mag classes and better figures.
Next up will be gmt figures, and decomposition of a function in Slepian functions.
release is set for friday.
Minor bugs in Tutorial 1 (https://github.com/SHTOOLS/SHTOOLS/blob/master/examples/notebooks/tutorial_1.ipynb): x and y labels of the first two power spectrum plots should be reversed.
I've added two new classes for dealing with Slepian basis functions, and the Slepian expansion coefficients of a function:
Slepian and SlepianCoeffs. If I get around to it, I will try to make a Jupyter notebook showing how they work (though they are not too different than SHWindow and SHCoeffs).
I'm fixing a few other minor things and will probably put out a new release before the end of the year.
and, two convenience functions for computing individual legendre and spherical harmonic functions: SHTOOLS/SHTOOLS#170
Thanks for the new release, Mark! Just an idea for a future one. I saw, that now you can get mass from
SHGravCoeffs, I think it would be great if we can calculate all sort of things from Stokes coeffs, i.e. coordinates of the center of mass from the first degree, moments of inertia from others (up to n=2 or 3 at least), etc. It is easy to implement, but useful in planetary studies. I can volunteer for that one, I believe
I just implemented saving SHGrid data to netcdf files. These can be used directly with GMT. SHTOOLS/SHTOOLS#178
Is it possible to retrieve the rotation matrix, i.e. if
sh is a SHCoeffs instance, sh.rotate(**args) gives the rotated coefficients, but is it possible to return the rotation matrix somehow ?
Hi @sekisushai , I think you can use the
djpi2 function for this: https://shtools.github.io/SHTOOLS/python-spherical-harmonic-rotations.html
So, I'm getting ready to put out a new release, but I'm having a small problem getting things to build on travis for linux (everything works fine for macOS). In short: I added a new module called
ftypes.f95, and the module SHTOOLS.f95 makes use or ftypes. The problem is that on Linux, setup.py is trying to compile SHTOOLS.f95 before ftypes.f95, and obviously, that doesn't work.
Does anyone know if there is a way in the
setup.py file to force ftypes.f95 to compile before SHTOOLS.f95?
I'm putting out a new release on monday. Here's the pull request:
SHTOOLS/SHTOOLS#190
There are a bunch of changes to the fortran code and how it calls FFTW, so if anyone wants to check it out beforehand, let me know if you run into any problems.
SHTOOLS/SHTOOLS#190
There are a bunch of changes to the fortran code and how it calls FFTW, so if anyone wants to check it out beforehand, let me know if you run into any problems.