[cig-commits] commit: process input in separate module and account for command-line arguments
Mercurial
hg at geodynamics.org
Tue Sep 20 12:13:05 PDT 2011
changeset: 9:7596917b8776
user: Sylvain Barbot <sylbar.vainbot at gmail.com>
date: Tue Apr 26 15:00:49 2011 -0700
files: .bzrignore .doxygen elastic3d.f90 examples/run1.sh makefile_imkl relax.f90 types.f90
description:
process input in separate module and account for command-line arguments
diff -r 9992112ff05a -r 7596917b8776 .bzrignore
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/.bzrignore Tue Apr 26 15:00:49 2011 -0700
@@ -0,0 +1,5 @@
+*.o
+*.mod
+*~
+doc
+./relax
diff -r 9992112ff05a -r 7596917b8776 .doxygen
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/.doxygen Tue Apr 26 15:00:49 2011 -0700
@@ -0,0 +1,1417 @@
+# Doxyfile 1.5.6
+
+# This file describes the settings to be used by the documentation system
+# doxygen (www.doxygen.org) for a project
+#
+# All text after a hash (#) is considered a comment and will be ignored
+# The format is:
+# TAG = value [value, ...]
+# For lists items can also be appended using:
+# TAG += value [value, ...]
+# Values that contain spaces should be placed between quotes (" ")
+
+#---------------------------------------------------------------------------
+# Project related configuration options
+#---------------------------------------------------------------------------
+
+# This tag specifies the encoding used for all characters in the config file
+# that follow. The default is UTF-8 which is also the encoding used for all
+# text before the first occurrence of this tag. Doxygen uses libiconv (or the
+# iconv built into libc) for the transcoding. See
+# http://www.gnu.org/software/libiconv for the list of possible encodings.
+
+DOXYFILE_ENCODING = UTF-8
+
+# The PROJECT_NAME tag is a single word (or a sequence of words surrounded
+# by quotes) that should identify the project.
+
+PROJECT_NAME = RELAX
+
+# The PROJECT_NUMBER tag can be used to enter a project or revision number.
+# This could be handy for archiving the generated documentation or
+# if some version control system is used.
+
+PROJECT_NUMBER =
+
+# The OUTPUT_DIRECTORY tag is used to specify the (relative or absolute)
+# base path where the generated documentation will be put.
+# If a relative path is entered, it will be relative to the location
+# where doxygen was started. If left blank the current directory will be used.
+
+OUTPUT_DIRECTORY = doc
+
+# If the CREATE_SUBDIRS tag is set to YES, then doxygen will create
+# 4096 sub-directories (in 2 levels) under the output directory of each output
+# format and will distribute the generated files over these directories.
+# Enabling this option can be useful when feeding doxygen a huge amount of
+# source files, where putting all generated files in the same directory would
+# otherwise cause performance problems for the file system.
+
+CREATE_SUBDIRS = NO
+
+# The OUTPUT_LANGUAGE tag is used to specify the language in which all
+# documentation generated by doxygen is written. Doxygen will use this
+# information to generate all constant output in the proper language.
+# The default language is English, other supported languages are:
+# Afrikaans, Arabic, Brazilian, Catalan, Chinese, Chinese-Traditional,
+# Croatian, Czech, Danish, Dutch, Farsi, Finnish, French, German, Greek,
+# Hungarian, Italian, Japanese, Japanese-en (Japanese with English messages),
+# Korean, Korean-en, Lithuanian, Norwegian, Macedonian, Persian, Polish,
+# Portuguese, Romanian, Russian, Serbian, Slovak, Slovene, Spanish, Swedish,
+# and Ukrainian.
+
+OUTPUT_LANGUAGE = English
+
+# If the BRIEF_MEMBER_DESC tag is set to YES (the default) Doxygen will
+# include brief member descriptions after the members that are listed in
+# the file and class documentation (similar to JavaDoc).
+# Set to NO to disable this.
+
+BRIEF_MEMBER_DESC = YES
+
+# If the REPEAT_BRIEF tag is set to YES (the default) Doxygen will prepend
+# the brief description of a member or function before the detailed description.
+# Note: if both HIDE_UNDOC_MEMBERS and BRIEF_MEMBER_DESC are set to NO, the
+# brief descriptions will be completely suppressed.
+
+REPEAT_BRIEF = YES
+
+# This tag implements a quasi-intelligent brief description abbreviator
+# that is used to form the text in various listings. Each string
+# in this list, if found as the leading text of the brief description, will be
+# stripped from the text and the result after processing the whole list, is
+# used as the annotated text. Otherwise, the brief description is used as-is.
+# If left blank, the following values are used ("$name" is automatically
+# replaced with the name of the entity): "The $name class" "The $name widget"
+# "The $name file" "is" "provides" "specifies" "contains"
+# "represents" "a" "an" "the"
+
+ABBREVIATE_BRIEF =
+
+# If the ALWAYS_DETAILED_SEC and REPEAT_BRIEF tags are both set to YES then
+# Doxygen will generate a detailed section even if there is only a brief
+# description.
+
+ALWAYS_DETAILED_SEC = NO
+
+# If the INLINE_INHERITED_MEMB tag is set to YES, doxygen will show all
+# inherited members of a class in the documentation of that class as if those
+# members were ordinary class members. Constructors, destructors and assignment
+# operators of the base classes will not be shown.
+
+INLINE_INHERITED_MEMB = NO
+
+# If the FULL_PATH_NAMES tag is set to YES then Doxygen will prepend the full
+# path before files name in the file list and in the header files. If set
+# to NO the shortest path that makes the file name unique will be used.
+
+FULL_PATH_NAMES = YES
+
+# If the FULL_PATH_NAMES tag is set to YES then the STRIP_FROM_PATH tag
+# can be used to strip a user-defined part of the path. Stripping is
+# only done if one of the specified strings matches the left-hand part of
+# the path. The tag can be used to show relative paths in the file list.
+# If left blank the directory from which doxygen is run is used as the
+# path to strip.
+
+STRIP_FROM_PATH =
+
+# The STRIP_FROM_INC_PATH tag can be used to strip a user-defined part of
+# the path mentioned in the documentation of a class, which tells
+# the reader which header file to include in order to use a class.
+# If left blank only the name of the header file containing the class
+# definition is used. Otherwise one should specify the include paths that
+# are normally passed to the compiler using the -I flag.
+
+STRIP_FROM_INC_PATH =
+
+# If the SHORT_NAMES tag is set to YES, doxygen will generate much shorter
+# (but less readable) file names. This can be useful is your file systems
+# doesn't support long names like on DOS, Mac, or CD-ROM.
+
+SHORT_NAMES = NO
+
+# If the JAVADOC_AUTOBRIEF tag is set to YES then Doxygen
+# will interpret the first line (until the first dot) of a JavaDoc-style
+# comment as the brief description. If set to NO, the JavaDoc
+# comments will behave just like regular Qt-style comments
+# (thus requiring an explicit @brief command for a brief description.)
+
+JAVADOC_AUTOBRIEF = NO
+
+# If the QT_AUTOBRIEF tag is set to YES then Doxygen will
+# interpret the first line (until the first dot) of a Qt-style
+# comment as the brief description. If set to NO, the comments
+# will behave just like regular Qt-style comments (thus requiring
+# an explicit \brief command for a brief description.)
+
+QT_AUTOBRIEF = NO
+
+# The MULTILINE_CPP_IS_BRIEF tag can be set to YES to make Doxygen
+# treat a multi-line C++ special comment block (i.e. a block of //! or ///
+# comments) as a brief description. This used to be the default behaviour.
+# The new default is to treat a multi-line C++ comment block as a detailed
+# description. Set this tag to YES if you prefer the old behaviour instead.
+
+MULTILINE_CPP_IS_BRIEF = NO
+
+# If the DETAILS_AT_TOP tag is set to YES then Doxygen
+# will output the detailed description near the top, like JavaDoc.
+# If set to NO, the detailed description appears after the member
+# documentation.
+
+DETAILS_AT_TOP = NO
+
+# If the INHERIT_DOCS tag is set to YES (the default) then an undocumented
+# member inherits the documentation from any documented member that it
+# re-implements.
+
+INHERIT_DOCS = YES
+
+# If the SEPARATE_MEMBER_PAGES tag is set to YES, then doxygen will produce
+# a new page for each member. If set to NO, the documentation of a member will
+# be part of the file/class/namespace that contains it.
+
+SEPARATE_MEMBER_PAGES = NO
+
+# The TAB_SIZE tag can be used to set the number of spaces in a tab.
+# Doxygen uses this value to replace tabs by spaces in code fragments.
+
+TAB_SIZE = 8
+
+# This tag can be used to specify a number of aliases that acts
+# as commands in the documentation. An alias has the form "name=value".
+# For example adding "sideeffect=\par Side Effects:\n" will allow you to
+# put the command \sideeffect (or @sideeffect) in the documentation, which
+# will result in a user-defined paragraph with heading "Side Effects:".
+# You can put \n's in the value part of an alias to insert newlines.
+
+ALIASES =
+
+# Set the OPTIMIZE_OUTPUT_FOR_C tag to YES if your project consists of C
+# sources only. Doxygen will then generate output that is more tailored for C.
+# For instance, some of the names that are used will be different. The list
+# of all members will be omitted, etc.
+
+OPTIMIZE_OUTPUT_FOR_C = NO
+
+# Set the OPTIMIZE_OUTPUT_JAVA tag to YES if your project consists of Java
+# sources only. Doxygen will then generate output that is more tailored for
+# Java. For instance, namespaces will be presented as packages, qualified
+# scopes will look different, etc.
+
+OPTIMIZE_OUTPUT_JAVA = NO
+
+# Set the OPTIMIZE_FOR_FORTRAN tag to YES if your project consists of Fortran
+# sources only. Doxygen will then generate output that is more tailored for
+# Fortran.
+
+OPTIMIZE_FOR_FORTRAN = YES
+
+# Set the OPTIMIZE_OUTPUT_VHDL tag to YES if your project consists of VHDL
+# sources. Doxygen will then generate output that is tailored for
+# VHDL.
+
+OPTIMIZE_OUTPUT_VHDL = NO
+
+# If you use STL classes (i.e. std::string, std::vector, etc.) but do not want
+# to include (a tag file for) the STL sources as input, then you should
+# set this tag to YES in order to let doxygen match functions declarations and
+# definitions whose arguments contain STL classes (e.g. func(std::string); v.s.
+# func(std::string) {}). This also make the inheritance and collaboration
+# diagrams that involve STL classes more complete and accurate.
+
+BUILTIN_STL_SUPPORT = NO
+
+# If you use Microsoft's C++/CLI language, you should set this option to YES to
+# enable parsing support.
+
+CPP_CLI_SUPPORT = NO
+
+# Set the SIP_SUPPORT tag to YES if your project consists of sip sources only.
+# Doxygen will parse them like normal C++ but will assume all classes use public
+# instead of private inheritance when no explicit protection keyword is present.
+
+SIP_SUPPORT = NO
+
+# For Microsoft's IDL there are propget and propput attributes to indicate getter
+# and setter methods for a property. Setting this option to YES (the default)
+# will make doxygen to replace the get and set methods by a property in the
+# documentation. This will only work if the methods are indeed getting or
+# setting a simple type. If this is not the case, or you want to show the
+# methods anyway, you should set this option to NO.
+
+IDL_PROPERTY_SUPPORT = YES
+
+# If member grouping is used in the documentation and the DISTRIBUTE_GROUP_DOC
+# tag is set to YES, then doxygen will reuse the documentation of the first
+# member in the group (if any) for the other members of the group. By default
+# all members of a group must be documented explicitly.
+
+DISTRIBUTE_GROUP_DOC = NO
+
+# Set the SUBGROUPING tag to YES (the default) to allow class member groups of
+# the same type (for instance a group of public functions) to be put as a
+# subgroup of that type (e.g. under the Public Functions section). Set it to
+# NO to prevent subgrouping. Alternatively, this can be done per class using
+# the \nosubgrouping command.
+
+SUBGROUPING = YES
+
+# When TYPEDEF_HIDES_STRUCT is enabled, a typedef of a struct, union, or enum
+# is documented as struct, union, or enum with the name of the typedef. So
+# typedef struct TypeS {} TypeT, will appear in the documentation as a struct
+# with name TypeT. When disabled the typedef will appear as a member of a file,
+# namespace, or class. And the struct will be named TypeS. This can typically
+# be useful for C code in case the coding convention dictates that all compound
+# types are typedef'ed and only the typedef is referenced, never the tag name.
+
+TYPEDEF_HIDES_STRUCT = NO
+
+#---------------------------------------------------------------------------
+# Build related configuration options
+#---------------------------------------------------------------------------
+
+# If the EXTRACT_ALL tag is set to YES doxygen will assume all entities in
+# documentation are documented, even if no documentation was available.
+# Private class members and static file members will be hidden unless
+# the EXTRACT_PRIVATE and EXTRACT_STATIC tags are set to YES
+
+EXTRACT_ALL = YES
+
+# If the EXTRACT_PRIVATE tag is set to YES all private members of a class
+# will be included in the documentation.
+
+EXTRACT_PRIVATE = NO
+
+# If the EXTRACT_STATIC tag is set to YES all static members of a file
+# will be included in the documentation.
+
+EXTRACT_STATIC = NO
+
+# If the EXTRACT_LOCAL_CLASSES tag is set to YES classes (and structs)
+# defined locally in source files will be included in the documentation.
+# If set to NO only classes defined in header files are included.
+
+EXTRACT_LOCAL_CLASSES = YES
+
+# This flag is only useful for Objective-C code. When set to YES local
+# methods, which are defined in the implementation section but not in
+# the interface are included in the documentation.
+# If set to NO (the default) only methods in the interface are included.
+
+EXTRACT_LOCAL_METHODS = NO
+
+# If this flag is set to YES, the members of anonymous namespaces will be
+# extracted and appear in the documentation as a namespace called
+# 'anonymous_namespace{file}', where file will be replaced with the base
+# name of the file that contains the anonymous namespace. By default
+# anonymous namespace are hidden.
+
+EXTRACT_ANON_NSPACES = NO
+
+# If the HIDE_UNDOC_MEMBERS tag is set to YES, Doxygen will hide all
+# undocumented members of documented classes, files or namespaces.
+# If set to NO (the default) these members will be included in the
+# various overviews, but no documentation section is generated.
+# This option has no effect if EXTRACT_ALL is enabled.
+
+HIDE_UNDOC_MEMBERS = NO
+
+# If the HIDE_UNDOC_CLASSES tag is set to YES, Doxygen will hide all
+# undocumented classes that are normally visible in the class hierarchy.
+# If set to NO (the default) these classes will be included in the various
+# overviews. This option has no effect if EXTRACT_ALL is enabled.
+
+HIDE_UNDOC_CLASSES = NO
+
+# If the HIDE_FRIEND_COMPOUNDS tag is set to YES, Doxygen will hide all
+# friend (class|struct|union) declarations.
+# If set to NO (the default) these declarations will be included in the
+# documentation.
+
+HIDE_FRIEND_COMPOUNDS = NO
+
+# If the HIDE_IN_BODY_DOCS tag is set to YES, Doxygen will hide any
+# documentation blocks found inside the body of a function.
+# If set to NO (the default) these blocks will be appended to the
+# function's detailed documentation block.
+
+HIDE_IN_BODY_DOCS = NO
+
+# The INTERNAL_DOCS tag determines if documentation
+# that is typed after a \internal command is included. If the tag is set
+# to NO (the default) then the documentation will be excluded.
+# Set it to YES to include the internal documentation.
+
+INTERNAL_DOCS = NO
+
+# If the CASE_SENSE_NAMES tag is set to NO then Doxygen will only generate
+# file names in lower-case letters. If set to YES upper-case letters are also
+# allowed. This is useful if you have classes or files whose names only differ
+# in case and if your file system supports case sensitive file names. Windows
+# and Mac users are advised to set this option to NO.
+
+CASE_SENSE_NAMES = NO
+
+# If the HIDE_SCOPE_NAMES tag is set to NO (the default) then Doxygen
+# will show members with their full class and namespace scopes in the
+# documentation. If set to YES the scope will be hidden.
+
+HIDE_SCOPE_NAMES = NO
+
+# If the SHOW_INCLUDE_FILES tag is set to YES (the default) then Doxygen
+# will put a list of the files that are included by a file in the documentation
+# of that file.
+
+SHOW_INCLUDE_FILES = YES
+
+# If the INLINE_INFO tag is set to YES (the default) then a tag [inline]
+# is inserted in the documentation for inline members.
+
+INLINE_INFO = YES
+
+# If the SORT_MEMBER_DOCS tag is set to YES (the default) then doxygen
+# will sort the (detailed) documentation of file and class members
+# alphabetically by member name. If set to NO the members will appear in
+# declaration order.
+
+SORT_MEMBER_DOCS = YES
+
+# If the SORT_BRIEF_DOCS tag is set to YES then doxygen will sort the
+# brief documentation of file, namespace and class members alphabetically
+# by member name. If set to NO (the default) the members will appear in
+# declaration order.
+
+SORT_BRIEF_DOCS = NO
+
+# If the SORT_GROUP_NAMES tag is set to YES then doxygen will sort the
+# hierarchy of group names into alphabetical order. If set to NO (the default)
+# the group names will appear in their defined order.
+
+SORT_GROUP_NAMES = NO
+
+# If the SORT_BY_SCOPE_NAME tag is set to YES, the class list will be
+# sorted by fully-qualified names, including namespaces. If set to
+# NO (the default), the class list will be sorted only by class name,
+# not including the namespace part.
+# Note: This option is not very useful if HIDE_SCOPE_NAMES is set to YES.
+# Note: This option applies only to the class list, not to the
+# alphabetical list.
+
+SORT_BY_SCOPE_NAME = NO
+
+# The GENERATE_TODOLIST tag can be used to enable (YES) or
+# disable (NO) the todo list. This list is created by putting \todo
+# commands in the documentation.
+
+GENERATE_TODOLIST = YES
+
+# The GENERATE_TESTLIST tag can be used to enable (YES) or
+# disable (NO) the test list. This list is created by putting \test
+# commands in the documentation.
+
+GENERATE_TESTLIST = YES
+
+# The GENERATE_BUGLIST tag can be used to enable (YES) or
+# disable (NO) the bug list. This list is created by putting \bug
+# commands in the documentation.
+
+GENERATE_BUGLIST = YES
+
+# The GENERATE_DEPRECATEDLIST tag can be used to enable (YES) or
+# disable (NO) the deprecated list. This list is created by putting
+# \deprecated commands in the documentation.
+
+GENERATE_DEPRECATEDLIST= YES
+
+# The ENABLED_SECTIONS tag can be used to enable conditional
+# documentation sections, marked by \if sectionname ... \endif.
+
+ENABLED_SECTIONS =
+
+# The MAX_INITIALIZER_LINES tag determines the maximum number of lines
+# the initial value of a variable or define consists of for it to appear in
+# the documentation. If the initializer consists of more lines than specified
+# here it will be hidden. Use a value of 0 to hide initializers completely.
+# The appearance of the initializer of individual variables and defines in the
+# documentation can be controlled using \showinitializer or \hideinitializer
+# command in the documentation regardless of this setting.
+
+MAX_INITIALIZER_LINES = 30
+
+# Set the SHOW_USED_FILES tag to NO to disable the list of files generated
+# at the bottom of the documentation of classes and structs. If set to YES the
+# list will mention the files that were used to generate the documentation.
+
+SHOW_USED_FILES = YES
+
+# If the sources in your project are distributed over multiple directories
+# then setting the SHOW_DIRECTORIES tag to YES will show the directory hierarchy
+# in the documentation. The default is NO.
+
+SHOW_DIRECTORIES = NO
+
+# Set the SHOW_FILES tag to NO to disable the generation of the Files page.
+# This will remove the Files entry from the Quick Index and from the
+# Folder Tree View (if specified). The default is YES.
+
+SHOW_FILES = YES
+
+# Set the SHOW_NAMESPACES tag to NO to disable the generation of the
+# Namespaces page. This will remove the Namespaces entry from the Quick Index
+# and from the Folder Tree View (if specified). The default is YES.
+
+SHOW_NAMESPACES = YES
+
+# The FILE_VERSION_FILTER tag can be used to specify a program or script that
+# doxygen should invoke to get the current version for each file (typically from
+# the version control system). Doxygen will invoke the program by executing (via
+# popen()) the command <command> <input-file>, where <command> is the value of
+# the FILE_VERSION_FILTER tag, and <input-file> is the name of an input file
+# provided by doxygen. Whatever the program writes to standard output
+# is used as the file version. See the manual for examples.
+
+FILE_VERSION_FILTER =
+
+#---------------------------------------------------------------------------
+# configuration options related to warning and progress messages
+#---------------------------------------------------------------------------
+
+# The QUIET tag can be used to turn on/off the messages that are generated
+# by doxygen. Possible values are YES and NO. If left blank NO is used.
+
+QUIET = NO
+
+# The WARNINGS tag can be used to turn on/off the warning messages that are
+# generated by doxygen. Possible values are YES and NO. If left blank
+# NO is used.
+
+WARNINGS = YES
+
+# If WARN_IF_UNDOCUMENTED is set to YES, then doxygen will generate warnings
+# for undocumented members. If EXTRACT_ALL is set to YES then this flag will
+# automatically be disabled.
+
+WARN_IF_UNDOCUMENTED = YES
+
+# If WARN_IF_DOC_ERROR is set to YES, doxygen will generate warnings for
+# potential errors in the documentation, such as not documenting some
+# parameters in a documented function, or documenting parameters that
+# don't exist or using markup commands wrongly.
+
+WARN_IF_DOC_ERROR = YES
+
+# This WARN_NO_PARAMDOC option can be abled to get warnings for
+# functions that are documented, but have no documentation for their parameters
+# or return value. If set to NO (the default) doxygen will only warn about
+# wrong or incomplete parameter documentation, but not about the absence of
+# documentation.
+
+WARN_NO_PARAMDOC = NO
+
+# The WARN_FORMAT tag determines the format of the warning messages that
+# doxygen can produce. The string should contain the $file, $line, and $text
+# tags, which will be replaced by the file and line number from which the
+# warning originated and the warning text. Optionally the format may contain
+# $version, which will be replaced by the version of the file (if it could
+# be obtained via FILE_VERSION_FILTER)
+
+WARN_FORMAT = "$file:$line: $text"
+
+# The WARN_LOGFILE tag can be used to specify a file to which warning
+# and error messages should be written. If left blank the output is written
+# to stderr.
+
+WARN_LOGFILE =
+
+#---------------------------------------------------------------------------
+# configuration options related to the input files
+#---------------------------------------------------------------------------
+
+# The INPUT tag can be used to specify the files and/or directories that contain
+# documented source files. You may enter file names like "myfile.cpp" or
+# directories like "/usr/src/myproject". Separate the files or directories
+# with spaces.
+
+INPUT = .
+
+# This tag can be used to specify the character encoding of the source files
+# that doxygen parses. Internally doxygen uses the UTF-8 encoding, which is
+# also the default input encoding. Doxygen uses libiconv (or the iconv built
+# into libc) for the transcoding. See http://www.gnu.org/software/libiconv for
+# the list of possible encodings.
+
+INPUT_ENCODING = UTF-8
+
+# If the value of the INPUT tag contains directories, you can use the
+# FILE_PATTERNS tag to specify one or more wildcard pattern (like *.cpp
+# and *.h) to filter out the source-files in the directories. If left
+# blank the following patterns are tested:
+# *.c *.cc *.cxx *.cpp *.c++ *.java *.ii *.ixx *.ipp *.i++ *.inl *.h *.hh *.hxx
+# *.hpp *.h++ *.idl *.odl *.cs *.php *.php3 *.inc *.m *.mm *.py *.f90
+
+FILE_PATTERNS = *.f90
+
+# The RECURSIVE tag can be used to turn specify whether or not subdirectories
+# should be searched for input files as well. Possible values are YES and NO.
+# If left blank NO is used.
+
+RECURSIVE = NO
+
+# The EXCLUDE tag can be used to specify files and/or directories that should
+# excluded from the INPUT source files. This way you can easily exclude a
+# subdirectory from a directory tree whose root is specified with the INPUT tag.
+
+EXCLUDE = ctfft.f mkl_dfti.f90
+
+# The EXCLUDE_SYMLINKS tag can be used select whether or not files or
+# directories that are symbolic links (a Unix filesystem feature) are excluded
+# from the input.
+
+EXCLUDE_SYMLINKS = NO
+
+# If the value of the INPUT tag contains directories, you can use the
+# EXCLUDE_PATTERNS tag to specify one or more wildcard patterns to exclude
+# certain files from those directories. Note that the wildcards are matched
+# against the file with absolute path, so to exclude all test directories
+# for example use the pattern */test/*
+
+EXCLUDE_PATTERNS =
+
+# The EXCLUDE_SYMBOLS tag can be used to specify one or more symbol names
+# (namespaces, classes, functions, etc.) that should be excluded from the
+# output. The symbol name can be a fully qualified name, a word, or if the
+# wildcard * is used, a substring. Examples: ANamespace, AClass,
+# AClass::ANamespace, ANamespace::*Test
+
+EXCLUDE_SYMBOLS =
+
+# The EXAMPLE_PATH tag can be used to specify one or more files or
+# directories that contain example code fragments that are included (see
+# the \include command).
+
+EXAMPLE_PATH =
+
+# If the value of the EXAMPLE_PATH tag contains directories, you can use the
+# EXAMPLE_PATTERNS tag to specify one or more wildcard pattern (like *.cpp
+# and *.h) to filter out the source-files in the directories. If left
+# blank all files are included.
+
+EXAMPLE_PATTERNS =
+
+# If the EXAMPLE_RECURSIVE tag is set to YES then subdirectories will be
+# searched for input files to be used with the \include or \dontinclude
+# commands irrespective of the value of the RECURSIVE tag.
+# Possible values are YES and NO. If left blank NO is used.
+
+EXAMPLE_RECURSIVE = NO
+
+# The IMAGE_PATH tag can be used to specify one or more files or
+# directories that contain image that are included in the documentation (see
+# the \image command).
+
+IMAGE_PATH =
+
+# The INPUT_FILTER tag can be used to specify a program that doxygen should
+# invoke to filter for each input file. Doxygen will invoke the filter program
+# by executing (via popen()) the command <filter> <input-file>, where <filter>
+# is the value of the INPUT_FILTER tag, and <input-file> is the name of an
+# input file. Doxygen will then use the output that the filter program writes
+# to standard output. If FILTER_PATTERNS is specified, this tag will be
+# ignored.
+
+INPUT_FILTER =
+
+# The FILTER_PATTERNS tag can be used to specify filters on a per file pattern
+# basis. Doxygen will compare the file name with each pattern and apply the
+# filter if there is a match. The filters are a list of the form:
+# pattern=filter (like *.cpp=my_cpp_filter). See INPUT_FILTER for further
+# info on how filters are used. If FILTER_PATTERNS is empty, INPUT_FILTER
+# is applied to all files.
+
+FILTER_PATTERNS =
+
+# If the FILTER_SOURCE_FILES tag is set to YES, the input filter (if set using
+# INPUT_FILTER) will be used to filter the input files when producing source
+# files to browse (i.e. when SOURCE_BROWSER is set to YES).
+
+FILTER_SOURCE_FILES = NO
+
+#---------------------------------------------------------------------------
+# configuration options related to source browsing
+#---------------------------------------------------------------------------
+
+# If the SOURCE_BROWSER tag is set to YES then a list of source files will
+# be generated. Documented entities will be cross-referenced with these sources.
+# Note: To get rid of all source code in the generated output, make sure also
+# VERBATIM_HEADERS is set to NO.
+
+SOURCE_BROWSER = NO
+
+# Setting the INLINE_SOURCES tag to YES will include the body
+# of functions and classes directly in the documentation.
+
+INLINE_SOURCES = NO
+
+# Setting the STRIP_CODE_COMMENTS tag to YES (the default) will instruct
+# doxygen to hide any special comment blocks from generated source code
+# fragments. Normal C and C++ comments will always remain visible.
+
+STRIP_CODE_COMMENTS = YES
+
+# If the REFERENCED_BY_RELATION tag is set to YES
+# then for each documented function all documented
+# functions referencing it will be listed.
+
+REFERENCED_BY_RELATION = NO
+
+# If the REFERENCES_RELATION tag is set to YES
+# then for each documented function all documented entities
+# called/used by that function will be listed.
+
+REFERENCES_RELATION = NO
+
+# If the REFERENCES_LINK_SOURCE tag is set to YES (the default)
+# and SOURCE_BROWSER tag is set to YES, then the hyperlinks from
+# functions in REFERENCES_RELATION and REFERENCED_BY_RELATION lists will
+# link to the source code. Otherwise they will link to the documentstion.
+
+REFERENCES_LINK_SOURCE = YES
+
+# If the USE_HTAGS tag is set to YES then the references to source code
+# will point to the HTML generated by the htags(1) tool instead of doxygen
+# built-in source browser. The htags tool is part of GNU's global source
+# tagging system (see http://www.gnu.org/software/global/global.html). You
+# will need version 4.8.6 or higher.
+
+USE_HTAGS = NO
+
+# If the VERBATIM_HEADERS tag is set to YES (the default) then Doxygen
+# will generate a verbatim copy of the header file for each class for
+# which an include is specified. Set to NO to disable this.
+
+VERBATIM_HEADERS = YES
+
+#---------------------------------------------------------------------------
+# configuration options related to the alphabetical class index
+#---------------------------------------------------------------------------
+
+# If the ALPHABETICAL_INDEX tag is set to YES, an alphabetical index
+# of all compounds will be generated. Enable this if the project
+# contains a lot of classes, structs, unions or interfaces.
+
+ALPHABETICAL_INDEX = NO
+
+# If the alphabetical index is enabled (see ALPHABETICAL_INDEX) then
+# the COLS_IN_ALPHA_INDEX tag can be used to specify the number of columns
+# in which this list will be split (can be a number in the range [1..20])
+
+COLS_IN_ALPHA_INDEX = 5
+
+# In case all classes in a project start with a common prefix, all
+# classes will be put under the same header in the alphabetical index.
+# The IGNORE_PREFIX tag can be used to specify one or more prefixes that
+# should be ignored while generating the index headers.
+
+IGNORE_PREFIX =
+
+#---------------------------------------------------------------------------
+# configuration options related to the HTML output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_HTML tag is set to YES (the default) Doxygen will
+# generate HTML output.
+
+GENERATE_HTML = YES
+
+# The HTML_OUTPUT tag is used to specify where the HTML docs will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `html' will be used as the default path.
+
+HTML_OUTPUT = html
+
+# The HTML_FILE_EXTENSION tag can be used to specify the file extension for
+# each generated HTML page (for example: .htm,.php,.asp). If it is left blank
+# doxygen will generate files with .html extension.
+
+HTML_FILE_EXTENSION = .html
+
+# The HTML_HEADER tag can be used to specify a personal HTML header for
+# each generated HTML page. If it is left blank doxygen will generate a
+# standard header.
+
+HTML_HEADER =
+
+# The HTML_FOOTER tag can be used to specify a personal HTML footer for
+# each generated HTML page. If it is left blank doxygen will generate a
+# standard footer.
+
+HTML_FOOTER =
+
+# The HTML_STYLESHEET tag can be used to specify a user-defined cascading
+# style sheet that is used by each HTML page. It can be used to
+# fine-tune the look of the HTML output. If the tag is left blank doxygen
+# will generate a default style sheet. Note that doxygen will try to copy
+# the style sheet file to the HTML output directory, so don't put your own
+# stylesheet in the HTML output directory as well, or it will be erased!
+
+HTML_STYLESHEET =
+
+# If the HTML_ALIGN_MEMBERS tag is set to YES, the members of classes,
+# files or namespaces will be aligned in HTML using tables. If set to
+# NO a bullet list will be used.
+
+HTML_ALIGN_MEMBERS = YES
+
+# If the GENERATE_HTMLHELP tag is set to YES, additional index files
+# will be generated that can be used as input for tools like the
+# Microsoft HTML help workshop to generate a compiled HTML help file (.chm)
+# of the generated HTML documentation.
+
+GENERATE_HTMLHELP = NO
+
+# If the GENERATE_DOCSET tag is set to YES, additional index files
+# will be generated that can be used as input for Apple's Xcode 3
+# integrated development environment, introduced with OSX 10.5 (Leopard).
+# To create a documentation set, doxygen will generate a Makefile in the
+# HTML output directory. Running make will produce the docset in that
+# directory and running "make install" will install the docset in
+# ~/Library/Developer/Shared/Documentation/DocSets so that Xcode will find
+# it at startup.
+
+GENERATE_DOCSET = NO
+
+# When GENERATE_DOCSET tag is set to YES, this tag determines the name of the
+# feed. A documentation feed provides an umbrella under which multiple
+# documentation sets from a single provider (such as a company or product suite)
+# can be grouped.
+
+DOCSET_FEEDNAME = "Doxygen generated docs"
+
+# When GENERATE_DOCSET tag is set to YES, this tag specifies a string that
+# should uniquely identify the documentation set bundle. This should be a
+# reverse domain-name style string, e.g. com.mycompany.MyDocSet. Doxygen
+# will append .docset to the name.
+
+DOCSET_BUNDLE_ID = org.doxygen.Project
+
+# If the HTML_DYNAMIC_SECTIONS tag is set to YES then the generated HTML
+# documentation will contain sections that can be hidden and shown after the
+# page has loaded. For this to work a browser that supports
+# JavaScript and DHTML is required (for instance Mozilla 1.0+, Firefox
+# Netscape 6.0+, Internet explorer 5.0+, Konqueror, or Safari).
+
+HTML_DYNAMIC_SECTIONS = NO
+
+# If the GENERATE_HTMLHELP tag is set to YES, the CHM_FILE tag can
+# be used to specify the file name of the resulting .chm file. You
+# can add a path in front of the file if the result should not be
+# written to the html output directory.
+
+CHM_FILE =
+
+# If the GENERATE_HTMLHELP tag is set to YES, the HHC_LOCATION tag can
+# be used to specify the location (absolute path including file name) of
+# the HTML help compiler (hhc.exe). If non-empty doxygen will try to run
+# the HTML help compiler on the generated index.hhp.
+
+HHC_LOCATION =
+
+# If the GENERATE_HTMLHELP tag is set to YES, the GENERATE_CHI flag
+# controls if a separate .chi index file is generated (YES) or that
+# it should be included in the master .chm file (NO).
+
+GENERATE_CHI = NO
+
+# If the GENERATE_HTMLHELP tag is set to YES, the CHM_INDEX_ENCODING
+# is used to encode HtmlHelp index (hhk), content (hhc) and project file
+# content.
+
+CHM_INDEX_ENCODING =
+
+# If the GENERATE_HTMLHELP tag is set to YES, the BINARY_TOC flag
+# controls whether a binary table of contents is generated (YES) or a
+# normal table of contents (NO) in the .chm file.
+
+BINARY_TOC = NO
+
+# The TOC_EXPAND flag can be set to YES to add extra items for group members
+# to the contents of the HTML help documentation and to the tree view.
+
+TOC_EXPAND = NO
+
+# The DISABLE_INDEX tag can be used to turn on/off the condensed index at
+# top of each HTML page. The value NO (the default) enables the index and
+# the value YES disables it.
+
+DISABLE_INDEX = NO
+
+# This tag can be used to set the number of enum values (range [1..20])
+# that doxygen will group on one line in the generated HTML documentation.
+
+ENUM_VALUES_PER_LINE = 4
+
+# The GENERATE_TREEVIEW tag is used to specify whether a tree-like index
+# structure should be generated to display hierarchical information.
+# If the tag value is set to FRAME, a side panel will be generated
+# containing a tree-like index structure (just like the one that
+# is generated for HTML Help). For this to work a browser that supports
+# JavaScript, DHTML, CSS and frames is required (for instance Mozilla 1.0+,
+# Netscape 6.0+, Internet explorer 5.0+, or Konqueror). Windows users are
+# probably better off using the HTML help feature. Other possible values
+# for this tag are: HIERARCHIES, which will generate the Groups, Directories,
+# and Class Hiererachy pages using a tree view instead of an ordered list;
+# ALL, which combines the behavior of FRAME and HIERARCHIES; and NONE, which
+# disables this behavior completely. For backwards compatibility with previous
+# releases of Doxygen, the values YES and NO are equivalent to FRAME and NONE
+# respectively.
+
+GENERATE_TREEVIEW = NONE
+
+# If the treeview is enabled (see GENERATE_TREEVIEW) then this tag can be
+# used to set the initial width (in pixels) of the frame in which the tree
+# is shown.
+
+TREEVIEW_WIDTH = 250
+
+# Use this tag to change the font size of Latex formulas included
+# as images in the HTML documentation. The default is 10. Note that
+# when you change the font size after a successful doxygen run you need
+# to manually remove any form_*.png images from the HTML output directory
+# to force them to be regenerated.
+
+FORMULA_FONTSIZE = 10
+
+#---------------------------------------------------------------------------
+# configuration options related to the LaTeX output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_LATEX tag is set to YES (the default) Doxygen will
+# generate Latex output.
+
+GENERATE_LATEX = NO
+
+# The LATEX_OUTPUT tag is used to specify where the LaTeX docs will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `latex' will be used as the default path.
+
+LATEX_OUTPUT = latex
+
+# The LATEX_CMD_NAME tag can be used to specify the LaTeX command name to be
+# invoked. If left blank `latex' will be used as the default command name.
+
+LATEX_CMD_NAME = latex
+
+# The MAKEINDEX_CMD_NAME tag can be used to specify the command name to
+# generate index for LaTeX. If left blank `makeindex' will be used as the
+# default command name.
+
+MAKEINDEX_CMD_NAME = makeindex
+
+# If the COMPACT_LATEX tag is set to YES Doxygen generates more compact
+# LaTeX documents. This may be useful for small projects and may help to
+# save some trees in general.
+
+COMPACT_LATEX = NO
+
+# The PAPER_TYPE tag can be used to set the paper type that is used
+# by the printer. Possible values are: a4, a4wide, letter, legal and
+# executive. If left blank a4wide will be used.
+
+PAPER_TYPE = a4wide
+
+# The EXTRA_PACKAGES tag can be to specify one or more names of LaTeX
+# packages that should be included in the LaTeX output.
+
+EXTRA_PACKAGES =
+
+# The LATEX_HEADER tag can be used to specify a personal LaTeX header for
+# the generated latex document. The header should contain everything until
+# the first chapter. If it is left blank doxygen will generate a
+# standard header. Notice: only use this tag if you know what you are doing!
+
+LATEX_HEADER =
+
+# If the PDF_HYPERLINKS tag is set to YES, the LaTeX that is generated
+# is prepared for conversion to pdf (using ps2pdf). The pdf file will
+# contain links (just like the HTML output) instead of page references
+# This makes the output suitable for online browsing using a pdf viewer.
+
+PDF_HYPERLINKS = YES
+
+# If the USE_PDFLATEX tag is set to YES, pdflatex will be used instead of
+# plain latex in the generated Makefile. Set this option to YES to get a
+# higher quality PDF documentation.
+
+USE_PDFLATEX = YES
+
+# If the LATEX_BATCHMODE tag is set to YES, doxygen will add the \\batchmode.
+# command to the generated LaTeX files. This will instruct LaTeX to keep
+# running if errors occur, instead of asking the user for help.
+# This option is also used when generating formulas in HTML.
+
+LATEX_BATCHMODE = NO
+
+# If LATEX_HIDE_INDICES is set to YES then doxygen will not
+# include the index chapters (such as File Index, Compound Index, etc.)
+# in the output.
+
+LATEX_HIDE_INDICES = NO
+
+#---------------------------------------------------------------------------
+# configuration options related to the RTF output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_RTF tag is set to YES Doxygen will generate RTF output
+# The RTF output is optimized for Word 97 and may not look very pretty with
+# other RTF readers or editors.
+
+GENERATE_RTF = NO
+
+# The RTF_OUTPUT tag is used to specify where the RTF docs will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `rtf' will be used as the default path.
+
+RTF_OUTPUT = rtf
+
+# If the COMPACT_RTF tag is set to YES Doxygen generates more compact
+# RTF documents. This may be useful for small projects and may help to
+# save some trees in general.
+
+COMPACT_RTF = NO
+
+# If the RTF_HYPERLINKS tag is set to YES, the RTF that is generated
+# will contain hyperlink fields. The RTF file will
+# contain links (just like the HTML output) instead of page references.
+# This makes the output suitable for online browsing using WORD or other
+# programs which support those fields.
+# Note: wordpad (write) and others do not support links.
+
+RTF_HYPERLINKS = NO
+
+# Load stylesheet definitions from file. Syntax is similar to doxygen's
+# config file, i.e. a series of assignments. You only have to provide
+# replacements, missing definitions are set to their default value.
+
+RTF_STYLESHEET_FILE =
+
+# Set optional variables used in the generation of an rtf document.
+# Syntax is similar to doxygen's config file.
+
+RTF_EXTENSIONS_FILE =
+
+#---------------------------------------------------------------------------
+# configuration options related to the man page output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_MAN tag is set to YES (the default) Doxygen will
+# generate man pages
+
+GENERATE_MAN = NO
+
+# The MAN_OUTPUT tag is used to specify where the man pages will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `man' will be used as the default path.
+
+MAN_OUTPUT = man
+
+# The MAN_EXTENSION tag determines the extension that is added to
+# the generated man pages (default is the subroutine's section .3)
+
+MAN_EXTENSION = .3
+
+# If the MAN_LINKS tag is set to YES and Doxygen generates man output,
+# then it will generate one additional man file for each entity
+# documented in the real man page(s). These additional files
+# only source the real man page, but without them the man command
+# would be unable to find the correct page. The default is NO.
+
+MAN_LINKS = NO
+
+#---------------------------------------------------------------------------
+# configuration options related to the XML output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_XML tag is set to YES Doxygen will
+# generate an XML file that captures the structure of
+# the code including all documentation.
+
+GENERATE_XML = NO
+
+# The XML_OUTPUT tag is used to specify where the XML pages will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `xml' will be used as the default path.
+
+XML_OUTPUT = xml
+
+# The XML_SCHEMA tag can be used to specify an XML schema,
+# which can be used by a validating XML parser to check the
+# syntax of the XML files.
+
+XML_SCHEMA =
+
+# The XML_DTD tag can be used to specify an XML DTD,
+# which can be used by a validating XML parser to check the
+# syntax of the XML files.
+
+XML_DTD =
+
+# If the XML_PROGRAMLISTING tag is set to YES Doxygen will
+# dump the program listings (including syntax highlighting
+# and cross-referencing information) to the XML output. Note that
+# enabling this will significantly increase the size of the XML output.
+
+XML_PROGRAMLISTING = YES
+
+#---------------------------------------------------------------------------
+# configuration options for the AutoGen Definitions output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_AUTOGEN_DEF tag is set to YES Doxygen will
+# generate an AutoGen Definitions (see autogen.sf.net) file
+# that captures the structure of the code including all
+# documentation. Note that this feature is still experimental
+# and incomplete at the moment.
+
+GENERATE_AUTOGEN_DEF = NO
+
+#---------------------------------------------------------------------------
+# configuration options related to the Perl module output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_PERLMOD tag is set to YES Doxygen will
+# generate a Perl module file that captures the structure of
+# the code including all documentation. Note that this
+# feature is still experimental and incomplete at the
+# moment.
+
+GENERATE_PERLMOD = NO
+
+# If the PERLMOD_LATEX tag is set to YES Doxygen will generate
+# the necessary Makefile rules, Perl scripts and LaTeX code to be able
+# to generate PDF and DVI output from the Perl module output.
+
+PERLMOD_LATEX = NO
+
+# If the PERLMOD_PRETTY tag is set to YES the Perl module output will be
+# nicely formatted so it can be parsed by a human reader. This is useful
+# if you want to understand what is going on. On the other hand, if this
+# tag is set to NO the size of the Perl module output will be much smaller
+# and Perl will parse it just the same.
+
+PERLMOD_PRETTY = YES
+
+# The names of the make variables in the generated doxyrules.make file
+# are prefixed with the string contained in PERLMOD_MAKEVAR_PREFIX.
+# This is useful so different doxyrules.make files included by the same
+# Makefile don't overwrite each other's variables.
+
+PERLMOD_MAKEVAR_PREFIX =
+
+#---------------------------------------------------------------------------
+# Configuration options related to the preprocessor
+#---------------------------------------------------------------------------
+
+# If the ENABLE_PREPROCESSING tag is set to YES (the default) Doxygen will
+# evaluate all C-preprocessor directives found in the sources and include
+# files.
+
+ENABLE_PREPROCESSING = YES
+
+# If the MACRO_EXPANSION tag is set to YES Doxygen will expand all macro
+# names in the source code. If set to NO (the default) only conditional
+# compilation will be performed. Macro expansion can be done in a controlled
+# way by setting EXPAND_ONLY_PREDEF to YES.
+
+MACRO_EXPANSION = NO
+
+# If the EXPAND_ONLY_PREDEF and MACRO_EXPANSION tags are both set to YES
+# then the macro expansion is limited to the macros specified with the
+# PREDEFINED and EXPAND_AS_DEFINED tags.
+
+EXPAND_ONLY_PREDEF = NO
+
+# If the SEARCH_INCLUDES tag is set to YES (the default) the includes files
+# in the INCLUDE_PATH (see below) will be search if a #include is found.
+
+SEARCH_INCLUDES = YES
+
+# The INCLUDE_PATH tag can be used to specify one or more directories that
+# contain include files that are not input files but should be processed by
+# the preprocessor.
+
+INCLUDE_PATH =
+
+# You can use the INCLUDE_FILE_PATTERNS tag to specify one or more wildcard
+# patterns (like *.h and *.hpp) to filter out the header-files in the
+# directories. If left blank, the patterns specified with FILE_PATTERNS will
+# be used.
+
+INCLUDE_FILE_PATTERNS =
+
+# The PREDEFINED tag can be used to specify one or more macro names that
+# are defined before the preprocessor is started (similar to the -D option of
+# gcc). The argument of the tag is a list of macros of the form: name
+# or name=definition (no spaces). If the definition and the = are
+# omitted =1 is assumed. To prevent a macro definition from being
+# undefined via #undef or recursively expanded use the := operator
+# instead of the = operator.
+
+PREDEFINED =
+
+# If the MACRO_EXPANSION and EXPAND_ONLY_PREDEF tags are set to YES then
+# this tag can be used to specify a list of macro names that should be expanded.
+# The macro definition that is found in the sources will be used.
+# Use the PREDEFINED tag if you want to use a different macro definition.
+
+EXPAND_AS_DEFINED =
+
+# If the SKIP_FUNCTION_MACROS tag is set to YES (the default) then
+# doxygen's preprocessor will remove all function-like macros that are alone
+# on a line, have an all uppercase name, and do not end with a semicolon. Such
+# function macros are typically used for boiler-plate code, and will confuse
+# the parser if not removed.
+
+SKIP_FUNCTION_MACROS = YES
+
+#---------------------------------------------------------------------------
+# Configuration::additions related to external references
+#---------------------------------------------------------------------------
+
+# The TAGFILES option can be used to specify one or more tagfiles.
+# Optionally an initial location of the external documentation
+# can be added for each tagfile. The format of a tag file without
+# this location is as follows:
+# TAGFILES = file1 file2 ...
+# Adding location for the tag files is done as follows:
+# TAGFILES = file1=loc1 "file2 = loc2" ...
+# where "loc1" and "loc2" can be relative or absolute paths or
+# URLs. If a location is present for each tag, the installdox tool
+# does not have to be run to correct the links.
+# Note that each tag file must have a unique name
+# (where the name does NOT include the path)
+# If a tag file is not located in the directory in which doxygen
+# is run, you must also specify the path to the tagfile here.
+
+TAGFILES =
+
+# When a file name is specified after GENERATE_TAGFILE, doxygen will create
+# a tag file that is based on the input files it reads.
+
+GENERATE_TAGFILE =
+
+# If the ALLEXTERNALS tag is set to YES all external classes will be listed
+# in the class index. If set to NO only the inherited external classes
+# will be listed.
+
+ALLEXTERNALS = NO
+
+# If the EXTERNAL_GROUPS tag is set to YES all external groups will be listed
+# in the modules index. If set to NO, only the current project's groups will
+# be listed.
+
+EXTERNAL_GROUPS = YES
+
+# The PERL_PATH should be the absolute path and name of the perl script
+# interpreter (i.e. the result of `which perl').
+
+PERL_PATH = /usr/bin/perl
+
+#---------------------------------------------------------------------------
+# Configuration options related to the dot tool
+#---------------------------------------------------------------------------
+
+# If the CLASS_DIAGRAMS tag is set to YES (the default) Doxygen will
+# generate a inheritance diagram (in HTML, RTF and LaTeX) for classes with base
+# or super classes. Setting the tag to NO turns the diagrams off. Note that
+# this option is superseded by the HAVE_DOT option below. This is only a
+# fallback. It is recommended to install and use dot, since it yields more
+# powerful graphs.
+
+CLASS_DIAGRAMS = YES
+
+# You can define message sequence charts within doxygen comments using the \msc
+# command. Doxygen will then run the mscgen tool (see
+# http://www.mcternan.me.uk/mscgen/) to produce the chart and insert it in the
+# documentation. The MSCGEN_PATH tag allows you to specify the directory where
+# the mscgen tool resides. If left empty the tool is assumed to be found in the
+# default search path.
+
+MSCGEN_PATH =
+
+# If set to YES, the inheritance and collaboration graphs will hide
+# inheritance and usage relations if the target is undocumented
+# or is not a class.
+
+HIDE_UNDOC_RELATIONS = YES
+
+# If you set the HAVE_DOT tag to YES then doxygen will assume the dot tool is
+# available from the path. This tool is part of Graphviz, a graph visualization
+# toolkit from AT&T and Lucent Bell Labs. The other options in this section
+# have no effect if this option is set to NO (the default)
+
+HAVE_DOT = NO
+
+# By default doxygen will write a font called FreeSans.ttf to the output
+# directory and reference it in all dot files that doxygen generates. This
+# font does not include all possible unicode characters however, so when you need
+# these (or just want a differently looking font) you can specify the font name
+# using DOT_FONTNAME. You need need to make sure dot is able to find the font,
+# which can be done by putting it in a standard location or by setting the
+# DOTFONTPATH environment variable or by setting DOT_FONTPATH to the directory
+# containing the font.
+
+DOT_FONTNAME = FreeSans
+
+# By default doxygen will tell dot to use the output directory to look for the
+# FreeSans.ttf font (which doxygen will put there itself). If you specify a
+# different font using DOT_FONTNAME you can set the path where dot
+# can find it using this tag.
+
+DOT_FONTPATH =
+
+# If the CLASS_GRAPH and HAVE_DOT tags are set to YES then doxygen
+# will generate a graph for each documented class showing the direct and
+# indirect inheritance relations. Setting this tag to YES will force the
+# the CLASS_DIAGRAMS tag to NO.
+
+CLASS_GRAPH = YES
+
+# If the COLLABORATION_GRAPH and HAVE_DOT tags are set to YES then doxygen
+# will generate a graph for each documented class showing the direct and
+# indirect implementation dependencies (inheritance, containment, and
+# class references variables) of the class with other documented classes.
+
+COLLABORATION_GRAPH = YES
+
+# If the GROUP_GRAPHS and HAVE_DOT tags are set to YES then doxygen
+# will generate a graph for groups, showing the direct groups dependencies
+
+GROUP_GRAPHS = YES
+
+# If the UML_LOOK tag is set to YES doxygen will generate inheritance and
+# collaboration diagrams in a style similar to the OMG's Unified Modeling
+# Language.
+
+UML_LOOK = NO
+
+# If set to YES, the inheritance and collaboration graphs will show the
+# relations between templates and their instances.
+
+TEMPLATE_RELATIONS = NO
+
+# If the ENABLE_PREPROCESSING, SEARCH_INCLUDES, INCLUDE_GRAPH, and HAVE_DOT
+# tags are set to YES then doxygen will generate a graph for each documented
+# file showing the direct and indirect include dependencies of the file with
+# other documented files.
+
+INCLUDE_GRAPH = YES
+
+# If the ENABLE_PREPROCESSING, SEARCH_INCLUDES, INCLUDED_BY_GRAPH, and
+# HAVE_DOT tags are set to YES then doxygen will generate a graph for each
+# documented header file showing the documented files that directly or
+# indirectly include this file.
+
+INCLUDED_BY_GRAPH = YES
+
+# If the CALL_GRAPH and HAVE_DOT options are set to YES then
+# doxygen will generate a call dependency graph for every global function
+# or class method. Note that enabling this option will significantly increase
+# the time of a run. So in most cases it will be better to enable call graphs
+# for selected functions only using the \callgraph command.
+
+CALL_GRAPH = NO
+
+# If the CALLER_GRAPH and HAVE_DOT tags are set to YES then
+# doxygen will generate a caller dependency graph for every global function
+# or class method. Note that enabling this option will significantly increase
+# the time of a run. So in most cases it will be better to enable caller
+# graphs for selected functions only using the \callergraph command.
+
+CALLER_GRAPH = NO
+
+# If the GRAPHICAL_HIERARCHY and HAVE_DOT tags are set to YES then doxygen
+# will graphical hierarchy of all classes instead of a textual one.
+
+GRAPHICAL_HIERARCHY = YES
+
+# If the DIRECTORY_GRAPH, SHOW_DIRECTORIES and HAVE_DOT tags are set to YES
+# then doxygen will show the dependencies a directory has on other directories
+# in a graphical way. The dependency relations are determined by the #include
+# relations between the files in the directories.
+
+DIRECTORY_GRAPH = YES
+
+# The DOT_IMAGE_FORMAT tag can be used to set the image format of the images
+# generated by dot. Possible values are png, jpg, or gif
+# If left blank png will be used.
+
+DOT_IMAGE_FORMAT = png
+
+# The tag DOT_PATH can be used to specify the path where the dot tool can be
+# found. If left blank, it is assumed the dot tool can be found in the path.
+
+DOT_PATH =
+
+# The DOTFILE_DIRS tag can be used to specify one or more directories that
+# contain dot files that are included in the documentation (see the
+# \dotfile command).
+
+DOTFILE_DIRS =
+
+# The DOT_GRAPH_MAX_NODES tag can be used to set the maximum number of
+# nodes that will be shown in the graph. If the number of nodes in a graph
+# becomes larger than this value, doxygen will truncate the graph, which is
+# visualized by representing a node as a red box. Note that doxygen if the
+# number of direct children of the root node in a graph is already larger than
+# DOT_GRAPH_MAX_NODES then the graph will not be shown at all. Also note
+# that the size of a graph can be further restricted by MAX_DOT_GRAPH_DEPTH.
+
+DOT_GRAPH_MAX_NODES = 50
+
+# The MAX_DOT_GRAPH_DEPTH tag can be used to set the maximum depth of the
+# graphs generated by dot. A depth value of 3 means that only nodes reachable
+# from the root by following a path via at most 3 edges will be shown. Nodes
+# that lay further from the root node will be omitted. Note that setting this
+# option to 1 or 2 may greatly reduce the computation time needed for large
+# code bases. Also note that the size of a graph can be further restricted by
+# DOT_GRAPH_MAX_NODES. Using a depth of 0 means no depth restriction.
+
+MAX_DOT_GRAPH_DEPTH = 0
+
+# Set the DOT_TRANSPARENT tag to YES to generate images with a transparent
+# background. This is enabled by default, which results in a transparent
+# background. Warning: Depending on the platform used, enabling this option
+# may lead to badly anti-aliased labels on the edges of a graph (i.e. they
+# become hard to read).
+
+DOT_TRANSPARENT = YES
+
+# Set the DOT_MULTI_TARGETS tag to YES allow dot to generate multiple output
+# files in one run (i.e. multiple -o and -T options on the command line). This
+# makes dot run faster, but since only newer versions of dot (>1.8.10)
+# support this, this feature is disabled by default.
+
+DOT_MULTI_TARGETS = NO
+
+# If the GENERATE_LEGEND tag is set to YES (the default) Doxygen will
+# generate a legend page explaining the meaning of the various boxes and
+# arrows in the dot generated graphs.
+
+GENERATE_LEGEND = YES
+
+# If the DOT_CLEANUP tag is set to YES (the default) Doxygen will
+# remove the intermediate dot files that are used to generate
+# the various graphs.
+
+DOT_CLEANUP = YES
+
+#---------------------------------------------------------------------------
+# Configuration::additions related to the search engine
+#---------------------------------------------------------------------------
+
+# The SEARCHENGINE tag specifies whether or not a search engine should be
+# used. If set to NO the values of all tags below this one will be ignored.
+
+SEARCHENGINE = NO
diff -r 9992112ff05a -r 7596917b8776 elastic3d.f90
--- a/elastic3d.f90 Mon Apr 11 08:59:18 2011 -0700
+++ b/elastic3d.f90 Tue Apr 26 15:00:49 2011 -0700
@@ -19,6 +19,7 @@
MODULE elastic3d
+ USE types
USE fourier
IMPLICIT NONE
@@ -30,53 +31,6 @@ MODULE elastic3d
REAL*8, PRIVATE, PARAMETER :: pid2 = 1.57079632679489655799898173427209258079_8
REAL*8, PRIVATE, PARAMETER :: DEG2RAD = 0.01745329251994329547437168059786927_8
- TYPE SOURCE_STRUCT
- SEQUENCE
- REAL*8 :: slip,x,y,z,width,length,strike,dip,rake
- END TYPE SOURCE_STRUCT
-
- TYPE PLANE_STRUCT
- SEQUENCE
- REAL*8 :: x,y,z,width,length,strike,dip,rake
- END TYPE PLANE_STRUCT
-
- TYPE LAYER_STRUCT
- SEQUENCE
- REAL*8 :: z,gammadot0,stressexponent,cohesion,friction
- END TYPE LAYER_STRUCT
-
- TYPE WEAK_STRUCT
- SEQUENCE
- REAL*8 :: dgammadot0,x,y,z,width,length,thickness,strike,dip
- END TYPE WEAK_STRUCT
-
- TYPE VECTOR_STRUCT
- SEQUENCE
- REAL*8 :: v1,v2,v3
- END TYPE VECTOR_STRUCT
-
- TYPE TENSOR
- SEQUENCE
- REAL*4 :: s11,s12,s13,s22,s23,s33
- END TYPE TENSOR
-
- TYPE TENSOR_LAYER_STRUCT
- SEQUENCE
- REAL*4 :: z,dum
- TYPE(TENSOR) :: t
- END TYPE TENSOR_LAYER_STRUCT
-
- TYPE SLIPPATCH_STRUCT
- SEQUENCE
- REAL*8 :: x1,x2,x3,lx,lz,slip,ss,ds
- END TYPE SLIPPATCH_STRUCT
-
- TYPE EVENT_STRUC
- REAL*8 :: time
- INTEGER*4 :: i,ns,nt,nm,nl
- TYPE(SOURCE_STRUCT), DIMENSION(:), ALLOCATABLE :: s,sc,ts,tsc,m,mc,l,lc
- END TYPE EVENT_STRUC
-
INTERFACE OPERATOR (.times.)
MODULE PROCEDURE tensorscalarprod
END INTERFACE
diff -r 9992112ff05a -r 7596917b8776 examples/run1.sh
--- a/examples/run1.sh Mon Apr 11 08:59:18 2011 -0700
+++ b/examples/run1.sh Tue Apr 26 15:00:49 2011 -0700
@@ -20,21 +20,23 @@
# type map.sh for a description of command-line options.
-time ../relax <<EOF
+time ../relax <<EOF | tee output1/in.param
# grid size (sx1,sx2,sx3)
256 256 256
# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
0.05 0.05 0.05 0.2 2
# origin position & rotation
0 0 0
+# geographic origin, UTM zone and length unit
+120 22 51 1e3
# observation depth (displacement and stress) (stress in only exported in GRD)
0 0.5
# output directory
./output1
# elastic parameters and gamma = (1-nu) rho g / mu = 8.33e-7 /m = 8.33e-4 /km
1 1 8.33e-4
-# integration time (t1)
-20 -1
+# integration time (t1), step and scaling
+20 -1 1
# number of observation planes
0
# number of observation points
@@ -65,6 +67,8 @@ 0
0
# number of dilatation sources
0
+# number of surface traction
+0
# time of second event
1
# number of shear dislocations
@@ -75,4 +79,6 @@ 0
0
# number of dilatation sources
0
+# number of surface traction
+0
EOF
diff -r 9992112ff05a -r 7596917b8776 makefile_imkl
--- a/makefile_imkl Mon Apr 11 08:59:18 2011 -0700
+++ b/makefile_imkl Tue Apr 26 15:00:49 2011 -0700
@@ -22,7 +22,7 @@ F77FLAGS=-zero
F77FLAGS=-zero
CFLAGS=-I/sw/include
-OBJ = mkl_dfti.o fourier.o green.o elastic3d.o friction3d.o viscoelastic3d.o writegrd4.2.o proj.o export.o getdata.o getopt.o relax.o
+OBJ = types.o mkl_dfti.o fourier.o green.o elastic3d.o friction3d.o viscoelastic3d.o writegrd4.2.o proj.o export.o getdata.o getopt.o input.o relax.o
%.o : %.c
$(CC) $(CFLAGS) -c $^
diff -r 9992112ff05a -r 7596917b8776 relax.f90
--- a/relax.f90 Mon Apr 11 08:59:18 2011 -0700
+++ b/relax.f90 Tue Apr 26 15:00:49 2011 -0700
@@ -158,6 +158,8 @@ PROGRAM relax
! and allow scaling of computed time steps.
!-----------------------------------------------------------------------
+ USE types
+ USE input
USE green
USE elastic3d
USE viscoelastic3d
@@ -168,51 +170,27 @@ PROGRAM relax
IMPLICIT NONE
- REAL*8, PARAMETER :: DEG2RAD = 0.01745329251994329547437168059786927_8
INTEGER, PARAMETER :: ITERATION_MAX = 9900
REAL*8, PARAMETER :: STEP_MAX = 1e7
- INTEGER :: i,k,sx1,sx2,sx3,e,ne,nv,np,nop,npl,nps,oi,nfc, &
- unit,iostatus,iargc,npts,skip=0,mech(3),nlwz,nnlwz
+ INTEGER :: i,k,e,oi,iostatus,mech(3)
#ifdef FFTW3_THREADS
INTEGER :: iret
!$ INTEGER :: omp_get_max_threads
#endif
- REAL*8 :: beta,lambda,mu,gam,x0,y0,interval, &
- minlength,minwidth,rot,maxwell(3),nyquist
-#ifdef PROJ
- REAL*8 :: lon0,lat0,umult
- INTEGER :: zone
-#endif
- CHARACTER(80) :: wdir,reporttimefilename,reportfilename, &
- inputfile,logfilename,inputfilename
+ REAL*8 :: maxwell(3)
+ TYPE(SIMULATION_STRUC) :: in
#ifdef VTK
- INTEGER :: j
- CHARACTER(80) :: rffilename,vcfilename,cgfilename
+ CHARACTER(80) :: vcfilename
CHARACTER(3) :: digit
#endif
- REAL*8 :: dx1,dx2,dx3,oz,ozs,t,Dt,tm,odt,tscale
- ! coseismic events
- TYPE(EVENT_STRUC), DIMENSION(:), ALLOCATABLE :: events
- TYPE(EVENT_STRUC) :: inter
-
- ! input dislocation (shear and tensile cracks)
- TYPE(PLANE_STRUCT), DIMENSION(:), ALLOCATABLE :: n, op
- TYPE(LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: linearlayer,nonlinearlayer
- TYPE(LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: faultcreeplayer
- TYPE(LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: linearstruc,nonlinearstruc
- TYPE(LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: faultcreepstruc
- TYPE(TENSOR_LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: stresslayer,stressstruc
- TYPE(WEAK_STRUCT), DIMENSION(:), ALLOCATABLE :: linearweakzone,linearweakzonec, &
- nonlinearweakzone,nonlinearweakzonec
+ REAL*8 :: t,Dt,tm
! arrays
REAL*4, DIMENSION(:,:,:), ALLOCATABLE :: v1,v2,v3,u1,u2,u3,gamma
REAL*4, DIMENSION(:,:), ALLOCATABLE :: t1,t2,t3
REAL*4, DIMENSION(:,:,:), ALLOCATABLE :: inter1,inter2,inter3
TYPE(TENSOR), DIMENSION(:,:,:), ALLOCATABLE :: tau,sig,moment
- TYPE(VECTOR_STRUCT), DIMENSION(:), ALLOCATABLE :: opts
- CHARACTER(LEN=4), DIMENSION(:), ALLOCATABLE :: ptsname
#ifdef FFTW3_THREADS
CALL sfftw_init_threads(iret)
@@ -223,74 +201,63 @@ PROGRAM relax
#endif
#endif
- ! read standard input or filename given in argument
- IF (0 .EQ. iargc()) THEN
- ! standard input
- unit=5
- ELSE
- ! open input file
- CALL getarg(1,inputfile)
+ ! read input parameters
+ CALL init(in)
- OPEN (UNIT=15,FILE=inputfile,IOSTAT=iostatus,FORM="FORMATTED")
- IF (iostatus .GT. 0) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("unable to access input file ",a)') inputfile
- STOP 1
- END IF
- ! input file
- unit=15
+ IF (in%isdryrun) THEN
+ PRINT '("abort calculation")'
+ END IF
+ IF (in%isdryrun .OR. in%ishelp) THEN
+ ! exit program
+ GOTO 100
END IF
- CALL init(UNIT=unit)
-
- ! close input file
- IF (iargc() .GT. 0) CLOSE(15)
-
- ALLOCATE (v1(sx1+2,sx2,sx3),v2(sx1+2,sx2,sx3),v3(sx1+2,sx2,sx3), &
- u1(sx1+2,sx2,sx3/2),u2(sx1+2,sx2,sx3/2),u3(sx1+2,sx2,sx3/2), &
- inter1(sx1+2,sx2,2),inter2(sx1+2,sx2,2),inter3(sx1+2,sx2,2), &
- tau(sx1,sx2,sx3/2),gamma(sx1+2,sx2,sx3/2), &
- t1(sx1+2,sx2),t2(sx1+2,sx2),t3(sx1+2,sx2), &
+ ! allocate memory
+ ALLOCATE (v1(in%sx1+2,in%sx2,in%sx3),v2(in%sx1+2,in%sx2,in%sx3),v3(in%sx1+2,in%sx2,in%sx3), &
+ u1(in%sx1+2,in%sx2,in%sx3/2),u2(in%sx1+2,in%sx2,in%sx3/2),u3(in%sx1+2,in%sx2,in%sx3/2), &
+ inter1(in%sx1+2,in%sx2,2),inter2(in%sx1+2,in%sx2,2),inter3(in%sx1+2,in%sx2,2), &
+ tau(in%sx1,in%sx2,in%sx3/2),gamma(in%sx1+2,in%sx2,in%sx3/2), &
+ t1(in%sx1+2,in%sx2),t2(in%sx1+2,in%sx2),t3(in%sx1+2,in%sx2), &
STAT=iostatus)
IF (iostatus>0) STOP "could not allocate memory"
v1=0;v2=0;v3=0;u1=0;u2=0;u3=0;gamma=0;t1=0;t2=0;t3=0
- CALL tensorfieldadd(tau,tau,sx1,sx2,sx3/2,c1=0._4,c2=0._4)
+ CALL tensorfieldadd(tau,tau,in%sx1,in%sx2,in%sx3/2,c1=0._4,c2=0._4)
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! - construct pre-stress structure
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- IF (ALLOCATED(stresslayer)) THEN
- CALL tensorstructure(stressstruc,stresslayer,dx3)
- DEALLOCATE(stresslayer)
+ IF (ALLOCATED(in%stresslayer)) THEN
+ CALL tensorstructure(in%stressstruc,in%stresslayer,in%dx3)
+ DEALLOCATE(in%stresslayer)
- DO k=1,sx3/2
- tau(:,:,k)=(-1._4) .times. stressstruc(k)%t
+ DO k=1,in%sx3/2
+ tau(:,:,k)=(-1._4) .times. in%stressstruc(k)%t
END DO
- DEALLOCATE(stressstruc)
+ DEALLOCATE(in%stressstruc)
END IF
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! - construct linear viscoelastic structure
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- IF (ALLOCATED(linearlayer)) THEN
- CALL viscoelasticstructure(linearstruc,linearlayer,dx3)
- DEALLOCATE(linearlayer)
+ IF (ALLOCATED(in%linearlayer)) THEN
+ CALL viscoelasticstructure(in%linearstruc,in%linearlayer,in%dx3)
+ DEALLOCATE(in%linearlayer)
END IF
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! - construct nonlinear viscoelastic structure
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- IF (ALLOCATED(nonlinearlayer)) THEN
- CALL viscoelasticstructure(nonlinearstruc,nonlinearlayer,dx3)
- DEALLOCATE(nonlinearlayer)
+ IF (ALLOCATED(in%nonlinearlayer)) THEN
+ CALL viscoelasticstructure(in%nonlinearstruc,in%nonlinearlayer,in%dx3)
+ DEALLOCATE(in%nonlinearlayer)
END IF
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! - construct nonlinear fault creep structure (rate-strenghtening)
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- IF (ALLOCATED(faultcreeplayer)) THEN
- CALL viscoelasticstructure(faultcreepstruc,faultcreeplayer,dx3)
- DEALLOCATE(faultcreeplayer)
+ IF (ALLOCATED(in%faultcreeplayer)) THEN
+ CALL viscoelasticstructure(in%faultcreepstruc,in%faultcreeplayer,in%dx3)
+ DEALLOCATE(in%faultcreeplayer)
END IF
! first event
@@ -299,94 +266,114 @@ PROGRAM relax
oi=1;
! sources
- CALL dislocations(events(e),lambda,mu,beta,sx1,sx2,sx3, &
- dx1,dx2,dx3,v1,v2,v3,t1,t2,t3,tau)
- CALL traction(mu,events(e),sx1,sx2,dx1,dx2,t3)
+ CALL dislocations(in%events(e),in%lambda,in%mu,in%beta,in%sx1,in%sx2,in%sx3, &
+ in%dx1,in%dx2,in%dx3,v1,v2,v3,t1,t2,t3,tau)
+ CALL traction(in%mu,in%events(e),in%sx1,in%sx2,in%dx1,in%dx2,t3)
PRINT '("coseismic event ",I3.3)', e
PRINT 0990
! export the amplitude of eigenstrain
- CALL exporteigenstrain(gamma,nop,op,x0,y0,dx1,dx2,dx3,sx1,sx2,sx3/2,wdir,0)
+ CALL exporteigenstrain(gamma,in%nop,in%op,in%x0,in%y0, &
+ in%dx1,in%dx2,in%dx3,in%sx1,in%sx2,in%sx3/2,in%wdir,0)
! export equivalent body forces
- IF (isoutput(skip,t,i,odt,oi,events(e)%time)) THEN
+ IF (isoutput(in%skip,t,i,in%odt,oi,in%events(e)%time)) THEN
#ifdef GRD_EQBF
- CALL exportgrd(v1,v2,v3,sx1,sx2,sx3/2,dx1,dx2,dx3,0.7_8,x0,y0,wdir,0,convention=3)
+ IF (in%isoutputgrd) THEN
+ CALL exportgrd(v1,v2,v3,in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,0.7_8,in%x0,in%y0,in%wdir,0,convention=3)
+ END IF
#endif
END IF
! test the presence of dislocations for coseismic calculation
- IF ((events(e)%nt .NE. 0) .OR. &
- (events(e)%ns .NE. 0) .OR. &
- (events(e)%nm .NE. 0)) THEN
+ IF ((in%events(e)%nt .NE. 0) .OR. &
+ (in%events(e)%ns .NE. 0) .OR. &
+ (in%events(e)%nm .NE. 0)) THEN
! apply the 3d elastic transfer function
- CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3,dx1,dx2,dx3,lambda,mu,gam)
+ CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3, &
+ in%dx1,in%dx2,in%dx3,in%lambda,in%mu,in%gam)
END IF
! transfer solution
- CALL fieldrep(u1,v1,sx1+2,sx2,sx3/2)
- CALL fieldrep(u2,v2,sx1+2,sx2,sx3/2)
- CALL fieldrep(u3,v3,sx1+2,sx2,sx3/2)
+ CALL fieldrep(u1,v1,in%sx1+2,in%sx2,in%sx3/2)
+ CALL fieldrep(u2,v2,in%sx1+2,in%sx2,in%sx3/2)
+ CALL fieldrep(u3,v3,in%sx1+2,in%sx2,in%sx3/2)
! export
#ifdef TXT
- CALL exporttxt(u1,u2,u3,sx1,sx2,sx3/2,oz,dx3,0,0._8,wdir,reportfilename)
+ IF (in%isoutputtxt) THEN
+ CALL exporttxt(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%oz,in%dx3,0,0._8,in%wdir,in%reportfilename)
+ END IF
#endif
#ifdef XYZ
- CALL exportxyz(u1,u2,u3,sx1,sx2,sx3/2,oz,dx1,dx2,dx3,0,wdir)
+ IF (in%isoutputxyz) THEN
+ CALL exportxyz(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%oz,in%dx1,in%dx2,in%dx3,0,in%wdir)
+ END IF
#endif
#ifdef GRD
- CALL exportgrd(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3,oz,x0,y0,wdir,0)
- CALL exportgrd(inter1,inter2,inter3,sx1,sx2,sx3/2, &
- dx1,dx2,dx3,0._8,x0,y0,wdir,0,convention=2)
+ IF (in%isoutputgrd) THEN
+ CALL exportgrd(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%oz,in%x0,in%y0,in%wdir,0)
+ CALL exportgrd(inter1,inter2,inter3,in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,0._8,in%x0,in%y0,in%wdir,0,convention=2)
+ END IF
#endif
#ifdef PROJ
- CALL exportproj(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3,oz, &
- x0,y0,lon0,lat0,zone,umult,wdir,0)
+ IF (in%isoutputproj) THEN
+ CALL exportproj(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%oz, &
+ in%x0,in%y0,in%lon0,in%lat0,in%zone,in%umult,in%wdir,0)
+ END IF
#endif
#ifdef VTK
- j=INDEX(wdir," ")
- vcfilename=wdir(1:j-1)//"/disp-000.vtr"
- CALL exportvtk_vectors(u1,u2,u3,sx1,sx2,sx3/4,dx1,dx2,dx3,8,8,8,vcfilename)
- !CALL exportvtk_vectors_slice(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3,oz,8,8,vcfilename)
+ IF (in%isoutputvtk) THEN
+ vcfilename=trim(in%wdir)//"/disp-000.vtr"
+ CALL exportvtk_vectors(u1,u2,u3,in%sx1,in%sx2,in%sx3/4,in%dx1,in%dx2,in%dx3,8,8,8,vcfilename)
+ !CALL exportvtk_vectors_slice(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%oz,8,8,vcfilename)
+ END IF
#endif
- IF (ALLOCATED(ptsname)) THEN
- CALL exportpoints(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3, &
- opts,ptsname,0._8,wdir,.true.,x0,y0,rot)
+ IF (ALLOCATED(in%ptsname)) THEN
+ CALL exportpoints(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3, &
+ in%opts,in%ptsname,0._8,in%wdir,.true.,in%x0,in%y0,in%rot)
END IF
- CALL reporttime(0,0._8,reporttimefilename)
+ CALL reporttime(0,0._8,in%reporttimefilename)
- PRINT 1101,0,0._8,0._8,0._8,0._8,0._8,interval,0._8,tensoramplitude(tau,dx1,dx2,dx3)
- IF (interval .LE. 0) THEN
+ PRINT 1101,0,0._8,0._8,0._8,0._8,0._8,in%interval,0._8,tensoramplitude(tau,in%dx1,in%dx2,in%dx3)
+ IF (in%interval .LE. 0) THEN
GOTO 100 ! no time integration
END IF
- ALLOCATE(moment(sx1,sx2,sx3/2),sig(sx1,sx2,sx3/2),STAT=iostatus)
+ ALLOCATE(moment(in%sx1,in%sx2,in%sx3/2),sig(in%sx1,in%sx2,in%sx3/2),STAT=iostatus)
IF (iostatus>0) STOP "could not allocate the mechanical structure"
- CALL tensorfieldadd(sig,sig,sx1,sx2,sx3/2,c1=0._4,c2=0._4)
- CALL tensorfieldadd(moment,moment,sx1,sx2,sx3/2,c1=0._4,c2=0._4)
+ CALL tensorfieldadd(sig,sig,in%sx1,in%sx2,in%sx3/2,c1=0._4,c2=0._4)
+ CALL tensorfieldadd(moment,moment,in%sx1,in%sx2,in%sx3/2,c1=0._4,c2=0._4)
t=0
DO i=1,ITERATION_MAX
- IF (t > (interval+1e-6)) GOTO 100 ! proper exit
+ IF (t > (in%interval+1e-6)) GOTO 100 ! proper exit
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! predictor
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- CALL tensorfieldadd(sig,tau,sx1,sx2,sx3/2,c1=0._4,c2=-1._4)
- CALL stressupdate(u1,u2,u3,lambda,mu,dx1,dx2,dx3,sx1,sx2,sx3/2,sig)
+ CALL tensorfieldadd(sig,tau,in%sx1,in%sx2,in%sx3/2,c1=0._4,c2=-1._4)
+ CALL stressupdate(u1,u2,u3,in%lambda,in%mu, &
+ in%dx1,in%dx2,in%dx3,in%sx1,in%sx2,in%sx3/2,sig)
- IF (isoutput(skip,t,i,odt,oi,events(e)%time)) THEN
+ IF (isoutput(in%skip,t,i,in%odt,oi,in%events(e)%time)) THEN
! export stress
#ifdef GRD
- CALL exportstressgrd(sig,sx1,sx2,sx3/2,dx1,dx2,dx3,ozs,x0,y0,wdir,i-1)
+ IF (in%isoutputgrd) THEN
+ CALL exportstressgrd(sig,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3, &
+ in%ozs,in%x0,in%y0,in%wdir,i-1)
+ END IF
#endif
#ifdef PROJ
- CALL exportstressproj(sig,sx1,sx2,sx3/2,dx1,dx2,dx3,ozs, &
- x0,y0,lon0,lat0,zone,umult,wdir,i-1)
+ IF (in%isoutputproj) THEN
+ CALL exportstressproj(sig,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%ozs, &
+ in%x0,in%y0,in%lon0,in%lat0,in%zone,in%umult,in%wdir,i-1)
+ END IF
#endif
END IF
@@ -398,31 +385,35 @@ PROGRAM relax
mech(:)=0
! initialize no forcing term in tensor space
- CALL tensorfieldadd(moment,moment,sx1,sx2,sx3/2,0._4,0._4)
+ CALL tensorfieldadd(moment,moment,in%sx1,in%sx2,in%sx3/2,0._4,0._4)
! power density from three mechanisms (linear and power-law viscosity
! and fault creep)
! 1- linear viscosity
- IF (ALLOCATED(linearstruc)) THEN
- CALL viscouseigenstress(mu,linearstruc,linearweakzone,sig,sx1,sx2,sx3/2, &
- dx1,dx2,dx3,moment,0.01_8,MAXWELLTIME=maxwell(1))
+ IF (ALLOCATED(in%linearstruc)) THEN
+ CALL viscouseigenstress(in%mu,in%linearstruc,in%linearweakzone, &
+ sig,in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,moment,0.01_8,MAXWELLTIME=maxwell(1))
mech(1)=1
END IF
! 2- powerlaw viscosity
- IF (ALLOCATED(nonlinearstruc)) THEN
- CALL viscouseigenstress(mu,nonlinearstruc,nonlinearweakzone,sig,sx1,sx2,sx3/2, &
- dx1,dx2,dx3,moment,0.01_8,MAXWELLTIME=maxwell(2))
+ IF (ALLOCATED(in%nonlinearstruc)) THEN
+ CALL viscouseigenstress(in%mu,in%nonlinearstruc,in%nonlinearweakzone, &
+ sig,in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,moment,0.01_8,MAXWELLTIME=maxwell(2))
mech(2)=1
END IF
! 3- nonlinear fault creep with rate-strengthening friction
- IF (ALLOCATED(faultcreepstruc)) THEN
- DO k=1,np
- CALL frictioneigenstress(n(k)%x,n(k)%y,n(k)%z, &
- n(k)%width,n(k)%length,n(k)%strike,n(k)%dip,n(k)%rake,beta, &
- sig,mu,faultcreepstruc,sx1,sx2,sx3/2,dx1,dx2,dx3,moment, &
- maxwelltime=maxwell(3))
+ IF (ALLOCATED(in%faultcreepstruc)) THEN
+ DO k=1,in%np
+ CALL frictioneigenstress(in%n(k)%x,in%n(k)%y,in%n(k)%z, &
+ in%n(k)%width,in%n(k)%length, &
+ in%n(k)%strike,in%n(k)%dip,in%n(k)%rake,in%beta, &
+ sig,in%mu,in%faultcreepstruc, &
+ in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3, &
+ moment,maxwelltime=maxwell(3))
END DO
mech(3)=1
END IF
@@ -435,230 +426,250 @@ PROGRAM relax
tm=MIN(tm,STEP_MAX)
! modify
- IF ((inter%ns .GT. 0) .OR. (inter%nt .GT. 0)) THEN
+ IF ((in%inter%ns .GT. 0) .OR. (in%inter%nt .GT. 0)) THEN
IF (tm .EQ. STEP_MAX) THEN
! no relaxation occurs, pick a small integration time
- tm=interval/20._8
+ tm=in%interval/20._8
END IF
END IF
! choose an integration time step
- CALL integrationstep(tm,Dt,t,oi,odt,skip,tscale,events,e,ne)
+ CALL integrationstep(tm,Dt,t,oi,in%odt,in%skip,in%tscale,in%events,e,in%ne)
- CALL tensorfieldadd(sig,moment,sx1,sx2,sx3/2,c1=0.0_4,c2=1._4)
+ CALL tensorfieldadd(sig,moment,in%sx1,in%sx2,in%sx3/2,c1=0.0_4,c2=1._4)
v1=0;v2=0;v3=0;t1=0;t2=0;t3=0;
- CALL equivalentbodyforce(sig,dx1,dx2,dx3,sx1,sx2,sx3/2,v1,v2,v3,t1,t2,t3)
- CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3,dx1,dx2,dx3,lambda,mu,gam)
+ CALL equivalentbodyforce(sig,in%dx1,in%dx2,in%dx3,in%sx1,in%sx2,in%sx3/2,v1,v2,v3,t1,t2,t3)
+ CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3,in%dx1,in%dx2,in%dx3,in%lambda,in%mu,in%gam)
! v1,v2,v3 contain the predictor displacement
- CALL fieldadd(v1,u1,sx1+2,sx2,sx3/2,c1=REAL(Dt/2))
- CALL fieldadd(v2,u2,sx1+2,sx2,sx3/2,c1=REAL(Dt/2))
- CALL fieldadd(v3,u3,sx1+2,sx2,sx3/2,c1=REAL(Dt/2))
- CALL tensorfieldadd(sig,tau,sx1,sx2,sx3/2,c1=-REAL(Dt/2),c2=-1._4)
+ CALL fieldadd(v1,u1,in%sx1+2,in%sx2,in%sx3/2,c1=REAL(Dt/2))
+ CALL fieldadd(v2,u2,in%sx1+2,in%sx2,in%sx3/2,c1=REAL(Dt/2))
+ CALL fieldadd(v3,u3,in%sx1+2,in%sx2,in%sx3/2,c1=REAL(Dt/2))
+ CALL tensorfieldadd(sig,tau,in%sx1,in%sx2,in%sx3/2,c1=-REAL(Dt/2),c2=-1._4)
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! corrector
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- CALL stressupdate(v1,v2,v3,lambda,mu,dx1,dx2,dx3,sx1,sx2,sx3/2,sig)
+ CALL stressupdate(v1,v2,v3,in%lambda,in%mu, &
+ in%dx1,in%dx2,in%dx3,in%sx1,in%sx2,in%sx3/2,sig)
! reinitialize moment density tensor
- CALL tensorfieldadd(moment,moment,sx1,sx2,sx3/2,0._4,0._4)
+ CALL tensorfieldadd(moment,moment,in%sx1,in%sx2,in%sx3/2,0._4,0._4)
- IF (ALLOCATED(linearstruc)) THEN
+ IF (ALLOCATED(in%linearstruc)) THEN
! linear viscosity
v1=0
- CALL viscouseigenstress(mu,linearstruc,linearweakzone,sig,sx1,sx2,sx3/2, &
- dx1,dx2,dx3,moment,0.01_8,GAMMA=v1)
+ CALL viscouseigenstress(in%mu,in%linearstruc,in%linearweakzone,sig, &
+ in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,moment,0.01_8,GAMMA=v1)
! update slip history
- CALL fieldadd(gamma,v1,sx1+2,sx2,sx3/2,c2=REAL(Dt))
+ CALL fieldadd(gamma,v1,in%sx1+2,in%sx2,in%sx3/2,c2=REAL(Dt))
END IF
- IF (ALLOCATED(nonlinearstruc)) THEN
+ IF (ALLOCATED(in%nonlinearstruc)) THEN
! powerlaw viscosity
v1=0
- CALL viscouseigenstress(mu,nonlinearstruc,nonlinearweakzone,sig,sx1,sx2,sx3/2, &
- dx1,dx2,dx3,moment,0.01_8,GAMMA=v1)
+ CALL viscouseigenstress(in%mu,in%nonlinearstruc,in%nonlinearweakzone,sig, &
+ in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,moment,0.01_8,GAMMA=v1)
! update slip history
- CALL fieldadd(gamma,v1,sx1+2,sx2,sx3/2,c2=REAL(Dt))
+ CALL fieldadd(gamma,v1,in%sx1+2,in%sx2,in%sx3/2,c2=REAL(Dt))
END IF
! nonlinear fault creep with rate-strengthening friction
- IF (ALLOCATED(faultcreepstruc)) THEN
+ IF (ALLOCATED(in%faultcreepstruc)) THEN
! use v1 as placeholders for the afterslip planes
v1=0
- DO k=1,np
+ DO k=1,in%np
! one may use optional arguments ...,VEL=v1) to convert
! fault slip to eigenstrain (scalar)
- CALL frictioneigenstress(n(k)%x,n(k)%y,n(k)%z, &
- n(k)%width,n(k)%length,n(k)%strike,n(k)%dip,n(k)%rake,beta, &
- sig,mu,faultcreepstruc,sx1,sx2,sx3/2,dx1,dx2,dx3,moment)
+ CALL frictioneigenstress(in%n(k)%x,in%n(k)%y,in%n(k)%z, &
+ in%n(k)%width,in%n(k)%length, &
+ in%n(k)%strike,in%n(k)%dip,in%n(k)%rake,in%beta, &
+ sig,in%mu,in%faultcreepstruc,in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,moment)
END DO
! update slip history
- CALL fieldadd(gamma,v1,sx1+2,sx2,sx3/2,c2=REAL(Dt))
+ CALL fieldadd(gamma,v1,in%sx1+2,in%sx2,in%sx3/2,c2=REAL(Dt))
! export strike and dip creep velocity
- IF (isoutput(skip,t,i,odt,oi,events(e)%time)) THEN
- CALL exportcreep(np,n,beta,sig,faultcreepstruc, &
- sx1,sx2,sx3/2,dx1,dx2,dx3,x0,y0,wdir,oi)
+ IF (isoutput(in%skip,t,i,in%odt,oi,in%events(e)%time)) THEN
+ CALL exportcreep(in%np,in%n,in%beta,sig,in%faultcreepstruc, &
+ in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%x0,in%y0,in%wdir,oi)
END IF
END IF
! interseismic loading
- IF ((inter%ns .GT. 0) .OR. (inter%nt .GT. 0)) THEN
+ IF ((in%inter%ns .GT. 0) .OR. (in%inter%nt .GT. 0)) THEN
! vectors v1,v2,v3 are not affected.
- CALL dislocations(inter,lambda,mu,beta,sx1,sx2,sx3, &
- dx1,dx2,dx3,v1,v2,v3,t1,t2,t3,tau,factor=Dt,eigenstress=moment)
+ CALL dislocations(in%inter,in%lambda,in%mu,in%beta,in%sx1,in%sx2,in%sx3, &
+ in%dx1,in%dx2,in%dx3,v1,v2,v3,t1,t2,t3,tau,factor=Dt,eigenstress=moment)
END IF
v1=0;v2=0;v3=0;t1=0;t2=0;t3=0;
- CALL equivalentbodyforce(moment,dx1,dx2,dx3,sx1,sx2,sx3/2,v1,v2,v3,t1,t2,t3)
+ CALL equivalentbodyforce(moment,in%dx1,in%dx2,in%dx3,in%sx1,in%sx2,in%sx3/2,v1,v2,v3,t1,t2,t3)
! export equivalent body forces
- IF (isoutput(skip,t,i,odt,oi,events(e)%time)) THEN
+ IF (isoutput(in%skip,t,i,in%odt,oi,in%events(e)%time)) THEN
#ifdef VTK_EQBF
- WRITE (digit,'(I3.3)') oi
- j=INDEX(wdir," ")
- vcfilename=wdir(1:j-1)//"/eqbf-"//digit//".vtr"
- CALL exportvtk_vectors(v1,v2,v3,sx1,sx2,sx3/4,dx1,dx2,dx3,8,8,8,vcfilename)
+ IF (in%isoutputvtk) THEN
+ WRITE (digit,'(I3.3)') oi
+ vcfilename=trim(in%wdir)//"/eqbf-"//digit//".vtr"
+ CALL exportvtk_vectors(v1,v2,v3,in%sx1,in%sx2,in%sx3/4,in%dx1,in%dx2,in%dx3,8,8,8,vcfilename)
+ END IF
#endif
#ifdef GRD_EQBF
- CALL exportgrd(v1,v2,v3,sx1,sx2,sx3/2,dx1,dx2,dx3,30.7_8,x0,y0,wdir,oi,convention=3)
+ IF (in%isoutputgrd) THEN
+ CALL exportgrd(v1,v2,v3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3, &
+ in%oz,in%x0,y0,in%wdir,oi,convention=3)
+ END IF
#endif
END IF
- CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3,dx1,dx2,dx3,lambda,mu,gam)
+ CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3,in%dx1,in%dx2,in%dx3,in%lambda,in%mu,in%gam)
! update deformation field
- CALL fieldadd(u1,v1,sx1+2,sx2,sx3/2,c2=REAL(Dt))
- CALL fieldadd(u2,v2,sx1+2,sx2,sx3/2,c2=REAL(Dt))
- CALL fieldadd(u3,v3,sx1+2,sx2,sx3/2,c2=REAL(Dt))
- CALL tensorfieldadd(tau,moment,sx1,sx2,sx3/2,c2=REAL(Dt))
+ CALL fieldadd(u1,v1,in%sx1+2,in%sx2,in%sx3/2,c2=REAL(Dt))
+ CALL fieldadd(u2,v2,in%sx1+2,in%sx2,in%sx3/2,c2=REAL(Dt))
+ CALL fieldadd(u3,v3,in%sx1+2,in%sx2,in%sx3/2,c2=REAL(Dt))
+ CALL tensorfieldadd(tau,moment,in%sx1,in%sx2,in%sx3/2,c2=REAL(Dt))
! keep track of the viscoelastic contribution alone
- CALL sliceadd(inter1(:,:,1),v1,sx1+2,sx2,sx3,int(oz/dx3)+1,c2=REAL(Dt))
- CALL sliceadd(inter2(:,:,1),v2,sx1+2,sx2,sx3,int(oz/dx3)+1,c2=REAL(Dt))
- CALL sliceadd(inter3(:,:,1),v3,sx1+2,sx2,sx3,int(oz/dx3)+1,c2=REAL(Dt))
+ CALL sliceadd(inter1(:,:,1),v1,in%sx1+2,in%sx2,in%sx3,int(in%oz/in%dx3)+1,c2=REAL(Dt))
+ CALL sliceadd(inter2(:,:,1),v2,in%sx1+2,in%sx2,in%sx3,int(in%oz/in%dx3)+1,c2=REAL(Dt))
+ CALL sliceadd(inter3(:,:,1),v3,in%sx1+2,in%sx2,in%sx3,int(in%oz/in%dx3)+1,c2=REAL(Dt))
! time increment
t=t+Dt
! next event
- IF (e .LT. ne) THEN
- IF (abs(t-events(e+1)%time) .LT. 1e-6) THEN
+ IF (e .LT. in%ne) THEN
+ IF (abs(t-in%events(e+1)%time) .LT. 1e-6) THEN
e=e+1
- events(e)%i=i
+ in%events(e)%i=i
PRINT '("coseismic event ",I3.3)', e
PRINT 0990
v1=0;v2=0;v3=0;t1=0;t2=0;t3=0;
- CALL dislocations(events(e),lambda,mu,beta,sx1,sx2,sx3, &
- dx1,dx2,dx3,v1,v2,v3,t1,t2,t3,tau)
- CALL traction(mu,events(e),sx1,sx2,dx1,dx2,t3)
+ CALL dislocations(in%events(e),in%lambda,in%mu, &
+ in%beta,in%sx1,in%sx2,in%sx3, &
+ in%dx1,in%dx2,in%dx3,v1,v2,v3,t1,t2,t3,tau)
+ CALL traction(in%mu,in%events(e),in%sx1,in%sx2,in%dx1,in%dx2,t3)
! apply the 3d elastic transfert function
- CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3,dx1,dx2,dx3,lambda,mu,gam)
+ CALL greenfunctioncowling(v1,v2,v3,t1,t2,t3, &
+ in%dx1,in%dx2,in%dx3,in%lambda,in%mu,in%gam)
! transfer solution
- CALL fieldadd(u1,v1,sx1+2,sx2,sx3/2)
- CALL fieldadd(u2,v2,sx1+2,sx2,sx3/2)
- CALL fieldadd(u3,v3,sx1+2,sx2,sx3/2)
+ CALL fieldadd(u1,v1,in%sx1+2,in%sx2,in%sx3/2)
+ CALL fieldadd(u2,v2,in%sx1+2,in%sx2,in%sx3/2)
+ CALL fieldadd(u3,v3,in%sx1+2,in%sx2,in%sx3/2)
END IF
END IF
! points are exported at all time steps
- IF (ALLOCATED(ptsname)) THEN
- CALL exportpoints(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3, &
- opts,ptsname,t,wdir,.false.,x0,y0,rot)
+ IF (ALLOCATED(in%ptsname)) THEN
+ CALL exportpoints(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3, &
+ in%opts,in%ptsname,t,in%wdir,.false.,in%x0,in%y0,in%rot)
END IF
! output only at discrete intervals (skip=0, odt>0),
! or every "skip" computational steps (skip>0, odt<0),
! or anytime a coseismic event occurs
- IF (isoutput(skip,t,i,odt,oi,events(e)%time)) THEN
+ IF (isoutput(in%skip,t,i,in%odt,oi,in%events(e)%time)) THEN
- CALL reporttime(1,t,reporttimefilename)
+ CALL reporttime(1,t,in%reporttimefilename)
! export
#ifdef TXT
- CALL exporttxt(u1,u2,u3,sx1,sx2,sx3/2,oz,dx3,oi,t,wdir,reportfilename)
+ IF (in%isoutputtxt) THEN
+ CALL exporttxt(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%oz,in%dx3,oi,t,in%wdir,in%reportfilename)
+ END IF
#endif
#ifdef XYZ
- CALL exportxyz(u1,u2,u3,sx1,sx2,sx3/2,oz,dx1,dx2,dx3,i,wdir)
- !CALL exportxyz(inter1,inter2,inter3,sx1,sx2,sx3/2,0.0_8,dx1,dx2,dx3,i,wdir)
+ IF (in%isoutputxyz) THEN
+ CALL exportxyz(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%oz,in%dx1,in%dx2,in%dx3,i,in%wdir)
+ !CALL exportxyz(inter1,inter2,inter3,in%sx1,in%sx2,in%sx3/2,0.0_8,in%dx1,in%dx2,in%dx3,i,in%wdir)
+ END IF
#endif
- CALL exporteigenstrain(gamma,nop,op,x0,y0,dx1,dx2,dx3,sx1,sx2,sx3/2,wdir,oi)
+ CALL exporteigenstrain(gamma,in%nop,in%op,in%x0,in%y0,in%dx1,in%dx2,in%dx3,in%sx1,in%sx2,in%sx3/2,in%wdir,oi)
#ifdef GRD
- CALL exportgrd(inter1,inter2,inter3,sx1,sx2,sx3/2, &
- dx1,dx2,dx3,0._8,x0,y0,wdir,oi,convention=2)
- CALL exportgrd(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3,oz,x0,y0,wdir,oi)
+ IF (in%isoutputgrd) THEN
+ CALL exportgrd(inter1,inter2,inter3,in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,0._8,in%x0,in%y0,in%wdir,oi,convention=2)
+ CALL exportgrd(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%oz,in%x0,in%y0,in%wdir,oi)
+ END IF
#endif
#ifdef PROJ
- CALL exportproj(inter1,inter2,inter3,sx1,sx2,sx3/2, &
- dx1,dx2,dx3,oz,x0,y0, &
- lon0,lat0,zone,umult,wdir,oi,convention=2)
- CALL exportproj(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3,oz,x0,y0, &
- lon0,lat0,zone,umult,wdir,oi)
+ IF (in%isoutputproj) THEN
+ CALL exportproj(inter1,inter2,inter3,in%sx1,in%sx2,in%sx3/2, &
+ in%dx1,in%dx2,in%dx3,in%oz,in%x0,in%y0, &
+ in%lon0,in%lat0,in%zone,in%umult,in%wdir,oi,convention=2)
+ CALL exportproj(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%oz,in%x0,in%y0, &
+ in%lon0,in%lat0,in%zone,in%umult,in%wdir,oi)
+ END IF
#endif
#ifdef VTK
- WRITE (digit,'(I3.3)') oi
- j=INDEX(wdir," ")
- ! export total displacement in VTK XML format
- vcfilename=wdir(1:j-1)//"/disp-"//digit//".vtr"
- CALL exportvtk_vectors(u1,u2,u3,sx1,sx2,sx3/4,dx1,dx2,dx3,8,8,8,vcfilename)
- !CALL exportvtk_vectors_slice(u1,u2,u3,sx1,sx2,sx3/2,dx1,dx2,dx3,oz,8,8,vcfilename)
+ IF (in%isoutputvtk) THEN
+ WRITE (digit,'(I3.3)') oi
+ ! export total displacement in VTK XML format
+ vcfilename=trim(in%wdir)//"/disp-"//digit//".vtr"
+ CALL exportvtk_vectors(u1,u2,u3,in%sx1,in%sx2,in%sx3/4,in%dx1,in%dx2,in%dx3,8,8,8,vcfilename)
+ !CALL exportvtk_vectors_slice(u1,u2,u3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%oz,8,8,vcfilename)
- ! export instantaneous velocity in VTK XML format
- vcfilename=wdir(1:j-1)//"/vel-"//digit//".vtr"
- CALL exportvtk_vectors(v1,v2,v3,sx1,sx2,sx3/4,dx1,dx2,dx3,8,8,8,vcfilename)
- !CALL exportvtk_vectors_slice(v1,v2,v3,sx1,sx2,sx3/2,dx1,dx2,dx3,oz,8,8,vcfilename)
+ ! export instantaneous velocity in VTK XML format
+ vcfilename=trim(in%wdir)//"/vel-"//digit//".vtr"
+ CALL exportvtk_vectors(v1,v2,v3,in%sx1,in%sx2,in%sx3/4,in%dx1,in%dx2,in%dx3,8,8,8,vcfilename)
+ !CALL exportvtk_vectors_slice(v1,v2,v3,in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%oz,8,8,vcfilename)
+ END IF
#endif
- PRINT 1101,i,Dt,maxwell,t,interval, &
- tensoramplitude(moment,dx1,dx2,dx3), &
- tensoramplitude(tau,dx1,dx2,dx3)
+ PRINT 1101,i,Dt,maxwell,t,in%interval, &
+ tensoramplitude(moment,in%dx1,in%dx2,in%dx3), &
+ tensoramplitude(tau,in%dx1,in%dx2,in%dx3)
! update output counter
oi=oi+1
ELSE
- PRINT 1100,i,Dt,maxwell,t,interval, &
- tensoramplitude(moment,dx1,dx2,dx3), &
- tensoramplitude(tau,dx1,dx2,dx3)
+ PRINT 1100,i,Dt,maxwell,t,in%interval, &
+ tensoramplitude(moment,in%dx1,in%dx2,in%dx3), &
+ tensoramplitude(tau,in%dx1,in%dx2,in%dx3)
END IF
END DO
100 CONTINUE
- DO i=1,ne
- IF (ALLOCATED(events(i)%s)) DEALLOCATE(events(i)%s,events(i)%sc)
- IF (ALLOCATED(events(i)%ts)) DEALLOCATE(events(i)%ts,events(i)%tsc)
+ DO i=1,in%ne
+ IF (ALLOCATED(in%events(i)%s)) DEALLOCATE(in%events(i)%s,in%events(i)%sc)
+ IF (ALLOCATED(in%events(i)%ts)) DEALLOCATE(in%events(i)%ts,in%events(i)%tsc)
END DO
- IF (ALLOCATED(events)) DEALLOCATE(events)
+ IF (ALLOCATED(in%events)) DEALLOCATE(in%events)
! free memory
IF (ALLOCATED(gamma)) DEALLOCATE(gamma)
- IF (ALLOCATED(opts)) DEALLOCATE(opts)
- IF (ALLOCATED(op)) DEALLOCATE(op)
- IF (ALLOCATED(n)) DEALLOCATE(n)
- IF (ALLOCATED(stressstruc)) DEALLOCATE(stressstruc)
- IF (ALLOCATED(linearstruc)) DEALLOCATE(linearstruc)
- IF (ALLOCATED(nonlinearstruc)) DEALLOCATE(nonlinearstruc)
- IF (ALLOCATED(faultcreepstruc)) DEALLOCATE(faultcreepstruc)
+ IF (ALLOCATED(in%opts)) DEALLOCATE(in%opts)
+ IF (ALLOCATED(in%op)) DEALLOCATE(in%op)
+ IF (ALLOCATED(in%n)) DEALLOCATE(in%n)
+ IF (ALLOCATED(in%stressstruc)) DEALLOCATE(in%stressstruc)
+ IF (ALLOCATED(in%linearstruc)) DEALLOCATE(in%linearstruc)
+ IF (ALLOCATED(in%nonlinearstruc)) DEALLOCATE(in%nonlinearstruc)
+ IF (ALLOCATED(in%faultcreepstruc)) DEALLOCATE(in%faultcreepstruc)
IF (ALLOCATED(sig)) DEALLOCATE(sig)
IF (ALLOCATED(tau)) DEALLOCATE(tau)
IF (ALLOCATED(moment)) DEALLOCATE(moment)
- IF (ALLOCATED(stresslayer)) DEALLOCATE(stresslayer)
- IF (ALLOCATED(linearlayer)) DEALLOCATE(linearlayer)
- IF (ALLOCATED(nonlinearlayer)) DEALLOCATE(nonlinearlayer)
- IF (ALLOCATED(faultcreeplayer)) DEALLOCATE(faultcreeplayer)
- DEALLOCATE(v1,v2,v3,t1,t2,t3)
- DEALLOCATE(u1,u2,u3)
- DEALLOCATE(inter1,inter2,inter3)
+ IF (ALLOCATED(in%stresslayer)) DEALLOCATE(in%stresslayer)
+ IF (ALLOCATED(in%linearlayer)) DEALLOCATE(in%linearlayer)
+ IF (ALLOCATED(in%nonlinearlayer)) DEALLOCATE(in%nonlinearlayer)
+ IF (ALLOCATED(in%faultcreeplayer)) DEALLOCATE(in%faultcreeplayer)
+ IF (ALLOCATED(v1)) DEALLOCATE(v1,v2,v3,t1,t2,t3)
+ IF (ALLOCATED(u1)) DEALLOCATE(u1,u2,u3)
+ IF (ALLOCATED(inter1)) DEALLOCATE(inter1,inter2,inter3)
#ifdef FFTW3_THREADS
@@ -672,39 +683,6 @@ 1200 FORMAT ("--------------------------
1200 FORMAT ("----------------------------------------------------------------------------")
CONTAINS
-
- !--------------------------------------------------------------------
- ! subroutine eqbf_mask
- ! fills an array with positive values if some linear/nonlinear/creep
- ! is expected at the corresponding depth, zero otherwise.
- !
- ! the mask can be given to the routine "equivalentBodyForce" to skip
- ! these depths where no creep happens.
- !--------------------------------------------------------------------
- SUBROUTINE eqbf_mask(mask,sx)
- INTEGER, INTENT(IN) :: sx
- REAL*4, DIMENSION(sx), INTENT(OUT) :: mask
-
- IF (ALLOCATED(linearstruc)) THEN
- DO k=1,sx
- mask(k)=MAX(mask(k),REAL(linearstruc(k)%gammadot0,4))
- END DO
- END IF
- IF (ALLOCATED(nonlinearstruc)) THEN
- DO k=1,sx
- mask(k)=MAX(mask(k),REAL(nonlinearstruc(k)%gammadot0,4))
- END DO
- END IF
- IF (ALLOCATED(faultcreepstruc)) THEN
- DO k=1,sx
- mask(k)=MAX(mask(k),REAL(faultcreepstruc(k)%gammadot0,4))
- END DO
- END IF
-
- ! smooth the mask in the depth direction
- mask(1:sx-2)=(mask(1:sx-2)+mask(2:sx-1)+mask(3:sx))/3._4
-
- END SUBROUTINE eqbf_mask
!---------------------------------------------------------------------
! subroutine Traction
@@ -741,7 +719,7 @@ CONTAINS
! not reverse in the postseismic time.
!--------------------------------------------------------------------
SUBROUTINE dislocations(event,lambda,mu,beta,sx1,sx2,sx3,dx1,dx2,dx3, &
- v1,v2,v3,t1,t2,t3,tau,factor,eigenstress)
+ v1,v2,v3,t1,t2,t3,tau,factor,eigenstress)
TYPE(EVENT_STRUC), INTENT(IN) :: event
INTEGER, INTENT(IN) :: sx1,sx2,sx3
REAL*8, INTENT(IN) :: lambda,mu,beta,dx1,dx2,dx3
@@ -855,946 +833,6 @@ CONTAINS
END SUBROUTINE dislocations
- SUBROUTINE init(unit)
- INTEGER, OPTIONAL, INTENT(INOUT) :: unit
-
- INTEGER :: k,iostatus,i,e
- CHARACTER(180) :: dataline
-#ifdef VTK
- INTEGER :: j
- CHARACTER(3) :: digit
-#endif
- INTEGER :: iunit
-!$ INTEGER :: omp_get_num_procs,omp_get_max_threads
- REAL*8 :: dummy,dum1,dum2
-
- ! default is standard input
- IF (.NOT. PRESENT(unit)) THEN
- iunit=5
- ELSE
- iunit=unit
- END IF
-
- PRINT 2000
- PRINT '(" nonlinear viscoelastic postseismic relaxation")'
-#ifdef FFTW3
-#ifdef FFTW3_THREADS
- PRINT '(" * FFTW3 (multi-threaded) implementation of the FFT")'
-#else
- PRINT '(" * FFTW3 implementation of the FFT")'
-#endif
-#else
-#ifdef SGI_FFT
- PRINT '(" * SGI_FFT implementation of the FFT")'
-#else
-#ifdef IMKL_FFT
- PRINT '(" * Intel MKL implementation of the FFT")'
-#else
- PRINT '(" * fourt implementation of the FFT")'
-#endif
-#endif
-#endif
-!$ PRINT '(" * parallel OpenMP implementation with ",I3.3,"/",I3.3," threads")', &
-!$ omp_get_max_threads(),omp_get_num_procs()
-#ifdef GRD
- PRINT '(" * export to GRD format")'
-#endif
-#ifdef TXT
- PRINT '(" * export to TXT format")'
-#endif
-#ifdef VTK
- PRINT '(" * export to VTK format")'
-#endif
-#ifdef PROJ
- PRINT '(" * export to longitude/latitude text format")'
-#endif
- PRINT 2000
-
- PRINT '(a)', "grid dimension (sx1,sx2,sx3)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) sx1,sx2,sx3
- PRINT '(3I5)', sx1,sx2,sx3
-
- PRINT '(a)', "sampling (dx1,dx2,dx3), smoothing (beta, nyquist)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) dx1,dx2,dx3,beta,nyquist
- PRINT '(5ES9.2E1)', dx1,dx2,dx3,beta,nyquist
-
- PRINT '(a)', "origin position (x0,y0) and rotation"
- CALL getdata(iunit,dataline)
- READ (dataline,*) x0, y0, rot
- PRINT '(3ES9.2E1)', x0, y0, rot
-
-#ifdef PROJ
- PRINT '(a)', "geographic origin (longitude, latitude, UTM zone, unit)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) lon0,lat0,zone,umult
- PRINT '(2ES9.2E1,I3.2,ES9.2E1)',lon0,lat0,zone,umult
- IF (zone.GT.60 .OR. zone.LT.1) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("invalid UTM zone ",I," (1<=zone<=60. exiting.)")') zone
- STOP 1
- END IF
-#endif
-
- PRINT '(a)', "observation depth (displacement and stress)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) oz,ozs
- PRINT '(2ES9.2E1)', oz,ozs
-
- PRINT '(a)', "output directory"
- CALL getdata(iunit,dataline)
- READ (dataline,'(a)') wdir
- i=INDEX(wdir," ")
- reporttimefilename=wdir(1:i-1)//"/time.txt"
- reportfilename=wdir(1:i-1)//"/report.txt"
- logfilename=wdir(1:i-1)//"/relax.log"
- inputfilename=wdir(1:i-1)//"/relax.inp"
-#ifdef TXT
- PRINT '(" ",a," (report: ",a,")")', wdir(1:i-1),reportfilename(1:i+10)
-#else
- PRINT '(" ",a," (time report: ",a,")")', wdir(1:i-1),reporttimefilename(1:i+8)
-#endif
-
- ! test write permissions on output directory
- OPEN (UNIT=14,FILE=reportfilename,POSITION="APPEND",&
- IOSTAT=iostatus,FORM="FORMATTED")
- IF (iostatus>0) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("unable to access ",a)') reporttimefilename(1:i+10)
- STOP 1
- END IF
- CLOSE(14)
- ! end test
-
-#ifdef VTK
- cgfilename=wdir(1:i-1)//"/cgrid.vtp"
- CALL exportvtk_grid(sx1,sx2,sx3,dx1,dx2,dx3,x0,y0,cgfilename)
-#endif
-
- PRINT '(a)', "lambda, mu, gamma (gamma = (1 - nu) rho g / mu)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) lambda,mu,gam
- PRINT '(3ES10.2E2)',lambda,mu,gam
-
- PRINT '(a)', "time interval, (positive time step) or (negative skip, scaling)"
- CALL getdata(unit,dataline)
- READ (dataline,*) interval, odt
- IF (odt .LT. 0.) THEN
- READ (dataline,*) dum1, dum2, tscale
- skip=ceiling(-odt)
- PRINT '(ES9.2E1," (output every ",I3.3," steps, dt scaled by ",ES7.2E1,")")', &
- interval,skip,tscale
- ELSE
- PRINT '(ES9.2E1," (output every ",ES9.2E1," time unit)")', interval,odt
- END IF
-
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! O B S E R V A T I O N P L A N E S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of observation planes"
- CALL getdata(unit,dataline)
- READ (dataline,*) nop
- PRINT '(I5)', nop
- IF (nop .gt. 0) THEN
- ALLOCATE(op(nop),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the observation plane list"
- PRINT 2000
- PRINT 2100
- PRINT 2000
- DO k=1,nop
- CALL getdata(unit,dataline)
- READ (dataline,*) i,op(k)%x,op(k)%y,op(k)%z,&
- op(k)%length,op(k)%width,op(k)%strike,op(k)%dip
-
- PRINT '(I3.3," ",5ES9.2E1,2f7.1)', &
- k,op(k)%x,op(k)%y,op(k)%z, &
- op(k)%length,op(k)%width,op(k)%strike,op(k)%dip
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,*) "error in input file: plane index misfit", k,"<>",i
- WRITE (0,*) op(k)
- STOP 1
- END IF
-
- ! comply to Wang's convention
- CALL wangconvention(dummy,op(k)%x,op(k)%y,op(k)%z,&
- op(k)%length,op(k)%width,op(k)%strike,op(k)%dip,dummy,rot)
-
- END DO
- END IF
-
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! O B S E R V A T I O N P O I N T S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of observation points"
- CALL getdata(iunit,dataline)
- READ (dataline,*) npts
- PRINT '(I5)', npts
- IF (npts .gt. 0) THEN
- ALLOCATE(opts(npts),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the observation point list"
- ALLOCATE(ptsname(npts),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the list of point name"
-
- PRINT 2000
- PRINT 2300
- PRINT 2000
- DO k=1,npts
- CALL getdata(iunit,dataline)
- READ (dataline,*) i,ptsname(k),opts(k)%v1,opts(k)%v2,opts(k)%v3
-
- PRINT '(I3.3," ",A4,3ES9.2E1)', i,ptsname(k), &
- opts(k)%v1,opts(k)%v2,opts(k)%v3
-
- ! shift and rotate coordinates
- opts(k)%v1=opts(k)%v1-x0
- opts(k)%v2=opts(k)%v2-y0
- CALL rotation(opts(k)%v1,opts(k)%v2,rot)
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: points index misfit")')
- STOP 1
- END IF
- END DO
-
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! P R E S T R E S S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of prestress interfaces"
- CALL getdata(unit,dataline)
- READ (dataline,*) nps
- PRINT '(I5)', nps
-
- IF (nps .GT. 0) THEN
- ALLOCATE(stresslayer(nps),stressstruc(sx3/2),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the stress layer structure"
-
- PRINT 2000
- PRINT '(a)', "no. depth sigma11 sigma12 sigma13 sigma22 sigma23 sigma33"
- PRINT 2000
- DO k=1,nps
- CALL getdata(unit,dataline)
- READ (dataline,*) i,stresslayer(k)%z, &
- stresslayer(k)%t%s11, stresslayer(k)%t%s12, &
- stresslayer(k)%t%s13, stresslayer(k)%t%s22, &
- stresslayer(k)%t%s23, stresslayer(k)%t%s33
-
- PRINT '(I3.3,7ES9.2E1)', i, &
- stresslayer(k)%z, &
- stresslayer(k)%t%s11, stresslayer(k)%t%s12, &
- stresslayer(k)%t%s13, stresslayer(k)%t%s22, &
- stresslayer(k)%t%s23, stresslayer(k)%t%s33
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: index misfit")')
- STOP 1
- END IF
- END DO
- END IF
-
-
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! L I N E A R V I S C O U S I N T E R F A C E
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of linear viscous interfaces"
- CALL getdata(unit,dataline)
- READ (dataline,*) nv
- PRINT '(I5)', nv
-
- IF (nv .GT. 0) THEN
- ALLOCATE(linearlayer(nv),linearstruc(sx3/2),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the layer structure"
-
- PRINT 2000
- PRINT '(a)', "no. depth gamma0 cohesion"
- PRINT 2000
- DO k=1,nv
- CALL getdata(unit,dataline)
- READ (dataline,*) i,linearlayer(k)%z, &
- linearlayer(k)%gammadot0, linearlayer(k)%cohesion
-
- linearlayer(k)%stressexponent=1
-
- PRINT '(I3.3,3ES10.2E2)', i, &
- linearlayer(k)%z, &
- linearlayer(k)%gammadot0, &
- linearlayer(k)%cohesion
-
- ! check positive strain rates
- IF (linearlayer(k)%gammadot0 .LT. 0) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: strain rates must be positive")')
- STOP 1
- END IF
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: index misfit")')
- STOP 1
- END IF
-#ifdef VTK
- ! export the viscous layer in VTK format
- j=INDEX(wdir," ")
- WRITE (digit,'(I3.3)') k
-
- rffilename=wdir(1:j-1)//"/linearlayer-"//digit//".vtp"
- CALL exportvtk_rectangle(0.d0,0.d0,linearlayer(k)%z, &
- DBLE(sx1)*dx1,DBLE(sx2)*dx2, &
- 0._8,1.57d0,rffilename)
-#endif
- END DO
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! L I N E A R W E A K Z O N E S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of linear weak zones (nlwz)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) nlwz
- PRINT '(I5)', nlwz
- IF (nlwz .GT. 0) THEN
- ALLOCATE(linearweakzone(nlwz),linearweakzonec(nlwz),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the linear weak zones"
- PRINT 2000
- PRINT '(a)', "no. dgammadot0 x1 x2 x3 length width thickn. strike dip"
- PRINT 2000
- DO k=1,nlwz
- CALL getdata(iunit,dataline)
- READ (dataline,*) i, &
- linearweakzone(k)%dgammadot0, &
- linearweakzone(k)%x,linearweakzone(k)%y,linearweakzone(k)%z,&
- linearweakzone(k)%length,linearweakzone(k)%width,linearweakzone(k)%thickness, &
- linearweakzone(k)%strike,linearweakzone(k)%dip
-
- linearweakzonec(k)=linearweakzone(k)
-
- PRINT '(I3.3,4ES9.2E1,3ES8.2E1,f7.1,f6.1)',k,&
- linearweakzone(k)%dgammadot0, &
- linearweakzone(k)%x,linearweakzone(k)%y,linearweakzone(k)%z, &
- linearweakzone(k)%length,linearweakzone(k)%width, &
- linearweakzone(k)%thickness, &
- linearweakzone(k)%strike,linearweakzone(k)%dip
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: source index misfit")')
- STOP 1
- END IF
- ! comply to Wang's convention
- CALL wangconvention( &
- dummy, &
- linearweakzone(k)%x,linearweakzone(k)%y,linearweakzone(k)%z, &
- linearweakzone(k)%length,linearweakzone(k)%width, &
- linearweakzone(k)%strike,linearweakzone(k)%dip,dummy,rot)
-#ifdef VTK
- ! export the ductile zone in VTK format
- j=INDEX(wdir," ")-1
- WRITE (digit,'(I3.3)') k
-
- rffilename=wdir(1:j)//"/weakzone-"//digit//".vtp"
- CALL exportvtk_brick(linearweakzone(k)%x,linearweakzone(k)%y,linearweakzone(k)%z, &
- linearweakzone(k)%length,linearweakzone(k)%width,linearweakzone(k)%thickness, &
- linearweakzone(k)%strike,linearweakzone(k)%dip,rffilename)
-#endif
- END DO
- END IF
- END IF ! end linear viscous
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! N O N L I N E A R V I S C O U S I N T E R F A C E
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of nonlinear viscous interfaces"
- CALL getdata(unit,dataline)
- READ (dataline,*) npl
- PRINT '(I5)', npl
-
- IF (npl .GT. 0) THEN
- ALLOCATE(nonlinearlayer(npl),nonlinearstruc(sx3/2),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the layer structure"
-
- PRINT 2000
- PRINT '(a)', "no. depth gamma0 power cohesion"
- PRINT 2000
- DO k=1,npl
- CALL getdata(unit,dataline)
-
- READ (dataline,*) i,nonlinearlayer(k)%z, &
- nonlinearlayer(k)%gammadot0, &
- nonlinearlayer(k)%stressexponent, &
- nonlinearlayer(k)%cohesion
-
- PRINT '(I3.3,4ES10.2E2)', i, &
- nonlinearlayer(k)%z, &
- nonlinearlayer(k)%gammadot0, &
- nonlinearlayer(k)%stressexponent, &
- nonlinearlayer(k)%cohesion
-
- ! check positive strain rates
- IF (nonlinearlayer(k)%gammadot0 .LT. 0) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: strain rates must be positive")')
- STOP 1
- END IF
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: index misfit")')
- STOP 1
- END IF
-
- END DO
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! N O N L I N E A R W E A K Z O N E S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of nonlinear weak zones (nnlwz)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) nnlwz
- PRINT '(I5)', nnlwz
- IF (nnlwz .GT. 0) THEN
- ALLOCATE(nonlinearweakzone(nnlwz),nonlinearweakzonec(nnlwz),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the nonlinear weak zones"
- PRINT 2000
- PRINT '(a)', "no. dgammadot0 x1 x2 x3 length width thickn. strike dip"
- PRINT 2000
- DO k=1,nnlwz
- CALL getdata(iunit,dataline)
- READ (dataline,*) i, &
- nonlinearweakzone(k)%dgammadot0, &
- nonlinearweakzone(k)%x,nonlinearweakzone(k)%y,nonlinearweakzone(k)%z,&
- nonlinearweakzone(k)%length,nonlinearweakzone(k)%width,nonlinearweakzone(k)%thickness, &
- nonlinearweakzone(k)%strike,nonlinearweakzone(k)%dip
-
- nonlinearweakzonec(k)=nonlinearweakzone(k)
-
- PRINT '(I3.3,4ES9.2E1,3ES8.2E1,f7.1,f6.1)',k,&
- nonlinearweakzone(k)%dgammadot0, &
- nonlinearweakzone(k)%x,nonlinearweakzone(k)%y,nonlinearweakzone(k)%z, &
- nonlinearweakzone(k)%length,nonlinearweakzone(k)%width, &
- nonlinearweakzone(k)%thickness, &
- nonlinearweakzone(k)%strike,nonlinearweakzone(k)%dip
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: source index misfit")')
- STOP 1
- END IF
- ! comply to Wang's convention
- CALL wangconvention( &
- dummy, &
- nonlinearweakzone(k)%x,nonlinearweakzone(k)%y,nonlinearweakzone(k)%z, &
- nonlinearweakzone(k)%length,nonlinearweakzone(k)%width, &
- nonlinearweakzone(k)%strike,nonlinearweakzone(k)%dip,dummy,rot)
- END DO
- END IF
- END IF ! end nonlinear viscous
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! F A U L T C R E E P
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of fault creep interfaces"
- CALL getdata(unit,dataline)
- READ (dataline,*) nfc
- PRINT '(I5)', nfc
-
- IF (nfc .GT. 0) THEN
- ALLOCATE(faultcreeplayer(nfc),faultcreepstruc(sx3/2),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the layer structure"
-
- PRINT 2000
- PRINT '(a)', "no. depth gamma0 (a-b)sig friction cohesion"
- PRINT 2000
- DO k=1,nfc
- CALL getdata(unit,dataline)
- READ (dataline,*) i,faultcreeplayer(k)%z, &
- faultcreeplayer(k)%gammadot0, &
- faultcreeplayer(k)%stressexponent, &
- faultcreeplayer(k)%friction, &
- faultcreeplayer(k)%cohesion
-
- PRINT '(I3.3,5ES9.2E1)', i, &
- faultcreeplayer(k)%z, &
- faultcreeplayer(k)%gammadot0, &
- faultcreeplayer(k)%stressexponent, &
- faultcreeplayer(k)%friction, &
- faultcreeplayer(k)%cohesion
-
- ! check positive strain rates
- IF (faultcreeplayer(k)%gammadot0 .LT. 0) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: slip rates must be positive")')
- STOP 1
- END IF
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: index misfit")')
- STOP 1
- END IF
-
- END DO
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! A F T E R S L I P P L A N E S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of afterslip planes"
- CALL getdata(unit,dataline)
- READ (dataline,*) np
- PRINT '(I5)', np
-
- IF (np .gt. 0) THEN
- ALLOCATE(n(np),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the plane list"
-
- PRINT 2000
- PRINT 2500
- PRINT 2000
-
- DO k=1,np
- CALL getdata(unit,dataline)
- READ (dataline,*) i,n(k)%x,n(k)%y,n(k)%z,&
- n(k)%length,n(k)%width,n(k)%strike,n(k)%dip,n(k)%rake
-
- PRINT '(I3.3," ",5ES9.2E1,3f7.1)',i, &
- n(k)%x,n(k)%y,n(k)%z, &
- n(k)%length,n(k)%width,n(k)%strike,n(k)%dip,n(k)%rake
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: plane index misfit")')
- STOP 1
- END IF
-
- ! modify rake for consistency with slip model
- IF (n(k)%rake .GE. 0.d0) THEN
- n(k)%rake=n(k)%rake-180.d0
- ELSE
- n(k)%rake=n(k)%rake+180.d0
- END IF
-
- ! comply to Wang's convention
- CALL wangconvention(dummy,n(k)%x,n(k)%y,n(k)%z,&
- n(k)%length,n(k)%width,n(k)%strike,n(k)%dip,n(k)%rake,rot)
-
-#ifdef VTK
- ! export the afterslip segment in VTK format
- j=INDEX(wdir," ")
- WRITE (digit,'(I3.3)') k
-
- rffilename=wdir(1:j-1)//"/aplane-"//digit//".vtp"
- CALL exportvtk_rectangle(n(k)%x,n(k)%y,n(k)%z,n(k)%length,n(k)%width, &
- n(k)%strike,n(k)%dip,rffilename)
-#endif
-
- END DO
- END IF
-
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! I N T E R - S E I S M I C L O A D I N G
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- minlength=sx1*dx1+sx2*dx2
- minwidth=sx3*dx3
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! S H E A R S O U R C E S R A T E
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of inter-seismic strike-slip segments"
- CALL getdata(iunit,dataline)
- READ (dataline,*) inter%ns
- PRINT '(I5)', inter%ns
- IF (inter%ns .GT. 0) THEN
- ALLOCATE(inter%s(inter%ns),inter%sc(inter%ns),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the source list"
- PRINT 2000
- PRINT '(a)',"no. slip xs ys zs length width strike dip rake"
- PRINT 2000
- DO k=1,inter%ns
- CALL getdata(iunit,dataline)
- READ (dataline,*) i,inter%s(k)%slip, &
- inter%s(k)%x,inter%s(k)%y,inter%s(k)%z, &
- inter%s(k)%length,inter%s(k)%width, &
- inter%s(k)%strike,inter%s(k)%dip,inter%s(k)%rake
- ! copy the input format for display
- inter%sc(k)=inter%s(k)
-
- PRINT '(I3.3,4ES9.2E1,2ES8.2E1,f7.1,f6.1,f7.1)',i, &
- inter%sc(k)%slip,&
- inter%sc(k)%x,inter%sc(k)%y,inter%sc(k)%z, &
- inter%sc(k)%length,inter%sc(k)%width, &
- inter%sc(k)%strike,inter%sc(k)%dip, &
- inter%sc(k)%rake
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: source index misfit")')
- STOP 1
- END IF
- IF (MAX(inter%s(k)%length,inter%s(k)%width) .LE. 0._8) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: lengths must be positive.")')
- STOP 1
- END IF
- IF (inter%s(k)%length .lt. minlength) THEN
- minlength=inter%s(k)%length
- END IF
- IF (inter%s(k)%width .lt. minwidth ) THEN
- minwidth =inter%s(k)%width
- END IF
-
- ! smooth out the slip distribution
- CALL antialiasingfilter(inter%s(k)%slip, &
- inter%s(k)%length,inter%s(k)%width, &
- dx1,dx2,dx3,nyquist)
-
- ! comply to Wang's convention
- CALL wangconvention(inter%s(k)%slip, &
- inter%s(k)%x,inter%s(k)%y,inter%s(k)%z, &
- inter%s(k)%length,inter%s(k)%width, &
- inter%s(k)%strike,inter%s(k)%dip, &
- inter%s(k)%rake,rot)
-
- END DO
- PRINT 2000
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! T E N S I L E S O U R C E S R A T E
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of inter-seismic tensile segments"
- CALL getdata(iunit,dataline)
- READ (dataline,*) inter%nt
- PRINT '(I5)', inter%nt
- IF (inter%nt .GT. 0) THEN
- ALLOCATE(inter%ts(inter%nt),inter%tsc(inter%nt),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the tensile source list"
- PRINT 2000
- PRINT '(a)',"no. opening xs ys zs length width strike dip"
- PRINT 2000
- DO k=1,inter%nt
- CALL getdata(iunit,dataline)
- READ (dataline,*) i,inter%ts(k)%slip, &
- inter%ts(k)%x,inter%ts(k)%y,inter%ts(k)%z, &
- inter%ts(k)%length,inter%ts(k)%width, &
- inter%ts(k)%strike,inter%ts(k)%dip
- ! copy the input format for display
- inter%tsc(k)=inter%ts(k)
-
- PRINT '(I3.3,4ES9.2E1,2ES8.2E1,f7.1,f6.1)', i, &
- inter%tsc(k)%slip,&
- inter%tsc(k)%x,inter%tsc(k)%y,inter%tsc(k)%z, &
- inter%tsc(k)%length,inter%tsc(k)%width, &
- inter%tsc(k)%strike,inter%tsc(k)%dip
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: tensile source index misfit")')
- STOP 1
- END IF
- IF (MAX(inter%ts(k)%length,inter%ts(k)%width) .LE. 0._8) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: lengths must be positive.")')
- STOP 1
- END IF
- IF (inter%ts(k)%length .lt. minlength) THEN
- minlength=inter%ts(k)%length
- END IF
- IF (inter%ts(k)%width .lt. minwidth) THEN
- minwidth =inter%ts(k)%width
- END IF
-
- ! smooth out the slip distribution
- CALL antialiasingfilter(inter%ts(k)%slip, &
- inter%ts(k)%length,inter%ts(k)%width, &
- dx1,dx2,dx3,nyquist)
-
- ! comply to Wang's convention
- CALL wangconvention(inter%ts(k)%slip, &
- inter%ts(k)%x,inter%ts(k)%y,inter%ts(k)%z, &
- inter%ts(k)%length,inter%ts(k)%width, &
- inter%ts(k)%strike,inter%ts(k)%dip,dummy,rot)
-
- END DO
- PRINT 2000
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! C 0 - S E I S M I C E V E N T S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of events"
- CALL getdata(iunit,dataline)
- READ (dataline,*) ne
- PRINT '(I5)', ne
- IF (ne .GT. 0) ALLOCATE(events(ne),STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the event list"
-
- DO e=1,ne
- IF (1 .NE. e) THEN
- PRINT '("time of next coseismic event")'
- CALL getdata(iunit,dataline)
- READ (dataline,*) events(e)%time
-
- IF (0 .EQ. skip) THEN
- ! change event time to multiples of output time step
- events(e)%time=fix(events(e)%time/odt)*odt
- END IF
-
- PRINT '(ES9.2E1," (multiple of ",ES9.2E1,")")', &
- events(e)%time,odt
-
- IF (events(e)%time .LE. events(e-1)%time) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'(a,a)') "input file error. ", &
- "coseismic source time must increase. interrupting."
- STOP 1
- END IF
- ELSE
- events(1)%time=0._8
- events(1)%i=0
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! S H E A R S O U R C E S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of coseismic strike-slip segments (ns)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) events(e)%ns
- PRINT '(I5)', events(e)%ns
- IF (events(e)%ns .GT. 0) THEN
- ALLOCATE(events(e)%s(events(e)%ns),events(e)%sc(events(e)%ns), &
- STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the source list"
- PRINT 2000
- PRINT '(a)',"no. slip xs ys zs length width strike dip rake"
- PRINT 2000
- DO k=1,events(e)%ns
- CALL getdata(iunit,dataline)
- READ (dataline,*) i,events(e)%s(k)%slip, &
- events(e)%s(k)%x,events(e)%s(k)%y,events(e)%s(k)%z, &
- events(e)%s(k)%length,events(e)%s(k)%width, &
- events(e)%s(k)%strike,events(e)%s(k)%dip,events(e)%s(k)%rake
- ! copy the input format for display
- events(e)%sc(k)=events(e)%s(k)
-
- PRINT '(I3.3,4ES9.2E1,2ES8.2E1,f7.1,f6.1,f7.1)',i, &
- events(e)%sc(k)%slip,&
- events(e)%sc(k)%x,events(e)%sc(k)%y,events(e)%sc(k)%z, &
- events(e)%sc(k)%length,events(e)%sc(k)%width, &
- events(e)%sc(k)%strike,events(e)%sc(k)%dip, &
- events(e)%sc(k)%rake
-
- IF (i .ne. k) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("invalid shear source definition ")')
- WRITE (0,'("error in input file: source index misfit")')
- STOP 1
- END IF
- IF (MAX(events(e)%s(k)%length,events(e)%s(k)%width) .LE. 0._8) THEN
- WRITE_DEBUG_INFO
- WRITE (0,'("error in input file: lengths must be positive.")')
- STOP 1
- END IF
- IF (events(e)%s(k)%length .lt. minlength) THEN
- minlength=events(e)%s(k)%length
- END IF
- IF (events(e)%s(k)%width .lt. minwidth ) THEN
- minwidth =events(e)%s(k)%width
- END IF
-
- ! smooth out the slip distribution
- CALL antialiasingfilter(events(e)%s(k)%slip, &
- events(e)%s(k)%length,events(e)%s(k)%width, &
- dx1,dx2,dx3,nyquist)
-
- ! comply to Wang's convention
- CALL wangconvention(events(e)%s(k)%slip, &
- events(e)%s(k)%x,events(e)%s(k)%y,events(e)%s(k)%z, &
- events(e)%s(k)%length,events(e)%s(k)%width, &
- events(e)%s(k)%strike,events(e)%s(k)%dip, &
- events(e)%s(k)%rake,rot)
-
- END DO
-
-#ifdef VTK
- ! export the fault segments in VTK format for the current event
- j=INDEX(wdir," ")
- WRITE (digit,'(I3.3)') e
-
- rffilename=wdir(1:j-1)//"/rfaults-"//digit//".vtp"
- CALL exportvtk_rfaults(events(e),rffilename)
-#endif
- rffilename=wdir(1:j-1)//"/rfaults-"//digit//".xy"
- CALL exportxy_rfaults(events(e),rffilename)
-
- PRINT 2000
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! T E N S I L E S O U R C E S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of coseismic tensile segments (nt)"
- CALL getdata(iunit,dataline)
- READ (dataline,*) events(e)%nt
- PRINT '(I5)', events(e)%nt
- IF (events(e)%nt .GT. 0) THEN
- ALLOCATE(events(e)%ts(events(e)%nt),events(e)%tsc(events(e)%nt), &
- STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the tensile source list"
- PRINT 2000
- PRINT '(a)',"no. opening xs ys zs length width strike dip"
- PRINT 2000
- DO k=1,events(e)%nt
- CALL getdata(iunit,dataline)
- READ (dataline,*) i,events(e)%ts(k)%slip, &
- events(e)%ts(k)%x,events(e)%ts(k)%y,events(e)%ts(k)%z, &
- events(e)%ts(k)%length,events(e)%ts(k)%width, &
- events(e)%ts(k)%strike,events(e)%ts(k)%dip
- ! copy the input format for display
- events(e)%tsc(k)=events(e)%ts(k)
-
- PRINT '(I3.3,4ES9.2E1,2ES8.2E1,f7.1,f6.1)',k, &
- events(e)%tsc(k)%slip,&
- events(e)%tsc(k)%x,events(e)%tsc(k)%y,events(e)%tsc(k)%z, &
- events(e)%tsc(k)%length,events(e)%tsc(k)%width, &
- events(e)%tsc(k)%strike,events(e)%tsc(k)%dip
-
- IF (i .ne. k) THEN
- PRINT *, "error in input file: source index misfit"
- STOP 1
- END IF
- IF (events(e)%ts(k)%length .lt. minlength) THEN
- minlength=events(e)%ts(k)%length
- END IF
- IF (events(e)%ts(k)%width .lt. minwidth) THEN
- minwidth =events(e)%ts(k)%width
- END IF
-
- ! smooth out the slip distribution
- CALL antialiasingfilter(events(e)%ts(k)%slip, &
- events(e)%ts(k)%length,events(e)%ts(k)%width, &
- dx1,dx2,dx3,nyquist)
-
- ! comply to Wang's convention
- CALL wangconvention(events(e)%ts(k)%slip, &
- events(e)%ts(k)%x,events(e)%ts(k)%y,events(e)%ts(k)%z, &
- events(e)%ts(k)%length,events(e)%ts(k)%width, &
- events(e)%ts(k)%strike,events(e)%ts(k)%dip,dummy,rot)
-
- END DO
- PRINT 2000
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! M O G I S O U R C E S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of coseismic dilatation point sources"
- CALL getdata(iunit,dataline)
- READ (dataline,*) events(e)%nm
- PRINT '(I5)', events(e)%nm
- IF (events(e)%nm .GT. 0) THEN
- ALLOCATE(events(e)%m(events(e)%nm),events(e)%mc(events(e)%nm), &
- STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the tensile source list"
- PRINT 2000
- PRINT '(a)',"no. strain (positive for extension) xs ys zs"
- PRINT 2000
- DO k=1,events(e)%nm
- CALL getdata(iunit,dataline)
- READ (dataline,*) i,events(e)%m(k)%slip, &
- events(e)%m(k)%x,events(e)%m(k)%y,events(e)%m(k)%z
- ! copy the input format for display
- events(e)%mc(k)=events(e)%m(k)
-
- PRINT '(I3.3,4ES9.2E1)',k, &
- events(e)%mc(k)%slip,&
- events(e)%mc(k)%x,events(e)%mc(k)%y,events(e)%mc(k)%z
-
- IF (i .ne. k) THEN
- PRINT *, "error in input file: source index misfit"
- STOP 1
- END IF
-
- ! rotate the source in the computational reference frame
- CALL rotation(events(e)%m(k)%x,events(e)%m(k)%y,rot)
- END DO
- PRINT 2000
- END IF
-
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- ! S U R F A C E L O A D S
- ! - - - - - - - - - - - - - - - - - - - - - - - - - -
- PRINT '(a)', "number of surface loads"
- CALL getdata(iunit,dataline)
- READ (dataline,*) events(e)%nl
- PRINT '(I5)', events(e)%nl
- IF (events(e)%nl .GT. 0) THEN
- ALLOCATE(events(e)%l(events(e)%nl),events(e)%lc(events(e)%nl), &
- STAT=iostatus)
- IF (iostatus>0) STOP "could not allocate the load list"
- PRINT 2000
- PRINT '(a)',"no. xs ys t3 (force/surface/rigidity, positive down)"
- PRINT 2000
- DO k=1,events(e)%nl
- CALL getdata(iunit,dataline)
- READ (dataline,*) i, &
- events(e)%l(k)%x,events(e)%l(k)%y,events(e)%l(k)%slip
- ! copy the input format for display
- events(e)%lc(k)=events(e)%l(k)
-
- PRINT '(I3.3,4ES9.2E1)',k, &
- events(e)%lc(k)%x,events(e)%lc(k)%y,events(e)%lc(k)%slip
-
- IF (i .NE. k) THEN
- PRINT *, "error in input file: source index misfit"
- STOP 1
- END IF
-
- ! rotate the source in the computational reference frame
- CALL rotation(events(e)%l(k)%x,events(e)%l(k)%y,rot)
- END DO
- PRINT 2000
- END IF
-
- END DO
-
- ! test the presence of dislocations for coseismic calculation
- IF ((events(1)%nt .EQ. 0) .AND. &
- (events(1)%ns .EQ. 0) .AND. &
- (events(1)%nm .EQ. 0) .AND. &
- (events(1)%nl .EQ. 0) .AND. &
- (interval .LE. 0._8)) THEN
-
- WRITE_DEBUG_INFO
- WRITE (0,'("**** error **** ")')
- WRITE (0,'("no input dislocations or dilatation point sources")')
- WRITE (0,'("or surface tractions for first event . exiting.")')
- STOP 1
- END IF
-
- ! maximum recommended sampling size
- PRINT '(a,2ES8.2E1)', &
- "max sampling size (hor.,vert.):", minlength/2.5_8,minwidth/2.5_8
-
- PRINT 2000
-
-2000 FORMAT ("----------------------------------------------------------------------------")
-2100 FORMAT ("no. x1 x2 x3 length width strike dip")
-2200 FORMAT ("no. slip x1 x2 x3 length width strike dip rake")
-2300 FORMAT ("no. name x1 x2 x3 (name is a 4-character string)")
-2400 FORMAT ("no. strain x1 x2 x3 (positive for extension)")
-2500 FORMAT ("no. x1 x2 x3 length width strike dip rake")
-
- END SUBROUTINE init
-
!--------------------------------------------------------------------
! function IsOutput
! checks if output should be written based on user choices: if output
@@ -1898,93 +936,4 @@ 2500 FORMAT ("no. x1 x2
END SUBROUTINE integrationstep
- !------------------------------------------------------------------
- ! subroutine Rotation
- ! rotates a point coordinate into the computational reference
- ! system.
- !
- ! sylvain barbot (04/16/09) - original form
- !------------------------------------------------------------------
- SUBROUTINE rotation(x,y,rot)
- REAL*8, INTENT(INOUT) :: x,y
- REAL*8, INTENT(IN) :: rot
-
- REAL*8 :: alpha,xx,yy
-
- alpha=rot*DEG2RAD
- xx=x
- yy=y
-
- x=+xx*cos(alpha)+yy*sin(alpha)
- y=-xx*sin(alpha)+yy*cos(alpha)
-
- END SUBROUTINE rotation
-
- !-------------------------------------------------------------------
- ! subroutine AntiAliasingFilter
- ! smoothes a slip distribution model to avoid aliasing of
- ! the source geometry. Aliasing occurs is a slip patch has
- ! dimensions (width or length) smaller than the grid sampling.
- !
- ! if a patch length is smaller than a critical size L=dx*nyquist, it
- ! is increased to L and the slip (or opening) is scaled accordingly
- ! so that the moment M = s*L*W is conserved.
- !
- ! sylvain barbot (12/08/09) - original form
- !-------------------------------------------------------------------
- SUBROUTINE antialiasingfilter(slip,length,width,dx1,dx2,dx3,nyquist)
- REAL*8, INTENT(INOUT) :: slip,length,width
- REAL*8, INTENT(IN) :: dx1,dx2,dx3,nyquist
-
- REAL*8 :: dx
-
- ! minimum slip patch dimension
- dx=MIN(dx1,dx2,dx3)*nyquist
-
- ! update length
- IF (length .LT. dx) THEN
- slip=slip*length/dx
- length=dx
- END IF
- ! update width
- IF (width .LT. dx) THEN
- slip=slip*width/dx
- width=dx
- END IF
-
- END SUBROUTINE antialiasingfilter
-
- !------------------------------------------------------------------
- ! subroutine WangConvention
- ! converts a fault slip model from a geologic description including
- ! fault length, width, strike, dip and rake into a description
- ! compatible with internal convention of the program.
- !
- ! Internal convention describes a fault patch by the location of
- ! its center, instead of an upper corner and its orientation by
- ! the deviation from the vertical, instead of the angle from the
- ! horizontal and by the angle from the x2 axis (East-West)
- !------------------------------------------------------------------
- SUBROUTINE wangconvention(slip,x,y,z,length,width,strike,dip,rake,rot)
- REAL*8, INTENT(OUT) :: slip, x,y,z,strike,dip,rake
- REAL*8, INTENT(IN) :: length,width,rot
-
- slip=-slip
- strike=-90._8-strike
- dip = 90._8-dip
-
- strike=strike*DEG2RAD
- dip=dip*DEG2RAD
- rake=rake*DEG2RAD
-
- x=x-x0-length/2._8*sin(strike)+width /2._8*sin(dip)*cos(strike)
- y=y-y0-length/2._8*cos(strike)-width /2._8*sin(dip)*sin(strike)
- z=z+width /2._8*cos(dip)
-
- CALL rotation(x,y,rot)
-
- strike=strike+rot*DEG2RAD
-
- END SUBROUTINE wangconvention
-
END PROGRAM relax
diff -r 9992112ff05a -r 7596917b8776 types.f90
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/types.f90 Tue Apr 26 15:00:49 2011 -0700
@@ -0,0 +1,189 @@
+!-----------------------------------------------------------------------
+! Copyright 2007, 2008, 2009 Sylvain Barbot
+!
+! This file is part of RELAX
+!
+! RELAX is free software: you can redistribute it and/or modify
+! it under the terms of the GNU General Public License as published by
+! the Free Software Foundation, either version 3 of the License, or
+! (at your option) any later version.
+!
+! RELAX is distributed in the hope that it will be useful,
+! but WITHOUT ANY WARRANTY; without even the implied warranty of
+! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+! GNU General Public License for more details.
+!
+! You should have received a copy of the GNU General Public License
+! along with RELAX. If not, see <http://www.gnu.org/licenses/>.
+!-----------------------------------------------------------------------
+
+#include 'include.f90'
+
+MODULE types
+
+ TYPE SOURCE_STRUCT
+ SEQUENCE
+ REAL*8 :: slip,x,y,z,width,length,strike,dip,rake
+ END TYPE SOURCE_STRUCT
+
+ TYPE PLANE_STRUCT
+ SEQUENCE
+ REAL*8 :: x,y,z,width,length,strike,dip,rake
+ END TYPE PLANE_STRUCT
+
+ TYPE LAYER_STRUCT
+ SEQUENCE
+ REAL*8 :: z,gammadot0,stressexponent,cohesion,friction
+ END TYPE LAYER_STRUCT
+
+ TYPE WEAK_STRUCT
+ SEQUENCE
+ REAL*8 :: dgammadot0,x,y,z,width,length,thickness,strike,dip
+ END TYPE WEAK_STRUCT
+
+ TYPE VECTOR_STRUCT
+ SEQUENCE
+ REAL*8 :: v1,v2,v3
+ END TYPE VECTOR_STRUCT
+
+ TYPE TENSOR
+ SEQUENCE
+ REAL*4 :: s11,s12,s13,s22,s23,s33
+ END TYPE TENSOR
+
+ TYPE TENSOR_LAYER_STRUCT
+ SEQUENCE
+ REAL*4 :: z,dum
+ TYPE(TENSOR) :: t
+ END TYPE TENSOR_LAYER_STRUCT
+
+ TYPE SLIPPATCH_STRUCT
+ SEQUENCE
+ REAL*8 :: x1,x2,x3,lx,lz,slip,ss,ds
+ END TYPE SLIPPATCH_STRUCT
+
+ TYPE EVENT_STRUC
+ REAL*8 :: time
+ INTEGER*4 :: i,ns,nt,nm,nl
+ TYPE(SOURCE_STRUCT), DIMENSION(:), ALLOCATABLE :: s,sc,ts,tsc,m,mc,l,lc
+ END TYPE EVENT_STRUC
+
+ TYPE, PUBLIC :: SIMULATION_STRUC
+ ! grid dimension
+ INTEGER :: sx1,sx2,sx3
+
+ ! sampling
+ REAL*8 :: dx1,dx2,dx3
+
+ ! smoothing factor
+ REAL*8 :: beta
+
+ ! filter parameter for slip models
+ REAL*8 :: nyquist
+
+ ! center coordinates and rotation
+ REAL*8 :: x0,y0,rot
+
+#ifdef PROJ
+ ! geographic coordinates of center, UTM zone, length unit
+ REAL*8 :: lon0,lat0,umult
+ INTEGER :: zone
+#endif
+
+ ! observation depths
+ REAL*8 :: oz,ozs
+
+ ! output directory
+ CHARACTER(80) :: wdir
+
+ ! filenames
+ CHARACTER(80) :: reportfilename,reporttimefilename
+
+ ! elastic moduli and gravity parameter
+ REAL*8 :: lambda,mu,gam
+
+ ! time step parameters
+ REAL*8 :: interval
+ REAL*8 :: odt,tscale
+ INTEGER :: skip=0
+
+ ! number of observation planes
+ INTEGER :: nop
+
+ ! observation planes
+ TYPE(PLANE_STRUCT), DIMENSION(:), ALLOCATABLE :: op
+
+ ! number of observation points
+ INTEGER :: npts
+
+ ! observation points
+ TYPE(VECTOR_STRUCT), DIMENSION(:), ALLOCATABLE :: opts
+
+ ! observation points name
+ CHARACTER(LEN=4), DIMENSION(:), ALLOCATABLE :: ptsname
+
+ ! number of prestress interfaces
+ INTEGER :: nps
+
+ ! stress layers and stress structure
+ TYPE(TENSOR_LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: stresslayer,stressstruc
+
+ ! number of linear viscous interfaces
+ INTEGER :: nv
+
+ ! linear viscous layers and structure
+ TYPE(LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: linearlayer,linearstruc
+
+ ! number of linear weak zones
+ INTEGER :: nlwz
+
+ ! linear weak zones
+ TYPE(WEAK_STRUCT), DIMENSION(:), ALLOCATABLE :: linearweakzone,linearweakzonec
+
+ ! number of nonlinear viscous interfaces
+ INTEGER :: npl
+
+ ! nonlinear viscous layers and structure
+ TYPE(LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: nonlinearlayer,nonlinearstruc
+
+ ! number of nonlinear weak zones
+ INTEGER :: nnlwz
+
+ ! nonlinear viscous layers and structure
+ TYPE(WEAK_STRUCT), DIMENSION(:), ALLOCATABLE :: nonlinearweakzone,nonlinearweakzonec
+
+ ! number of fault creep interfaces
+ INTEGER :: nfc
+
+ ! fault creep interfaces
+ TYPE(LAYER_STRUCT), DIMENSION(:), ALLOCATABLE :: faultcreeplayer,faultcreepstruc
+
+ ! number of afterslip planes
+ INTEGER :: np
+
+ ! afterslip planes
+ TYPE(PLANE_STRUCT), DIMENSION(:), ALLOCATABLE :: n
+
+ ! interseismic event
+ TYPE(EVENT_STRUC) :: inter
+
+ ! number of coseismic events
+ INTEGER :: ne
+
+ ! coseismic events
+ TYPE(EVENT_STRUC), DIMENSION(:), ALLOCATABLE :: events
+
+ ! overrides output to formats
+ LOGICAL :: isoutputproj=.TRUE.
+ LOGICAL :: isoutputtxt=.TRUE.
+ LOGICAL :: isoutputvtk=.TRUE.
+ LOGICAL :: isoutputgrd=.TRUE.
+ LOGICAL :: isoutputxyz=.TRUE.
+
+ ! other options
+ LOGICAL :: isdryrun=.FALSE.
+ LOGICAL :: ishelp=.FALSE.
+
+ END TYPE SIMULATION_STRUC
+
+END MODULE types
More information about the CIG-COMMITS
mailing list