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

Eleonora Ficini eleonora.ficini at uniroma1.it
Mon Oct 30 02:37:49 PDT 2017


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:

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?

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.*

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)*


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


Thank you in advance.


All the best,

Eleonora

2017-10-26 18:54 GMT+02:00 <aspect-devel-request at geodynamics.org>:

> 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. 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>
> To: aspect-devel at geodynamics.org
> Subject: [aspect-devel] Plate motion reconstrucion - Layered rheology
> Message-ID:
>         <CACT2GZZCwO52XL_Bvq1Cd42GAKgJ0CCJ+qUxTLA_
> 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>
>
> ------------------------------
>
> Message: 2
> Date: Wed, 25 Oct 2017 14:20:55 -0700
> From: Maxwell Rudolph <maxrudolph at ucdavis.edu>
> To: aspect-devel at geodynamics.org
> Subject: Re: [aspect-devel] Plate motion reconstrucion - Layered
>         rheology
> Message-ID:
>         <CA+pU5678ZkCzKi9VKBLeT+4=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> 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
> > 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>
>
> ------------------------------
>
> Message: 3
> Date: Wed, 25 Oct 2017 23:33:07 +0000
> From: "Glerum, A.C. (Anne)" <A.C.Glerum at uu.nl>
> To: "aspect-devel at geodynamics.org" <aspect-devel at geodynamics.org>
> Subject: Re: [aspect-devel] Plate motion reconstrucion - Layered
>         rheology
> Message-ID: <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>> 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>> 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
>
> _______________________________________________
> 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
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <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>
> To: aspect-devel at geodynamics.org
> Subject: Re: [aspect-devel] Plate motion reconstrucion - Layered
>         rheology
> Message-ID: <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>> 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>>
> 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
> >
> >>
> >>
> >> _______________________________________________
> >> 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
> >
> >
> >
> > _______________________________________________
> > 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/20171026/21377303/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 71, Issue 12
> ********************************************
>



-- 
Eleonora Ficini, PhD student
Department of Earth Sciences
Sapienza University of Rome
P.le A. Moro, 5 - 00185, Rome, Italy
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20171030/7084cfd8/attachment-0001.html>


More information about the Aspect-devel mailing list