[aspect-devel] Aspect-devel Digest, Vol 71, Issue 14 - Plate reconstruction
Eleonora Ficini
eleonora.ficini at uniroma1.it
Thu Nov 2 09:16:20 PDT 2017
Dear Rene,
thanks for the answers and the explanations!
I have been successful in creating my velocity files and import them in
ASPECT, finally! But now, when I try to run my velocity files (I am trying
with a 200 Myr reconstruction by Seton et al., 2012), ASPECT reads
correctly the first two files but then it suddenly jumps to the 200 and 201
files (my files go from 0 to 200 every 1Myr). It assumes that the first
timestep is the entire time of the computation (for instance, the End Time,
set to 2e8).
Can this depend from the input velocity files? I might have done something
wrong during the exportation.
Thanks again!
Eleonora
2017-10-30 23:44 GMT+01: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. Re: Aspect-devel Digest, Vol 71, Issue 12 - Plate motion
> reconstrucion - Layered rheology - (Rene Gassmoeller)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 30 Oct 2017 16:48:58 -0600
> From: Rene Gassmoeller <rene.gassmoeller at mailbox.org>
> To: aspect-devel at geodynamics.org
> Subject: Re: [aspect-devel] Aspect-devel Digest, Vol 71, Issue 12 -
> Plate motion reconstrucion - Layered rheology -
> Message-ID: <251df413-0a47-db0d-681f-0d4ad034a1d0 at mailbox.org>
> Content-Type: text/plain; charset="utf-8"; Format="flowed"
>
> 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/gplate
> s/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 geodynami
> cs.org>
> > Subject: [aspect-devel] Plate motion reconstrucion - Layered rheology
> > Message-ID:
> >
> > <CACT2GZZCwO52XL_Bvq1Cd42GAKgJ0CCJ+qUxTLA_eGuHUxAHqw at mail.gmail.com
> > <mailto:CACT2GZZCwO52XL_Bvq1Cd42GAKgJ0CCJ%2BqUxTLA_eGuHUxAH
> qw 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/att
> achments/20171025/43f26a72/attachment-0001.html
> > <http://lists.geodynamics.org/pipermail/aspect-devel/att
> achments/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 geodynami
> cs.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/att
> achments/20171025/7507efa7/attachment-0001.html
> > <http://lists.geodynamics.org/pipermail/aspect-devel/att
> achments/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
> @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@
> 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@
> 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/att
> achments/20171025/bf6742ba/attachment-0001.html
> > <http://lists.geodynamics.org/pipermail/aspect-devel/att
> achments/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 geodynami
> cs.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/att
> achments/20171026/21377303/attachment.html
> > <http://lists.geodynamics.org/pipermail/aspect-devel/att
> achments/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/attachm
> ents/20171030/95a51917/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 14
> ********************************************
>
--
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/20171102/6a2da353/attachment-0001.html>
More information about the Aspect-devel
mailing list