[cig-commits] r22378 - in short/3D/PyLith/trunk: . examples/meshing/surface_nurbs examples/meshing/surface_nurbs/subduction
brad at geodynamics.org
brad at geodynamics.org
Wed Jun 19 14:39:59 PDT 2013
Author: brad
Date: 2013-06-19 14:39:59 -0700 (Wed, 19 Jun 2013)
New Revision: 22378
Added:
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.am
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.in
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/README
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/bc.jou
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/geometry.jou
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_netsurf.py
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_points.txt
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/mesh.jou
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_points.txt
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_skinsurf.py
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_netsurf.py
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_points.txt
Modified:
short/3D/PyLith/trunk/configure.ac
short/3D/PyLith/trunk/examples/meshing/surface_nurbs/Makefile.am
Log:
Added subduction meshing examples from v1.7-trunk.
Modified: short/3D/PyLith/trunk/configure.ac
===================================================================
--- short/3D/PyLith/trunk/configure.ac 2013-06-19 20:41:17 UTC (rev 22377)
+++ short/3D/PyLith/trunk/configure.ac 2013-06-19 21:39:59 UTC (rev 22378)
@@ -396,6 +396,7 @@
examples/meshing/surface_nurbs/dem/vlines/Makefile
examples/meshing/surface_nurbs/triangles/Makefile
examples/meshing/surface_nurbs/merge_surfs/Makefile
+ examples/meshing/surface_nurbs/subduction/Makefile
examples/meshing/cubit_cellsize/Makefile
templates/Makefile
share/Makefile
Modified: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/Makefile.am 2013-06-19 20:41:17 UTC (rev 22377)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/Makefile.am 2013-06-19 21:39:59 UTC (rev 22378)
@@ -20,7 +20,8 @@
contours \
dem \
triangles \
- merge_surfs
+ merge_surfs \
+ subduction
dist_noinst_DATA = \
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.am (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.am 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,32 @@
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010-2013 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+dist_noinst_DATA = \
+ README \
+ geometry.jou \
+ mesh.jou \
+ bc.jou \
+ topobath_netsurf.py \
+ topobath_points.txt \
+ interface_netsurf.py \
+ interface_points.txt \
+ splay_skinsurf.py \
+ splay_points.txt
+
+
+# End of file
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.in
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.in (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/Makefile.in 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,451 @@
+# Makefile.in generated by automake 1.11.1 from Makefile.am.
+# @configure_input@
+
+# Copyright (C) 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002,
+# 2003, 2004, 2005, 2006, 2007, 2008, 2009 Free Software Foundation,
+# Inc.
+# This Makefile.in is free software; the Free Software Foundation
+# gives unlimited permission to copy and/or distribute it,
+# with or without modifications, as long as this notice is preserved.
+
+# This program is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY, to the extent permitted by law; without
+# even the implied warranty of MERCHANTABILITY or FITNESS FOR A
+# PARTICULAR PURPOSE.
+
+ at SET_MAKE@
+
+# -*- Makefile -*-
+#
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010-2013 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+
+VPATH = @srcdir@
+pkgdatadir = $(datadir)/@PACKAGE@
+pkgincludedir = $(includedir)/@PACKAGE@
+pkglibdir = $(libdir)/@PACKAGE@
+pkglibexecdir = $(libexecdir)/@PACKAGE@
+am__cd = CDPATH="$${ZSH_VERSION+.}$(PATH_SEPARATOR)" && cd
+install_sh_DATA = $(install_sh) -c -m 644
+install_sh_PROGRAM = $(install_sh) -c
+install_sh_SCRIPT = $(install_sh) -c
+INSTALL_HEADER = $(INSTALL_DATA)
+transform = $(program_transform_name)
+NORMAL_INSTALL = :
+PRE_INSTALL = :
+POST_INSTALL = :
+NORMAL_UNINSTALL = :
+PRE_UNINSTALL = :
+POST_UNINSTALL = :
+build_triplet = @build@
+host_triplet = @host@
+subdir = examples/meshing/surface_nurbs/subduction
+DIST_COMMON = README $(dist_noinst_DATA) $(srcdir)/Makefile.am \
+ $(srcdir)/Makefile.in
+ACLOCAL_M4 = $(top_srcdir)/aclocal.m4
+am__aclocal_m4_deps = $(top_srcdir)/m4/ac_pkg_swig.m4 \
+ $(top_srcdir)/m4/cit_cppunit.m4 $(top_srcdir)/m4/cit_cuda.m4 \
+ $(top_srcdir)/m4/cit_funcstring.m4 \
+ $(top_srcdir)/m4/cit_hdf5.m4 $(top_srcdir)/m4/cit_mpi.m4 \
+ $(top_srcdir)/m4/cit_netcdf.m4 $(top_srcdir)/m4/cit_numpy.m4 \
+ $(top_srcdir)/m4/cit_petsc.m4 $(top_srcdir)/m4/cit_proj4.m4 \
+ $(top_srcdir)/m4/cit_python.m4 \
+ $(top_srcdir)/m4/cit_spatialdata.m4 \
+ $(top_srcdir)/m4/libtool.m4 $(top_srcdir)/m4/ltoptions.m4 \
+ $(top_srcdir)/m4/ltsugar.m4 $(top_srcdir)/m4/ltversion.m4 \
+ $(top_srcdir)/m4/lt~obsolete.m4 $(top_srcdir)/configure.ac
+am__configure_deps = $(am__aclocal_m4_deps) $(CONFIGURE_DEPENDENCIES) \
+ $(ACLOCAL_M4)
+mkinstalldirs = $(install_sh) -d
+CONFIG_HEADER = $(top_builddir)/portinfo
+CONFIG_CLEAN_FILES =
+CONFIG_CLEAN_VPATH_FILES =
+SOURCES =
+DIST_SOURCES =
+DATA = $(dist_noinst_DATA)
+DISTFILES = $(DIST_COMMON) $(DIST_SOURCES) $(TEXINFOS) $(EXTRA_DIST)
+ACLOCAL = @ACLOCAL@
+AMTAR = @AMTAR@
+AM_LDFLAGS = @AM_LDFLAGS@
+AR = @AR@
+AUTOCONF = @AUTOCONF@
+AUTOHEADER = @AUTOHEADER@
+AUTOMAKE = @AUTOMAKE@
+AWK = @AWK@
+CC = @CC@
+CCDEPMODE = @CCDEPMODE@
+CFLAGS = @CFLAGS@
+CPP = @CPP@
+CPPFLAGS = @CPPFLAGS@
+CXX = @CXX@
+CXXCPP = @CXXCPP@
+CXXDEPMODE = @CXXDEPMODE@
+CXXFLAGS = @CXXFLAGS@
+CYGPATH_W = @CYGPATH_W@
+DEFS = @DEFS@
+DEPDIR = @DEPDIR@
+DLLTOOL = @DLLTOOL@
+DSYMUTIL = @DSYMUTIL@
+DUMPBIN = @DUMPBIN@
+ECHO_C = @ECHO_C@
+ECHO_N = @ECHO_N@
+ECHO_T = @ECHO_T@
+EGREP = @EGREP@
+EXEEXT = @EXEEXT@
+FGREP = @FGREP@
+GREP = @GREP@
+INSTALL = @INSTALL@
+INSTALL_DATA = @INSTALL_DATA@
+INSTALL_PROGRAM = @INSTALL_PROGRAM@
+INSTALL_SCRIPT = @INSTALL_SCRIPT@
+INSTALL_STRIP_PROGRAM = @INSTALL_STRIP_PROGRAM@
+LD = @LD@
+LDFLAGS = @LDFLAGS@
+LIBOBJS = @LIBOBJS@
+LIBS = @LIBS@
+LIBTOOL = @LIBTOOL@
+LIPO = @LIPO@
+LN_S = @LN_S@
+LTLIBOBJS = @LTLIBOBJS@
+MAKEINFO = @MAKEINFO@
+MANIFEST_TOOL = @MANIFEST_TOOL@
+MKDIR_P = @MKDIR_P@
+NM = @NM@
+NMEDIT = @NMEDIT@
+NUMPY_INCDIR = @NUMPY_INCDIR@
+NVCC = @NVCC@
+OBJDUMP = @OBJDUMP@
+OBJEXT = @OBJEXT@
+OTOOL = @OTOOL@
+OTOOL64 = @OTOOL64@
+PACKAGE = @PACKAGE@
+PACKAGE_BUGREPORT = @PACKAGE_BUGREPORT@
+PACKAGE_NAME = @PACKAGE_NAME@
+PACKAGE_STRING = @PACKAGE_STRING@
+PACKAGE_TARNAME = @PACKAGE_TARNAME@
+PACKAGE_URL = @PACKAGE_URL@
+PACKAGE_VERSION = @PACKAGE_VERSION@
+PATH_SEPARATOR = @PATH_SEPARATOR@
+PETSC_ARCH = @PETSC_ARCH@
+PETSC_CC = @PETSC_CC@
+PETSC_CC_INCLUDES = @PETSC_CC_INCLUDES@
+PETSC_CXX = @PETSC_CXX@
+PETSC_CXX_INCLUDE = @PETSC_CXX_INCLUDE@
+PETSC_CXX_LIB = @PETSC_CXX_LIB@
+PETSC_DIR = @PETSC_DIR@
+PETSC_FC = @PETSC_FC@
+PETSC_FC_INCLUDES = @PETSC_FC_INCLUDES@
+PETSC_FORTRAN_LIB = @PETSC_FORTRAN_LIB@
+PETSC_LIB = @PETSC_LIB@
+PETSC_MPI_INCLUDE = @PETSC_MPI_INCLUDE@
+PETSC_MPI_LIB = @PETSC_MPI_LIB@
+PETSC_SIEVE_FLAGS = @PETSC_SIEVE_FLAGS@
+PETSC_VERSION_MAJOR = @PETSC_VERSION_MAJOR@
+PETSC_VERSION_MINOR = @PETSC_VERSION_MINOR@
+PETSC_VERSION_SUBMINOR = @PETSC_VERSION_SUBMINOR@
+PYLITH_SWIG_CPPFLAGS = @PYLITH_SWIG_CPPFLAGS@
+PYTHON = @PYTHON@
+PYTHONPATH = @PYTHONPATH@
+PYTHON_BLDLIBRARY = @PYTHON_BLDLIBRARY@
+PYTHON_EGG_CFLAGS = @PYTHON_EGG_CFLAGS@
+PYTHON_EGG_CPPFLAGS = @PYTHON_EGG_CPPFLAGS@
+PYTHON_EGG_LDFLAGS = @PYTHON_EGG_LDFLAGS@
+PYTHON_EGG_LIBS = @PYTHON_EGG_LIBS@
+PYTHON_EGG_PYXFLAGS = @PYTHON_EGG_PYXFLAGS@
+PYTHON_EXEC_PREFIX = @PYTHON_EXEC_PREFIX@
+PYTHON_INCDIR = @PYTHON_INCDIR@
+PYTHON_LA_LDFLAGS = @PYTHON_LA_LDFLAGS@
+PYTHON_LDFLAGS = @PYTHON_LDFLAGS@
+PYTHON_LDLAST = @PYTHON_LDLAST@
+PYTHON_LDLIBRARY = @PYTHON_LDLIBRARY@
+PYTHON_LIBDIR = @PYTHON_LIBDIR@
+PYTHON_LIBP = @PYTHON_LIBP@
+PYTHON_LIBPL = @PYTHON_LIBPL@
+PYTHON_LIBS = @PYTHON_LIBS@
+PYTHON_LINKFORSHARED = @PYTHON_LINKFORSHARED@
+PYTHON_MODLIBS = @PYTHON_MODLIBS@
+PYTHON_PLATFORM = @PYTHON_PLATFORM@
+PYTHON_PREFIX = @PYTHON_PREFIX@
+PYTHON_SYSLIBS = @PYTHON_SYSLIBS@
+PYTHON_VERSION = @PYTHON_VERSION@
+RANLIB = @RANLIB@
+SED = @SED@
+SET_MAKE = @SET_MAKE@
+SHELL = @SHELL@
+STRIP = @STRIP@
+SWIG = @SWIG@
+SWIG_LIB = @SWIG_LIB@
+VERSION = @VERSION@
+abs_builddir = @abs_builddir@
+abs_srcdir = @abs_srcdir@
+abs_top_builddir = @abs_top_builddir@
+abs_top_srcdir = @abs_top_srcdir@
+ac_ct_AR = @ac_ct_AR@
+ac_ct_CC = @ac_ct_CC@
+ac_ct_CXX = @ac_ct_CXX@
+ac_ct_DUMPBIN = @ac_ct_DUMPBIN@
+am__include = @am__include@
+am__leading_dot = @am__leading_dot@
+am__quote = @am__quote@
+am__tar = @am__tar@
+am__untar = @am__untar@
+bindir = @bindir@
+build = @build@
+build_alias = @build_alias@
+build_cpu = @build_cpu@
+build_os = @build_os@
+build_vendor = @build_vendor@
+builddir = @builddir@
+datadir = @datadir@
+datarootdir = @datarootdir@
+docdir = @docdir@
+dvidir = @dvidir@
+exec_prefix = @exec_prefix@
+have_doxygen = @have_doxygen@
+host = @host@
+host_alias = @host_alias@
+host_cpu = @host_cpu@
+host_os = @host_os@
+host_vendor = @host_vendor@
+htmldir = @htmldir@
+includedir = @includedir@
+infodir = @infodir@
+install_sh = @install_sh@
+libdir = @libdir@
+libexecdir = @libexecdir@
+localedir = @localedir@
+localstatedir = @localstatedir@
+mandir = @mandir@
+mkdir_p = @mkdir_p@
+oldincludedir = @oldincludedir@
+pdfdir = @pdfdir@
+pkgpyexecdir = @pkgpyexecdir@
+pkgpythondir = @pkgpythondir@
+prefix = @prefix@
+program_transform_name = @program_transform_name@
+psdir = @psdir@
+pyexecdir = @pyexecdir@
+pythondir = @pythondir@
+sbindir = @sbindir@
+sharedstatedir = @sharedstatedir@
+srcdir = @srcdir@
+sysconfdir = @sysconfdir@
+target_alias = @target_alias@
+top_build_prefix = @top_build_prefix@
+top_builddir = @top_builddir@
+top_srcdir = @top_srcdir@
+dist_noinst_DATA = \
+ README \
+ geometry.jou \
+ mesh.jou \
+ bc.jou \
+ topobath_netsurf.py \
+ topobath_points.txt \
+ interface_netsurf.py \
+ interface_points.txt \
+ splay_skinsurf.py \
+ splay_points.txt
+
+all: all-am
+
+.SUFFIXES:
+$(srcdir)/Makefile.in: $(srcdir)/Makefile.am $(am__configure_deps)
+ @for dep in $?; do \
+ case '$(am__configure_deps)' in \
+ *$$dep*) \
+ ( cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh ) \
+ && { if test -f $@; then exit 0; else break; fi; }; \
+ exit 1;; \
+ esac; \
+ done; \
+ echo ' cd $(top_srcdir) && $(AUTOMAKE) --foreign examples/meshing/surface_nurbs/subduction/Makefile'; \
+ $(am__cd) $(top_srcdir) && \
+ $(AUTOMAKE) --foreign examples/meshing/surface_nurbs/subduction/Makefile
+.PRECIOUS: Makefile
+Makefile: $(srcdir)/Makefile.in $(top_builddir)/config.status
+ @case '$?' in \
+ *config.status*) \
+ cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh;; \
+ *) \
+ echo ' cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe)'; \
+ cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe);; \
+ esac;
+
+$(top_builddir)/config.status: $(top_srcdir)/configure $(CONFIG_STATUS_DEPENDENCIES)
+ cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
+
+$(top_srcdir)/configure: $(am__configure_deps)
+ cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
+$(ACLOCAL_M4): $(am__aclocal_m4_deps)
+ cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
+$(am__aclocal_m4_deps):
+
+mostlyclean-libtool:
+ -rm -f *.lo
+
+clean-libtool:
+ -rm -rf .libs _libs
+tags: TAGS
+TAGS:
+
+ctags: CTAGS
+CTAGS:
+
+
+distdir: $(DISTFILES)
+ @srcdirstrip=`echo "$(srcdir)" | sed 's/[].[^$$\\*]/\\\\&/g'`; \
+ topsrcdirstrip=`echo "$(top_srcdir)" | sed 's/[].[^$$\\*]/\\\\&/g'`; \
+ list='$(DISTFILES)'; \
+ dist_files=`for file in $$list; do echo $$file; done | \
+ sed -e "s|^$$srcdirstrip/||;t" \
+ -e "s|^$$topsrcdirstrip/|$(top_builddir)/|;t"`; \
+ case $$dist_files in \
+ */*) $(MKDIR_P) `echo "$$dist_files" | \
+ sed '/\//!d;s|^|$(distdir)/|;s,/[^/]*$$,,' | \
+ sort -u` ;; \
+ esac; \
+ for file in $$dist_files; do \
+ if test -f $$file || test -d $$file; then d=.; else d=$(srcdir); fi; \
+ if test -d $$d/$$file; then \
+ dir=`echo "/$$file" | sed -e 's,/[^/]*$$,,'`; \
+ if test -d "$(distdir)/$$file"; then \
+ find "$(distdir)/$$file" -type d ! -perm -700 -exec chmod u+rwx {} \;; \
+ fi; \
+ if test -d $(srcdir)/$$file && test $$d != $(srcdir); then \
+ cp -fpR $(srcdir)/$$file "$(distdir)$$dir" || exit 1; \
+ find "$(distdir)/$$file" -type d ! -perm -700 -exec chmod u+rwx {} \;; \
+ fi; \
+ cp -fpR $$d/$$file "$(distdir)$$dir" || exit 1; \
+ else \
+ test -f "$(distdir)/$$file" \
+ || cp -p $$d/$$file "$(distdir)/$$file" \
+ || exit 1; \
+ fi; \
+ done
+check-am: all-am
+check: check-am
+all-am: Makefile $(DATA)
+installdirs:
+install: install-am
+install-exec: install-exec-am
+install-data: install-data-am
+uninstall: uninstall-am
+
+install-am: all-am
+ @$(MAKE) $(AM_MAKEFLAGS) install-exec-am install-data-am
+
+installcheck: installcheck-am
+install-strip:
+ $(MAKE) $(AM_MAKEFLAGS) INSTALL_PROGRAM="$(INSTALL_STRIP_PROGRAM)" \
+ install_sh_PROGRAM="$(INSTALL_STRIP_PROGRAM)" INSTALL_STRIP_FLAG=-s \
+ `test -z '$(STRIP)' || \
+ echo "INSTALL_PROGRAM_ENV=STRIPPROG='$(STRIP)'"` install
+mostlyclean-generic:
+
+clean-generic:
+
+distclean-generic:
+ -test -z "$(CONFIG_CLEAN_FILES)" || rm -f $(CONFIG_CLEAN_FILES)
+ -test . = "$(srcdir)" || test -z "$(CONFIG_CLEAN_VPATH_FILES)" || rm -f $(CONFIG_CLEAN_VPATH_FILES)
+
+maintainer-clean-generic:
+ @echo "This command is intended for maintainers to use"
+ @echo "it deletes files that may require special tools to rebuild."
+clean: clean-am
+
+clean-am: clean-generic clean-libtool mostlyclean-am
+
+distclean: distclean-am
+ -rm -f Makefile
+distclean-am: clean-am distclean-generic
+
+dvi: dvi-am
+
+dvi-am:
+
+html: html-am
+
+html-am:
+
+info: info-am
+
+info-am:
+
+install-data-am:
+
+install-dvi: install-dvi-am
+
+install-dvi-am:
+
+install-exec-am:
+
+install-html: install-html-am
+
+install-html-am:
+
+install-info: install-info-am
+
+install-info-am:
+
+install-man:
+
+install-pdf: install-pdf-am
+
+install-pdf-am:
+
+install-ps: install-ps-am
+
+install-ps-am:
+
+installcheck-am:
+
+maintainer-clean: maintainer-clean-am
+ -rm -f Makefile
+maintainer-clean-am: distclean-am maintainer-clean-generic
+
+mostlyclean: mostlyclean-am
+
+mostlyclean-am: mostlyclean-generic mostlyclean-libtool
+
+pdf: pdf-am
+
+pdf-am:
+
+ps: ps-am
+
+ps-am:
+
+uninstall-am:
+
+.MAKE: install-am install-strip
+
+.PHONY: all all-am check check-am clean clean-generic clean-libtool \
+ distclean distclean-generic distclean-libtool distdir dvi \
+ dvi-am html html-am info info-am install install-am \
+ install-data install-data-am install-dvi install-dvi-am \
+ install-exec install-exec-am install-html install-html-am \
+ install-info install-info-am install-man install-pdf \
+ install-pdf-am install-ps install-ps-am install-strip \
+ installcheck installcheck-am installdirs maintainer-clean \
+ maintainer-clean-generic mostlyclean mostlyclean-generic \
+ mostlyclean-libtool pdf pdf-am ps ps-am uninstall uninstall-am
+
+
+# End of file
+
+# Tell versions [3.59,3.63) of GNU make to not export all variables.
+# Otherwise a system limit (for SysV at least) may be exceeded.
+.NOEXPORT:
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/README
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/README (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/README 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,51 @@
+CUBIT example illustrating methods for creating nonplanar surfaces in
+meshing a subduction interface with a steeply-dipping splay fault and
+topography/bathymetry.
+
+To generate the mesh for this example, we perform the following steps:
+
+1. Generate CUBIT journal files for creation of the surfaces:
+
+ a. Generate the subduction interface journal file by first
+ running the Python script interface_netsurf.py. This reads points
+ from interface_points.txt, which defines the interface geometry,
+ and creates a CUBIT journal file to create the surface using a net
+ surface.
+
+ python interface_netsurf.py
+
+ b. Generate the splay fault journal file by running the Python
+ script splay_skinsurf.py. This reads points from the file
+ splay_points.txt, which defines the splay fault geometry, and
+ creates a CUBIT journal file to create the surface using a skin
+ surface.
+
+ python splay_skinsurf.py
+
+ c. Generate the topography/bathymetry journal file by running the
+ Python script topobath_netsurf.py script. This reads points from
+ topobath_points.txt, which defines the topography/bathymetry, and
+ creates a CUBIT journal file to create the surface using a net
+ surface.
+
+ python splay_skinsurf.py
+
+2. Create the surfaces by running each of the three journal files
+produced by the Python scripts within CUBIT.
+
+ a. interface_netsurf.jou creates the interface surface.
+
+ b. splay_skinsurf.jou creates the splay fault surface.
+
+ c. topobath_netsurf.jou creates the DEM surface.
+
+
+3. Generate the mesh using these surfaces by running the
+ mesh.jou journal file within CUBIT. Note that this file also
+ makes use of the following journal files:
+
+ geometry.jou: Imports the surfaces created in step 2 and creates the
+ volumes needed for meshing.
+
+ bc.jou: Creates material blocks and nodesets needed by the
+ finite element model.
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/bc.jou
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/bc.jou (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/bc.jou 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,108 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+#
+# Charles A. Williams, GNS Science
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file with boundary condition info for M7.1 Darfield earthquake modeling.
+#
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1
+block 1 name "plate"
+block 2 volume 6
+block 2 name "slab"
+block 3 volume 7
+block 3 name "sediments"
+
+# ----------------------------------------------------------------------
+# Create nodeset for subduction interface
+# ----------------------------------------------------------------------
+group "subd_fault" add node in surface 33 38
+nodeset 100 group subd_fault
+nodeset 100 name "subd_fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for splay fault
+# ----------------------------------------------------------------------
+group "splay_fault" add node in surface 30
+
+# We must remove the nodes that are already contained in the subduction interface.
+group "splay_fault" remove node in group subd_fault
+
+nodeset 101 group splay_fault
+nodeset 101 name "splay_fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 12
+nodeset 10 group face_xpos
+nodeset 10 name "face xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 21 26
+nodeset 11 group face_xneg
+nodeset 11 name "face xneg"
+
+# To apply BC on face_xneg, we also create a nodeset that has fault nodes removed.
+group "face_xneg_nofault" add node in group face_xneg
+group "face_xneg_nofault" remove node in group subd_fault
+nodeset 16 group face_xneg_nofault
+nodeset 16 name "face_xneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 32 27 39
+nodeset 12 group face_ypos
+nodeset 12 name "face ypos"
+
+# To apply BC on face_ypos, we also create a nodeset that has fault nodes removed.
+group "face_ypos_nofault" add node in group face_ypos
+group "face_ypos_nofault" remove node in group subd_fault
+group "face_ypos_nofault" remove node in group splay_fault
+nodeset 17 group face_ypos_nofault
+nodeset 17 name "face_ypos_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 29 34 37
+nodeset 13 group face_yneg
+nodeset 13 name "face yneg"
+
+# To apply BC on face_yneg, we also create a nodeset that has fault nodes removed.
+group "face_yneg_nofault" add node in group face_yneg
+group "face_yneg_nofault" remove node in group subd_fault
+group "face_yneg_nofault" remove node in group splay_fault
+nodeset 18 group face_yneg_nofault
+nodeset 18 name "face_yneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 2
+nodeset 14 group face_zneg
+nodeset 14 name "face zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 28 36 31
+nodeset 15 group face_zpos
+nodeset 15 name "face zpos"
+
+
+# End of file
+
+
+
+
+
+
+
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/geometry.jou
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/geometry.jou (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/geometry.jou 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,89 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010-2013 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file with geometry for example showing how to import
+# several surfaces (topography, subduction interface, splay fault), and
+# mesh the resulting volumes.
+#
+# We also use the CUBIT support for APREPRO (the expressions in
+# braces), which includes units and mathematical functions. See the
+# APREPRO section in the appendix of the CUBIT documentation for more
+# information.
+#
+# ----------------------------------------------------------------------
+# Set units to SI.
+# ----------------------------------------------------------------------
+# {Units('si')}
+#
+# ----------------------------------------------------------------------
+# Reset geometry.
+# ----------------------------------------------------------------------
+reset
+
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+# Block is 250 km x 500 km x 62 km
+#{blockLength=250.0*km}
+#{blockWidth=500.0*km}
+#{blockHeight=80.0*km}
+
+brick x {blockLength} y {blockWidth} z {blockHeight}
+#{idVol=Id("volume")}
+
+# Translate block so top can be chopped with topography and so faults
+# intersect upper surface rather than side.
+#{moveX=50.0*km}
+#{moveY=0.0*km}
+
+volume {idVol} move x {moveX} y {moveY} z {-0.5*blockHeight+10.0*km}
+
+# ----------------------------------------------------------------------
+# Import topography, subduction interface, and splay fault.
+# ----------------------------------------------------------------------
+import Acis "topobath_surf.sat"
+import Acis "interface_surf.sat"
+import Acis "splay_surf.sat"
+
+# ----------------------------------------------------------------------
+# Webcut block with topography and delete volume lying above topography.
+# ----------------------------------------------------------------------
+webcut volume {idVol} with sheet surface 7
+delete volume 5
+delete body 2
+
+# ----------------------------------------------------------------------
+# Webcut remaining portion of block with subduction interface.
+# ----------------------------------------------------------------------
+webcut volume 1 with sheet surface 8
+delete body 3
+
+# ----------------------------------------------------------------------
+# Webcut volume above subduction interface with splay fault.
+# ----------------------------------------------------------------------
+webcut volume 1 with sheet surface 9
+delete body 4
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, then merge.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# End of file
+
+
+
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_netsurf.py
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_netsurf.py (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_netsurf.py 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,82 @@
+#!/usr/bin/env python
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010-2013 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+# PREREQUISITES: numpy
+
+# ======================================================================
+import numpy
+
+# Define parameters.
+interfaceFile = "interface_points.txt"
+numContours = 4
+pointsPerContour = 3
+journalFile = "interface_netsurf.jou"
+acisFile = "interface_surf.sat"
+
+# Journal file formatting, etc.
+separator = "# ----------------------------------------------------------\n"
+journalBeg = \
+ "# CUBIT journal file generated by interface_netsurf.py.\n" + \
+ "#\n" + \
+ "# Create an ACIS NURBS surface from intersecting lines.\n" + \
+ "#\n" + separator + \
+ "reset\n"
+lineBeg = "create curve spline"
+splineFmt = " location %10.2e %10.2e %10.2e"
+netFmt = "create surface net u curve %d to %d v curve %d to %d\n"
+delCmd = "delete curve all\n"
+expCmd = "export Acis '" + acisFile + "'\n"
+
+
+# Read coordinates and reshape them.
+intCoords = numpy.loadtxt(interfaceFile, dtype=numpy.float64).reshape(numContours, pointsPerContour, 3)
+
+j = open(journalFile, 'w')
+j.write(journalBeg)
+
+# Loop over contours (u-lines).
+for contour in range(numContours):
+ points = intCoords[contour,:,:]
+ j.write(lineBeg)
+ for pointNum in range(pointsPerContour):
+ point = points[pointNum,:]
+ j.write(splineFmt % (point[0], point[1], point[2]))
+
+ j.write("\n")
+
+# Loop over profiles (v-lines).
+for profile in range(pointsPerContour):
+ points = intCoords[:,profile,:]
+ j.write(lineBeg)
+ for pointNum in range(numContours):
+ point = points[pointNum,:]
+ j.write(splineFmt % (point[0], point[1], point[2]))
+
+ j.write("\n")
+
+# Create net surface.
+uline1 = 1
+uline2 = numContours
+vline1 = numContours + 1
+vline2 = vline1 + pointsPerContour - 1
+j.write(netFmt % (uline1, uline2, vline1, vline2))
+
+# Delete spline curves and export Acis file.
+j.write(delCmd)
+j.write(expCmd)
+j.close()
+
+# End of file
Property changes on: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_netsurf.py
___________________________________________________________________
Added: svn:executable
+ *
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_points.txt
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_points.txt (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/interface_points.txt 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,17 @@
+# Curvilinear grid of points organized as horizontal contours.
+#
+-158363.49995 -320000.00000 -80000.0
+-105810.54254 23.046240000 -80000.0
+-151349.45888 320000.00000 -80000.0
+#
+-125972.41361 -320000.00000 -60000.0
+-73054.343245 232.36802000 -60000.0
+-118727.49631 320000.00000 -60000.0
+#
+-36204.490870 -320000.00000 -20000.0
+18146.082654 -1571.7179170 -20000.0
+-28188.208958 320000.00000 -20000.0
+#
+34030.592403 -320000.00000 5000.0
+154341.36394 740.47119400 5000.0
+76062.567758 320000.00000 5000.0
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/mesh.jou
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/mesh.jou (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/mesh.jou 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,61 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010-2013 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file to create an ExodusII file with a mesh at uniform
+# resolution.
+#
+# ----------------------------------------------------------------------
+# Generate geometry.
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Create tet4 mesh at 6.0 km resolution.
+# ----------------------------------------------------------------------
+
+#{dx=6.0*km}
+volume all size {dx}
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Generate the mesh.
+# ----------------------------------------------------------------------
+mesh surface all
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth mesh to imporve quality.
+# ----------------------------------------------------------------------
+cleanup volume all
+volume all smooth scheme condition number beta 2.0 cpu 4
+smooth volume all
+#
+# ----------------------------------------------------------------------
+# Boundary conditions.
+# ----------------------------------------------------------------------
+playback 'bc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file.
+# ----------------------------------------------------------------------
+set large exodus off
+export mesh "mesh.exo" dimension 3 overwrite
+
+
+# End of file
+
+
+
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_points.txt
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_points.txt (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_points.txt 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,26 @@
+# Vertical profiles for splay fault geometry.
+#
+-60000.0 -320000.0 -30000.0
+-45000.0 -320000.0 -20000.0
+-20000.0 -320000.0 0.0
+ 5000.0 -320000.0 25000.0
+#
+-10000.0 -150000.0 -30000.0
+ 5000.0 -150000.0 -20000.0
+ 30000.0 -150000.0 0.0
+ 55000.0 -150000.0 25000.0
+#
+ 15000.0 0.0 -30000.0
+ 30000.0 0.0 -20000.0
+ 55000.0 0.0 0.0
+ 80000.0 0.0 25000.0
+#
+-10000.0 150000.0 -30000.0
+ 5000.0 150000.0 -20000.0
+ 30000.0 150000.0 0.0
+ 55000.0 150000.0 25000.0
+#
+-35000.0 320000.0 -30000.0
+-20000.0 320000.0 -20000.0
+ 5000.0 320000.0 0.0
+ 30000.0 320000.0 25000.0
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_skinsurf.py
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_skinsurf.py (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_skinsurf.py 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,68 @@
+#!/usr/bin/env python
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010-2013 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+# PREREQUISITES: numpy
+
+# ======================================================================
+import numpy
+
+# Define parameters.
+splayFile = "splay_points.txt"
+numProfiles = 5
+pointsPerProfile = 4
+journalFile = "splay_skinsurf.jou"
+acisFile = "splay_surf.sat"
+
+# Journal file formatting, etc.
+separator = "# ----------------------------------------------------------\n"
+journalBeg = \
+ "# CUBIT journal file generated by splay_skinsurf.py.\n" + \
+ "#\n" + \
+ "# Create an ACIS NURBS surface from a set of profiles.\n" + \
+ "#\n" + separator + \
+ "reset\n"
+lineBeg = "create curve spline"
+splineFmt = " location %10.2e %10.2e %10.2e"
+skinCmd = "create surface skin curve all\n"
+delCmd = "delete curve all\n"
+expCmd = "export Acis '" + acisFile + "'\n"
+
+
+# Read coordinates and reshape them.
+splayCoords = numpy.loadtxt(splayFile, dtype=numpy.float64).reshape(numProfiles, pointsPerProfile, 3)
+
+j = open(journalFile, 'w')
+j.write(journalBeg)
+
+# Loop over profiles.
+for profile in range(numProfiles):
+ points = splayCoords[profile,:,:]
+ j.write(lineBeg)
+ for pointNum in range(pointsPerProfile):
+ point = points[pointNum,:]
+ j.write(splineFmt % (point[0], point[1], point[2]))
+
+ j.write("\n")
+
+# Create skin surface.
+j.write(skinCmd)
+
+# Delete spline curves and export Acis file.
+j.write(delCmd)
+j.write(expCmd)
+j.close()
+
+# End of file
Property changes on: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/splay_skinsurf.py
___________________________________________________________________
Added: svn:executable
+ *
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_netsurf.py
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_netsurf.py (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_netsurf.py 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,82 @@
+#!/usr/bin/env python
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+# Charles A. Williams, GNS Science
+# Matthew G. Knepley, University of Chicago
+#
+# This code was developed as part of the Computational Infrastructure
+# for Geodynamics (http://geodynamics.org).
+#
+# Copyright (c) 2010-2013 University of California, Davis
+#
+# See COPYING for license information.
+#
+# ----------------------------------------------------------------------
+#
+# PREREQUISITES: numpy
+
+# ======================================================================
+import numpy
+
+# Define parameters.
+demFile = "topobath_points.txt"
+numProfiles = 4
+pointsPerProfile = 3
+journalFile = "topobath_netsurf.jou"
+acisFile = "topobath_surf.sat"
+
+# Journal file formatting, etc.
+separator = "# ----------------------------------------------------------\n"
+journalBeg = \
+ "# CUBIT journal file generated by topobath_netsurf.py.\n" + \
+ "#\n" + \
+ "# Create an ACIS NURBS surface from intersecting lines.\n" + \
+ "#\n" + separator + \
+ "reset\n"
+lineBeg = "create curve spline"
+splineFmt = " location %10.2e %10.2e %10.2e"
+netFmt = "create surface net u curve %d to %d v curve %d to %d\n"
+delCmd = "delete curve all\n"
+expCmd = "export Acis '" + acisFile + "'\n"
+
+
+# Read coordinates and reshape them.
+demCoords = numpy.loadtxt(demFile, dtype=numpy.float64).reshape(numProfiles, pointsPerProfile, 3)
+
+j = open(journalFile, 'w')
+j.write(journalBeg)
+
+# Loop over profiles (u-lines).
+for profile in range(numProfiles):
+ points = demCoords[profile,:,:]
+ j.write(lineBeg)
+ for pointNum in range(pointsPerProfile):
+ point = points[pointNum,:]
+ j.write(splineFmt % (point[0], point[1], point[2]))
+
+ j.write("\n")
+
+# Loop over contours (v-lines).
+for contour in range(pointsPerProfile):
+ points = demCoords[:,contour,:]
+ j.write(lineBeg)
+ for pointNum in range(numProfiles):
+ point = points[pointNum,:]
+ j.write(splineFmt % (point[0], point[1], point[2]))
+
+ j.write("\n")
+
+# Create net surface.
+uline1 = 1
+uline2 = numProfiles
+vline1 = numProfiles + 1
+vline2 = vline1 + pointsPerProfile - 1
+j.write(netFmt % (uline1, uline2, vline1, vline2))
+
+# Delete spline curves and export Acis file.
+j.write(delCmd)
+j.write(expCmd)
+j.close()
+
+# End of file
Property changes on: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_netsurf.py
___________________________________________________________________
Added: svn:executable
+ *
Added: short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_points.txt
===================================================================
--- short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_points.txt (rev 0)
+++ short/3D/PyLith/trunk/examples/meshing/surface_nurbs/subduction/topobath_points.txt 2013-06-19 21:39:59 UTC (rev 22378)
@@ -0,0 +1,17 @@
+# Grid of points to define DEM of topography/bathymetry.
+#
+-200000.0 -320000.0 3000.0
+ 0.0 -320000.0 2000.0
+ 200000.0 -320000.0 -8000.0
+#
+-200000.0 -110000.0 2000.0
+ 0.0 -110000.0 1000.0
+ 200000.0 -110000.0 -6000.0
+#
+-200000.0 110000.0 3000.0
+ 0.0 110000.0 -2000.0
+ 200000.0 110000.0 -9000.0
+#
+-200000.0 320000.0 1000.0
+ 0.0 320000.0 -3000.0
+ 200000.0 320000.0 -8000.0
More information about the CIG-COMMITS
mailing list