Merge pull request #1828 from nqdu/devel

CMT + point force simulation at the same time.

Author: danielpeter

DATA/CMTSOLUTION

@@ -0,0 +1,13 @@
+PDE  1999 01 01 00 00 00.00  67000 67000 -25000 4.2 4.2 homog_test
+event name:       homog_test
+time shift:       0.0000
+half duration:    5.0
+latorUTM:       67000.0
+longorUTM:      67000.0
+depth:          30.0
+Mrr:      -7.600000e+27
+Mtt:       7.700000e+27
+Mpp:      -2.000000e+26
+Mrt:      -2.500000e+28
+Mrp:       4.000000e+26
+Mtp:      -2.500000e+27

DATA/CMTSOLUTION

@@ -0,0 +1,4 @@
+X20 DB 67000.00 22732.14 0.0 0.0
+X30 DB 67000.00 34696.43 0.0 0.0
+X40 DB 67000.00 46660.71 0.0 0.0
+X50 DB 67000.00 58625.00 0.0 0.0

DATA/Par_file

@@ -0,0 +1,10 @@
+A0 XA 556.079455 813.961930 0. 0.
+A1 XA 1360.180039 -1437.744152 0. 0.
+A2 XA -1488.155201 400.944705 0. 0.
+A3 XA 36.576790 746.411648 0. 0.
+A4 XA 937.862885 -4.478963 0. 0.
+A5 XA 337.578200 -825.610063 0. 0.
+A6 XA 665.265952 -905.811406 0. 0.
+A7 XA -624.371795 781.592137 0. 0.
+A8 XA 1253.322368 -992.667490 0. 0.
+A9 XA 643.727350 -1234.980557 0. 0.

DATA/Par_file

@@ -0,0 +1,102 @@
+#-----------------------------------------------------------
+#
+# Meshing input parameters
+#
+#-----------------------------------------------------------
+
+# coordinates of mesh block in latitude/longitude and depth in km
+LATITUDE_MIN                    = -2500.0
+LATITUDE_MAX                    = 2500.0
+LONGITUDE_MIN                   = -2500.0
+LONGITUDE_MAX                   = 2500.0
+DEPTH_BLOCK_KM                  = 5.d0
+UTM_PROJECTION_ZONE             = 0
+SUPPRESS_UTM_PROJECTION         = .true.
+
+# file that contains the interfaces of the model / mesh
+INTERFACES_FILE                 = interfaces.dat
+
+# file that contains the cavity
+CAVITY_FILE                     = no_cavity.dat
+
+# number of elements at the surface along edges of the mesh at the surface
+# (must be 8 * multiple of NPROC below if mesh is not regular and contains mesh doublings)
+# (must be multiple of NPROC below if mesh is regular)
+NEX_XI                          = 32
+NEX_ETA                         = 32
+
+# number of MPI processors along xi and eta (can be different)
+NPROC_XI                        = 2
+NPROC_ETA                       = 2
+
+#-----------------------------------------------------------
+#
+# Doubling layers
+#
+#-----------------------------------------------------------
+
+# Regular/irregular mesh
+USE_REGULAR_MESH                = .true.
+# Only for irregular meshes, number of doubling layers and their position
+NDOUBLINGS                      = 0
+# NZ_DOUBLING_1 is the parameter to set up if there is only one doubling layer
+# (more doubling entries can be added if needed to match NDOUBLINGS value)
+NZ_DOUBLING_1                   = 0
+NZ_DOUBLING_2                   = 0
+
+#-----------------------------------------------------------
+#
+# Visualization
+#
+#-----------------------------------------------------------
+
+# create mesh files for visualisation or further checking
+CREATE_ABAQUS_FILES             = .false.
+CREATE_DX_FILES                 = .false.
+CREATE_VTK_FILES                = .true.
+
+# stores mesh files as cubit-exported files into directory MESH/ (for single process run)
+SAVE_MESH_AS_CUBIT              = .false.
+
+# path to store the databases files
+LOCAL_PATH                      = ./DATABASES_MPI
+
+#-----------------------------------------------------------
+#
+# CPML
+#
+#-----------------------------------------------------------
+
+# CPML perfectly matched absorbing layers
+THICKNESS_OF_X_PML              = 0.0d0
+THICKNESS_OF_Y_PML              = 0.0d0
+THICKNESS_OF_Z_PML              = 0.0d0
+
+#-----------------------------------------------------------
+#
+# Domain materials
+#
+#-----------------------------------------------------------
+
+# number of materials
+NMATERIALS                      = 1
+# define the different materials in the model as:
+# #material_id  #rho  #vp  #vs  #Q_Kappa  #Q_mu  #anisotropy_flag  #domain_id
+#     Q_Kappa          : Q_Kappa attenuation quality factor
+#     Q_mu             : Q_mu attenuation quality factor
+#     anisotropy_flag  : 0 = no anisotropy / 1,2,... check the implementation in file aniso_model.f90
+#     domain_id        : 1 = acoustic / 2 = elastic
+1  2600  5175  3025  9999. 9999.  0  2
+
+#-----------------------------------------------------------
+#
+# Domain regions
+#
+#-----------------------------------------------------------
+
+# number of regions
+NREGIONS                        = 1
+# define the different regions of the model as :
+#NEX_XI_BEGIN  #NEX_XI_END  #NEX_ETA_BEGIN  #NEX_ETA_END  #NZ_BEGIN #NZ_END  #material_id
+1              32           1              32           1         32      1
+