local drex model paths; added drex_model case

bvanderbeek · Jul 1, 2024 · 36babdb · 36babdb
1 parent f68e481
commit 36babdb
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 2 deletions.
diff --git a/src/meshfem3D/model_aniso_mantle_drex.f90 b/src/meshfem3D/model_aniso_mantle_drex.f90
@@ -623,8 +623,8 @@ subroutine read_aniso_mantle_model_drex()
   integer :: ier
   double precision :: Mcur,Gcur,Lcur
   double precision, DIMENSION(21) :: XE 
-  character(len=*), parameter :: drex_model = '~/geodynamic_models/plume/drex_model.xyz'
-  character(len=*), parameter :: iso_prem = '~/specfem3d_globe/DATA/DREX/iso_prem.dat'
+  character(len=*), parameter :: drex_model = 'DATA/DREX/drex_model.xyz'
+  character(len=*), parameter :: iso_prem = 'DATA/DREX/iso_prem.dat'
 
 ! read the model iso_prem model
 

diff --git a/src/shared/get_model_parameters.F90 b/src/shared/get_model_parameters.F90
@@ -466,6 +466,17 @@ subroutine get_model_parameters_flags()
     ONE_CRUST = .true.
     TRANSVERSE_ISOTROPY = .true.
 
+  case('drex_model')
+    HONOR_1D_SPHERICAL_MOHO = .true.
+    ! CRUSTAL = .true. ! with 3D crust: depends on 3D mantle reference model
+    CASE_3D = .false. ! crustal moho stretching
+    ONE_CRUST = .false. ! if true 1 element layer in top crust region
+    REFERENCE_1D_MODEL = REFERENCE_MODEL_PREM
+    TRANSVERSE_ISOTROPY = .false.
+    MODEL_3D_MANTLE_PERTUBATIONS = .true.
+    ANISOTROPIC_3D_MANTLE = .true. ! treats mantle elements as fully anisotropic
+    THREE_D_MODEL = THREE_D_MODEL_ANISO_MANTLE_DREX
+
   case ('s20rts','s20rts_paper')
     CASE_3D = .true.
     CRUSTAL = .true.