[cig-commits] r18216 - in seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D: . utm_example
carltape at geodynamics.org
carltape at geodynamics.org
Sat Apr 9 22:14:33 PDT 2011
Author: carltape
Date: 2011-04-09 22:14:33 -0700 (Sat, 09 Apr 2011)
New Revision: 18216
Added:
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/CMTSOLUTION_utm
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Par_file_utm
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/README
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/STATIONS_utm
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/interfaces_utm.dat
Modified:
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/process.sh
Log:
added two sub-examples to handle utm suppression for the socal1D internal mesh example
Modified: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/process.sh
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/process.sh 2011-04-10 00:33:12 UTC (rev 18215)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/process.sh 2011-04-10 05:14:33 UTC (rev 18216)
@@ -64,6 +64,7 @@
mv in_out_files/OUTPUT_FILES/output_mesher.txt in_out_files/OUTPUT_FILES/output_meshfem3D.txt
# stores setup
+cp -r in_data_files/meshfem3D_files/ in_out_files/OUTPUT_FILES/
cp in_data_files/Par_file in_out_files/OUTPUT_FILES/
cp in_data_files/CMTSOLUTION in_out_files/OUTPUT_FILES/
cp in_data_files/STATIONS in_out_files/OUTPUT_FILES/
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/CMTSOLUTION_utm
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/CMTSOLUTION_utm (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/CMTSOLUTION_utm 2011-04-10 05:14:33 UTC (rev 18216)
@@ -0,0 +1,13 @@
+PDE 2001 9 9 23 59 17.78 34.0745 -118.3792 6.4 4.2 4.2 HOLLYWOOD POLET SEM
+event name: 9703873
+time shift: 0.0000
+half duration: 0.3114
+latitude: 3771080.0
+longitude: 372738.0
+depth: 5.4000
+Mrr: -0.002000e+23
+Mtt: -0.064000e+23
+Mpp: 0.066000e+23
+Mrt: -0.090000e+23
+Mrp: -0.002000e+23
+Mtp: 0.188000e+23
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Mesh_Par_file_001_utm 2011-04-10 05:14:33 UTC (rev 18216)
@@ -0,0 +1,58 @@
+
+# coordinates of mesh block in latitude/longitude and depth in km
+LATITUDE_MIN = 3739440.d0
+LATITUDE_MAX = 3781380.d0
+LONGITUDE_MIN = 340080.d0
+LONGITUDE_MAX = 385720.d0
+DEPTH_BLOCK_KM = 60.d0
+UTM_PROJECTION_ZONE = 11
+SUPPRESS_UTM_PROJECTION = .true.
+
+# 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 = 24
+NEX_ETA = 24
+
+# number of MPI processors along xi and eta (can be different)
+NPROC_XI = 1
+NPROC_ETA = 1
+
+# Regular/irregular mesh
+USE_REGULAR_MESH = .false.
+# Only for irregular meshes, number of doubling layers (1 or 2) and their position
+NDOUBLINGS = 2
+# NZ_DOUGLING_1 is the parameter to set up if there is only one doubling layer
+NZ_DOUGLING_1 = 6
+NZ_DOUGLING_2 = 9
+
+# 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 = 4
+# 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 3000 7800 4500 0.0 0 2
+2 2800 6700 3870 0.0 0 2
+3 2670 6300 3640 0.0 0 2
+4 2400 5500 3180 0.0 0 2
+# number of regions
+NREGIONS = 4
+# 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 24 1 24 1 3 1
+1 24 1 24 4 6 2
+1 24 1 24 7 9 3
+1 24 1 24 10 12 4
+
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Par_file_utm
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Par_file_utm (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/Par_file_utm 2011-04-10 05:14:33 UTC (rev 18216)
@@ -0,0 +1,55 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
+NOISE_TOMOGRAPHY = 0 # 0 = earthquake simulation, 1/2/3 = three steps in noise simulation
+SAVE_FORWARD = .false.
+
+# UTM projection parameters
+UTM_PROJECTION_ZONE = 11
+SUPPRESS_UTM_PROJECTION = .true.
+
+# number of MPI processors
+NPROC = 1
+
+# time step parameters
+NSTEP = 840
+DT = 0.012
+
+# parameters describing the model
+OCEANS = .false.
+TOPOGRAPHY = .false.
+ATTENUATION = .false.
+USE_OLSEN_ATTENUATION = .false.
+ANISOTROPY = .false.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS = .true.
+
+# save AVS or OpenDX movies
+MOVIE_SURFACE = .false.
+MOVIE_VOLUME = .false.
+NTSTEP_BETWEEN_FRAMES = 200
+CREATE_SHAKEMAP = .false.
+SAVE_DISPLACEMENT = .false.
+USE_HIGHRES_FOR_MOVIES = .false.
+HDUR_MOVIE = 0.0
+
+# save AVS or OpenDX mesh files to check the mesh
+SAVE_MESH_FILES = .true.
+
+# path to store the local database file on each node
+LOCAL_PATH = ../in_out_files/DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 100
+
+# interval in time steps for writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
+
+# interval in time steps for reading adjoint traces
+NTSTEP_BETWEEN_READ_ADJSRC = 0 # 0 = read the whole adjoint sources at the same time
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .true.
+
+
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/README
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/README (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/README 2011-04-10 05:14:33 UTC (rev 18216)
@@ -0,0 +1,27 @@
+README
+(SPECFEM3D/examples/meshfem3D_examples/socal1D/utm/README)
+
+This example is mainly to show the ease of using either UTM or lon-lat coordinates.
+The utm example is better suited for direct comparison with external meshes for
+"flat models" that do not need to take into account the sphericity of Earth
+(for example, the 3D socal model used in Komatitsch et al. 2004).
+
+First, run the default socal1D example with 1 core.
+
+EXAMPLE A: Input geographic bounds in UTM, everything else in lon-lat.
+
+ > cp Mesh_Par_file_001_utm ../Mesh_Par_file
+ > cd ..
+ > ./process.sh
+
+EXAMPLE B: All input values in UTM.
+
+ > cp CMTSOLUTION_utm ../CMTSOLUTION
+ > cp STATIONS_utm ../STATIONS
+ > cp interfaces_utm.dat ../interfaces.dat
+ > cp Mesh_Par_file_001_utm ../Mesh_Par_file
+ > cp Par_file_utm ../Par_file
+ > cd ..
+ > ./process.sh
+
+=======================
\ No newline at end of file
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/STATIONS_utm
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/STATIONS_utm (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/STATIONS_utm 2011-04-10 05:14:33 UTC (rev 18216)
@@ -0,0 +1,30 @@
+ 24851 CE 3770540.194093 375794.778375 0.00 0.00
+ 24853 CE 3771245.142806 381341.534900 0.00 0.00
+ K851 CE 3770540.194093 375794.778375 0.00 0.00
+ K853 CE 3771245.142806 381341.534900 0.00 0.00
+ BVH CI 3771287.223592 371199.161741 0.00 0.00
+ CAB CI 3780444.660252 348729.637807 0.00 0.00
+ DJJ CI 3774690.021889 365788.275263 0.00 0.00
+ DJJB CI 3774690.021889 365788.275263 0.00 0.00
+ GR2 CI 3775840.420798 380163.280642 0.00 0.00
+ GRF CI 3775941.393569 380072.320783 0.00 0.00
+ LAF CI 3748222.391683 376849.747924 0.00 0.00
+ LAX CI 3756632.948247 369335.533371 0.00 0.00
+ LCG CI 3762853.919174 372743.220601 0.00 0.00
+ PDR CI 3758759.514975 367228.127192 0.00 0.00
+ SMF CI 3765299.026424 366421.222454 0.00 0.00
+ SMS CI 3764511.850697 365540.030956 0.00 0.00
+ SPF CI 3769762.226999 348077.027970 0.00 0.00
+ TCF CI 3772403.408916 352470.414102 0.00 0.00
+ USC CI 3764835.171833 381232.340025 0.00 0.00
+ WNS CI 3775572.000924 372724.507168 0.00 0.00
+ WTT CI 3756982.464007 383984.272469 0.00 0.00
+ FABS FA 3770467.255447 355014.474196 0.00 0.00
+ 5081 NP 3772430.358145 352378.540706 0.00 0.00
+ 5399 NP 3756538.331182 369602.304617 0.00 0.00
+ 5402 NP 3771266.130124 371281.939298 0.00 0.00
+ 5404 NP 3780363.811004 348721.882056 0.00 0.00
+ BVH NP 3771287.223592 371199.161741 0.00 0.00
+ GRF NP 3775941.393569 380072.320783 0.00 0.00
+ LAX NP 3756632.948247 369335.533371 0.00 0.00
+ TCF NP 3772403.408916 352470.414102 0.00 0.00
Added: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/interfaces_utm.dat
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/interfaces_utm.dat (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/utm_example/interfaces_utm.dat 2011-04-10 05:14:33 UTC (rev 18216)
@@ -0,0 +1,30 @@
+# number of interfaces
+ 4
+#
+# We describe each interface below, structured as a 2D-grid, with several parameters :
+# number of points along XI and ETA, minimal XI ETA coordinates
+# and spacing between points which must be constant.
+# Then the records contain the Z coordinates of the NXI x NETA points.
+#
+# interface number 1 (moho)
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
+ interface1.dat
+# interface number 2
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
+ interface2.dat
+# interface number 3
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
+ interface3.dat
+# interface number 4 (topography, top of the mesh)
+ .true. 2 2 122000.d0 3547000.d0 600000.d0 600000.d0
+ interface4.dat
+#
+# for each layer, we give the number of spectral elements in the vertical direction
+# layer number 1 (bottom layer)
+ 3
+# layer number 2
+ 3
+# layer number 4
+ 3
+# layer number 5 (top layer)
+ 3
\ No newline at end of file
More information about the CIG-COMMITS
mailing list