<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="color:#1F497D">Dear Emanuele,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D">Please find the scripts and the volume in attachment. I use block_mesh.py after making some modification to it. I am using Cubit 13.0.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D">Thanks<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D">Best regards,<br>
Saad Khan<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> cig-seismo-bounces@geodynamics.org [mailto:cig-seismo-bounces@geodynamics.org]
<b>On Behalf Of </b>emanuele casarotti<br>
<b>Sent:</b> Monday, January 13, 2014 8:51 PM<br>
<b>To:</b> cig-seismo@geodynamics.org<br>
<b>Subject:</b> Re: [CIG-SEISMO] cubit2specfem3d export<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hi<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Can you share the scripta and the volumes?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">And the version of cubit...<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><br>
Emanuele<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
Il giorno 13/gen/2014, alle ore 19:47, "<a href="mailto:s.khan-1@utwente.nl">s.khan-1@utwente.nl</a>" <<a href="mailto:s.khan-1@utwente.nl">s.khan-1@utwente.nl</a>> ha scritto:<o:p></o:p></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal">Dear All,<o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal">The boundary_definition.py script does not catch my topographic surface. I have to define it manually after running the script (boundary_definition.py). Then I export the mesh with run_cubit2specfem3d.py, which says
<o:p></o:p></p>
<p class="MsoNormal"><span lang="EN">:</span><o:p></o:p></p>
<p class="MsoNormal">Writing MESH/free_surface_file.....<o:p></o:p></p>
<p class="MsoNormal"> block name: face_topo id: 2<o:p></o:p></p>
<p class="MsoNormal"> number of faces = 0<o:p></o:p></p>
<p class="MsoNormal">:<o:p></o:p></p>
<p class="MsoNormal">:<o:p></o:p></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN"># all files needed by SCOTCH are now in directory MESH</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN">Rest of the process goes smoothly and I also get result in the OUTPUT. The question is, what does “number of face=0” mean any does it affect the results?</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN">Another thing, when I use very fine mesh (e.g. mesh interval 3000) the result (in seismograms) produces a lot of NaNs and at some places Infinities. Why is it so?</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN">Note: The topo surface has been modified while the rest of geo model is the same as given in homogeneous_halfspace example. The topo surface is made of several surfaces which are combined to one using simplify option in
cubit. I am using SPECFEM3D_Cartesian 2.1.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN">Thanks in advance.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN">Best regards,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN">Saad Khan</span><o:p></o:p></p>
</div>
</blockquote>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman","serif"">_______________________________________________<br>
CIG-SEISMO mailing list<br>
<a href="mailto:CIG-SEISMO@geodynamics.org">CIG-SEISMO@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo</a><o:p></o:p></span></p>
</div>
</blockquote>
</div>
</body>
</html>