# coordinates of mesh block in latitude/longitude and depth in km
LATITUDE_MIN                    = 34.0d0
LATITUDE_MAX                    = 38.0d0
LONGITUDE_MIN                   = -91.0d0
LONGITUDE_MAX                   = -87.0d0
DEPTH_BLOCK_KM                  = 40.d0
UTM_PROJECTION_ZONE             = 16
SUPPRESS_UTM_PROJECTION         = .false.

# file that contains the interfaces of the model / mesh
INTERFACES_FILE                 = interfaces.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                          = 128
NEX_ETA                         = 128

# number of MPI processors along xi and eta (can be different)
NPROC_XI                        = 2
NPROC_ETA                       = 2

# Regular/irregular mesh
USE_REGULAR_MESH                = .false.
# Only for irregular meshes, number of doubling layers (1 or 2) and their position
NDOUBLINGS                      = 1
# NZ_DOUGLING_1 is the parameter to set up if there is only one doubling layer
NZ_DOUGLING_1                   = 35
NZ_DOUGLING_2                   = 0

# create mesh files for visualisation or further checking
CREATE_ABAQUS_FILES             = .true.
CREATE_DX_FILES                 = .true.

# path to store the databases files
LOCAL_PATH                      = ../in_out_files/DATABASES_MPI

# number of materials
NMATERIALS                      = 6
# define the different materials in the model as :
# #material_id  #rho  #vp  #vs  #Q  #anisotropy_flag #domain_id
#     Q                : quality factor
#     anisotropy_flag  : 0=no anisotropy/ 1,2,.. check with implementation in aniso_model.f90
#     domain_id        : 1=acoustic / 2=elastic
1  3134  7300  4200  40.0  0  2
2  2925  6600  3800  50.0  0  2
3  2800  6170  3570  0     0  2
4  2543  4830  3200  50.0  0  2
5  2772  6020  3560  50.0  0  2
6  1774  1800  600   0     0  2
# number of regions
NREGIONS                        = 6
# 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              128           1              128             1          13       1
1              128           1              128             14         24       2
1              128           1              128             25         36       3
1              128           1              128             37         39       4
1              128           1              128             40         41       5
1              128           1              128             42         42       6