[aspect-devel] 3 composition/instability drip

[Payman Janbakhsh]

Thank you John
I tried to run the prm file you sent, but gives an error message that I have attached. I tried different things but kept getting the same message. Not sure why.

Thanks

payman
From: Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org] On Behalf Of John Naliboff
Sent: October 3, 2016 12:15 PM
To: aspect-devel at geodynamics.org
Subject: Re: [aspect-devel] 3 composition/instability drip

Hi Payman,

Sorry for the delay.  How are things going with the models? If it helps, I’ve attached an input that is designed to replicate the numerical model design in Schmeling et al.  2008.

The mesh resolution is fixed (10 km everywhere), so it will run much slower than models with AMR.  If the other input file is still giving odd behavior, this should be a reasonable starting point for varying AMR parameters, etc.  Behavior looks as I would expect it to after 10 Myr.

Hope this helps and let me know if you have any questions about the input file, etc.

Cheers,
John

*************************************************
John Naliboff
Assistant Project Scientist, CIG
Earth & Planetary Sciences Dept., UC Davis




On Sep 28, 2016, at 1:52 PM, Payman Janbakhsh <payman.janbakhsh at mail.utoronto.ca<mailto:payman.janbakhsh at mail.utoronto.ca>> wrote:

Hi John
I got it working without any error. Right now the lithosphere start to look strange when I add even one adaptive refinement. I'll keep working on it .
Thanks again
Payman



Sent from my Samsung Galaxy smartphone.


-------- Original message --------
From: John Naliboff <jbnaliboff at ucdavis.edu<mailto:jbnaliboff at ucdavis.edu>>
Date: 2016-09-28 12:43 PM (GMT-05:00)
To: aspect-devel at geodynamics.org<mailto:aspect-devel at geodynamics.org>
Subject: Re: [aspect-devel] 3 composition/instability drip

Hi Payman,

Thanks for sending this over.  To further assist, can you point to a specific model setup/publication this is based on (if it is based on one?). Also, can you send over an image of your initial compositional field?

The solver is not happy about something and it might be easier to diagnose with a few images.  I should have time later today or tomorrow morning to take a close look.  Currently playing catchup after traveling.

Cheers,
John



*************************************************

John Naliboff

Assistant Project Scientist, CIG

Earth & Planetary Sciences Dept., UC Davis
On 09/28/2016 04:20 AM, Payman Janbakhsh wrote:
Good morning John
I have attached the prm and error files. I’m using version 1.5.
Hope you can point out the issue.
Many thanks


payman


From: Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org] On Behalf Of John Naliboff
Sent: September 27, 2016 2:51 PM
To: aspect-devel at geodynamics.org<mailto:aspect-devel at geodynamics.org>
Subject: Re: [aspect-devel] 3 composition/instability drip


Hi Payman,


I’m guessing the “scheming_empty.prm” model refers to the Scheling et al 2008 (Phys. Earth. Planet Inter.) subduction benchchmark?


To help diagnose the problem, can you send over the .prm you used and a file with the error messages you received.  If the input file is from 2014, a few things may have changed in ASPECT, but it is good to hear the model still runs.


Also, what version of ASPECT are you currently using (version 1.4, master branch from github, etc).


Material model - I think one the existing models should be able to accommodate the parameters used in the Schmeling et al. 2008 paper, especially the linear viscous case?


Cheers,
John


*************************************************
John Naliboff
Assistant Project Scientist, CIG
Earth & Planetary Sciences Dept., UC Davis









On Sep 27, 2016, at 11:06 AM, Payman Janbakhsh <payman.janbakhsh at mail.utoronto.ca<mailto:payman.janbakhsh at mail.utoronto.ca>> wrote:


Hi everyone,
I’m fairly new to ASPECT and trying to create an instability drip or a subduction model similar to what was used in scheming_empty.prm file in CIDER 2014.
I have mantle, lithosphere and crust each with their own density and viscosity.
I tried just to set up the parameter file with multicomponent material model without the use of a plugin but get errors related to convergence. Not sure if is doable without a plug in! let me know if so.
the virtualbox tutorial for CIDER 2014 which contained the related files is no longer available on wiki except the GeoMod2014 tutorial iv pdf file.
Unfortunately I’m stuck in creating the material model plug-in from the appropriate .cc and .h files.
Would you be able to direct me to the right resources so I can get this model running?


Thank you in advance


Payman


_______________________________________________
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

_______________________________________________
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/20161003/7ceefee9/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Error-schmeling2008.odt
Type: application/octet-stream
Size: 16268 bytes
Desc: Error-schmeling2008.odt
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161003/7ceefee9/attachment-0001.obj>