[cig-commits] r22878 - in seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace: . MESH

carltape at geodynamics.org carltape at geodynamics.org
Fri Sep 27 00:17:42 PDT 2013 

Author: carltape
Date: 2013-09-27 00:17:41 -0700 (Fri, 27 Sep 2013)
New Revision: 22878

Added:
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference.verticalboundary
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH_GEOCUBIT/
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlayer.sh
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/waterlayer_only.cfg
Removed:
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/boundary_definition.py
   seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/cubit2specfem3d.py
Log:
updated waterlayered_halfspace example; did not delete the old .py scripts


Added: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference	                        (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference	2013-09-27 07:17:41 UTC (rev 22878)
@@ -0,0 +1,3 @@
+1   1 1028 1480 0 9999.0 9999.0 0  
+1   2 1028 1480 0 9999.0 9999.0 0 
+1   3 1028 1480 0 9999.0 9999.0 0
\ No newline at end of file

Added: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference.verticalboundary
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference.verticalboundary	                        (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/MESH/nummaterial_velocity_file.reference.verticalboundary	2013-09-27 07:17:41 UTC (rev 22878)
@@ -0,0 +1,2 @@
+1   1 2300 3000 0 9999.0 9999.0 0
+2   2 2300 3000 0.1 9999.0 9999.0 0
\ No newline at end of file

Deleted: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/boundary_definition.py
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/boundary_definition.py	2013-09-27 06:46:29 UTC (rev 22877)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/boundary_definition.py	2013-09-27 07:17:41 UTC (rev 22878)
@@ -1 +0,0 @@
-link ../../CUBIT_GEOCUBIT/geocubitlib/boundary_definition.py
\ No newline at end of file

Deleted: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/cubit2specfem3d.py
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/cubit2specfem3d.py	2013-09-27 06:46:29 UTC (rev 22877)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/cubit2specfem3d.py	2013-09-27 07:17:41 UTC (rev 22878)
@@ -1 +0,0 @@
-link ../../CUBIT_GEOCUBIT/geocubitlib/cubit2specfem3d.py
\ No newline at end of file

Added: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py	                        (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py	2013-09-27 07:17:41 UTC (rev 22878)
@@ -0,0 +1,33 @@
+#!/usr/bin/env python
+
+SEMoutput='MESH'
+CUBIToutput='MESH_GEOCUBIT'
+
+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=CUBIToutput)
+exportlib.e2SEM(outdir=SEMoutput)


Property changes on: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlater_verticalboundary.py
___________________________________________________________________
Added: svn:executable
   + *

Added: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlayer.sh
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlayer.sh	                        (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlayer.sh	2013-09-27 07:17:41 UTC (rev 22878)
@@ -0,0 +1,19 @@
+#!/bin/bash
+
+#
+# Temporary instructions
+#
+# 1. set path to GEOCUBIT base directory (../../CUBIT_GEOCUBIT/), for example:
+#    export PYTHONPATH=$PYTHONPATH:/import/c/w/tape/3D/SPECFEM3D/CUBIT_GEOCUBIT
+#    export PATH=$PATH:/import/c/w/tape/3D/SPECFEM3D/CUBIT_GEOCUBIT
+#
+#    check path:
+#    which GEOCUBIT.py
+#
+# 2. run this script to generate mesh
+#    ./make_mesh_waterlayer.sh
+#
+
+GEOCUBIT.py --build_volume --mesh --cfg=waterlayer_only.cfg
+GEOCUBIT.py --collect --meshfiles=MESH_GEOCUBIT/mesh_vol_0.e --export2SPECFEM3D --SEMoutput=MESH
+cp MESH/nummaterial_velocity_file.reference MESH/nummaterial_velocity_file
\ No newline at end of file


Property changes on: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/make_mesh_waterlayer.sh
___________________________________________________________________
Added: svn:executable
   + *

Added: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/waterlayer_only.cfg
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/waterlayer_only.cfg	                        (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/waterlayer_only.cfg	2013-09-27 07:17:41 UTC (rev 22878)
@@ -0,0 +1,35 @@
+[cubit.options]
+cubit_info=on
+echo_info=on
+jou_info=on
+jer_info=on
+working_dir=tmp
+output_dir=MESH_GEOCUBIT
+save_geometry_cubit = True
+save_surface_cubit = False
+monitored_cpu=0
+
+[simulation.cpu_parameters]
+nodes=1
+#
+[geometry.volumes]
+volume_type                     = layercake_volume_ascii_regulargrid_regularmap
+latitude_min                    = 0.
+latitude_max                    = 134000.
+longitude_min                   = 0.
+longitude_max                   = 134000.
+unit                            = utm
+
+[geometry.volumes.layercake]
+nz = 4
+#included the bottom
+bottomflat = True
+depth_bottom = -60000
+geometry_format=regmesh
+zdepth=-60000,-9000,-3000,0
+
+[meshing]
+map_meshing_type=regularmap
+iv_interval=17,2,1
+size=3000
+or_mesh_scheme=map


Property changes on: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/waterlayer_only.cfg
___________________________________________________________________
Added: svn:executable
   + *

More information about the CIG-COMMITS mailing list