[cig-commits] r15634 - in seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta: . setup
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sat Aug 29 09:27:32 PDT 2009
Author: dkomati1
Date: 2009-08-29 09:27:32 -0700 (Sat, 29 Aug 2009)
New Revision: 15634
Added:
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/setup/constants_topo.h
Removed:
seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/constants_topo.h
Log:
moved include file to "setup" directory
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/constants_topo.h
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/constants_topo.h 2009-08-29 16:14:01 UTC (rev 15633)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/constants_topo.h 2009-08-29 16:27:32 UTC (rev 15634)
@@ -1,515 +0,0 @@
-!=====================================================================
-!
-! S p e c f e m 3 D G l o b e V e r s i o n 4 . 1
-! --------------------------------------------------
-!
-! Main authors: Dimitri Komatitsch and Jeroen Tromp
-! Seismological Laboratory, California Institute of Technology, USA
-! and University of Pau / CNRS / INRIA, France
-! (c) California Institute of Technology and University of Pau / CNRS / INRIA
-! August 2008
-!
-! This program 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 2 of the License, or
-! (at your option) any later version.
-!
-! This program 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 this program; if not, write to the Free Software Foundation, Inc.,
-! 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
-!
-!=====================================================================
-
-! constants.h. Generated from constants.h.in by configure.
-
-!
-!--- user can modify parameters below
-!
-
-! exaggeration factor to display topography in OpenDX
- double precision, parameter :: EXAGGERATION_FACTOR_TOPO = 22.d0 ! 35.d0
-
-! this for non blocking assembly
- logical, parameter :: USE_NONBLOCKING_COMMS = .true.
- integer, parameter :: ELEMENTS_NONBLOCKING_CM_IC = 1500
- integer, parameter :: ELEMENTS_NONBLOCKING_OC = 3000
-
- logical, parameter :: DEBUG_NONBLOCKING_COMMS = .false.
- logical, parameter :: DEBUG_USING_OPENDX = .false.
-
-!! DK DK temporary patch for the large Gordon Bell runs: set RECEIVERS_CAN_BE_BURIED
-!! DK DK to false in all cases etc
- logical, parameter :: PATCH_FOR_GORDON_BELL = .true.
-
-! save partial seismograms every 10,000 time steps or not
- logical, parameter :: SAVE_PARTIAL_SEISMOGRAMS = .false.
-
-! (much) faster detection of receivers at high resolution: use grid points only
- logical, parameter :: FASTER_RECEIVERS_POINTS_ONLY = .true.
-
-! suppress calculation and storage of seismograms if needed
- logical, parameter :: COMPUTE_STORE_SEISMOGRAMS = .true.
-
-!! DK DK for Gordon Bell
-! integer, parameter :: SEA99_VS_DIM1 = 100, SEA99_VS_DIM2 = 100, SEA99_VS_DIM3 = 100
- integer, parameter :: SEA99_VS_DIM1 = 1, SEA99_VS_DIM2 = 1, SEA99_VS_DIM3 = 1
-
-!! DK DK for Gordon Bell
-! integer, parameter :: AMM_V_DIM1 = 14, AMM_V_DIM2 = 34, AMM_V_DIM3 = 37, AMM_V_DIM4 = 73
-! use 6 for first index to avoid a warning when compiling
- integer, parameter :: AMM_V_DIM1 = 6, AMM_V_DIM2 = 1, AMM_V_DIM3 = 1, AMM_V_DIM4 = 1
-
-!
-! solver in single or double precision depending on the machine (4 or 8 bytes)
-!
-! ALSO CHANGE FILE precision.h ACCORDINGLY
-!
- integer, parameter :: SIZE_REAL = 4, SIZE_DOUBLE = 8
-
-! usually the size of integer and logical variables is the same as regular single-precision real variable
- integer, parameter :: SIZE_INTEGER = SIZE_REAL
- integer, parameter :: SIZE_LOGICAL = SIZE_REAL
-
-! set to SIZE_REAL to run in single precision
-! set to SIZE_DOUBLE to run in double precision (increases memory size by 2)
- integer, parameter :: CUSTOM_REAL = SIZE_REAL
-
-! input, output and main MPI I/O files
- integer, parameter :: ISTANDARD_OUTPUT = 6
- integer, parameter :: IIN = 40,IOUT = 41,IOUT_SAC = 903
-! local file unit for output of buffers
- integer, parameter :: IOUT_BUFFERS = 35
-! uncomment this to write messages to a text file
- integer, parameter :: IMAIN = 42
-! uncomment this to write messages to the screen (slows down the code)
-! integer, parameter :: IMAIN = ISTANDARD_OUTPUT
-
-! number of values read in Par_file that we need to broadcast
- integer, parameter :: NVALUES_bcast_integer = 41
- integer, parameter :: NVALUES_bcast_double_precision = 30
- integer, parameter :: NVALUES_bcast_logical = 33
-
-! R_EARTH is the radius of the bottom of the oceans (radius of Earth in m)
- double precision, parameter :: R_EARTH = 6371000.d0
-! uncomment line below for PREM with oceans
-! double precision, parameter :: R_EARTH = 6368000.d0
-
-! average density in the full Earth to normalize equation
- double precision, parameter :: RHOAV = 5514.3d0
-
-! for topography/bathymetry model
-
-!!--- ETOPO5 5-minute model, smoothed Harvard version
-!! size of topography and bathymetry file
-! integer, parameter :: NX_BATHY = 4320,NY_BATHY = 2160
-!! resolution of topography file in minutes
-! integer, parameter :: RESOLUTION_TOPO_FILE = 5
-!! pathname of the topography file
-! character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo5_smoothed_Harvard.dat'
-
-!--- ETOPO4 4-minute model created by subsampling and smoothing etopo-2
-! size of topography and bathymetry file
- integer, parameter :: NX_BATHY = 5400,NY_BATHY = 2700
-!!!!!!!!!! integer, parameter :: NX_BATHY = 1,NY_BATHY = 1
-! resolution of topography file in minutes
- integer, parameter :: RESOLUTION_TOPO_FILE = 4
-! pathname of the topography file
-! character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo4_smoothed_window_7.dat'
- character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo4_smoothed_window_11_minmax_3500.dat'
-
-!!--- ETOPO2 2-minute model, not implemented yet
-!! size of topography and bathymetry file
-! integer, parameter :: NX_BATHY = 10800,NY_BATHY = 5400
-!! resolution of topography file in minutes
-! integer, parameter :: RESOLUTION_TOPO_FILE = 2
-!! pathname of the topography file
-! character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo2_smoothed_window7.dat'
-
-! minimum thickness in meters to include the effect of the oceans and topo
- double precision, parameter :: MINIMUM_THICKNESS_3D_OCEANS = 100.d0
-
-! number of GLL points in each direction of an element (degree plus one)
- integer, parameter :: NGLLX = 5
- integer, parameter :: NGLLY = NGLLX
- integer, parameter :: NGLLZ = NGLLX
-
-! flag to exclude elements that are too far from target in source detection
- logical, parameter :: USE_DISTANCE_CRITERION = .true.
-
-! flag to display detailed information about location of stations
- logical, parameter :: DISPLAY_DETAILS_STATIONS = .false.
-
-! maximum length of station and network name for receivers
- integer, parameter :: MAX_LENGTH_STATION_NAME = 32
- integer, parameter :: MAX_LENGTH_NETWORK_NAME = 8
-
-! we mimic a triangle of half duration equal to half_duration_triangle
-! using a Gaussian having a very close shape, as explained in Figure 4.2
-! of the manual. This source decay rate to mimic an equivalent triangle
-! was found by trial and error
- double precision, parameter :: SOURCE_DECAY_MIMIC_TRIANGLE = 1.628d0
-
-! maximum number of sources to locate simultaneously (do not set it too high otherwise you will use a lot of memory)
- integer, parameter :: NSOURCES_SUBSET_MAX = 100
-
-! distance threshold (in km) above which we consider that a receiver
-! is located outside the mesh and therefore excluded from the station list
- double precision, parameter :: THRESHOLD_EXCLUDE_STATION = 50.d0
-
-! the first doubling is implemented right below the Moho
-! it seems optimal to implement the three other doublings at these depths
-! in the mantle
- double precision, parameter :: DEPTH_SECOND_DOUBLING_OPTIMAL = 1650000.d0
-! in the outer core
- double precision, parameter :: DEPTH_THIRD_DOUBLING_OPTIMAL = 3860000.d0
-! in the outer core
- double precision, parameter :: DEPTH_FOURTH_DOUBLING_OPTIMAL = 5000000.d0
-
-! Boundary Mesh -- save Moho, 400, 670 km discontinuity topology files (in
-! the mesher) and use them for the computation of boundary kernel (in the solver)
- logical, parameter :: SAVE_BOUNDARY_MESH = .false.
-
-! this parameter must be set to .true. to compute anisotropic kernels
-! in crust and mantle (related to the 21 Cij in geographical coordinates)
-! default is .false. to compute isotropic kernels (related to alpha and beta)
- logical, parameter :: ANISOTROPIC_KL = .false.
-
-! print date and time estimate of end of run in another country,
-! in addition to local time.
-! For instance: the code runs at Caltech in California but the person
-! running the code is connected remotely from France, which has 9 hours more.
-! The time difference with that remote location can be positive or negative
- logical, parameter :: ADD_TIME_ESTIMATE_ELSEWHERE = .false.
- integer, parameter :: HOURS_TIME_DIFFERENCE = +9
- integer, parameter :: MINUTES_TIME_DIFFERENCE = +0
-
-!
-!--- debugging flags
-!
-
-! flags to actually assemble with MPI or not
-! and to actually match fluid and solid regions of the Earth or not
-! should always be set to true except when debugging code
- logical, parameter :: ACTUALLY_ASSEMBLE_MPI_SLICES = .true.
- logical, parameter :: ACTUALLY_ASSEMBLE_MPI_CHUNKS = .true.
- logical, parameter :: ACTUALLY_COUPLE_FLUID_CMB = .true.
- logical, parameter :: ACTUALLY_COUPLE_FLUID_ICB = .true.
-
-! flag to turn off the conversion of geographic to geocentric coordinates for
-! the seismic source and the stations; i.e. assume a perfect sphere, which
-! can be useful for benchmarks of a spherical Earth with fictitious sources and stations
- logical, parameter :: ASSUME_PERFECT_SPHERE = .false.
-
-! flag to only create the mesh but not start the solver (for instance to check the mesh obtained)
- logical, parameter :: MESHER_ONLY = .false.
-
-! flag to put a fictitious source in each region in the case of a serial test
- logical, parameter :: PUT_SOURCE_IN_EACH_REGION = .false.
-
-!------------------------------------------------------
-!----------- do not modify anything below -------------
-!------------------------------------------------------
-
-! on some processors (e.g. Pentiums) it is necessary to suppress underflows
-! by using a small initial field instead of zero
- logical, parameter :: FIX_UNDERFLOW_PROBLEM = .true.
-
-! some useful constants
- double precision, parameter :: PI = 3.141592653589793d0
- double precision, parameter :: TWO_PI = 2.d0 * PI
- double precision, parameter :: PI_OVER_FOUR = PI / 4.d0
-
-! to convert angles from degrees to radians
- double precision, parameter :: DEGREES_TO_RADIANS = PI / 180.d0
-
-! 3-D simulation
- integer, parameter :: NDIM = 3
-
-! dimension of the boundaries of the slices
- integer, parameter :: NDIM2D = 2
-
-! number of nodes for 2D and 3D shape functions for hexahedra with 27 nodes
- integer, parameter :: NGNOD = 27, NGNOD2D = 9
-
-! gravitational constant
- double precision, parameter :: GRAV = 6.6723d-11
-
-! a few useful constants
- double precision, parameter :: ZERO = 0.d0,ONE = 1.d0,TWO = 2.d0,HALF = 0.5d0
-
- real(kind=CUSTOM_REAL), parameter :: &
- ONE_THIRD = 1._CUSTOM_REAL/3._CUSTOM_REAL, &
- TWO_THIRDS = 2._CUSTOM_REAL/3._CUSTOM_REAL, &
- FOUR_THIRDS = 4._CUSTOM_REAL/3._CUSTOM_REAL
-
-! very large and very small values
- double precision, parameter :: HUGEVAL = 1.d+30,TINYVAL = 1.d-9
-
-! very large real value declared independently of the machine
- real(kind=CUSTOM_REAL), parameter :: HUGEVAL_SNGL = 1.e+30_CUSTOM_REAL
-
-! very large integer value
- integer, parameter :: HUGEINT = 100000000
-
-! normalized radius of free surface
- double precision, parameter :: R_UNIT_SPHERE = ONE
-
-! same radius in km
- double precision, parameter :: R_EARTH_KM = R_EARTH / 1000.d0
-
-! fixed thickness of 3 km for PREM oceans
- double precision, parameter :: THICKNESS_OCEANS_PREM = 3000.d0 / R_EARTH
-
-! shortest radius at which crust is implemented (80 km depth)
-! to be constistent with the D80 discontinuity, we impose the crust only above it
- double precision, parameter :: R_DEEPEST_CRUST = (R_EARTH - 80000.d0) / R_EARTH
-
-! maximum number of chunks (full sphere)
- integer, parameter :: NCHUNKS_MAX = 6
-
-! define block type based upon chunk number (between 1 and 6)
-! do not change this numbering, chunk AB must be number 1 for central cube
- integer, parameter :: CHUNK_AB = 1
- integer, parameter :: CHUNK_AC = 2
- integer, parameter :: CHUNK_BC = 3
- integer, parameter :: CHUNK_AC_ANTIPODE = 4
- integer, parameter :: CHUNK_BC_ANTIPODE = 5
- integer, parameter :: CHUNK_AB_ANTIPODE = 6
-
-! maximum number of regions in the mesh
- integer, parameter :: MAX_NUM_REGIONS = 3
-
-! define flag for regions of the global Earth mesh
- integer, parameter :: IREGION_CRUST_MANTLE = 1
- integer, parameter :: IREGION_OUTER_CORE = 2
- integer, parameter :: IREGION_INNER_CORE = 3
-
-! define flag for elements
- integer, parameter :: IFLAG_CRUST = 1
-
- integer, parameter :: IFLAG_80_MOHO = 2
- integer, parameter :: IFLAG_220_80 = 3
- integer, parameter :: IFLAG_670_220 = 4
- integer, parameter :: IFLAG_MANTLE_NORMAL = 5
-
- integer, parameter :: IFLAG_OUTER_CORE_NORMAL = 6
-
- integer, parameter :: IFLAG_INNER_CORE_NORMAL = 7
- integer, parameter :: IFLAG_MIDDLE_CENTRAL_CUBE = 8
- integer, parameter :: IFLAG_BOTTOM_CENTRAL_CUBE = 9
- integer, parameter :: IFLAG_TOP_CENTRAL_CUBE = 10
- integer, parameter :: IFLAG_IN_FICTITIOUS_CUBE = 11
-
- integer, parameter :: NSPEC2D_XI_SUPERBRICK = 8
- integer, parameter :: NSPEC2D_ETA_SUPERBRICK = 8
- integer, parameter :: NSPEC2D_XI_SUPERBRICK_1L = 6
- integer, parameter :: NSPEC2D_ETA_SUPERBRICK_1L = 6
-
-! dummy flag used for mesh display purposes only
- integer, parameter :: IFLAG_DUMMY = 100
-
-! max number of layers that are used in the radial direction to build the full mesh
- integer, parameter :: MAX_NUMBER_OF_MESH_LAYERS = 15
-
-! define number of spectral elements and points in basic symmetric mesh doubling superbrick
- integer, parameter :: NSPEC_DOUBLING_SUPERBRICK = 32
- integer, parameter :: NGLOB_DOUBLING_SUPERBRICK = 67
- integer, parameter :: NSPEC_SUPERBRICK_1L = 28
- integer, parameter :: NGLOB_SUPERBRICK_1L = 58
- integer, parameter :: NGNOD_EIGHT_CORNERS = 8
-
-! define flag for reference 1D Earth model
- integer, parameter :: REFERENCE_MODEL_PREM = 1
- integer, parameter :: REFERENCE_MODEL_IASP91 = 2
- integer, parameter :: REFERENCE_MODEL_1066A = 3
- integer, parameter :: REFERENCE_MODEL_AK135 = 4
- integer, parameter :: REFERENCE_MODEL_REF = 5
- integer, parameter :: REFERENCE_MODEL_JP1D = 6
- integer, parameter :: REFERENCE_MODEL_SEA1D = 7
-
-! define flag for 3D Earth model
- integer, parameter :: THREE_D_MODEL_S20RTS = 1
- integer, parameter :: THREE_D_MODEL_S362ANI = 2
- integer, parameter :: THREE_D_MODEL_S362WMANI = 3
- integer, parameter :: THREE_D_MODEL_S362ANI_PREM = 4
- integer, parameter :: THREE_D_MODEL_S29EA = 5
- integer, parameter :: THREE_D_MODEL_SEA99_JP3D = 6
- integer, parameter :: THREE_D_MODEL_SEA99 = 7
- integer, parameter :: THREE_D_MODEL_JP3D = 8
-
-! define flag for regions of the global Earth for attenuation
- integer, parameter :: NUM_REGIONS_ATTENUATION = 5
-
- integer, parameter :: IREGION_ATTENUATION_INNER_CORE = 1
- integer, parameter :: IREGION_ATTENUATION_CMB_670 = 2
- integer, parameter :: IREGION_ATTENUATION_670_220 = 3
- integer, parameter :: IREGION_ATTENUATION_220_80 = 4
- integer, parameter :: IREGION_ATTENUATION_80_SURFACE = 5
- integer, parameter :: IREGION_ATTENUATION_UNDEFINED = 6
-
-! number of standard linear solids for attenuation
- integer, parameter :: N_SLS = 3
-
-! computation of standard linear solids in meshfem3D
-! ATTENUATION_COMP_RESOLUTION: Number of Digits after decimal
-! ATTENUATION_COMP_MAXIMUM: Maximum Q Value
- integer, parameter :: ATTENUATION_COMP_RESOLUTION = 1
- integer, parameter :: ATTENUATION_COMP_MAXIMUM = 5000
-
-! for lookup table for attenuation every 100 m in radial direction of Earth model
- integer, parameter :: NRAD_ATTENUATION = 70000
- double precision, parameter :: TABLE_ATTENUATION = R_EARTH_KM * 10.0d0
-
-!!!!!!!!!!!!!! for determination of the attenuation period range
-!!!!!!!!!!!!!! if this is set to .true. then the hardcoded values will be used
-!!!!!!!!!!!!!! otherwise they are computed automatically from the number of elements
-!!!!!!!!!!!!!! This *may* be a useful parameter for benchmarking against older versions
-!!!!!!!!!!!!! logical, parameter :: ATTENUATION_RANGE_PREDEFINED = .false.
-
-! flag for the four edges of each slice and for the bottom edge
- integer, parameter :: XI_MIN = 1
- integer, parameter :: XI_MAX = 2
- integer, parameter :: ETA_MIN = 3
- integer, parameter :: ETA_MAX = 4
- integer, parameter :: BOTTOM = 5
-
-! flags to select the right corner in each slice
- integer, parameter :: ILOWERLOWER = 1
- integer, parameter :: ILOWERUPPER = 2
- integer, parameter :: IUPPERLOWER = 3
- integer, parameter :: IUPPERUPPER = 4
-
-! number of points in each AVS or OpenDX quadrangular cell for movies
- integer, parameter :: NGNOD2D_AVS_DX = 4
-
-! number of faces a given slice can share with other slices
-! this is at most 2, except when there is only once slice per chunk
-! in which case it is 4
- integer, parameter :: NUMFACES_SHARED = 2 !!!!! DK DK removed support for one slice only 4
-
-! number of corners a given slice can share with other slices
-! this is at most 1, except when there is only once slice per chunk
-! in which case it is 4
- integer, parameter :: NUMCORNERS_SHARED = 1 !!!!!! DK DK removed support for one slice only 4
-
-! number of layers in PREM
- integer, parameter :: NR = 640
-
-! smallest real number on many machines = 1.1754944E-38
-! largest real number on many machines = 3.4028235E+38
-! small negligible initial value to avoid very slow underflow trapping
-! but not too small to avoid trapping on velocity and acceleration in Newmark
- real(kind=CUSTOM_REAL), parameter :: VERYSMALLVAL = 1.E-24_CUSTOM_REAL
-
-! displacement threshold above which we consider that the code became unstable
- real(kind=CUSTOM_REAL), parameter :: STABILITY_THRESHOLD = 1.E+25_CUSTOM_REAL
-
-! geometrical tolerance for boundary detection
- double precision, parameter :: SMALLVAL = 0.00001d0
-
-! small tolerance for conversion from x y z to r theta phi
- double precision, parameter :: SMALL_VAL_ANGLE = 1.d-10
-
-! geometry tolerance parameter to calculate number of independent grid points
-! sensitive to actual size of model, assumes reference sphere of radius 1
-! this is an absolute value for normalized coordinates in the Earth
- double precision, parameter :: SMALLVALTOL = 1.d-10
-
-! do not use tags for MPI messages, use dummy tag instead
- integer, parameter :: itag = 0,itag2 = 0
-
-! for the Gauss-Lobatto-Legendre points and weights
- double precision, parameter :: GAUSSALPHA = 0.d0,GAUSSBETA = 0.d0
-
-! number of lines per source in CMTSOLUTION file
- integer, parameter :: NLINES_PER_CMTSOLUTION_SOURCE = 13
-
-! number of iterations to solve the non linear system for xi and eta
- integer, parameter :: NUM_ITER = 4
-
-! number of hours per day for rotation rate of the Earth
- double precision, parameter :: HOURS_PER_DAY = 24.d0
-
-! for lookup table for gravity every 100 m in radial direction of Earth model
- integer, parameter :: NRAD_GRAVITY = 70000
-
-! number of layers in DATA/1066a/1066a.dat
- integer, parameter :: NR_1066A = 160
-
-! number of layers in DATA/ak135/ak135.dat
- integer, parameter :: NR_AK135 = 144
-
-! number of layers in DATA/s362ani/REF
- integer, parameter :: NR_REF = 750
-
-! number of layers in DATA/Lebedev_sea99 1D model
- integer, parameter :: NR_SEA1D = 163
-
-! three_d_mantle_model_constants
- integer, parameter :: NK = 20,NS = 20,ND = 1
-
-! Japan 3D model (Zhao, 1994) constants
- integer, parameter :: MPA=42,MRA=48,MHA=21,MPB=42,MRB=48,MHB=18
- integer, parameter :: MKA=2101,MKB=2101
-
-! The meaningful range of Zhao et al.'s model (1994) is as follows:
-! latitude : 32 - 45 N
-! longitude: 130-145 E
-! depth : 0 - 500 km
-! The deepest Moho beneath Japan is 40 km
- double precision,parameter :: LAT_MAX = 45.d0
- double precision,parameter :: LAT_MIN = 32.d0
- double precision,parameter :: LON_MAX = 145.d0
- double precision,parameter :: LON_MIN = 130.d0
- double precision,parameter :: DEP_MAX = 500.d0
-
-! use sedimentary layers of crust 2.0
- logical, parameter :: INCLUDE_SEDIMENTS_CRUST = .true.
-
-! number of points per degree in smoothed crust2.0
- integer, parameter :: NFACTOR_CRUST = 2
- integer, parameter :: NLON_CRUST = 180 * NFACTOR_CRUST
- integer, parameter :: NLAT_CRUST = 90 * NFACTOR_CRUST
-
-! to inflate the central cube (set to 0.d0 for a non-inflated cube)
- double precision, parameter :: CENTRAL_CUBE_INFLATE_FACTOR = 0.41d0
-
-! for the stretching of crustal elements in the case of 3D models
- double precision, parameter :: MAX_RATIO_CRUST_STRETCHING = 0.6d0
-
-! to suppress the crustal layers (replaced by an extension of the mantle: R_EARTH is not modified, but no more crustal doubling)
- logical, parameter :: SUPPRESS_CRUSTAL_MESH = .false.
-
-! to add a fourth doubling at the bottom of the outer core
- logical, parameter :: ADD_4TH_DOUBLING = .false.
-
-! parameters to cut the doubling brick
-
-! this to cut the superbrick: 3 possibilities, 4 cases max / possibility
-! three possibilities: (cut in xi and eta) or (cut in xi) or (cut in eta)
-! case 1: (ximin and etamin) or ximin or etamin
-! case 2: (ximin and etamax) or ximax or etamax
-! case 3: ximax and etamin
-! case 4: ximax and etamax
- integer, parameter :: NB_CUT_CASE = 4
-
-! corner 1: ximin and etamin
-! corner 2: ximax and etamin
-! corner 3: ximax and etamax
-! corner 4: ximin and etamax
- integer, parameter :: NB_SQUARE_CORNERS = 4
-
-! two possibilities: xi or eta
-! face 1: ximin or etamin
-! face 2: ximax or etamax
- integer, parameter :: NB_SQUARE_EDGES_ONEDIR = 2
-
-! this for the geometry of the basic doubling brick
- integer, parameter :: NSPEC_DOUBLING_BASICBRICK = 8
- integer, parameter :: NGLOB_DOUBLING_BASICBRICK = 27
-
Copied: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/setup/constants_topo.h (from rev 15632, seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/constants_topo.h)
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/setup/constants_topo.h (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/setup/constants_topo.h 2009-08-29 16:27:32 UTC (rev 15634)
@@ -0,0 +1,515 @@
+!=====================================================================
+!
+! S p e c f e m 3 D G l o b e V e r s i o n 4 . 1
+! --------------------------------------------------
+!
+! Main authors: Dimitri Komatitsch and Jeroen Tromp
+! Seismological Laboratory, California Institute of Technology, USA
+! and University of Pau / CNRS / INRIA, France
+! (c) California Institute of Technology and University of Pau / CNRS / INRIA
+! August 2008
+!
+! This program 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 2 of the License, or
+! (at your option) any later version.
+!
+! This program 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 this program; if not, write to the Free Software Foundation, Inc.,
+! 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
+!
+!=====================================================================
+
+! constants.h. Generated from constants.h.in by configure.
+
+!
+!--- user can modify parameters below
+!
+
+! exaggeration factor to display topography in OpenDX
+ double precision, parameter :: EXAGGERATION_FACTOR_TOPO = 22.d0 ! 35.d0
+
+! this for non blocking assembly
+ logical, parameter :: USE_NONBLOCKING_COMMS = .true.
+ integer, parameter :: ELEMENTS_NONBLOCKING_CM_IC = 1500
+ integer, parameter :: ELEMENTS_NONBLOCKING_OC = 3000
+
+ logical, parameter :: DEBUG_NONBLOCKING_COMMS = .false.
+ logical, parameter :: DEBUG_USING_OPENDX = .false.
+
+!! DK DK temporary patch for the large Gordon Bell runs: set RECEIVERS_CAN_BE_BURIED
+!! DK DK to false in all cases etc
+ logical, parameter :: PATCH_FOR_GORDON_BELL = .true.
+
+! save partial seismograms every 10,000 time steps or not
+ logical, parameter :: SAVE_PARTIAL_SEISMOGRAMS = .false.
+
+! (much) faster detection of receivers at high resolution: use grid points only
+ logical, parameter :: FASTER_RECEIVERS_POINTS_ONLY = .true.
+
+! suppress calculation and storage of seismograms if needed
+ logical, parameter :: COMPUTE_STORE_SEISMOGRAMS = .true.
+
+!! DK DK for Gordon Bell
+! integer, parameter :: SEA99_VS_DIM1 = 100, SEA99_VS_DIM2 = 100, SEA99_VS_DIM3 = 100
+ integer, parameter :: SEA99_VS_DIM1 = 1, SEA99_VS_DIM2 = 1, SEA99_VS_DIM3 = 1
+
+!! DK DK for Gordon Bell
+! integer, parameter :: AMM_V_DIM1 = 14, AMM_V_DIM2 = 34, AMM_V_DIM3 = 37, AMM_V_DIM4 = 73
+! use 6 for first index to avoid a warning when compiling
+ integer, parameter :: AMM_V_DIM1 = 6, AMM_V_DIM2 = 1, AMM_V_DIM3 = 1, AMM_V_DIM4 = 1
+
+!
+! solver in single or double precision depending on the machine (4 or 8 bytes)
+!
+! ALSO CHANGE FILE precision.h ACCORDINGLY
+!
+ integer, parameter :: SIZE_REAL = 4, SIZE_DOUBLE = 8
+
+! usually the size of integer and logical variables is the same as regular single-precision real variable
+ integer, parameter :: SIZE_INTEGER = SIZE_REAL
+ integer, parameter :: SIZE_LOGICAL = SIZE_REAL
+
+! set to SIZE_REAL to run in single precision
+! set to SIZE_DOUBLE to run in double precision (increases memory size by 2)
+ integer, parameter :: CUSTOM_REAL = SIZE_REAL
+
+! input, output and main MPI I/O files
+ integer, parameter :: ISTANDARD_OUTPUT = 6
+ integer, parameter :: IIN = 40,IOUT = 41,IOUT_SAC = 903
+! local file unit for output of buffers
+ integer, parameter :: IOUT_BUFFERS = 35
+! uncomment this to write messages to a text file
+ integer, parameter :: IMAIN = 42
+! uncomment this to write messages to the screen (slows down the code)
+! integer, parameter :: IMAIN = ISTANDARD_OUTPUT
+
+! number of values read in Par_file that we need to broadcast
+ integer, parameter :: NVALUES_bcast_integer = 41
+ integer, parameter :: NVALUES_bcast_double_precision = 30
+ integer, parameter :: NVALUES_bcast_logical = 33
+
+! R_EARTH is the radius of the bottom of the oceans (radius of Earth in m)
+ double precision, parameter :: R_EARTH = 6371000.d0
+! uncomment line below for PREM with oceans
+! double precision, parameter :: R_EARTH = 6368000.d0
+
+! average density in the full Earth to normalize equation
+ double precision, parameter :: RHOAV = 5514.3d0
+
+! for topography/bathymetry model
+
+!!--- ETOPO5 5-minute model, smoothed Harvard version
+!! size of topography and bathymetry file
+! integer, parameter :: NX_BATHY = 4320,NY_BATHY = 2160
+!! resolution of topography file in minutes
+! integer, parameter :: RESOLUTION_TOPO_FILE = 5
+!! pathname of the topography file
+! character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo5_smoothed_Harvard.dat'
+
+!--- ETOPO4 4-minute model created by subsampling and smoothing etopo-2
+! size of topography and bathymetry file
+ integer, parameter :: NX_BATHY = 5400,NY_BATHY = 2700
+!!!!!!!!!! integer, parameter :: NX_BATHY = 1,NY_BATHY = 1
+! resolution of topography file in minutes
+ integer, parameter :: RESOLUTION_TOPO_FILE = 4
+! pathname of the topography file
+! character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo4_smoothed_window_7.dat'
+ character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo4_smoothed_window_11_minmax_3500.dat'
+
+!!--- ETOPO2 2-minute model, not implemented yet
+!! size of topography and bathymetry file
+! integer, parameter :: NX_BATHY = 10800,NY_BATHY = 5400
+!! resolution of topography file in minutes
+! integer, parameter :: RESOLUTION_TOPO_FILE = 2
+!! pathname of the topography file
+! character(len=*), parameter :: PATHNAME_TOPO_FILE = 'DATA/topo_bathy/topo_bathy_etopo2_smoothed_window7.dat'
+
+! minimum thickness in meters to include the effect of the oceans and topo
+ double precision, parameter :: MINIMUM_THICKNESS_3D_OCEANS = 100.d0
+
+! number of GLL points in each direction of an element (degree plus one)
+ integer, parameter :: NGLLX = 5
+ integer, parameter :: NGLLY = NGLLX
+ integer, parameter :: NGLLZ = NGLLX
+
+! flag to exclude elements that are too far from target in source detection
+ logical, parameter :: USE_DISTANCE_CRITERION = .true.
+
+! flag to display detailed information about location of stations
+ logical, parameter :: DISPLAY_DETAILS_STATIONS = .false.
+
+! maximum length of station and network name for receivers
+ integer, parameter :: MAX_LENGTH_STATION_NAME = 32
+ integer, parameter :: MAX_LENGTH_NETWORK_NAME = 8
+
+! we mimic a triangle of half duration equal to half_duration_triangle
+! using a Gaussian having a very close shape, as explained in Figure 4.2
+! of the manual. This source decay rate to mimic an equivalent triangle
+! was found by trial and error
+ double precision, parameter :: SOURCE_DECAY_MIMIC_TRIANGLE = 1.628d0
+
+! maximum number of sources to locate simultaneously (do not set it too high otherwise you will use a lot of memory)
+ integer, parameter :: NSOURCES_SUBSET_MAX = 100
+
+! distance threshold (in km) above which we consider that a receiver
+! is located outside the mesh and therefore excluded from the station list
+ double precision, parameter :: THRESHOLD_EXCLUDE_STATION = 50.d0
+
+! the first doubling is implemented right below the Moho
+! it seems optimal to implement the three other doublings at these depths
+! in the mantle
+ double precision, parameter :: DEPTH_SECOND_DOUBLING_OPTIMAL = 1650000.d0
+! in the outer core
+ double precision, parameter :: DEPTH_THIRD_DOUBLING_OPTIMAL = 3860000.d0
+! in the outer core
+ double precision, parameter :: DEPTH_FOURTH_DOUBLING_OPTIMAL = 5000000.d0
+
+! Boundary Mesh -- save Moho, 400, 670 km discontinuity topology files (in
+! the mesher) and use them for the computation of boundary kernel (in the solver)
+ logical, parameter :: SAVE_BOUNDARY_MESH = .false.
+
+! this parameter must be set to .true. to compute anisotropic kernels
+! in crust and mantle (related to the 21 Cij in geographical coordinates)
+! default is .false. to compute isotropic kernels (related to alpha and beta)
+ logical, parameter :: ANISOTROPIC_KL = .false.
+
+! print date and time estimate of end of run in another country,
+! in addition to local time.
+! For instance: the code runs at Caltech in California but the person
+! running the code is connected remotely from France, which has 9 hours more.
+! The time difference with that remote location can be positive or negative
+ logical, parameter :: ADD_TIME_ESTIMATE_ELSEWHERE = .false.
+ integer, parameter :: HOURS_TIME_DIFFERENCE = +9
+ integer, parameter :: MINUTES_TIME_DIFFERENCE = +0
+
+!
+!--- debugging flags
+!
+
+! flags to actually assemble with MPI or not
+! and to actually match fluid and solid regions of the Earth or not
+! should always be set to true except when debugging code
+ logical, parameter :: ACTUALLY_ASSEMBLE_MPI_SLICES = .true.
+ logical, parameter :: ACTUALLY_ASSEMBLE_MPI_CHUNKS = .true.
+ logical, parameter :: ACTUALLY_COUPLE_FLUID_CMB = .true.
+ logical, parameter :: ACTUALLY_COUPLE_FLUID_ICB = .true.
+
+! flag to turn off the conversion of geographic to geocentric coordinates for
+! the seismic source and the stations; i.e. assume a perfect sphere, which
+! can be useful for benchmarks of a spherical Earth with fictitious sources and stations
+ logical, parameter :: ASSUME_PERFECT_SPHERE = .false.
+
+! flag to only create the mesh but not start the solver (for instance to check the mesh obtained)
+ logical, parameter :: MESHER_ONLY = .false.
+
+! flag to put a fictitious source in each region in the case of a serial test
+ logical, parameter :: PUT_SOURCE_IN_EACH_REGION = .false.
+
+!------------------------------------------------------
+!----------- do not modify anything below -------------
+!------------------------------------------------------
+
+! on some processors (e.g. Pentiums) it is necessary to suppress underflows
+! by using a small initial field instead of zero
+ logical, parameter :: FIX_UNDERFLOW_PROBLEM = .true.
+
+! some useful constants
+ double precision, parameter :: PI = 3.141592653589793d0
+ double precision, parameter :: TWO_PI = 2.d0 * PI
+ double precision, parameter :: PI_OVER_FOUR = PI / 4.d0
+
+! to convert angles from degrees to radians
+ double precision, parameter :: DEGREES_TO_RADIANS = PI / 180.d0
+
+! 3-D simulation
+ integer, parameter :: NDIM = 3
+
+! dimension of the boundaries of the slices
+ integer, parameter :: NDIM2D = 2
+
+! number of nodes for 2D and 3D shape functions for hexahedra with 27 nodes
+ integer, parameter :: NGNOD = 27, NGNOD2D = 9
+
+! gravitational constant
+ double precision, parameter :: GRAV = 6.6723d-11
+
+! a few useful constants
+ double precision, parameter :: ZERO = 0.d0,ONE = 1.d0,TWO = 2.d0,HALF = 0.5d0
+
+ real(kind=CUSTOM_REAL), parameter :: &
+ ONE_THIRD = 1._CUSTOM_REAL/3._CUSTOM_REAL, &
+ TWO_THIRDS = 2._CUSTOM_REAL/3._CUSTOM_REAL, &
+ FOUR_THIRDS = 4._CUSTOM_REAL/3._CUSTOM_REAL
+
+! very large and very small values
+ double precision, parameter :: HUGEVAL = 1.d+30,TINYVAL = 1.d-9
+
+! very large real value declared independently of the machine
+ real(kind=CUSTOM_REAL), parameter :: HUGEVAL_SNGL = 1.e+30_CUSTOM_REAL
+
+! very large integer value
+ integer, parameter :: HUGEINT = 100000000
+
+! normalized radius of free surface
+ double precision, parameter :: R_UNIT_SPHERE = ONE
+
+! same radius in km
+ double precision, parameter :: R_EARTH_KM = R_EARTH / 1000.d0
+
+! fixed thickness of 3 km for PREM oceans
+ double precision, parameter :: THICKNESS_OCEANS_PREM = 3000.d0 / R_EARTH
+
+! shortest radius at which crust is implemented (80 km depth)
+! to be constistent with the D80 discontinuity, we impose the crust only above it
+ double precision, parameter :: R_DEEPEST_CRUST = (R_EARTH - 80000.d0) / R_EARTH
+
+! maximum number of chunks (full sphere)
+ integer, parameter :: NCHUNKS_MAX = 6
+
+! define block type based upon chunk number (between 1 and 6)
+! do not change this numbering, chunk AB must be number 1 for central cube
+ integer, parameter :: CHUNK_AB = 1
+ integer, parameter :: CHUNK_AC = 2
+ integer, parameter :: CHUNK_BC = 3
+ integer, parameter :: CHUNK_AC_ANTIPODE = 4
+ integer, parameter :: CHUNK_BC_ANTIPODE = 5
+ integer, parameter :: CHUNK_AB_ANTIPODE = 6
+
+! maximum number of regions in the mesh
+ integer, parameter :: MAX_NUM_REGIONS = 3
+
+! define flag for regions of the global Earth mesh
+ integer, parameter :: IREGION_CRUST_MANTLE = 1
+ integer, parameter :: IREGION_OUTER_CORE = 2
+ integer, parameter :: IREGION_INNER_CORE = 3
+
+! define flag for elements
+ integer, parameter :: IFLAG_CRUST = 1
+
+ integer, parameter :: IFLAG_80_MOHO = 2
+ integer, parameter :: IFLAG_220_80 = 3
+ integer, parameter :: IFLAG_670_220 = 4
+ integer, parameter :: IFLAG_MANTLE_NORMAL = 5
+
+ integer, parameter :: IFLAG_OUTER_CORE_NORMAL = 6
+
+ integer, parameter :: IFLAG_INNER_CORE_NORMAL = 7
+ integer, parameter :: IFLAG_MIDDLE_CENTRAL_CUBE = 8
+ integer, parameter :: IFLAG_BOTTOM_CENTRAL_CUBE = 9
+ integer, parameter :: IFLAG_TOP_CENTRAL_CUBE = 10
+ integer, parameter :: IFLAG_IN_FICTITIOUS_CUBE = 11
+
+ integer, parameter :: NSPEC2D_XI_SUPERBRICK = 8
+ integer, parameter :: NSPEC2D_ETA_SUPERBRICK = 8
+ integer, parameter :: NSPEC2D_XI_SUPERBRICK_1L = 6
+ integer, parameter :: NSPEC2D_ETA_SUPERBRICK_1L = 6
+
+! dummy flag used for mesh display purposes only
+ integer, parameter :: IFLAG_DUMMY = 100
+
+! max number of layers that are used in the radial direction to build the full mesh
+ integer, parameter :: MAX_NUMBER_OF_MESH_LAYERS = 15
+
+! define number of spectral elements and points in basic symmetric mesh doubling superbrick
+ integer, parameter :: NSPEC_DOUBLING_SUPERBRICK = 32
+ integer, parameter :: NGLOB_DOUBLING_SUPERBRICK = 67
+ integer, parameter :: NSPEC_SUPERBRICK_1L = 28
+ integer, parameter :: NGLOB_SUPERBRICK_1L = 58
+ integer, parameter :: NGNOD_EIGHT_CORNERS = 8
+
+! define flag for reference 1D Earth model
+ integer, parameter :: REFERENCE_MODEL_PREM = 1
+ integer, parameter :: REFERENCE_MODEL_IASP91 = 2
+ integer, parameter :: REFERENCE_MODEL_1066A = 3
+ integer, parameter :: REFERENCE_MODEL_AK135 = 4
+ integer, parameter :: REFERENCE_MODEL_REF = 5
+ integer, parameter :: REFERENCE_MODEL_JP1D = 6
+ integer, parameter :: REFERENCE_MODEL_SEA1D = 7
+
+! define flag for 3D Earth model
+ integer, parameter :: THREE_D_MODEL_S20RTS = 1
+ integer, parameter :: THREE_D_MODEL_S362ANI = 2
+ integer, parameter :: THREE_D_MODEL_S362WMANI = 3
+ integer, parameter :: THREE_D_MODEL_S362ANI_PREM = 4
+ integer, parameter :: THREE_D_MODEL_S29EA = 5
+ integer, parameter :: THREE_D_MODEL_SEA99_JP3D = 6
+ integer, parameter :: THREE_D_MODEL_SEA99 = 7
+ integer, parameter :: THREE_D_MODEL_JP3D = 8
+
+! define flag for regions of the global Earth for attenuation
+ integer, parameter :: NUM_REGIONS_ATTENUATION = 5
+
+ integer, parameter :: IREGION_ATTENUATION_INNER_CORE = 1
+ integer, parameter :: IREGION_ATTENUATION_CMB_670 = 2
+ integer, parameter :: IREGION_ATTENUATION_670_220 = 3
+ integer, parameter :: IREGION_ATTENUATION_220_80 = 4
+ integer, parameter :: IREGION_ATTENUATION_80_SURFACE = 5
+ integer, parameter :: IREGION_ATTENUATION_UNDEFINED = 6
+
+! number of standard linear solids for attenuation
+ integer, parameter :: N_SLS = 3
+
+! computation of standard linear solids in meshfem3D
+! ATTENUATION_COMP_RESOLUTION: Number of Digits after decimal
+! ATTENUATION_COMP_MAXIMUM: Maximum Q Value
+ integer, parameter :: ATTENUATION_COMP_RESOLUTION = 1
+ integer, parameter :: ATTENUATION_COMP_MAXIMUM = 5000
+
+! for lookup table for attenuation every 100 m in radial direction of Earth model
+ integer, parameter :: NRAD_ATTENUATION = 70000
+ double precision, parameter :: TABLE_ATTENUATION = R_EARTH_KM * 10.0d0
+
+!!!!!!!!!!!!!! for determination of the attenuation period range
+!!!!!!!!!!!!!! if this is set to .true. then the hardcoded values will be used
+!!!!!!!!!!!!!! otherwise they are computed automatically from the number of elements
+!!!!!!!!!!!!!! This *may* be a useful parameter for benchmarking against older versions
+!!!!!!!!!!!!! logical, parameter :: ATTENUATION_RANGE_PREDEFINED = .false.
+
+! flag for the four edges of each slice and for the bottom edge
+ integer, parameter :: XI_MIN = 1
+ integer, parameter :: XI_MAX = 2
+ integer, parameter :: ETA_MIN = 3
+ integer, parameter :: ETA_MAX = 4
+ integer, parameter :: BOTTOM = 5
+
+! flags to select the right corner in each slice
+ integer, parameter :: ILOWERLOWER = 1
+ integer, parameter :: ILOWERUPPER = 2
+ integer, parameter :: IUPPERLOWER = 3
+ integer, parameter :: IUPPERUPPER = 4
+
+! number of points in each AVS or OpenDX quadrangular cell for movies
+ integer, parameter :: NGNOD2D_AVS_DX = 4
+
+! number of faces a given slice can share with other slices
+! this is at most 2, except when there is only once slice per chunk
+! in which case it is 4
+ integer, parameter :: NUMFACES_SHARED = 2 !!!!! DK DK removed support for one slice only 4
+
+! number of corners a given slice can share with other slices
+! this is at most 1, except when there is only once slice per chunk
+! in which case it is 4
+ integer, parameter :: NUMCORNERS_SHARED = 1 !!!!!! DK DK removed support for one slice only 4
+
+! number of layers in PREM
+ integer, parameter :: NR = 640
+
+! smallest real number on many machines = 1.1754944E-38
+! largest real number on many machines = 3.4028235E+38
+! small negligible initial value to avoid very slow underflow trapping
+! but not too small to avoid trapping on velocity and acceleration in Newmark
+ real(kind=CUSTOM_REAL), parameter :: VERYSMALLVAL = 1.E-24_CUSTOM_REAL
+
+! displacement threshold above which we consider that the code became unstable
+ real(kind=CUSTOM_REAL), parameter :: STABILITY_THRESHOLD = 1.E+25_CUSTOM_REAL
+
+! geometrical tolerance for boundary detection
+ double precision, parameter :: SMALLVAL = 0.00001d0
+
+! small tolerance for conversion from x y z to r theta phi
+ double precision, parameter :: SMALL_VAL_ANGLE = 1.d-10
+
+! geometry tolerance parameter to calculate number of independent grid points
+! sensitive to actual size of model, assumes reference sphere of radius 1
+! this is an absolute value for normalized coordinates in the Earth
+ double precision, parameter :: SMALLVALTOL = 1.d-10
+
+! do not use tags for MPI messages, use dummy tag instead
+ integer, parameter :: itag = 0,itag2 = 0
+
+! for the Gauss-Lobatto-Legendre points and weights
+ double precision, parameter :: GAUSSALPHA = 0.d0,GAUSSBETA = 0.d0
+
+! number of lines per source in CMTSOLUTION file
+ integer, parameter :: NLINES_PER_CMTSOLUTION_SOURCE = 13
+
+! number of iterations to solve the non linear system for xi and eta
+ integer, parameter :: NUM_ITER = 4
+
+! number of hours per day for rotation rate of the Earth
+ double precision, parameter :: HOURS_PER_DAY = 24.d0
+
+! for lookup table for gravity every 100 m in radial direction of Earth model
+ integer, parameter :: NRAD_GRAVITY = 70000
+
+! number of layers in DATA/1066a/1066a.dat
+ integer, parameter :: NR_1066A = 160
+
+! number of layers in DATA/ak135/ak135.dat
+ integer, parameter :: NR_AK135 = 144
+
+! number of layers in DATA/s362ani/REF
+ integer, parameter :: NR_REF = 750
+
+! number of layers in DATA/Lebedev_sea99 1D model
+ integer, parameter :: NR_SEA1D = 163
+
+! three_d_mantle_model_constants
+ integer, parameter :: NK = 20,NS = 20,ND = 1
+
+! Japan 3D model (Zhao, 1994) constants
+ integer, parameter :: MPA=42,MRA=48,MHA=21,MPB=42,MRB=48,MHB=18
+ integer, parameter :: MKA=2101,MKB=2101
+
+! The meaningful range of Zhao et al.'s model (1994) is as follows:
+! latitude : 32 - 45 N
+! longitude: 130-145 E
+! depth : 0 - 500 km
+! The deepest Moho beneath Japan is 40 km
+ double precision,parameter :: LAT_MAX = 45.d0
+ double precision,parameter :: LAT_MIN = 32.d0
+ double precision,parameter :: LON_MAX = 145.d0
+ double precision,parameter :: LON_MIN = 130.d0
+ double precision,parameter :: DEP_MAX = 500.d0
+
+! use sedimentary layers of crust 2.0
+ logical, parameter :: INCLUDE_SEDIMENTS_CRUST = .true.
+
+! number of points per degree in smoothed crust2.0
+ integer, parameter :: NFACTOR_CRUST = 2
+ integer, parameter :: NLON_CRUST = 180 * NFACTOR_CRUST
+ integer, parameter :: NLAT_CRUST = 90 * NFACTOR_CRUST
+
+! to inflate the central cube (set to 0.d0 for a non-inflated cube)
+ double precision, parameter :: CENTRAL_CUBE_INFLATE_FACTOR = 0.41d0
+
+! for the stretching of crustal elements in the case of 3D models
+ double precision, parameter :: MAX_RATIO_CRUST_STRETCHING = 0.6d0
+
+! to suppress the crustal layers (replaced by an extension of the mantle: R_EARTH is not modified, but no more crustal doubling)
+ logical, parameter :: SUPPRESS_CRUSTAL_MESH = .false.
+
+! to add a fourth doubling at the bottom of the outer core
+ logical, parameter :: ADD_4TH_DOUBLING = .false.
+
+! parameters to cut the doubling brick
+
+! this to cut the superbrick: 3 possibilities, 4 cases max / possibility
+! three possibilities: (cut in xi and eta) or (cut in xi) or (cut in eta)
+! case 1: (ximin and etamin) or ximin or etamin
+! case 2: (ximin and etamax) or ximax or etamax
+! case 3: ximax and etamin
+! case 4: ximax and etamax
+ integer, parameter :: NB_CUT_CASE = 4
+
+! corner 1: ximin and etamin
+! corner 2: ximax and etamin
+! corner 3: ximax and etamax
+! corner 4: ximin and etamax
+ integer, parameter :: NB_SQUARE_CORNERS = 4
+
+! two possibilities: xi or eta
+! face 1: ximin or etamin
+! face 2: ximax or etamax
+ integer, parameter :: NB_SQUARE_EDGES_ONEDIR = 2
+
+! this for the geometry of the basic doubling brick
+ integer, parameter :: NSPEC_DOUBLING_BASICBRICK = 8
+ integer, parameter :: NGLOB_DOUBLING_BASICBRICK = 27
+
More information about the CIG-COMMITS
mailing list