|
| 1 | +#----------------------------------------------------------- |
| 2 | +# |
| 3 | +# Meshing input parameters |
| 4 | +# |
| 5 | +#----------------------------------------------------------- |
| 6 | + |
| 7 | +# coordinates of mesh block in latitude/longitude and depth in km |
| 8 | +LATITUDE_MIN = 40.7 |
| 9 | +LATITUDE_MAX = 41.0 |
| 10 | +LONGITUDE_MIN = 13.88 |
| 11 | +LONGITUDE_MAX = 14.48 |
| 12 | +DEPTH_BLOCK_KM = 40.0 |
| 13 | +UTM_PROJECTION_ZONE = 33 |
| 14 | +SUPPRESS_UTM_PROJECTION = .false. |
| 15 | + |
| 16 | +# file that contains the interfaces of the model / mesh |
| 17 | +INTERFACES_FILE = interfaces.dat |
| 18 | + |
| 19 | +# file that contains the cavity |
| 20 | +CAVITY_FILE = no_cavity.dat |
| 21 | + |
| 22 | +# number of elements at the surface along edges of the mesh at the surface |
| 23 | +# (must be 8 * multiple of NPROC below if mesh is not regular and contains mesh doublings) |
| 24 | +# (must be multiple of NPROC below if mesh is regular) |
| 25 | +NEX_XI = 256 |
| 26 | +NEX_ETA = 256 |
| 27 | + |
| 28 | +# number of MPI processors along xi and eta (can be different) |
| 29 | +NPROC_XI = 4 |
| 30 | +NPROC_ETA = 4 |
| 31 | + |
| 32 | +#----------------------------------------------------------- |
| 33 | +# |
| 34 | +# Doubling layers |
| 35 | +# |
| 36 | +#----------------------------------------------------------- |
| 37 | + |
| 38 | +# Regular/irregular mesh |
| 39 | +USE_REGULAR_MESH = .false. |
| 40 | +# Only for irregular meshes, number of doubling layers and their position |
| 41 | +NDOUBLINGS = 2 |
| 42 | +# NZ_DOUBLING_1 is the parameter to set up if there is only one doubling layer |
| 43 | +# (more doubling entries can be added if needed to match NDOUBLINGS value) |
| 44 | +NZ_DOUBLING_1 = 17 |
| 45 | +NZ_DOUBLING_1 = 25 |
| 46 | + |
| 47 | +#----------------------------------------------------------- |
| 48 | +# |
| 49 | +# Visualization |
| 50 | +# |
| 51 | +#----------------------------------------------------------- |
| 52 | + |
| 53 | +# create mesh files for visualisation or further checking |
| 54 | +CREATE_ABAQUS_FILES = .false. |
| 55 | +CREATE_DX_FILES = .false. |
| 56 | +CREATE_VTK_FILES = .true. |
| 57 | + |
| 58 | +# stores mesh files as cubit-exported files into directory MESH/ (for single process run) |
| 59 | +SAVE_MESH_AS_CUBIT = .false. |
| 60 | + |
| 61 | +# path to store the databases files |
| 62 | +LOCAL_PATH = ./DATABASES_MPI |
| 63 | + |
| 64 | +#----------------------------------------------------------- |
| 65 | +# |
| 66 | +# CPML |
| 67 | +# |
| 68 | +#----------------------------------------------------------- |
| 69 | + |
| 70 | +# CPML perfectly matched absorbing layers |
| 71 | +THICKNESS_OF_X_PML = 0.02 # with 0.005 deg ~ 550 m |
| 72 | +THICKNESS_OF_Y_PML = 0.02 |
| 73 | +THICKNESS_OF_Z_PML = 0.02 |
| 74 | + |
| 75 | +# add PML layers as extra outer mesh layers |
| 76 | +ADD_PML_AS_EXTRA_MESH_LAYERS = .true. |
| 77 | +NUMBER_OF_PML_LAYERS_TO_ADD = 4 |
| 78 | + |
| 79 | +#----------------------------------------------------------- |
| 80 | +# |
| 81 | +# Domain materials |
| 82 | +# |
| 83 | +#----------------------------------------------------------- |
| 84 | + |
| 85 | +## number of materials |
| 86 | +NMATERIALS = 8 |
| 87 | +## define the different materials in the model as: |
| 88 | +## #material_id #rho #vp #vs #Q_Kappa #Q_mu #anisotropy_flag #domain_id |
| 89 | +## Q_Kappa : Q_Kappa attenuation quality factor |
| 90 | +## Q_mu : Q_mu attenuation quality factor |
| 91 | +## anisotropy_flag : 0 = no anisotropy / 1,2,... check the implementation in file aniso_model.f90 |
| 92 | +## domain_id : 1 = acoustic / 2 = elastic / 3 = poroelastic |
| 93 | +1 1378 1900 840 9999. 9999.0 0 2 |
| 94 | +2 1378 1900 1000 9999. 9999.0 0 2 |
| 95 | +3 1666 2800 1500 9999. 9999.0 0 2 |
| 96 | +4 1954 3700 2000 9999. 9999.0 0 2 |
| 97 | +5 2130 4250 2350 9999. 9999.0 0 2 |
| 98 | +6 2370 5000 2840 9999. 9999.0 0 2 |
| 99 | +7 2693 6010 3400 9999. 9999.0 0 2 |
| 100 | +8 2850 6500 3400 9999. 9999.0 0 2 |
| 101 | + |
| 102 | +#----------------------------------------------------------- |
| 103 | +# |
| 104 | +# Domain regions |
| 105 | +# |
| 106 | +#----------------------------------------------------------- |
| 107 | + |
| 108 | +# number of regions |
| 109 | +NREGIONS = 8 |
| 110 | +# define the different regions of the model as : |
| 111 | +#NEX_XI_BEGIN #NEX_XI_END #NEX_ETA_BEGIN #NEX_ETA_END #NZ_BEGIN #NZ_END #material_id |
| 112 | +1 -1 1 -1 1 17 8 |
| 113 | +1 -1 1 -1 18 19 7 |
| 114 | +1 -1 1 -1 20 21 6 |
| 115 | +1 -1 1 -1 22 23 5 |
| 116 | +1 -1 1 -1 24 25 4 |
| 117 | +1 -1 1 -1 26 30 3 |
| 118 | +1 -1 1 -1 31 34 2 |
| 119 | +1 -1 1 -1 35 35 1 |
| 120 | + |
0 commit comments