<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi,
<div><br>
</div>
<div>I am a PhD scholar at ITC- The Netherlands and learning SPECFEM3D for my research. </div>
<div><br>
</div>
<div>I am using Cubit for mesh generation and following the instruction given the the example of homogeneous_halfspace. There it tells to play a journal file "block_mesh.py" but when it is open I am supposed to get some file in MESH sub directory, which infact
are not created. I receive the following errors (end of email) in the command line. I will be grateful if somebody can respond to the problem. </div>
<div><br>
</div>
<div>Thanks in advance.</div>
<div><br>
</div>
<div>Kind regards,</div>
<div>Saad Khan</div>
<div><br>
</div>
<div>
<div>********** Now playing /home/khan31332/specfem/examples/homogeneous_halfspace/block_mesh.py **********</div>
<div><br>
</div>
<div><br>
</div>
<div>CUBIT> </div>
<div>CUBIT> #!/usr/bin/env python</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>Traceback (most recent call last):</div>
<div> File "<string>", line 1, in <module></div>
<div>ImportError: No module named boundary_definition</div>
<div>CUBIT> </div>
<div>Traceback (most recent call last):</div>
<div> File "<string>", line 1, in <module></div>
<div>ImportError: No module named cubit2specfem3d</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> reset</div>
<div>Default CUBIT model file is '/home/khan31332/cubit/cubit01.cub'</div>
<div><br>
</div>
<div>Journaled Command: reset</div>
<div><br>
</div>
<div>CUBIT> brick x 134000 y 134000 z 60000</div>
<div>Successfully created brick volume 1 </div>
<div>Journaled Command: brick x 134000 y 134000 z 60000</div>
<div><br>
</div>
<div>CUBIT> volume 1 move x 67000 y 67000 z -30000</div>
<div>Journaled Command: volume 1 move x 67000 y 67000 z -30000</div>
<div><br>
</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> # Meshing the volumes</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> volume 1 size '+str(elementsize))</div>
<div>Journaled Command: volume 1 size 3750</div>
<div><br>
</div>
<div>CUBIT> mesh volume 1</div>
<div>Volume 1 (Volume 1) meshing scheme automatically set to submap.</div>
<div>Matching intervals successful.</div>
<div>Meshing Volume 1 (Volume 1)</div>
<div>Generated 1296 faces for Surface 1 (Surface 1).</div>
<div>Generated 1296 faces for Surface 2 (Surface 2).</div>
<div>Generated 576 faces for Surface 3 (Surface 3).</div>
<div>Generated 576 faces for Surface 4 (Surface 4).</div>
<div>Generated 576 faces for Surface 5 (Surface 5).</div>
<div>Generated 576 faces for Surface 6 (Surface 6).</div>
<div>Generated 20736 hexes for Volume 1 (Volume 1).</div>
<div>Volume 1 (Volume 1) meshing completed using scheme: submap</div>
<div>Journaled Command: mesh volume 1</div>
<div><br>
</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> #### End of meshing</div>
<div>CUBIT> </div>
<div>CUBIT> ###### This is boundary_definition.py of GEOCUBIT</div>
<div>CUBIT> #..... which extracts the bounding faces and defines them into blocks</div>
<div>CUBIT> </div>
<div>Traceback (most recent call last):</div>
<div> File "<string>", line 1, in <module></div>
<div>NameError: name 'boundary_definition' is not defined</div>
<div>CUBIT> </div>
<div>Traceback (most recent call last):</div>
<div> File "<string>", line 1, in <module></div>
<div>NameError: name 'boundary_definition' is not defined</div>
<div>CUBIT> </div>
<div>CUBIT> #### Define material properties for the 3 volumes ################</div>
<div>CUBIT> #### DEFINE MATERIAL PROPERTIES #######################</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> block 1 name "elastic 1" </div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> block 1 attribute count 6</div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> block 1 attribute index 1 1</div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> block 1 attribute index 2 2800</div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> block 1 attribute index 3 1500</div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> block 1 attribute index 4 2300</div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> block 1 attribute index 5 9000.0</div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> block 1 attribute index 6 0 </div>
<div>ERROR: No entity with ID 1 was found</div>
<div>ERROR: No entities specified.</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> #cubit.cmd('block 1 name "acoustic 1" ') # acoustic material region</div>
<div>CUBIT> #cubit.cmd('block 1 attribute count 4')</div>
<div>CUBIT> #cubit.cmd('block 1 attribute index 1 1 ') # material 1</div>
<div>CUBIT> #cubit.cmd('block 1 attribute index 2 1480 ') # vp</div>
<div>CUBIT> #cubit.cmd('block 1 attribute index 3 0 ') # vs</div>
<div>CUBIT> #cubit.cmd('block 1 attribute index 4 1028 ') # rho (ocean salt water density:</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> export mesh "top.e" dimension 3 overwrite</div>
<div><br>
</div>
<div>Using ExodusII API version 4.98 and DB version 4.98 to export "top.e".</div>
<div><br>
</div>
<div><br>
</div>
<div>Executive Exodus summary:</div>
<div> Number of dimensions = 3</div>
<div> Number of element blocks = 1</div>
<div> Number of sidesets = 0</div>
<div> Number of nodesets = 0</div>
<div> Number of bc sets = 1</div>
<div> Number of elements = 20736</div>
<div> Number of nodes = 23273</div>
<div><br>
</div>
<div>Finished writing top.e</div>
<div>Journaled Command: export mesh "top.e" dimension 3 overwrite</div>
<div><br>
</div>
<div>CUBIT> save as "meshing.cub" overwrite</div>
<div><br>
</div>
<div>Exported: 1 Volume</div>
<div><br>
</div>
<div>Successfully saved CUBIT file 'meshing.cub'</div>
<div>Journaled Command: save as "meshing.cub" overwrite</div>
<div><br>
</div>
<div>CUBIT> </div>
<div>CUBIT> #### Export to SESAME format using cubit2specfem3d.py of GEOCUBIT</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>CUBIT> </div>
<div>Traceback (most recent call last):</div>
<div> File "<string>", line 1, in <module></div>
<div>NameError: name 'cubit2specfem3d' is not defined</div>
<div>CUBIT> </div>
<div>CUBIT> # all files needed by SCOTCH are now in directory MESH</div>
<div>CUBIT> </div>
<div>CUBIT> </div>
</div>
</div>
<span id="skype_highlighting_settings" display="none" autoextractnumbers="1" hightlightingswitched="0"></span>
</body>
</html>