Hi Jun-wei,<br><br><div class="gmail_quote">On Mon, May 24, 2010 at 7:03 PM, Huang, Jun-Wei <span dir="ltr"><<a href="mailto:Jun-Wei.Huang@nrcan-rncan.gc.ca">Jun-Wei.Huang@nrcan-rncan.gc.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div>
<p><font face="Arial" size="2">Hello:</font>
<br><font face="Arial" size="2">I just started using specfem3d v1.4.3 and have a few basic questions. Could anyone in this email list offer some help? Thank you in advance. </font></p>
<p><font face="Arial" size="2">1). In the manual page 6, it says "MODEL must be one of the following: SoCal and Harvard_LA". Can I provide my own models, and how to create my own models?</font></p></div></blockquote>
<div>You have to know what you are doing when replacing the model, in terms of model dimension and how various components (topography, crystalline basement and Moho) are loaded. <br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div>
<p><font face="Arial" size="2">2). I plan to use this program for an exploration seismology problem, which involves 1000~3000 sources. In the specfem3d.f90 code, line 528, I notice that an error message will be issued if NSOURCES (number of sources?) is larger than 1000. Will this limit my application, in other words, can specfem3d allows more than 1000 shots?</font></p>
</div></blockquote><div>This theoretically can be modified in the code (which I don't think is problematic for the forward simulation you are doing).<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div>
<p><font face="Arial" size="2">3). I initially compiled specfem3d on a single PC (compiled --without-mpi, FC=Fortran) and followed the manual till page 19, where I cannot use UTILS/<a href="http://slice_number.pl" target="_blank">slice_number.pl</a> (actually UTILS/Vistualization/<a href="http://slice_number.pl" target="_blank">slice_number.pl</a>) to generate a slices_file. I have matlab installed. The machine is complaining no sph2utm command. Where and what is the sph2utm command? Is this the only way to load/view the binary kernel output files in matlab or paraview?</font></p>
</div></blockquote><div>You don't have to use the mesh2vtu tools to generate data files for visualization. Saving mesh in AVS format is also fine. sph2utm converts spherical coordinates to UTM coordinates for source and receivers. I'll make sure I supply a version to the repository. But if you need it, you can download it from here:<br>
<br><a href="http://www.physics.utoronto.ca/~liuqy/temp/sph2utm.tar.gz">http://www.physics.utoronto.ca/~liuqy/temp/sph2utm.tar.gz</a><br><br>Good luck with your application,<br></div></div>Qinya<br>