[CIG-SEISMO] Error in homogeneous_halfspace_example with cubit

[s.khan-1 at utwente.nl]

Hi,

I am a PhD scholar at ITC- The Netherlands and learning SPECFEM3D for my research.

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.

Thanks in advance.

Kind regards,
Saad Khan

********** Now playing /home/khan31332/specfem/examples/homogeneous_halfspace/block_mesh.py **********


CUBIT>
CUBIT> #!/usr/bin/env python
CUBIT>
CUBIT>
CUBIT>
Traceback (most recent call last):
  File "<string>", line 1, in <module>
ImportError: No module named boundary_definition
CUBIT>
Traceback (most recent call last):
  File "<string>", line 1, in <module>
ImportError: No module named cubit2specfem3d
CUBIT>
CUBIT>
CUBIT>
CUBIT>
CUBIT> reset
Default CUBIT model file is '/home/khan31332/cubit/cubit01.cub'

Journaled Command: reset

CUBIT> brick x 134000 y 134000 z 60000
Successfully created brick volume 1
Journaled Command: brick x 134000 y 134000 z 60000

CUBIT> volume 1 move x 67000 y 67000 z -30000
Journaled Command: volume 1 move x 67000 y 67000 z -30000

CUBIT>
CUBIT>
CUBIT> # Meshing the volumes
CUBIT>
CUBIT>
CUBIT> volume 1 size '+str(elementsize))
Journaled Command: volume 1 size 3750

CUBIT> mesh volume 1
Volume 1 (Volume 1) meshing scheme automatically set to submap.
Matching intervals successful.
Meshing Volume 1 (Volume 1)
Generated 1296 faces for Surface 1 (Surface 1).
Generated 1296 faces for Surface 2 (Surface 2).
Generated 576 faces for Surface 3 (Surface 3).
Generated 576 faces for Surface 4 (Surface 4).
Generated 576 faces for Surface 5 (Surface 5).
Generated 576 faces for Surface 6 (Surface 6).
Generated 20736 hexes for Volume 1 (Volume 1).
Volume 1 (Volume 1) meshing completed using scheme: submap
Journaled Command: mesh volume 1

CUBIT>
CUBIT>
CUBIT> #### End of meshing
CUBIT>
CUBIT> ###### This is boundary_definition.py of GEOCUBIT
CUBIT> #..... which extracts the bounding faces and defines them into blocks
CUBIT>
Traceback (most recent call last):
  File "<string>", line 1, in <module>
NameError: name 'boundary_definition' is not defined
CUBIT>
Traceback (most recent call last):
  File "<string>", line 1, in <module>
NameError: name 'boundary_definition' is not defined
CUBIT>
CUBIT> #### Define material properties for the 3 volumes ################
CUBIT> #### DEFINE MATERIAL PROPERTIES #######################
CUBIT>
CUBIT>
CUBIT> block 1 name "elastic 1"
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT> block 1 attribute count 6
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT> block 1 attribute index 1 1
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT> block 1 attribute index 2 2800
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT> block 1 attribute index 3 1500
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT> block 1 attribute index 4 2300
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT> block 1 attribute index 5 9000.0
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT> block 1 attribute index 6 0
ERROR: No entity with ID 1 was found
ERROR: No entities specified.
CUBIT>
CUBIT>
CUBIT> #cubit.cmd('block 1 name "acoustic 1" ')       # acoustic material region
CUBIT> #cubit.cmd('block 1 attribute count 4')
CUBIT> #cubit.cmd('block 1 attribute index 1 1  ')     # material 1
CUBIT> #cubit.cmd('block 1 attribute index 2 1480 ')  # vp
CUBIT> #cubit.cmd('block 1 attribute index 3 0 ')      # vs
CUBIT> #cubit.cmd('block 1 attribute index 4 1028 ')  # rho (ocean salt water density:
CUBIT>
CUBIT>
CUBIT> export mesh "top.e" dimension 3 overwrite

Using ExodusII API version 4.98 and DB version 4.98 to export "top.e".


Executive Exodus summary:
  Number of dimensions     = 3
  Number of element blocks = 1
  Number of sidesets       = 0
  Number of nodesets       = 0
  Number of bc sets        = 1
  Number of elements       = 20736
  Number of nodes          = 23273

Finished writing top.e
Journaled Command: export mesh "top.e" dimension 3 overwrite

CUBIT> save as "meshing.cub" overwrite

Exported:   1 Volume

Successfully saved CUBIT file 'meshing.cub'
Journaled Command: save as "meshing.cub" overwrite

CUBIT>
CUBIT> #### Export to SESAME format using cubit2specfem3d.py of GEOCUBIT
CUBIT>
CUBIT>
CUBIT>
Traceback (most recent call last):
  File "<string>", line 1, in <module>
NameError: name 'cubit2specfem3d' is not defined
CUBIT>
CUBIT> # all files needed by SCOTCH are now in directory MESH
CUBIT>
CUBIT>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-seismo/attachments/20130815/f4a10b3a/attachment-0001.html>