[aspect-devel] Aspect-devel Digest, Vol 71, Issue 12 - Plate motion reconstrucion - Layered rheology -

Rene Gassmoeller rene.gassmoeller at mailbox.org
Mon Oct 30 15:48:58 PDT 2017


Hi Eleonora,


On 10/30/2017 03:37 AM, Eleonora Ficini wrote:
> Dear all,
> thank you for your answers and suggestions.
> What I wanted to say in the previews mail (sorry if I didn't explain 
> it in a correct way, it probably sounded very rude!) was just that I 
> had a number of error messages in learning the first steps of how 
> include GPlates reconstruction in the ASPECT code:
>

Do not worry, we all started at some point, questions are normal :-)

> When I try to run the GPlates cookbook example what I get is:
>
> / Loading GPlates data boundary file 
> /home/ubuntu/aspect/data/velocity-boundary-conditions/gplates/current_day.gpml./
> /
> /
> /   Loading new velocity file did not succeed./
> /   Assuming constant boundary conditions for rest of model run./
>
> Does this mean that ASPECT use the current_day velocities for the 
> entire simulation? Or there's an error occuring when loading the 
> velocity file?
>

This message means ASPECT loaded a single file, and noticed the file is 
not part of a series of data files. Thus it will use this single file 
for the whole computation. This is not an error in itself. The example 
input file simply prescribes to use the same velocity file for the whole 
simulation time. I admit that the warning message could be better 
formulated. If you want to use time-dependent velocities, you will have 
to create a dataset in GPlates and use that one, because we do not 
include such large datasets inside of ASPECT (the process for this is 
described in the "Time-dependent boundary conditions" paragraph of the 
GPlates section in the manual).

> I tried then to look for the velocity file (to see how it was built) 
> but I wasn't able to find the source file, neither from the ASPECT 
> link or from the websites mentioned in the manual (in particular, the 
> "Global reconstructions with continuously closing plates from 140 Ma 
> to the present"). The error I got was:
>
> /404 Not Found The requested URL /~gurnis/GPlates/gplates.html was not 
> found on this server./

Thanks for this hint, apparently the website has moved. I have opened an 
issue here: https://github.com/geodynamics/aspect/issues/1979. If we 
find a solution it will be posted there.

> /
> /
> Lastly, with the manual instructions (that also Anne suggested), I 
> tried to build a new velocity file by myself using Seton et al., 2012 
> reconstruction (that comes with GPlates software 1.5.0)  but what I 
> got was:
>
> /An error occurred while saving the file 
> 'lat_lon_velocity_domain_9_18.gpml': /
>
> /ErrorOpeningFileForWritingException: Error opening file 
> 'lat_lon_velocity_domain_9_18.gpml' for writing/
>
> /Call stack trace:/
>
> /(/Users/jcannon/gplates/builds/1.5/src/file-io/GpmlOutputVisitor.cc, 
> 551)/
>

Hmm, this looks like a problem inside of GPlates. The Seton 2012 dataset 
should be fine to use (I used that before). Which version of Gplates did 
you use? And it seems GPlates could not write the file to your disk, are 
you sure you tried to save it in a location you have write-access to?

> /
> /
>
> I am sorry if any of the mistakes I am making is basic, but I just 
> started to deal with the code.
>
>

Dont worry, some of your points are genuine problems, and others are 
normal questions, that is what this mailing list is for.

Best,
Rene

> Thank you in advance.
>
>
> All the best,
>
> Eleonora
>
>
> 2017-10-26 18:54 GMT+02:00 <aspect-devel-request at geodynamics.org 
> <mailto:aspect-devel-request at geodynamics.org>>:
>
>     Send Aspect-devel mailing list submissions to
>     aspect-devel at geodynamics.org <mailto: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
>     <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
>     <mailto:aspect-devel-request at geodynamics.org>
>
>     You can reach the person managing the list at
>     aspect-devel-owner at geodynamics.org
>     <mailto: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. Plate motion reconstrucion - Layered rheology (Eleonora Ficini)
>        2. Re: Plate motion reconstrucion - Layered rheology
>           (Maxwell Rudolph)
>        3. Re: Plate motion reconstrucion - Layered rheology
>           (Glerum, A.C. (Anne))
>        4. Re: Plate motion reconstrucion - Layered rheology
>           (Rene Gassmoeller)
>
>
>     ----------------------------------------------------------------------
>
>     Message: 1
>     Date: Wed, 25 Oct 2017 15:41:46 +0200
>     From: Eleonora Ficini <eleonora.ficini at uniroma1.it
>     <mailto:eleonora.ficini at uniroma1.it>>
>     To: aspect-devel at geodynamics.org <mailto:aspect-devel at geodynamics.org>
>     Subject: [aspect-devel] Plate motion reconstrucion - Layered rheology
>     Message-ID:
>            
>     <CACT2GZZCwO52XL_Bvq1Cd42GAKgJ0CCJ+qUxTLA_eGuHUxAHqw at mail.gmail.com
>     <mailto:CACT2GZZCwO52XL_Bvq1Cd42GAKgJ0CCJ%2BqUxTLA_eGuHUxAHqw at mail.gmail.com>>
>     Content-Type: text/plain; charset="utf-8"
>
>     Hello everyone,
>     I am a PhD student, I recently started to use ASPECT and I have two
>     different questions:
>
>     1) I would like to know if someone tried to include plate motion
>     reconstruction (from GPlates) into ASPECT. If yes, what kind of
>     files does
>     ASPECT need? It seems that the manual doesn't mention the correct
>     way to
>     built the velocity/time file.
>
>     2) In a global spherical model, how is the layered rheological
>     structure
>     (especially regarding a 100 km thick lithosphere)?
>
>     Thank you in advance for your time and attention.
>
>     All the best,
>     Eleonora
>     -------------- next part --------------
>     An HTML attachment was scrubbed...
>     URL:
>     <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171025/43f26a72/attachment-0001.html
>     <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171025/43f26a72/attachment-0001.html>>
>
>     ------------------------------
>
>     Message: 2
>     Date: Wed, 25 Oct 2017 14:20:55 -0700
>     From: Maxwell Rudolph <maxrudolph at ucdavis.edu
>     <mailto:maxrudolph at ucdavis.edu>>
>     To: aspect-devel at geodynamics.org <mailto:aspect-devel at geodynamics.org>
>     Subject: Re: [aspect-devel] Plate motion reconstrucion - Layered
>             rheology
>     Message-ID:
>            
>     <CA+pU5678ZkCzKi9VKBLeT+4=E3OVK9oO5c31LgmKFRkATjwoeQ at mail.gmail.com
>     <mailto:E3OVK9oO5c31LgmKFRkATjwoeQ at mail.gmail.com>>
>     Content-Type: text/plain; charset="utf-8"
>
>     Eleonora,
>     There is a 'depth dependent' material model that allows you to
>     prescribe a
>     depth-dependent viscosity prefactor to any of the other material
>     models.
>
>     Max
>
>     On Wed, Oct 25, 2017 at 6:41 AM, Eleonora Ficini <
>     eleonora.ficini at uniroma1.it <mailto:eleonora.ficini at uniroma1.it>>
>     wrote:
>
>     > Hello everyone,
>     > I am a PhD student, I recently started to use ASPECT and I have two
>     > different questions:
>     >
>     > 1) I would like to know if someone tried to include plate motion
>     > reconstruction (from GPlates) into ASPECT. If yes, what kind of
>     files does
>     > ASPECT need? It seems that the manual doesn't mention the
>     correct way to
>     > built the velocity/time file.
>     >
>     > 2) In a global spherical model, how is the layered rheological
>     structure
>     > (especially regarding a 100 km thick lithosphere)?
>     >
>     > Thank you in advance for your time and attention.
>     >
>     > All the best,
>     > Eleonora
>     >
>     >
>     > _______________________________________________
>     > Aspect-devel mailing list
>     > Aspect-devel at geodynamics.org <mailto:Aspect-devel at geodynamics.org>
>     >
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <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/20171025/7507efa7/attachment-0001.html
>     <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171025/7507efa7/attachment-0001.html>>
>
>     ------------------------------
>
>     Message: 3
>     Date: Wed, 25 Oct 2017 23:33:07 +0000
>     From: "Glerum, A.C. (Anne)" <A.C.Glerum at uu.nl
>     <mailto:A.C.Glerum at uu.nl>>
>     To: "aspect-devel at geodynamics.org
>     <mailto:aspect-devel at geodynamics.org>"
>     <aspect-devel at geodynamics.org <mailto:aspect-devel at geodynamics.org>>
>     Subject: Re: [aspect-devel] Plate motion reconstrucion - Layered
>             rheology
>     Message-ID: <EC9309A6-D0E9-4CD6-A475-415B65A8377B at uu.nl
>     <mailto:EC9309A6-D0E9-4CD6-A475-415B65A8377B at uu.nl>>
>     Content-Type: text/plain; charset="utf-8"
>
>     Hi Eleonora,
>
>     GPlates files can be used as boundary conditions in ASPECT, see
>     for example Heister et al. 2017 GJI.
>
>     There is a section in the manual that describes how to create the
>     GPlates velocity files:
>     Section 5.3.4 Using reconstructed surface velocities by GPlates
>
>      What I did recently to create a file for present day velocities was:
>
>       1.  (1)  Load the Reconstruction Tree
>
>       2.  (2)  Load the Dynamic Polygons
>
>       3.  (3)  Generate a grid for the velocity output with Features →
>     Generate Velocity Domain Points
>
>       4.  (4)  This will show velocity vectors for the grid points.
>     Export through Reconstruction → Export
>
>     I hope that gets you started,
>     Anne
>
>
>
>     On 26 Oct 2017, at 08:20, Maxwell Rudolph <maxrudolph at ucdavis.edu
>     <mailto:maxrudolph at ucdavis.edu><mailto:maxrudolph at ucdavis.edu
>     <mailto:maxrudolph at ucdavis.edu>>> wrote:
>
>     Eleonora,
>     There is a 'depth dependent' material model that allows you to
>     prescribe a depth-dependent viscosity prefactor to any of the
>     other material models.
>
>     Max
>
>     On Wed, Oct 25, 2017 at 6:41 AM, Eleonora Ficini
>     <eleonora.ficini at uniroma1.it
>     <mailto:eleonora.ficini at uniroma1.it><mailto:eleonora.ficini at uniroma1.it
>     <mailto:eleonora.ficini at uniroma1.it>>> wrote:
>     Hello everyone,
>     I am a PhD student, I recently started to use ASPECT and I have
>     two different questions:
>
>     1) I would like to know if someone tried to include plate motion
>     reconstruction (from GPlates) into ASPECT. If yes, what kind of
>     files does ASPECT need? It seems that the manual doesn't mention
>     the correct way to built the velocity/time file.
>
>     2) In a global spherical model, how is the layered rheological
>     structure (especially regarding a 100 km thick lithosphere)?
>
>     Thank you in advance for your time and attention.
>
>     All the best,
>     Eleonora
>
>
>     _______________________________________________
>     Aspect-devel mailing list
>     Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org><mailto:Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org>>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>
>
>     _______________________________________________
>     Aspect-devel mailing list
>     Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org><mailto:Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org>>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <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/20171025/bf6742ba/attachment-0001.html
>     <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171025/bf6742ba/attachment-0001.html>>
>
>     ------------------------------
>
>     Message: 4
>     Date: Thu, 26 Oct 2017 10:58:54 -0600
>     From: Rene Gassmoeller <rene.gassmoeller at mailbox.org
>     <mailto:rene.gassmoeller at mailbox.org>>
>     To: aspect-devel at geodynamics.org <mailto:aspect-devel at geodynamics.org>
>     Subject: Re: [aspect-devel] Plate motion reconstrucion - Layered
>             rheology
>     Message-ID: <4f51844b-6959-7098-df33-99722b105448 at mailbox.org
>     <mailto:4f51844b-6959-7098-df33-99722b105448 at mailbox.org>>
>     Content-Type: text/plain; charset="utf-8"; Format="flowed"
>
>     Hi Eleonora,
>
>     you stated that
>
>     > It seems that the manual doesn't mention the correct way to
>     built the
>     > velocity/time file.
>
>     Could you elaborate on that? What is the error message that you see if
>     you follow the instructions that Anne send?
>
>     If this continues to fail that might mean GPlates has again
>     changed its
>     data format (they did that once in the past, and we currently support
>     both formats). If Anne's suggestions do not help, could you send me
>     which version of ASPECT and GPlates you are using (and ideally
>     also the
>     plate reconstruction you used if that is possible, though better send
>     that in a private email off the mailing list)? Then I can try to
>     reproduce the issue and look for a fix.
>
>     Best,
>
>     Rene
>
>
>
>     On 10/25/2017 05:33 PM, Glerum, A.C. (Anne) wrote:
>     > Hi Eleonora,
>     >
>     > GPlates files can be used as boundary conditions in ASPECT, see for
>     > example Heister et al. 2017 GJI.
>     >
>     > There is a section in the manual that describes how to create the
>     > GPlates velocity files:
>     > Section 5.3.4 Using reconstructed surface velocities by GPlates
>     >
>     >  What I did recently to create a file for present day velocities
>     was:
>     >
>     > 1.
>     >
>     >     (1)  Load the Reconstruction Tree
>     >
>     > 2.
>     >
>     >     (2)  Load the Dynamic Polygons
>     >
>     > 3.
>     >
>     >     (3)  Generate a grid for the velocity output with Features →
>     >     Generate Velocity Domain Points
>     >
>     > 4.
>     >
>     >     (4)  This will show velocity vectors for the grid points. Export
>     >     through Reconstruction → Export
>     >
>     > I hope that gets you started,
>     > Anne
>     >
>     >
>     >
>     >> On 26 Oct 2017, at 08:20, Maxwell Rudolph
>     <maxrudolph at ucdavis.edu <mailto:maxrudolph at ucdavis.edu>
>     >> <mailto:maxrudolph at ucdavis.edu
>     <mailto:maxrudolph at ucdavis.edu>>> wrote:
>     >>
>     >> Eleonora,
>     >> There is a 'depth dependent' material model that allows you to
>     >> prescribe a depth-dependent viscosity prefactor to any of the other
>     >> material models.
>     >>
>     >> Max
>     >>
>     >> On Wed, Oct 25, 2017 at 6:41 AM, Eleonora Ficini
>     >> <eleonora.ficini at uniroma1.it
>     <mailto:eleonora.ficini at uniroma1.it>
>     <mailto:eleonora.ficini at uniroma1.it
>     <mailto:eleonora.ficini at uniroma1.it>>> wrote:
>     >>
>     >>     Hello everyone,
>     >>     I am a PhD student, I recently started to use ASPECT and I have
>     >>     two different questions:
>     >>
>     >>     1) I would like to know if someone tried to include plate
>     motion
>     >>     reconstruction (from GPlates) into ASPECT. If yes, what kind of
>     >>     files does ASPECT need? It seems that the manual doesn't
>     mention
>     >>     the correct way to built the velocity/time file.
>     >>
>     >>     2) In a global spherical model, how is the layered rheological
>     >>     structure (especially regarding a 100 km thick lithosphere)?
>     >>
>     >>     Thank you in advance for your time and attention.
>     >>
>     >>     All the best,
>     >>     Eleonora
>     >>
>     >>
>     >>     _______________________________________________
>     >>     Aspect-devel mailing list
>     >> Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org>
>     <mailto:Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org>>
>     >>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>
>     >>   
>      <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>>
>     >>
>     >>
>     >> _______________________________________________
>     >> Aspect-devel mailing list
>     >> Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org>
>     <mailto:Aspect-devel at geodynamics.org
>     <mailto:Aspect-devel at geodynamics.org>>
>     >>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>
>     >
>     >
>     >
>     > _______________________________________________
>     > Aspect-devel mailing list
>     > Aspect-devel at geodynamics.org <mailto:Aspect-devel at geodynamics.org>
>     >
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>
>
>     --
>     Rene Gassmoeller
>     http://www.math.colostate.edu/~gassmoel/
>     <http://www.math.colostate.edu/%7Egassmoel/>
>
>     -------------- next part --------------
>     An HTML attachment was scrubbed...
>     URL:
>     <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171026/21377303/attachment.html
>     <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171026/21377303/attachment.html>>
>
>     ------------------------------
>
>     Subject: Digest Footer
>
>     _______________________________________________
>     Aspect-devel mailing list
>     Aspect-devel at geodynamics.org <mailto:Aspect-devel at geodynamics.org>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>
>
>     ------------------------------
>
>     End of Aspect-devel Digest, Vol 71, Issue 12
>     ********************************************
>
>
>
>
> -- 
> Eleonora Ficini, PhD student
> Department of Earth Sciences
> Sapienza University of Rome
> P.le A. Moro, 5 - 00185, Rome, Italy
>
>
>
>
> _______________________________________________
> 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/20171030/95a51917/attachment-0001.html>


More information about the Aspect-devel mailing list