[aspect-devel] Some helps about ASPECT

[Rene Gassmoeller]

Hi Elisa,
Your task certainly sounds doable with ASPECT, however you might need to
do some slight adjustments to the code. As John mentioned the 'Initial
condition/Ascii data' plugin is a powerful and flexible tool to
prescribe arbitrary temperature fields as initial condition (a very
simple example model is in tests/ascii_data_initial_condition_2d_box.prm
and the corresponding data is in
data/initial-conditions/ascii-data/test/box_2d.txt). To use this plugin
however, you will need to provide a temperature field for the whole
domain (not just the thermal boundary age), so you would need to prepare
a more complete file.
A relatively new addition to the code is the adiabatic-boundary plugin
(thanks to Tahiry Rajaonarison and D. Sarah Stamps for contributing),
which seems closer to your problem, but is intended for a ellipsoid
geometry - you will need to check if this fits your application. It
prescribes a depth (not an age) for the lithosphere-asthenosphere
temperature profile transition in a data file. A simple test model can
be found in tests/adiabatic_boundary.prm and the corresponding data is
in data/initial-conditions/adiabatic-boundary/adiabatic-boundary.txt.
Have a look if this fits your needs, I guess Tahiry can answer more in
depth questions than I can.
In general you find all available plugins in source/initial_conditions/,
or in the 'Parameters in section Initial conditions' section of the
manual. If you find no plugin that fits your need, we can also discuss
how to best add a new plugin using a boundary age, I once ran a model
with a similar setup and just wrote something crude myself without ever
generalizing it, but it seems a general enough application.

Let us know about your progress,
Best regards,
Rene

On 05/19/2016 11:22 PM, John Naliboff wrote:
> Hi Elisa,
>
> Welcome to ASPECT!
>
>> If I create a /lithosphere age file (a .txt file /with different ages
>> in) and I input it, is ASPECT able to read it? 
>
> Yes, you can certainly setup ASPECT to do this.  I know of one person
> who has done something similar to what you describe.  To get to know
> the code better, it would be a good exercise to work through this.  
>
> As guidance for reading in data for initial conditions, look at the
> section of the manual (and relevant source code) that describes
> reading in temperature boundary data ‘Boundary temperature model/Ascii
> data model’.
>
> Note that you can also create functions defined within the input file
> that can create lateral and depth-dependent temperature field.  See
> ‘Initial conditions/Function’ parameter section in the manual.  You
> could distinguish between continental and oceanic domains in this
> manner as well.
>
> Cheers,
> John
>
> *************************************************
> Assistant Project Scientist, CIG
> Earth & Planetary Sciences Dept., UC Davis
>
>
>
>
>
>> On May 19, 2016, at 2:20 AM, ELISA FIERRO <elisa.fierro at unicam.it
>> <mailto:elisa.fierro at unicam.it>> wrote:
>>
>> Hi All,
>>
>> I'm an italian PhD student in Geodynamics and the topic of my PhD is
>> studying mantle dynamics beneath the Red Sea during its various
>> opening phases. To do this I downloaded the latest version of Aspect
>> (1.5.0-pre) (after trying CitcomS) and I'm learning to use it.
>> I need to distinguish the oceanic lithosphere from the continental
>> domain and I would do it defining a thermal boundary layer and
>> different lithosphere ages. 
>> My questions are: 
>> If I create a /lithosphere age file (a .txt file /with different ages
>> in) and I input it, is ASPECT able to read it? 
>> if it so, how can I set up this in the input file? 
>> if it's not, is there a way to distinguish the two lithospheric domains?
>> Thanks!!
>>
>> Elisa Fierro
>> _______________________________________________
>> 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20160520/a2564090/attachment-0001.html>