[cig-commits] r19147 - in seismo/3D/SPECFEM3D_GLOBE/trunk: setup src/meshfem3D src/shared

[danielpeter at geodynamics.org]

Author: danielpeter
Date: 2011-11-03 10:44:33 -0700 (Thu, 03 Nov 2011)
New Revision: 19147

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D_models.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/get_model_parameters.f90
Log:
updates some source code comments and initializes crustal flag

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in	2011-11-03 05:52:54 UTC (rev 19146)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/setup/constants.h.in	2011-11-03 17:44:33 UTC (rev 19147)
@@ -551,9 +551,11 @@
 
 !!!!!!!!!!!!!! parameters added for the thread-safe version of the code
 
+!!-----------------------------------------------------------
 !!
 !! model constants
 !!
+!!-----------------------------------------------------------
 
 ! number of layers in DATA/1066a/1066a.dat
   integer, parameter :: NR_1066A = 160
@@ -581,7 +583,7 @@
   integer,parameter :: NKQ=8,MAXL_Q=12
   integer,parameter :: NSQ=(MAXL_Q+1)**2,NDEPTHS_REFQ=913
 
-! The meaningful range of Zhao et al.'s model (1994) is as follows:
+! The meaningful range of Zhao et al. (1994) model is as follows:
 !        latitude : 32 - 45 N
 !        longitude: 130-145 E
 !        depth    : 0  - 500 km
@@ -625,9 +627,11 @@
 
 !!!!!!!!!!!!!! end of parameters added for the thread-safe version of the code
 
+!!-----------------------------------------------------------
 !!
 !! crustal stretching
 !!
+!!-----------------------------------------------------------
 
 ! for the stretching of crustal elements in the case of 3D models
 ! (values are chosen for 3D models to have RMOHO_FICTICIOUS at 35 km
@@ -647,9 +651,11 @@
 !!-- uncomment for deep moho cases, e.g. Asia case (Himalayan moho)
 !  double precision, parameter :: RMOHO_STRETCH_ADJUSTEMENT = -20000.d0  ! moho mesh boundary down to 60km
 
+!!-----------------------------------------------------------
 !!
 !! mesh tweaking
 !!
+!!-----------------------------------------------------------
 
 ! to suppress the crustal layers
 ! (replaced by an extension of the mantle: R_EARTH is not modified, but no more crustal doubling)

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D_models.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D_models.f90	2011-11-03 05:52:54 UTC (rev 19146)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/meshfem3D_models.f90	2011-11-03 17:44:33 UTC (rev 19147)
@@ -1084,7 +1084,8 @@
 ! ADD YOUR MODEL HERE
 !
 !---
-
+  found_crust = .false.
+  
   ! crustal model can vary for different 3-D models
   select case (THREE_D_MODEL )
 

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/get_model_parameters.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/get_model_parameters.f90	2011-11-03 05:52:54 UTC (rev 19146)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/shared/get_model_parameters.f90	2011-11-03 17:44:33 UTC (rev 19147)
@@ -348,9 +348,13 @@
 
   else if(MODEL_ROOT == '3D_anisotropic') then
     ANISOTROPIC_3D_MANTLE = .true.
-    CASE_3D = .true.
-    ONE_CRUST = .true.
-    TRANSVERSE_ISOTROPY = .true.
+    CASE_3D = .true. ! crustal moho stretching
+    ONE_CRUST = .true. ! 1 element layer in top crust region
+    TRANSVERSE_ISOTROPY = .true. ! to use transverse isotropic PREM 1D ref model
+    ! CRUSTAL = .true. ! with 3D crust: depends on 3D mantle reference model
+    ! THREE_D_MODEL = 0 ! for default crustal model
+    ! REFERENCE_1D_MODEL = REFERENCE_MODEL_AK135 
+    ! TRANSVERSE_ISOTROPY = .false. ! for AK135 ref model
 
   else if(MODEL_ROOT == 'heterogen') then
     CASE_3D = .true.