[aspect-devel] lateral variations in viscosity & present-day seismic image as inputs

[Austermann, Jacqueline]

Hi Petar,

Yes, both is possible (Wolfgang and Timo correct me if I'm wrong). For the initial thermal structure you could take a look at S40RTS_perturbation.cc, which inputs a scaled S40RTS initial condition. This should give you an idea of how you could read in tomography data and use them as your initial condition. For the material model, you could look at the steinberger.cc material model on how to read in information on the viscosity.
I've actually started working on a material model that I call mitrovica_forte.cc, which is likely close to some of the things you are looking for and inbetween the simple model and the steinberger model in terms of complexity. I'm still testing it so it's not in aspect yet but if you want we could work on this together (maybe use it for a benchmark with your code?) and eventually push it to aspect. Feel free to e-mail me directly if you want to do that (jaustermann at fas.harvard.edu<mailto:jaustermann at fas.harvard.edu>) or keep it on the aspect-devel thread.

Cheers,
Jacky



On Oct 1, 2014, at 12:20 PM, Petar Glisovic wrote:

Hi Wolfgang & Timo,

I would like to run an instantaneous flow calculation in 3-D spherical shell using
1. initial thermal structure that includes a present-day seismic image, and
2. lateral variations in viscosity as a 3-D input, i.e., I do not want to use viscosity prefactors (set Lateral viscosity file name option).

Is that possible?

Regards,
Petar

--
Petar Glisovic, Dr.
postdoctoral fellow

address:  GEOTOP, Université du Québec à Montréal, CP 8888, succursale Centre-Ville, Montréal, Québec, Canada H3C 3P8  | email:  pglisovic at gmail.com<mailto:pglisovic at gmail.com> & glisovic.petar at courrier.uqam.ca<mailto:glisovic.petar at courrier.uqam.ca>  | phone:  1-514-987-3000 #3592 & #4600  | fax:  1-514-987-3635
_______________________________________________
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/20141001/b4127489/attachment.html>