[cig-commits] [commit] devel, master: bug fix for no gravity case and non-Deville routine (compute_forces_outer_core.f90) (692fd6a)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 6 08:04:13 PST 2014
Repository : https://github.com/geodynamics/specfem3d_globe
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f
>---------------------------------------------------------------
commit 692fd6aebff27f3ece6ae44e050545f81db99263
Author: Daniel Peter <peterda at ethz.ch>
Date: Tue Mar 23 21:28:16 2010 +0000
bug fix for no gravity case and non-Deville routine (compute_forces_outer_core.f90)
>---------------------------------------------------------------
692fd6aebff27f3ece6ae44e050545f81db99263
old/Par_file.default | 119 ---------------------------------------------------
old/Par_file.in | 119 ---------------------------------------------------
2 files changed, 238 deletions(-)
diff --git a/old/Par_file.default b/old/Par_file.default
deleted file mode 100644
index 22e165b..0000000
--- a/old/Par_file.default
+++ /dev/null
@@ -1,119 +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 = 6
-
-# angular width of the first chunk (not used if full sphere with six chunks)
-ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
-ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
-CENTER_LATITUDE_IN_DEGREES = 40.d0
-CENTER_LONGITUDE_IN_DEGREES = 10.d0
-GAMMA_ROTATION_AZIMUTH = 20.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 = 240
-NEX_ETA = 240
-
-# number of MPI processors along the two sides of the first chunk
-NPROC_XI = 5
-NPROC_ETA = 5
-
-# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
-#
-# 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
-#
-# fully 3D models:
-# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
-# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
-#
-# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
-# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
-# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
-MODEL = s362ani
-
-# parameters describing the Earth model
-OCEANS = .true.
-ELLIPTICITY = .true.
-TOPOGRAPHY = .true.
-GRAVITY = .true.
-ROTATION = .true.
-ATTENUATION = .true.
-
-# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS = .false.
-
-# record length in minutes
-RECORD_LENGTH_IN_MINUTES = 15.0d0
-
-# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
-#MOVIE_COARSE does not work with create_movie_AVS_DX
-MOVIE_SURFACE = .false.
-MOVIE_VOLUME = .false.
-MOVIE_COARSE = .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_TYPE = 2
-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 = ./DATABASES_MPI
-
-# 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 = .true.
-OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
-OUTPUT_SEISMOS_SAC_BINARY = .false.
-
-# 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.
-
diff --git a/old/Par_file.in b/old/Par_file.in
deleted file mode 100644
index 22e165b..0000000
--- a/old/Par_file.in
+++ /dev/null
@@ -1,119 +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 = 6
-
-# angular width of the first chunk (not used if full sphere with six chunks)
-ANGULAR_WIDTH_XI_IN_DEGREES = 90.d0 # angular size of a chunk
-ANGULAR_WIDTH_ETA_IN_DEGREES = 90.d0
-CENTER_LATITUDE_IN_DEGREES = 40.d0
-CENTER_LONGITUDE_IN_DEGREES = 10.d0
-GAMMA_ROTATION_AZIMUTH = 20.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 = 240
-NEX_ETA = 240
-
-# number of MPI processors along the two sides of the first chunk
-NPROC_XI = 5
-NPROC_ETA = 5
-
-# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
-#
-# 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
-#
-# fully 3D models:
-# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
-# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
-#
-# 3D models with 1D crust: append "_1Dcrust" the the 3D model name
-# to take the 1D crustal model from the
-# associated reference model rather than the default 3D crustal model
-# e.g. s20rts_1Dcrust, s362ani_1Dcrust, etc.
-MODEL = s362ani
-
-# parameters describing the Earth model
-OCEANS = .true.
-ELLIPTICITY = .true.
-TOPOGRAPHY = .true.
-GRAVITY = .true.
-ROTATION = .true.
-ATTENUATION = .true.
-
-# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS = .false.
-
-# record length in minutes
-RECORD_LENGTH_IN_MINUTES = 15.0d0
-
-# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME)
-#MOVIE_COARSE does not work with create_movie_AVS_DX
-MOVIE_SURFACE = .false.
-MOVIE_VOLUME = .false.
-MOVIE_COARSE = .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_TYPE = 2
-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 = ./DATABASES_MPI
-
-# 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 = .true.
-OUTPUT_SEISMOS_SAC_ALPHANUM = .false.
-OUTPUT_SEISMOS_SAC_BINARY = .false.
-
-# 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.
-
More information about the CIG-COMMITS
mailing list