<div dir="ltr"><div><div>Dear Yingzi,<br><br></div>Thank you for the comments. If I find/write a better solution, I will post it to the user group.<br>Cheers,<br><br></div>Dylan<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Wed, Oct 30, 2013 at 11:28 AM, Yingzi Ying <span dir="ltr"><<a href="mailto:yingzi.ying@me.com" target="_blank">yingzi.ying@me.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Dylan,<br>
<br>
I am in similar situation to do 2d simulations with some velocity structures, but with fluid-solid coupling...<br>
<br>
One solution to is to supply with "DATA/model_velocity.dat_<u></u>input" which defines value in each GLL points. I tried this method but no success. PS. I have reported my problem here but no response...<br>
<br>
Another way is to give value of each spectral element, that you can define "nbmodels" and "nbregions" in "DATA/Par_file" with your own.<br>
<br>
If you have any better solution, please also let me know.<br>
<br>
Good luck,<br>
Yingzi<br>
On 10/30/2013 03:14 PM, <a href="mailto:cig-seismo-request@geodynamics.org" target="_blank">cig-seismo-request@<u></u>geodynamics.org</a> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Send CIG-SEISMO mailing list submissions to<br>
<a href="mailto:cig-seismo@geodynamics.org" target="_blank">cig-seismo@geodynamics.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo" target="_blank">http://geodynamics.org/cgi-<u></u>bin/mailman/listinfo/cig-<u></u>seismo</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:cig-seismo-request@geodynamics.org" target="_blank">cig-seismo-request@<u></u>geodynamics.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:cig-seismo-owner@geodynamics.org" target="_blank">cig-seismo-owner@geodynamics.<u></u>org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of CIG-SEISMO digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. SPECFEM2D - random velocity models and meshing (Dylan Mikesell)<br>
<br>
<br>
------------------------------<u></u>------------------------------<u></u>----------<br>
<br>
Message: 1<br>
Date: Wed, 30 Oct 2013 10:42:51 -0400<br>
From: Dylan Mikesell <<a href="mailto:mikesell@MIT.EDU" target="_blank">mikesell@MIT.EDU</a>><br>
To: <a href="mailto:cig-seismo@geodynamics.org" target="_blank">cig-seismo@geodynamics.org</a><br>
Subject: [CIG-SEISMO] SPECFEM2D - random velocity models and meshing<br>
Message-ID:<br>
<CADgqVRb+<u></u>3JFKBvTF7c5VChkhPsv7HKiNOLk-<u></u>q78Gd=<a href="mailto:ULucOM8w@mail.gmail.com" target="_blank">ULucOM8w@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Dear SPECFEM2D group,<br>
<br>
I have a question regarding meshing and the EXAMPLES/INDUSTRIAL_FORMAT<br>
example.<br>
<br>
I would like to simulate p-wave propagation with SPECFEM2D in the attached<br>
velocity model.<br>
(the left plot shows the background velocity field, the right plot shows<br>
the field after draping a random perturbation field over the background. I<br>
am using the right plot for simulation).<br>
<br>
I was trying to use/understand the INDUSTRIAL_FORMAT example so that I<br>
could just drop this velocity model into place and go. I have a SeisUnix<br>
velocity model. As far as I could tell, that example created a mesh with<br>
xmeshfem2D and then mapped some given SEP velocity model onto that mesh<br>
using a nearest neighbor routine. I can modify things to get an SU model to<br>
work, but I am not sure that is the way I need to go.<br>
<br>
Since I do not know too much about the meshing, I thought I would ask the<br>
pros. Will the internal xmeshfem2D work for this velocity model or should I<br>
use something else? I see Gmsh in the 2D user manual and Cubit in the 3D<br>
Cartesian user manual.<br>
<br>
The boundary conditions are not a big deal at the moment. I see there is a<br>
little extra work to get PMLs using the external meshers.<br>
<br>
Thank you for any thoughts.<br>
<br>
Dylan<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://geodynamics.org/pipermail/cig-seismo/attachments/20131030/c7ff00fe/attachment.html" target="_blank">http://geodynamics.org/<u></u>pipermail/cig-seismo/<u></u>attachments/20131030/c7ff00fe/<u></u>attachment.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: rand_model.pdf<br>
Type: application/pdf<br>
Size: 2583634 bytes<br>
Desc: not available<br>
URL: <<a href="http://geodynamics.org/pipermail/cig-seismo/attachments/20131030/c7ff00fe/attachment.pdf" target="_blank">http://geodynamics.org/<u></u>pipermail/cig-seismo/<u></u>attachments/20131030/c7ff00fe/<u></u>attachment.pdf</a>><br>
<br>
------------------------------<br>
<br>
______________________________<u></u>_________________<br>
CIG-SEISMO mailing list<br>
<a href="mailto:CIG-SEISMO@geodynamics.org" target="_blank">CIG-SEISMO@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo" target="_blank">http://geodynamics.org/cgi-<u></u>bin/mailman/listinfo/cig-<u></u>seismo</a><br>
<br>
End of CIG-SEISMO Digest, Vol 69, Issue 6<br>
******************************<u></u>***********<br>
<br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
CIG-SEISMO mailing list<br>
<a href="mailto:CIG-SEISMO@geodynamics.org" target="_blank">CIG-SEISMO@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo" target="_blank">http://geodynamics.org/cgi-<u></u>bin/mailman/listinfo/cig-<u></u>seismo</a><br>
</blockquote></div><br></div>