Dear Pylith team,
I’m currently running viscoelastic simulations (for post-seismic study) using Burgers rheology. Burger rheology is represented by a generalized Maxwell model in PyLith. However, when setting μ₃ to zero, I encountered the following error: Error converting spatial database values for maxwell_time at ( 347588 72137.1 -71649.6) in spatial database ‘mantle Maxwell viscoelastic properties’. Found non-positive viscosity 3 (0).Found non-positive shear modulus ratio 3 (0).To resolve this, I modified the parameters to: viscosity = [0.9744×10¹⁸, 0.0256×10¹⁸, 1.0×10⁻¹⁰] Pa·s,
shear_modulus_ratio = [0.474, 0.525, 1.0×10⁻⁶]. Could you verify whether this parameter configuration properly represents Burgers rheology behavior?
Best,
Yali Shao