adds tiso and qmu support to EMC model

Author: danielpeter

src/meshfem3D/meshfem3D_models.F90

@@ -1776,6 +1776,17 @@ subroutine meshfem3D_models_getatten_val(idoubling,r_prem,theta,phi, &
 
   endif
 
+#ifdef USE_EMC
+  ! IRIS EMC model
+  ! over-imposes Qmu values
+  if (EMC_MODEL_QMU) then
+    ! EMC model w/ 3D Qmu parameterization
+    r_used = r_prem  ! takes actual position (between CMB and surface)
+    call model_EMC_attenuation(r_used,theta,phi,Qmu)
+  endif
+#endif
+
+
   ! sponge layer
   if (ABSORB_USING_GLOBAL_SPONGE) then
     ! get distance to chunk center

src/meshfem3D/meshfem3D_par.f90

@@ -90,7 +90,7 @@ module meshfem_models_par
     MODEL_3D_MANTLE_PERTUBATIONS,HETEROGEN_3D_MANTLE,ANISOTROPIC_3D_MANTLE,ANISOTROPIC_INNER_CORE, &
     ATTENUATION_3D, ATTENUATION_GLL, &
     CEM_REQUEST, CEM_ACCEPT, &
-    EMC_MODEL, &
+    EMC_MODEL,EMC_MODEL_QMU, &
     NX_BATHY,NY_BATHY, &
     R_PLANET,R_PLANET_KM,R_DEEPEST_CRUST,REGIONAL_MOHO_MESH