search

underworldcode / UWGeodynamics

Underworld Geodynamics
Other
81 stars 32 forks source link

Viscosity problem when replicating Tutorial_1 in Underworld2 #153

Closed alexstruc closed 4 years ago

alexstruc commented 4 years ago

Dear developers: I am trying to replicate the Tutorial_1 of UWGeodynamics in Underworld2. To test the reliability of the model, I input the identical viscosity and plasticity parameters with their units. However, my UW2 model gives me different vertical viscosity profile compared with the solved viscosity profile of Tutorial_1. Here is the comparison of the two results. The first picture comes from the Tutorial_1, and the second one is for my UW2 model. image image

I attach my code file as text here. viscosity test for Tutorial_1 in UW2.txt

I don't understand what causes the big difference. Hope to get your reply. Thanks a lot Alexstruc

PatriceFRey commented 4 years ago

Hi Alex Log scale ? Patrice

Sent from my iPhone

On 17 Dec 2019, at 8:59 pm, alexstruc notifications@github.com wrote:

 Dear developers: I am trying to replicate the Tutorial_1 of UWGeodynamics in Underworld2. To test the reliability of the model, I input the identical viscosity and plasticity parameters with their units. However, my UW2 model gives me different vertical viscosity profile compared with the solved viscosity profile of Tutorial_1. Here is the comparison of the two results. The first picture comes from the Tutorial_1, and the second one is for my UW2 model.

I attach my code file as text here. viscosity test for Tutorial_1 in UW2.txt

I don't understand what causes the big difference. Hope to get your reply. Thanks a lot Alexstruc

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or unsubscribe.

alexstruc commented 4 years ago

Hello, Patrice Yes. You are right. The problem is my viscosity scale. I get the same profile without using the log scale. Haha, It is a stupid mistake. Thanks a lot for your help. Alexstruc