[cig-commits] r13299 - in seismo/3D/SPECFEM3D_GLOBE/branches/pluggable: . MODELS/3D_mantle MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0
leif at geodynamics.org
leif at geodynamics.org
Wed Nov 12 13:19:08 PST 2008
Author: leif
Date: 2008-11-12 13:19:08 -0800 (Wed, 12 Nov 2008)
New Revision: 13299
Added:
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/build.mk
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/prem_plus_crust.f90
Modified:
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/Makefile.in
Log:
Resurrected the model 'transversely_isotropic_prem_plus_3D_crust_2.0',
because it is offered as an option by the portal (under the name
"crust2.0+prem"). From the plug-in perspective, this particular 3D
model is basically a no-op, since it sets both ISOTROPIC_3D_MANTLE and
ANISOTROPIC_3D_MANTLE to false.
Added: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/build.mk
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/build.mk (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/build.mk 2008-11-12 21:19:08 UTC (rev 13299)
@@ -0,0 +1,9 @@
+
+include $S/MODELS/1D_ref/prem/build.mk
+
+mantle_model_OBJECTS = \
+ $O/prem_plus_crust.o \
+ $(EMPTY_MACRO)
+
+$O/prem_plus_crust.o: constants.h $S/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/prem_plus_crust.f90
+ ${MPIFCCOMPILE_CHECK} -c -o $O/prem_plus_crust.o ${FCFLAGS_f90} $S/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/prem_plus_crust.f90
Added: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/prem_plus_crust.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/prem_plus_crust.f90 (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/prem_plus_crust.f90 2008-11-12 21:19:08 UTC (rev 13299)
@@ -0,0 +1,116 @@
+!=====================================================================
+!
+! S p e c f e m 3 D G l o b e V e r s i o n 4 . 0
+! --------------------------------------------------
+!
+! Main authors: Dimitri Komatitsch and Jeroen Tromp
+! Seismological Laboratory, California Institute of Technology, USA
+! and University of Pau / CNRS / INRIA, France
+! (c) California Institute of Technology and University of Pau / CNRS / INRIA
+! February 2008
+!
+! This program is free software; you can redistribute it and/or modify
+! it under the terms of the GNU General Public License as published by
+! the Free Software Foundation; either version 2 of the License, or
+! (at your option) any later version.
+!
+! This program is distributed in the hope that it will be useful,
+! but WITHOUT ANY WARRANTY; without even the implied warranty of
+! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+! GNU General Public License for more details.
+!
+! You should have received a copy of the GNU General Public License along
+! with this program; if not, write to the Free Software Foundation, Inc.,
+! 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
+!
+!=====================================================================
+
+
+subroutine get_model_properties(HONOR_1D_SPHERICAL_MOHO,ONE_CRUST, &
+ TRANSVERSE_ISOTROPY, &
+ ISOTROPIC_3D_MANTLE,ANISOTROPIC_3D_MANTLE,ANISOTROPIC_INNER_CORE, &
+ CRUSTAL,CASE_3D, &
+ ATTENUATION_3D)
+
+ implicit none
+
+ logical HONOR_1D_SPHERICAL_MOHO,ONE_CRUST,&
+ TRANSVERSE_ISOTROPY,&
+ ISOTROPIC_3D_MANTLE,ANISOTROPIC_3D_MANTLE,ANISOTROPIC_INNER_CORE,&
+ CRUSTAL,CASE_3D,&
+ ATTENUATION_3D
+
+ TRANSVERSE_ISOTROPY = .true.
+ ISOTROPIC_3D_MANTLE = .false.
+ ANISOTROPIC_3D_MANTLE = .false.
+ ANISOTROPIC_INNER_CORE = .false.
+ CRUSTAL = .true.
+ ATTENUATION_3D = .false.
+ ONE_CRUST = .true.
+ CASE_3D = .true.
+
+end subroutine get_model_properties
+
+
+subroutine read_3d_mantle_model(myrank)
+
+ implicit none
+
+ integer myrank
+
+end subroutine read_3d_mantle_model
+
+
+subroutine iso_mantle_model(radius,theta,phi, &
+ vpv,vph,vsv,vsh,rho,eta_aniso)
+
+ implicit none
+
+ double precision radius,theta,phi
+ double precision vpv,vph,vsv,vsh,rho,eta_aniso
+ double precision dvp,dvs,drho
+
+end subroutine iso_mantle_model
+
+
+subroutine aniso_mantle_model(r,theta,phi,rho, &
+ c11,c12,c13,c14,c15,c16,c22,c23,c24,c25,c26,c33,c34,c35,c36,c44,c45,c46,c55,c56,c66)
+
+ implicit none
+
+ double precision r,theta,phi
+ double precision rho
+ double precision c11,c12,c13,c14,c15,c16,c22,c23,c24,c25,c26, &
+ c33,c34,c35,c36,c44,c45,c46,c55,c56,c66
+
+end subroutine aniso_mantle_model
+
+
+subroutine add_moho_topography(myrank,xelm,yelm,zelm,RMOHO,R220)
+
+ implicit none
+
+ include "constants.h"
+
+ integer myrank
+ double precision xelm(NGNOD)
+ double precision yelm(NGNOD)
+ double precision zelm(NGNOD)
+ double precision RMOHO,R220
+
+end subroutine add_moho_topography
+
+
+subroutine add_topography_410_650(myrank,xelm,yelm,zelm,R220,R400,R670,R771)
+
+ implicit none
+
+ include "constants.h"
+
+ integer myrank
+ double precision xelm(NGNOD)
+ double precision yelm(NGNOD)
+ double precision zelm(NGNOD)
+ double precision R220,R400,R670,R771
+
+end subroutine add_topography_410_650
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/Makefile.in
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/Makefile.in 2008-11-12 20:38:05 UTC (rev 13298)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/Makefile.in 2008-11-12 21:19:08 UTC (rev 13299)
@@ -189,6 +189,7 @@
#include $S/MODELS/3D_mantle/s362ani_prem/build.mk
#include $S/MODELS/3D_mantle/s362iso/build.mk
#include $S/MODELS/3D_mantle/s362wmani/build.mk
+#include $S/MODELS/3D_mantle/transversely_isotropic_prem_plus_3D_crust_2.0/build.mk
#include $S/MODELS/3D_mantle/Montagner/build.mk
#include $S/MODELS/3D_mantle/none/build.mk
More information about the CIG-COMMITS
mailing list