[cig-commits] r11756 - in seismo/3D/SPECFEM3D_GLOBE/trunk: DATA UTILS
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sun Apr 6 13:27:32 PDT 2008
Author: dkomati1
Date: 2008-04-06 13:27:32 -0700 (Sun, 06 Apr 2008)
New Revision: 11756
Removed:
seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file_Japan_Tomo
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/extract_all_seismos_from_large_file.f90
Log:
declared a parameter (which has been moved to DATA/Par_file)
Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file_Japan_Tomo
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file_Japan_Tomo 2008-04-06 20:12:29 UTC (rev 11755)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file_Japan_Tomo 2008-04-06 20:27:32 UTC (rev 11756)
@@ -1,112 +0,0 @@
-
-# forward or adjoint simulation
-SIMULATION_TYPE = 1
-SAVE_FORWARD = .false. # save last frame of forward simulation or not
-
-# number of chunks (1,2,3 or 6)
-NCHUNKS = 1
-
-# angular width of the first chunk (not used if full sphere with six chunks)
-ANGULAR_WIDTH_XI_IN_DEGREES = 44.d0 # angular size of a chunk
-ANGULAR_WIDTH_ETA_IN_DEGREES = 33.d0
-CENTER_LATITUDE_IN_DEGREES = 34.5d0
-CENTER_LONGITUDE_IN_DEGREES = 130.0d0
-GAMMA_ROTATION_AZIMUTH = 0.d0
-
-# number of elements at the surface along the two sides of the first chunk
-# (must be multiple of 16 and 8 * multiple of NPROC below)
-NEX_XI = 384
-NEX_ETA = 288
-
-# number of MPI processors along the two sides of the first chunk
-NPROC_XI = 24
-NPROC_ETA = 18
-
-# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso
-#
-# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust, 1D_jp3d,1D_sea99
-#
-# fully 3D models:
-# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s29ea,sea99_jp3d1994,sea99,jp3d1994
-MODEL = s20rts
-
-# parameters describing the Earth model
-OCEANS = .false.
-ELLIPTICITY = .true.
-TOPOGRAPHY = .true.
-GRAVITY = .false.
-ROTATION = .false.
-ATTENUATION = .true.
-
-# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS = .true.
-
-# record length in minutes
-RECORD_LENGTH_IN_MINUTES = 35.0d0
-
-# save AVS or OpenDX movies
-MOVIE_SURFACE = .false.
-MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 100
-HDUR_MOVIE = 0.d0
-
-# save movie in volume. Will save element if center of element is in prescribed volume
-# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
-# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
-# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
-# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
-# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
-#MOVIE_VOLUME_COARSE saves movie only at corners of elements
-MOVIE_VOLUME_TYPE = 2
-MOVIE_VOLUME_COARSE = .true.
-MOVIE_TOP_KM = -100.0
-MOVIE_BOTTOM_KM = 1000.0
-MOVIE_WEST_DEG = -90.0
-MOVIE_EAST_DEG = 90.0
-MOVIE_NORTH_DEG = 90.0
-MOVIE_SOUTH_DEG = -90.0
-MOVIE_START = 0
-MOVIE_STOP = 40000
-
-# save mesh files to check the mesh
-SAVE_MESH_FILES = .false.
-
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
-NUMBER_OF_RUNS = 1
-NUMBER_OF_THIS_RUN = 1
-
-# path to store the local database files on each node
-LOCAL_PATH = /scratch/mchen/DATABASES_MPI.378137
-
-# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO = 1000
-
-# interval in time steps for temporary writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000
-NTSTEP_BETWEEN_READ_ADJSRC = 1000
-
-# output format for the seismograms (one can use either or all of the three formats)
-OUTPUT_SEISMOS_ASCII_TEXT = .false.
-OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
-OUTPUT_SEISMOS_SAC_BINARY = .true.
-
-# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
-ROTATE_SEISMOGRAMS_RT = .false.
-
-# decide if master process writes all the seismograms or if all processes do it in parallel
-WRITE_SEISMOGRAMS_BY_MASTER = .true.
-
-# save all seismograms in one large combined file instead of one file per seismogram
-# to avoid overloading shared non-local file systems such as GPFS for instance
-SAVE_ALL_SEISMOS_IN_ONE_FILE = .false.
-USE_BINARY_FOR_LARGE_FILE = .false.
-
-# flag to impose receivers at the surface or allow them to be buried
-RECEIVERS_CAN_BE_BURIED = .true.
-
-# print source time function
-PRINT_SOURCE_TIME_FUNCTION = .false.
-
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/extract_all_seismos_from_large_file.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/extract_all_seismos_from_large_file.f90 2008-04-06 20:12:29 UTC (rev 11755)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/extract_all_seismos_from_large_file.f90 2008-04-06 20:27:32 UTC (rev 11756)
@@ -34,12 +34,16 @@
include "../constants.h"
+! binary or ASCII storage format (see its value in DATA/Par_file instead)
+ logical, parameter :: USE_BINARY_FOR_LARGE_FILE = .true.
+
! number of seismogram files stored in the unique large file
- integer, parameter :: N_COMPONENTS = 1
- integer, parameter :: NREC = 451 * N_COMPONENTS
+ integer, parameter :: N_STATIONS = 647 ! see number of lines in DATA/STATIONS
+ integer, parameter :: N_COMPONENTS = 3
+ integer, parameter :: NREC = N_STATIONS * N_COMPONENTS
! number of time steps in each seismogram file
- integer, parameter :: NSTEP = 61600
+ integer, parameter :: NSTEP = 37200 ! see its value in OUTPUT_FILES/values_from_mesher.h
integer :: irec,istep,irepeat
real :: time,U_value
More information about the cig-commits
mailing list