Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • Jan 29 09:17
    mjziebarth starred SHTOOLS/SHTOOLS
  • Jan 23 16:48
    WX-78 starred SHTOOLS/SHTOOLS
  • Jan 23 02:27
    alameday starred SHTOOLS/SHTOOLS
  • Jan 02 11:11
    ashwinvis commented #147
  • Jan 02 11:10
    ashwinvis starred SHTOOLS/SHTOOLS
  • Dec 18 2018 20:18
    ajacquey starred SHTOOLS/SHTOOLS
  • Dec 17 2018 16:19

    MarkWieczorek on develop

    Start development branch for Ve… Merge pull request #175 from Ma… (compare)

  • Dec 17 2018 16:19
    MarkWieczorek closed #175
  • Dec 17 2018 11:34
    MarkWieczorek opened #175
  • Dec 17 2018 11:14

    MarkWieczorek on v4.4.1

    (compare)

  • Dec 17 2018 11:13

    MarkWieczorek on master

    V4.4.1 (compare)

  • Dec 17 2018 10:33

    MarkWieczorek on www

    Update web documentation for v4… Update web documentation for v4… (compare)

  • Dec 17 2018 10:25

    MarkWieczorek on v4.4

    (compare)

  • Dec 17 2018 10:21

    MarkWieczorek on master

    Set version number to 4.4 on de… Change plt.show() to fig.show() update references and 33 more (compare)

  • Dec 17 2018 10:21
    MarkWieczorek closed #173
  • Dec 17 2018 09:59
    MarkWieczorek synchronize #173
  • Dec 17 2018 09:59

    MarkWieczorek on develop

    Update Gemfile.lock Merge pull request #174 from Ma… (compare)

  • Dec 17 2018 09:59
    MarkWieczorek closed #174
  • Dec 17 2018 09:37
    MarkWieczorek opened #174
  • Dec 12 2018 16:38
    MarkWieczorek synchronize #171
MMesch
@MMesch
I installed with anaconda and used their fftw
Mark Wieczorek
@MarkWieczorek
That doesn't work in Travis either: https://travis-ci.org/SHTOOLS/SHTOOLS/jobs/388356217
Mark Wieczorek
@MarkWieczorek
@MMesch actually, this appears to work: conda install -n test-environment -c eumetsat fftw3. I forgot to install in into the test environment. Thanks.
MMesch
@MMesch
great!
btw: there are other fftw3 libs on anaconda that don't have fortran bindings and therefore don't work ...
Mark Wieczorek
@MarkWieczorek
well, this one works , for the time being.... I'll try to figure out what is going on get-apt. We are finding the blas and lapack libraries ok, so maybe something changed with libfftw3-dev.
Mark Wieczorek
@MarkWieczorek
My summer project will be to make new classes for working with gravity and magnetic field data: SHTOOLS/SHTOOLS#153
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.)
Ilya Oshchepkov
@ioshchepkov
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.
Mark Wieczorek
@MarkWieczorek
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...
Elliott Sales de Andrade
@QuLogic
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.
Mark Wieczorek
@MarkWieczorek
@QuLogic My guess is that Travis isn't setting the paths correctly. If other people run into the same problem, I'll look into this more...
Mark Wieczorek
@MarkWieczorek
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
Mark Wieczorek
@MarkWieczorek
The constant submodule is now done. Let me know if you have any comments: SHTOOLS/SHTOOLS#157
Mark Wieczorek
@MarkWieczorek
The gravity classes are done: SHTOOLS/SHTOOLS#162
MMesch
@MMesch
nice
I might move the tests to pytest at some point soon Mark
Mark Wieczorek
@MarkWieczorek
After I finalize this, I'll make a copy of this class for magnetic field data. I don't think subclassing is going to work because all the documentation needs to be changed, the plotting labels are different, and the methods all are slightly different.
Then I was going to think about make a class SHSlepianCoeffs, for the coefficients of a function expanded using a local Slepian basis.
Mark Wieczorek
@MarkWieczorek
The magnetic field classes are now done. If there aren't any comments, I am going to merge this all to develop tomorrow.
Mark Wieczorek
@MarkWieczorek
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.
Mark Wieczorek
@MarkWieczorek
Verison 4.3 pull request: SHTOOLS/SHTOOLS#165
release is set for friday.
cryocuhk
@cryocuhk
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.
MMesch
@MMesch
thanks for reporting!
Mark Wieczorek
@MarkWieczorek
Thanks, I just fixed this on my local branch. It should be fixed in the next release at the end of the year....
Mark Wieczorek
@MarkWieczorek
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.
Mark Wieczorek
@MarkWieczorek
weighted least squares spherical harmonic expansion: SHTOOLS/SHTOOLS#171
and, two convenience functions for computing individual legendre and spherical harmonic functions: SHTOOLS/SHTOOLS#170
Mark Wieczorek
@MarkWieczorek
MMesch
@MMesch
:+1:
Ilya Oshchepkov
@ioshchepkov
Thanks for the new release, Mark! Just an idea for a future one. I saw, that now you can get mass fromSHGravCoeffs, 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
Mark Wieczorek
@MarkWieczorek
sounds like a good idea. If you want to help out, that would be great. Now that I got this release out of the way, I need to get some real work done!
Mark Wieczorek
@MarkWieczorek
I just implemented saving SHGrid data to netcdf files. These can be used directly with GMT. SHTOOLS/SHTOOLS#178
sekisushai
@sekisushai
Hello, I've got a question regarding rotation matrix of SH coefficients:
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 ?
MMesch
@MMesch
Hi @sekisushai , I think you can use the djpi2 function for this: https://shtools.github.io/SHTOOLS/python-spherical-harmonic-rotations.html
sekisushai
@sekisushai
@MMesch thank you
Mark Wieczorek
@MarkWieczorek
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?
Elliott Sales de Andrade
@QuLogic
os.listdir is not deterministic, probably just sort the result, or don't glob
Mark Wieczorek
@MarkWieczorek
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.