[aspect-devel] natural strain / geoid / ParaView questions (Shangxin)

Tue Nov 14 11:25:20 PST 2017

Hi Ludovic,

The geoid in ASPECT is computed in spherical harmonic domain. The
traditional way to calculate geoid in spherical harmonics, as CitcomS does,
is to do the layered integral of spherical harmonics from CMB to surface
(as equation (22) shows in Zhong et al 2008). However, due to ASPECT's
adaptive mesh refinement features, this way is not appropriate for ASPECT.
So I did some equation transform to enable the integral of spherical
harmonics in AMR. I have not organized my methodology into a
formal technical notes (That's my plan in mind), but the benchmark of geoid
is already done. The geoid.prm file in ASPECT's tests directory also
includes a simple global input single degree-order harmonic perturbation
case.

The dynamic topography part in geoid directly comes from dynamic
topography post processor, which is implemented by Ian Rose using
consistent boundary flux method.

Best,
Shangxin

On Tue, Nov 14, 2017 at 1:40 PM, <aspect-devel-request at geodynamics.org>
wrote:

> Send Aspect-devel mailing list submissions to
>         aspect-devel at geodynamics.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> or, via email, send a message with subject or body 'help' to
>         aspect-devel-request at geodynamics.org
>
> You can reach the person managing the list at
>         aspect-devel-owner at geodynamics.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Aspect-devel digest..."
>
>
> Today's Topics:
>
>    1. Re: natural strain / geoid / ParaView questions (Rene Gassmoeller)
>    2. Re: Fwd: Aspect installation error on OS X El Capitan
>       (Hammond, Noah)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 14 Nov 2017 10:49:19 -0700
> From: Rene Gassmoeller <rene.gassmoeller at mailbox.org>
> To: aspect-devel at geodynamics.org
> Subject: Re: [aspect-devel] natural strain / geoid / ParaView
>         questions
> Message-ID: <cd4b083a-7072-b7b3-3ffa-3a13d1cbbec7 at mailbox.org>
> Content-Type: text/plain; charset="windows-1252"
>
> Hi Ludovic,
>
>
> About 1: Currently only the deformation gradient tensor F (in the
> language of the manual) is tracked and output by ASPECT. Everything else
> in that figure I created in VisIt (which has the nice option to compute
> eigenvectors of tensors). I am sure this is somehow possible in Paraview
> as well (probably using the python programmable filter?), but for this
> application Visit seemed easier (let me know if you want more advice on
> how to do it). On the other hand it might make sense to create a
> postprocessor for ASPECT that does the same thing internally, if you are
> interested we can talk about how to implement such a postprocessor.
>
>
> 2. Shangxin Liu would be the best contact for the geoid as he imlemented
> the methodology. He would know if it is already published, or which
> references are appropriate. The dynamic topography methodology has a
> reference in the manual section, I am not sure if there are publications
> by Jacky, Ian, or Shangxin (the authors of that section) that describe
> it in more detail than the general reference.
>
>
> 3. Timo already said it, but unfortunately the support for cartographic
> projections in Paraview is not really existent. If you calculate
> spherical coordinates (e.g. with Timo's script, or by setting up
> calculators that do that), you can create a programmable filter with my
> following script to generate a (crude) Mollweide projection for a global
> model (not sure that is what you want?):
>
>
> theta =inputs[0].PointData['Latitude']/1.01
>
> phi =inputs[0].PointData['Longitude']
>
> bigt =theta
>
> fori inrange(100):
>
> bigt =bigt -(2*bigt +sin(2*bigt) -3.1415*sin(theta))/(2+2*cos(2*bigt))
>
> x =2*sqrt(2)*phi*cos(bigt)/3.1415
>
> y =sqrt(2)*sin(bigt)
>
> output.PointData.append(x, "x")
>
> output.PointData.append(y, "y")
>
> output.PointData.append(bigt, "bigt")
>
> An alternative approach can be to preprocess your data in Paraview, write
> it into text files ("Save Data"),
> and then use GMT to create a map. Cumbersome. Again, a map projection of
> output data would
> be a nice to have feature, let us know if you are interested in adding
> something.
>
> Best,
> Rene
>
>
>
> On 11/13/2017 10:41 AM, Jeanniot, L. (Ludovic) wrote:
> > Hi!
> >
> > I am Ludovic Jeanniot and I do a postdoc with Cedric Thieulot.
> > I am currently learning ASPECT and there are several points where I
> > need some enlightenments, please:
> >
> > *1. "Tracking finite strain":  *
> > I got confused with Figure 31b. Does the natural strain distribution
> > should be a postprocess output of ASPECT or is it later derived from
> > the strain in ParaView?
> > Also it would be great to get to know how to add the eigenvectors.
> > *
> > 2. "Postprocessing spherical 3D convection":*
> > Is there any references about the methodology and calculation of the
> > geoid anomalies (and possibly with it the dynamic topography)?
> >
> > *3. Cartographic projections:*
> > How is it possible to represent the different cartographic projections
> > in ParaView? Does the unofficial Evenden (2008) library work or is it
> > too old or maybe someone else has another solution?
> >
> > Thanks!
> > Ludovic
> >
> >
> > _______________________________________________
> > Aspect-devel mailing list
> > Aspect-devel at geodynamics.org
> > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>
> --
> Rene Gassmoeller
> http://www.math.colostate.edu/~gassmoel/
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.geodynamics.org/pipermail/aspect-devel/
> attachments/20171114/ec061278/attachment-0001.html>
>
> ------------------------------
>
> Message: 2
> Date: Tue, 14 Nov 2017 13:44:24 -0500
> From: "Hammond, Noah" <noah_hammond at alumni.brown.edu>
> To: aspect-devel <aspect-devel at geodynamics.org>
> Subject: Re: [aspect-devel] Fwd: Aspect installation error on OS X El
>         Capitan
> Message-ID:
>         <CAL43Gjrm6TQABkgM-FLVfoy0REaM5572x5XoXcgZRUTorS0
> 28Q at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> After updating to 10.13, I am trying to use candi to install deal-ii. It
> fails with an error when building petsc 3.7.6:
>
> [ Configure does not support Python 3 yet, please run as:
> python2 './configure'
> '--prefix=/Users/noahhammond/deal.ii-candi/petsc-3.7.6'
> '--with-debugging=0' '--with-shared-libraries=1' '--with-mpi=1'
> '--with-x=0' '--with-64-bit-indices=0' '--download-hypre=1' 'CC=mpicc'
> 'CXX=mpicxx' 'FC=mpif90'
> '--with-parmetis-dir=/Users/noahhammond/deal.ii-candi/parmetis-4.0.3'
> '--with-metis-dir=/Users/noahhammond/deal.ii-candi/parmetis-4.0.3'
> '--download-scalapack=1' '--download-mumps=1']
>
> After running this command and make, petsc installs but with some warnings.
> However, when I rerun candi.sh, it does not skip over the petsc build, like
> it does with the other packages that it previously installed.
>
> Any suggestions? Thanks for everyone's help.
> -Noah Hammond
>
> On Fri, Nov 10, 2017 at 10:46 AM, Maxwell Rudolph <maxrudolph at ucdavis.edu>
> wrote:
>
> > I compiled and installed it with zero problems on 10.13.
> >
> > On Nov 10, 2017, at 5:26 AM, Timo Heister <heister at clemson.edu> wrote:
> >
> > >> I will try updating my mac to 10.12 and see how it goes. Will Deal.ii
> > 8.5
> > >> work on the latest release, that is High Seirra 10.13.1?
> > >
> > > 10.13 is still quite new so we don't have good experience with it yet.
> > > It is likely that you have to a) wait until we make a specific
> > > installer, b) switch to deal.II 9.0.0 (maybe a prerelease), or c) (I
> > > have heard about some success) need to install the 10.12 OSX SDK in
> > > xcode.
> > >
> > > --
> > > Timo Heister
> > > http://www.math.clemson.edu/~heister/
> > > _______________________________________________
> > > Aspect-devel mailing list
> > > Aspect-devel at geodynamics.org
> > > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.geodynamics.org/pipermail/aspect-devel/
> attachments/20171114/527f2389/attachment.html>
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>
> ------------------------------
>
> End of Aspect-devel Digest, Vol 72, Issue 10
> ********************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171114/d1530126/attachment-0001.html>