Guillaume’s animation “seismicity”

[gd-a]

Volcanic seismicity

• [x] create a static mock up
• [ ] polish labels and stuff
• [ ] animate

[gd-a]

I started to work on it in December, but halted for 2 reasons :

1. I wasn’t able to have a nice rotation around the volcano (see here)
2. Even though I have a nice workaround for the PostScript ”spheres” using -Se, the trick collapses for azimuth 0

Also, I position everything manually, which can be ugly code-wise.

• On the left : 3D map of La Soufriere (Guadeloupe) volcano edifice (waiting to be fixed)
• Top right : a more regional eye-bird view (done)
• Bottom right : a profile view with the same regional extent (need some embellishments)

[Esteban82]

• I wasn’t able to have a nice rotation around the volcano (see here)

How you try using +v for-p(as in the animation 03)?

[gd-a]

I haven't. But for demonstration material, it's a little weird to project from lon-lat to x-y while +w exists ... no?

[gd-a]

(for illustration) BTW, I couldn't make the topography semi-transparent (grdview).

[PaulWessel]

BTW, I couldn't make the topography semi-transparent (grdview).

Using -t would make the entire grdview layer transparent. Would that work?

[PaulWessel]

Sorry this is so slow. Might you be able to share the script/data that makes this and that chokes with -p+w ?

[gd-a]

Ouch, sorry yeah, I'll try to do come back to it this week. New job... so I expect to be worse than slow from now on.

[gd-a]

Archive.zip

It "was" work in progress ... From what I remember, -t was buggy too.

[Esteban82]

@gd-a is this file DATA/2014-2019_catalog.csv on github? If not, could you add it so I can help with your script.

[gd-a]

done

[Esteban82]

@gd-a I run your script and I got this figure:

I will work on the comments here that I copy below:

• Control the rotation of the left part with [-p+wlon0/lat0] (https://github.com/GenericMappingTools/gmt/issues/7066)
• avoid the whitespace between the topography and the "sea area" (coast not aligned with grdview)
• make grdview more transparent
• project the hypocenters on the top right figure (profile from track)

[gd-a]

Dude! You’ll make me feel bad if you’re doing the work that I should be doing :o

[Esteban82]

Uhh ok sorry. Let me know if do you need help with something.

[gd-a]

Well… depending on @PaulWessel, either the -p problem is fixed, or I completely switch to a 2x2 subplot with map (recycled), time-series(new), cross section (recycled but incomplete) and heatmap (recycled from forum).

but of you find a way, I’ll take your solution :)

[Esteban82]

Hi @gd-a !

I am checking this script. When I run it, I only get an image but I didn't got an animation. Do you get the same? Do you have a finish version of the script?

[gd-a]

It is "work in progress". I wanted to settle on the visual (static) first before going to animation.

I had several issues :

• Do I go for the 3D view despite the perspective bug ?
• Should I switch to a 2x2 subplot with a map, some time-series, two depth-distance projections with orthogonal angles
• How do I do this depth-distance projection (the angle is obvious, one orthogonal to the subduction, the second parallel)?
• How do I add the slab profile in the depth-distance projection ?

(-X,-Y) [ 3 [map] D [profile] view] [time-series]

or

(subplot) [map] [time-series] [profile X] [profile Y]