[cig-commits] r13414 - in seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS: 1D_ref 1D_ref/blank 3D 3D/blank crustal crustal/blank
leif at geodynamics.org
leif at geodynamics.org
Wed Nov 26 15:47:52 PST 2008
Author: leif
Date: 2008-11-26 15:47:52 -0800 (Wed, 26 Nov 2008)
New Revision: 13414
Added:
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/1D_ref/blank/
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/1D_ref/blank/my_1D_ref_model.f90
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/1D_ref
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/crust
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/crustal/blank/
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/crustal/blank/my_crustal_model.f90
Modified:
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/Makefile
Log:
Completed a blank 3D model which users can use as a template.
Added: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/1D_ref/blank/my_1D_ref_model.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/1D_ref/blank/my_1D_ref_model.f90 (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/1D_ref/blank/my_1D_ref_model.f90 2008-11-26 23:47:52 UTC (rev 13414)
@@ -0,0 +1,119 @@
+!=====================================================================
+!
+! 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_reference_1d_model_radii(ROCEAN,RMIDDLE_CRUST,RMOHO, &
+ R80,R120,R220,R400,R600,R670,R771, &
+ RTOPDDOUBLEPRIME,RCMB,RICB,RHO_TOP_OC,RHO_BOTTOM_OC)
+
+ implicit none
+
+ double precision ROCEAN,RMIDDLE_CRUST,RMOHO, &
+ R80,R120,R220,R400,R600,R670,R771, &
+ RTOPDDOUBLEPRIME,RCMB,RICB,RHO_TOP_OC,RHO_BOTTOM_OC
+
+ ! Define radii here.
+
+end subroutine get_reference_1d_model_radii
+
+
+subroutine reference_1d_model_adjust_time_step(DT)
+
+ double precision DT
+
+ ! If necessary, adjust time step here.
+
+end subroutine reference_1d_model_adjust_time_step
+
+
+subroutine reference_1d_model_aniso_inner_core_model_setup(x,vp,vp0,vs,vs0,rho,rho0)
+
+ implicit none
+
+ double precision x,vp,vp0,vs,vs0,rho,rho0
+
+ ! Leave this blank.
+
+ stop 'anisotropic inner core not implemented'
+
+end subroutine reference_1d_model_aniso_inner_core_model_setup
+
+
+subroutine define_reference_1d_model(USE_EXTERNAL_CRUSTAL_MODEL)
+
+ logical USE_EXTERNAL_CRUSTAL_MODEL
+
+ ! Initialize any data structures here.
+
+end subroutine define_reference_1d_model
+
+
+subroutine reference_1d_model_attenuation_model_setup(AM_V)
+
+ implicit none
+
+ include "constants.h"
+
+ ! attenuation_model_variables
+ type attenuation_model_variables
+ sequence
+ double precision min_period, max_period
+ double precision :: QT_c_source ! Source Frequency
+ double precision, dimension(:), pointer :: Qtau_s ! tau_sigma
+ double precision, dimension(:), pointer :: QrDisc ! Discontinutitues Defined
+ double precision, dimension(:), pointer :: Qr ! Radius
+ integer, dimension(:), pointer :: interval_Q ! Steps
+ double precision, dimension(:), pointer :: Qmu ! Shear Attenuation
+ double precision, dimension(:,:), pointer :: Qtau_e ! tau_epsilon
+ double precision, dimension(:), pointer :: Qomsb, Qomsb2 ! one_minus_sum_beta
+ double precision, dimension(:,:), pointer :: Qfc, Qfc2 ! factor_common
+ double precision, dimension(:), pointer :: Qsf, Qsf2 ! scale_factor
+ integer, dimension(:), pointer :: Qrmin ! Max and Mins of idoubling
+ integer, dimension(:), pointer :: Qrmax ! Max and Mins of idoubling
+ integer :: Qn ! Number of points
+ end type attenuation_model_variables
+
+ type (attenuation_model_variables) AM_V
+ ! attenuation_model_variables
+
+ ! Define attenuation model here.
+
+end subroutine reference_1d_model_attenuation_model_setup
+
+
+subroutine get_reference_1d_model(myrank,r_prem,rho,vpv,vph,vsv,vsh,eta_aniso, &
+ Qkappa,Qmu,idoubling,iregion_code,CRUSTAL,ONE_CRUST,TRANSVERSE_ISOTROPY, &
+ ISOTROPIC_3D_MANTLE)
+
+ implicit none
+
+ logical CRUSTAL,ONE_CRUST,TRANSVERSE_ISOTROPY,ISOTROPIC_3D_MANTLE
+ integer myrank,idoubling,iregion_code
+ double precision r_prem,rho,Qkappa,Qmu,vpv,vph,vsv,vsh,eta_aniso
+
+ ! Implement 1D reference model here.
+
+end subroutine get_reference_1d_model
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/Makefile
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/Makefile 2008-11-26 21:47:49 UTC (rev 13413)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/Makefile 2008-11-26 23:47:52 UTC (rev 13414)
@@ -2,6 +2,7 @@
tar_flags = --exclude=.svn --exclude=build.mk -czhf
portal_models = \
+ blank.tgz \
s20rts.tgz \
s29ea.tgz \
s362ani.tgz \
@@ -16,6 +17,9 @@
clean:
rm *.tgz
+blank.tgz:
+ tar $(tar_flags) blank.tgz blank
+
s20rts.tgz:
tar $(tar_flags) s20rts.tgz s20rts
Added: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/1D_ref
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/1D_ref (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/1D_ref 2008-11-26 23:47:52 UTC (rev 13414)
@@ -0,0 +1 @@
+link ../../1D_ref/blank
\ No newline at end of file
Property changes on: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/1D_ref
___________________________________________________________________
Name: svn:special
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/crust
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/crust (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/crust 2008-11-26 23:47:52 UTC (rev 13414)
@@ -0,0 +1 @@
+link ../../crustal/blank
\ No newline at end of file
Property changes on: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/3D/blank/crust
___________________________________________________________________
Name: svn:special
+ *
Added: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/crustal/blank/my_crustal_model.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/crustal/blank/my_crustal_model.f90 (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/MODELS/crustal/blank/my_crustal_model.f90 2008-11-26 23:47:52 UTC (rev 13414)
@@ -0,0 +1,46 @@
+!=====================================================================
+!
+! 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 read_crust()
+
+ implicit none
+
+ ! Insert code to read data files here.
+
+end subroutine read_crust
+
+
+subroutine crustal_model(lat,lon,x,vp,vs,rho,moho,found_crust)
+
+ implicit none
+
+ double precision lat,lon,x,vp,vs,rho,moho
+ logical found_crust
+
+ ! Define crustal model here.
+
+end subroutine crustal_model
More information about the CIG-COMMITS
mailing list