cubit.cmd('block ' + str(id_block) + ' attribute index 5 400.0') # Qkappa cubit.cmd('block ' + str(id_block) + ' attribute index 6 200.0') # Qmu cubit.cmd('block ' + str(id_block) + ' attribute index 7 0') # anisotropy_flag else: print("domain setup not supported yet, exiting") sys.exit(1) print("") print("exporting to SPECFEM3D-format:") print("") # Export to SPECFEM3D format # for hex27 elements if 1 == 1: cubit2specfem3d.export2SPECFEM3D('MESH/', hex27=True) else: cubit2specfem3d.export2SPECFEM3D('MESH/') # backup cubit cubit.cmd('save as "MESH/meshing.cub" overwrite') else: from geocubitlib import exportlib print("") print("exporting to SPECFEM3D-format:") print("") # Export to SPECFEM3D format # note: exportlib-commands will overwrite material properties # for hex27 elements if 1 == 1: exportlib.define_blocks(outdir='MESH/', save_cubfile=True,
cubit.cmd('block 2 attribute count 6') cubit.cmd('block 2 attribute index 1 2 ') # volume 2 cubit.cmd('block 2 attribute index 2 7500 ') cubit.cmd('block 2 attribute index 3 4300 ') cubit.cmd('block 2 attribute index 4 3200 ') cubit.cmd('block 2 attribute index 5 9999.0 ') cubit.cmd('block 2 attribute index 6 0 ') # anisotropy_flag cubit.cmd('block 3 name "elastic 2" ') # material region cubit.cmd('block 3 attribute count 6') cubit.cmd('block 3 attribute index 1 3 ') # same properties as for volume 2 cubit.cmd('block 3 attribute index 2 7500 ') cubit.cmd('block 3 attribute index 3 4300 ') cubit.cmd('block 3 attribute index 4 3200 ') cubit.cmd('block 3 attribute index 5 9999.0 ') cubit.cmd('block 3 attribute index 6 0 ') # anisotropy_flag cubit.cmd('export mesh "' + CUBIToutput + '/top.e" dimension 3 overwrite') cubit.cmd('save as "' + CUBIToutput + '/top.cub" overwrite') #### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT # file export cubit2specfem3d.export2SPECFEM3D(SEMoutput) # screen shot cubit.cmd('hardcopy "' + CUBIToutput + '/waterlayer_fig8.png" png') # all files needed by SCOTCH are now in directory MESH
# Define material properties print "#### DEFINE MATERIAL PROPERTIES #######################" # elastic model cubit.cmd('block '+str(id_block)+' name "elastic tomography_model.xyz 1" ') # elastic material region cubit.cmd('block '+str(id_block)+' attribute count 2') cubit.cmd('block '+str(id_block)+' attribute index 1 -1') # flag for material: -1 for 1. undefined material cubit.cmd('block '+str(id_block)+' attribute index 2 2') # flag for tomographic model # acoustic model #cubit.cmd('block '+str(id_block)+' name "acoustic tomographic 1" ') # acoustic material region #cubit.cmd('block '+str(id_block)+' attribute count 2') #cubit.cmd('block '+str(id_block)+' attribute index 1 -1') # material 1 #cubit.cmd('block '+str(id_block)+' attribute index 2 2') # tomographic model flag # creates MESH/ directory for file output os.system('mkdir -p MESH') cubit.cmd('export mesh "MESH/top.e" dimension 3 overwrite') cubit.cmd('save as "MESH/meshing.cub" overwrite') # Export to SPECFEM3D format #note: exportlib-commands will overwrite the above material properties # exportlib.collect(outdir='MESH/') # exportlib.e2SEM(outdir='MESH/') cubit2specfem3d.export2SPECFEM3D('MESH/') # all files needed by xdecompose_mesh are now in directory MESH/
cubit.cmd('block '+str(id_block)+' attribute index 4 2300') # rho cubit.cmd('block '+str(id_block)+' attribute index 5 9999.0') # Qkappa cubit.cmd('block '+str(id_block)+' attribute index 6 9000.0') # Qmu cubit.cmd('block '+str(id_block)+' attribute index 7 0') # anisotropy_flag # acoustic material #cubit.cmd('block '+str(id_block)+' name "acoustic 1" ') # acoustic material region #cubit.cmd('block '+str(id_block)+' attribute count 4') #cubit.cmd('block '+str(id_block)+' attribute index 1 1 ') # material 1 #cubit.cmd('block '+str(id_block)+' attribute index 2 1480 ') # vp #cubit.cmd('block '+str(id_block)+' attribute index 3 0 ') # vs #cubit.cmd('block '+str(id_block)+' attribute index 4 1028 ') # rho (ocean salt water density: print "" print "exporting to SPECFEM3D-format:" print "" # Export to SPECFEM3D format cubit2specfem3d.export2SPECFEM3D('MESH/') else: from geocubitlib import exportlib print "" print "exporting to SPECFEM3D-format:" print "" # Export to SPECFEM3D format # note: exportlib-commands will overwrite material properties exportlib.collect(outdir='MESH/') exportlib.e2SEM('MESH/') # Define material properties print "#### DEFINE MATERIAL PROPERTIES #######################" # elastic material material_cfg=[{'material region':'2','id_block':'1','vp':'2800','vs':'1500','rho':'2300','Qkappa':'9990.0','Qmu':'9000.0','anisotropy_flag':'0'}] # modifies material file nummaterial_velocity_file='MESH/nummaterial_velocity_file'
cubit.cmd('block 2 attribute index 6 0 ') # anisotropy_flag cubit.cmd('block 3 name "elastic 3" ') # elastic material region cubit.cmd('block 3 attribute count 6') cubit.cmd('block 3 attribute index 1 3 ') # same properties as for volume 2 cubit.cmd('block 3 attribute index 2 7500 ') cubit.cmd('block 3 attribute index 3 4300 ') cubit.cmd('block 3 attribute index 4 3200 ') cubit.cmd('block 3 attribute index 5 9000.0 ') cubit.cmd('block 3 attribute index 6 0 ') # anisotropy_flag cubit.cmd('export mesh "' + CUBIToutput + '/top.e" dimension 3 overwrite') cubit.cmd('save as "' + CUBIToutput + '/top.cub" overwrite') #### Export to SPECFEM3D format using cubit2specfem3d.py of GEOCUBIT cubit2specfem3d.export2SPECFEM3D(SEMoutput) # screen shot # (could crash version < 16.4) if version_major >= 16 and version_minor >= 4: cubit.cmd('view top') cubit.cmd('hardcopy "' + CUBIToutput + '/layered-fig8.png" png') # all files needed by SCOTCH are now in directory MESH