Python collect示例

示例#1

文件： make_mesh_waterlayer_verticalboundary.py 项目： geodynamics/specfem3d

except:
    pass

cubit.cmd('reset')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 1 move x 33500 y 67000 z -30000')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 2 move x 100500 y 67000 z -30000')
cubit.cmd('merge all')

# Meshing the volumes
elementsize = 3000.0

cubit.cmd('volume 1 size '+str(elementsize))
cubit.cmd('volume 2 size '+str(elementsize))
cubit.cmd('mesh volume 1 2')


from geocubitlib import boundary_definition,exportlib

# boundary
boundary_definition.define_bc(parallel=True)

# file export
exportlib.collect(outdir=CUBIToutput)
exportlib.e2SEM(outdir=SEMoutput)

cubit.cmd('save as "meshing.cub" overwrite')

示例#2

文件： GEOCUBIT.py 项目： hejunzhu/specfem3d

        except:
            output = 'totalmesh_merged'
        output = output.upper()
        #
        collect(menu.cpuxmin,
                menu.cpuxmax,
                menu.cpuymin,
                menu.cpuymax,
                menu.cpux,
                menu.cpuy,
                menu.cubfiles,
                menu.ckbound_method1,
                menu.ckbound_method2,
                menu.merge_tolerance,
                curverefining=menu.curverefining,
                outfilename=output,
                qlog=menu.qlog,
                export2SPECFEM3D=menu.export2SPECFEM3D,
                listblock=menu.listblock,
                listflag=menu.listflag,
                outdir=menu.SPECFEM3D_output_dir,
                add_sea=menu.add_sea,
                decimate=menu.decimate,
                cpml=menu.cpml,
                cpml_size=menu.cpml_size,
                top_absorbing=menu.top_absorbing,
                hex27=menu.hex27)

    if menu.export2SPECFEM3D and not menu.collect:
        from geocubitlib.exportlib import e2SEM
        print menu.cubfiles

示例#3

文件： create_mesh.py 项目： AlaaHadji/specfem3d

    #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'
    f=open(nummaterial_velocity_file,'w')
    for block in material_cfg:
        print block
        s=block['material region']+' '
        s=s+block['id_block']+' '
        s=s+block['rho']+' '
        s=s+block['vp']+' '
        s=s+block['vs']+' '

示例#4

文件： GEOCUBIT.py 项目： EtienneBachmann/specfem3d

    
    #GEOMETRY
    if menu.build_surface:              
        from geocubitlib import surfaces
        surfaces.surfaces()                 
    if menu.build_volume:                   
        from geocubitlib import volumes     
        volumes.volumes()                   
               
    #MESHING
    if menu.meshing:
        from geocubitlib import mesh_volume
        mesh_volume.mesh()             
    
    #MERGING and EXPORTING
    if menu.collect:
        from geocubitlib.exportlib import collect
        try:
            output=menu.output
        except:
            output='totalmesh_merged'
        output=output.upper()
        #
        collect(menu.cpuxmin,menu.cpuxmax,menu.cpuymin,menu.cpuymax,menu.cpux,menu.cpuy,menu.cubfiles,menu.ckbound_method1,menu.ckbound_method2,menu.merge_tolerance,curverefining=menu.curverefining,outfilename=output,qlog=menu.qlog,export2SPECFEM3D=menu.export2SPECFEM3D,listblock=menu.listblock,listflag=menu.listflag,outdir=menu.SPECFEM3D_output_dir,add_sea=menu.add_sea,decimate=menu.decimate,cpml=menu.cpml,cpml_size=menu.cpml_size,top_absorbing=menu.top_absorbing,hex27=menu.hex27)
        
    if menu.export2SPECFEM3D and not menu.collect:
        from geocubitlib.exportlib import e2SEM
        print menu.cubfiles
        print 'hex27 ',menu.hex27
        e2SEM(files=menu.cubfiles,listblock=menu.listblock,listflag=menu.listflag,outdir=menu.SPECFEM3D_output_dir,hex27=menu.hex27)

示例#5

文件： create_mesh.py 项目： geodynamics/specfem3d

#!/usr/bin/env python



import cubit
try:
    cubit.init([""])
except:
    pass

cfg='fullpath_myparameterfile.cfg'
o='fullpath_specfem3d_meshfiles_outputdir'
o1='fullpath_cubit_meshfiles_outputdir'

from geocubitlib import volumes, mesh_volume, exportlib
volumes.volumes(cfg)
mesh_volume.mesh(cfg)
exportlib.collect(outdir=o1)
#exportlib.e2SEM(outdir=o,listblock=[4],listflag=[-1])
exportlib.e2SEM(outdir=o)

示例#6

文件： create_mesh.py 项目： geodynamics/specfem3d

import cubit
try:
    cubit.init([""])
except:
    pass



cubit.cmd('reset')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 1 move x 33500 y 67000 z -30000')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 2 move x 100500 y 67000 z -30000')
cubit.cmd('merge all')

# Meshing the volumes
elementsize = 3000.0

cubit.cmd('volume 1 size '+str(elementsize))
cubit.cmd('volume 2 size '+str(elementsize))
cubit.cmd('mesh volume 1 2')


from geocubitlib import boundary_definition,exportlib

boundary_definition.define_bc(parallel=True)

exportlib.collect(outdir='.')
exportlib.e2SEM(outdir='.')

示例#7

#!/usr/bin/env python

import cubit
try:
    cubit.init([""])
except:
    pass

cubit.cmd('reset')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 1 move x 33500 y 67000 z -30000')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 2 move x 100500 y 67000 z -30000')
cubit.cmd('merge all')

# Meshing the volumes
elementsize = 3000.0

cubit.cmd('volume 1 size ' + str(elementsize))
cubit.cmd('volume 2 size ' + str(elementsize))
cubit.cmd('mesh volume 1 2')

from geocubitlib import boundary_definition, exportlib

boundary_definition.define_bc(parallel=True)

exportlib.collect(outdir='.')
exportlib.e2SEM(outdir='.')

示例#8

文件： create_mesh.py 项目： lixuekai2001/Specfem-3d

#!/usr/bin/env python

import cubit
try:
    cubit.init([""])
except:
    pass

cfg = './tomographic_model//tomographic_model.cfg'
o = './tomographic_model/specfem3d_mesh_files'
o1 = './tomographic_model/cubit_mesh_files'

from geocubitlib import volumes, mesh_volume, exportlib
volumes.volumes(cfg)
mesh_volume.mesh(cfg)
exportlib.collect(outdir=o1)
exportlib.e2SEM(outdir=o, listblock=[4], listflag=[-1])

示例#9

文件： create_mesh.py 项目： zyex1108/specfem3d

    #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

示例#10

import cubit
try:
    cubit.init([""])
except:
    pass

cubit.cmd('reset')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 1 move x 33500 y 67000 z -30000')
cubit.cmd('brick x 67000 y 134000 z 60000')
cubit.cmd('volume 2 move x 100500 y 67000 z -30000')
cubit.cmd('merge all')

# Meshing the volumes
elementsize = 3000.0

cubit.cmd('volume 1 size ' + str(elementsize))
cubit.cmd('volume 2 size ' + str(elementsize))
cubit.cmd('mesh volume 1 2')

from geocubitlib import boundary_definition, exportlib

# boundary
boundary_definition.define_bc(parallel=True)

# file export
exportlib.collect(outdir=CUBIToutput)
exportlib.e2SEM(outdir=SEMoutput)

cubit.cmd('save as "meshing.cub" overwrite')