[cig-commits] [commit] devel, master: Added mantle heterogeneity model 'heterogen', where the user can specify an arbitrary density anomaly field. See for more information on implementation and use, DATA/heterogen/README_heterogen. (e48f803)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d_globe

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

>---------------------------------------------------------------

commit e48f8036436615e85d59afaad93575779025c247
Author: Laura Alisic <alisic at gps.caltech.edu>
Date:   Tue Jun 17 22:47:37 2008 +0000

    Added mantle heterogeneity model 'heterogen', where the user can specify an arbitrary density anomaly field. See for more information on implementation and use, DATA/heterogen/README_heterogen.


>---------------------------------------------------------------

e48f8036436615e85d59afaad93575779025c247
 Par_file                   | 20 ++++++++++----------
 heterogen/README_heterogen | 46 ++++++++++++++++++++++++++++++++++++++++++++++
 2 files changed, 56 insertions(+), 10 deletions(-)

diff --git a/Par_file b/Par_file
index 749ca14..3ec27ba 100644
--- a/Par_file
+++ b/Par_file
@@ -19,8 +19,8 @@ NEX_XI                          = 96
 NEX_ETA                         = 96
 
 # number of MPI processors along the two sides of the first chunk
-NPROC_XI                        = 6
-NPROC_ETA                       = 6
+NPROC_XI                        = 2
+NPROC_ETA                       = 2
 
 # 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
@@ -30,8 +30,8 @@ NPROC_ETA                       = 6
 #
 # fully 3D models:
 # transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994
-MODEL                           = 1D_isotropic_prem
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+MODEL                           = heterogen
 
 # parameters describing the Earth model
 OCEANS                          = .false.
@@ -45,7 +45,7 @@ ATTENUATION                     = .false.
 ABSORBING_CONDITIONS            = .false.
 
 # record length in minutes
-RECORD_LENGTH_IN_MINUTES        = 75.0d0
+RECORD_LENGTH_IN_MINUTES        = 10.0d0
 
 # save AVS or OpenDX movies
 #MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME) 
@@ -73,26 +73,26 @@ MOVIE_START                     = 0
 MOVIE_STOP                      = 40000
 
 # save mesh files to check the mesh
-SAVE_MESH_FILES                 = .false.
+SAVE_MESH_FILES                 = .true.
 
 # 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/
+LOCAL_PATH                      = /scratch/alisic/DATABASES_MPI
 
 # interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO      = 1000
+NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 
 # 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_ASCII_TEXT       = .false.
 OUTPUT_SEISMOS_SAC_ALPHANUM     = .false.
-OUTPUT_SEISMOS_SAC_BINARY       = .false.
+OUTPUT_SEISMOS_SAC_BINARY       = .true.
 
 # rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
 ROTATE_SEISMOGRAMS_RT           = .false.
diff --git a/heterogen/README_heterogen b/heterogen/README_heterogen
new file mode 100644
index 0000000..3f74803
--- /dev/null
+++ b/heterogen/README_heterogen
@@ -0,0 +1,46 @@
+----------------------------------------------------------------
+README for the heterogen model
+
+06/11/08
+Laura Alisic, Caltech
+alisic [at] gps [dot] caltech [dot] edu
+----------------------------------------------------------------
+
+
+The option 'heterogen' for the parameter MODEL in the Par_file
+allows the user to read in a density model created by the user. 
+It should be located in this directory, 
+SPECFEM3D_GLOBE/DATA/heterogen/. It uses the subroutine 
+'heterogen_mantle_model.f90'.
+
+The file must be named 'heterogen.dat', and be in
+direct access, formatted form with record length 20. The 
+number of points in radial, longitudinal and latitudinal
+direction has to be specified in 'constants.h.in' (N_R,
+N_THETA, N_PHI).
+NOTE: currently only N_R = N_THETA = N_PHI is supported.
+
+The subroutine then interpolates the value of the density
+anomaly field (which should be the fractional perturbation
+to the background model) onto the locations of the gll points.
+
+In the subroutine 'heterogen_mantle_model.f90', one can 
+specify the radial domain where the anomalous density model
+is applied, using 'r_inner' and 'r_outer' in meters.
+NOTE: currently, r_outer HAS to be < R_EARTH due to the used
+interpolation scheme. r_outer = R_EARTH-1.0 meter works fine.
+
+In 'read_compute_parameters.f90', one can modify the model that
+the heterogeneity model is applied to by changing the flags.
+
+For copying .bin files from the nodes to a local machine, use
+the script 'copy_m_globe_database_dvp.pl' instead of the normal
+script 'copy_m_globe_database.pl'. One can then visualize the
+applied density perturbations with 'xcombine_vol_data' and 
+Paraview. If the .bin output is unwanted, the user can modify
+'save_arrays_solver.f90'.
+
+
+Please email me with any problems or questions!
+
+----------------------------------------------------------------