[cig-commits] r21764 - in seismo/3D/SPECFEM3D/branches/coupling_vadim/src: shared specfem3D

[vadim at geodynamics.org]

Author: vadim
Date: 2013-04-08 09:54:30 -0700 (Mon, 08 Apr 2013)
New Revision: 21764

Modified:
   seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/constants.h.in
   seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/read_parameter_file.f90
   seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/initialize_simulation.f90
   seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/read_mesh_databases.f90
   seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/specfem3D_par.f90
Log:
open dsm tract file


Modified: seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/constants.h.in
===================================================================
--- seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/constants.h.in	2013-04-08 16:39:25 UTC (rev 21763)
+++ seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/constants.h.in	2013-04-08 16:54:30 UTC (rev 21764)
@@ -430,4 +430,4 @@
 
 ! vadim tests
    logical, parameter :: USE_VADIM=.true.
-
+   integer, parameter :: IIN_veloc_dsm = 51, IIN_tract_dsm = 52, Ntime_step_dsm = 100 

Modified: seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/read_parameter_file.f90
===================================================================
--- seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/read_parameter_file.f90	2013-04-08 16:39:25 UTC (rev 21763)
+++ seismo/3D/SPECFEM3D/branches/coupling_vadim/src/shared/read_parameter_file.f90	2013-04-08 16:54:30 UTC (rev 21764)
@@ -171,7 +171,7 @@
   call read_value_logical(PRINT_SOURCE_TIME_FUNCTION, 'solver.PRINT_SOURCE_TIME_FUNCTION')
   if(err_occurred() /= 0) return
   !! VM VM read the traction path directory
-  if (USE_DSM) then
+  if (USE_VADIM) then
      call read_value_string(TRAC_PATH, 'TRAC_PATH')
      if(err_occurred() /= 0) return
   end if

Modified: seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/initialize_simulation.f90
===================================================================
--- seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/initialize_simulation.f90	2013-04-08 16:39:25 UTC (rev 21763)
+++ seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/initialize_simulation.f90	2013-04-08 16:54:30 UTC (rev 21764)
@@ -135,6 +135,7 @@
 
   ! reads in numbers of spectral elements and points for this process' domain
   call create_name_database(prname,myrank,LOCAL_PATH)
+  if (USE_VADIM) call create_name_database(dsmname,myrank,TRAC_PATH)  !! VM VM 
   open(unit=27,file=prname(1:len_trim(prname))//'external_mesh.bin',status='old',&
         action='read',form='unformatted',iostat=ier)
   if( ier /= 0 ) then

Modified: seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/read_mesh_databases.f90
===================================================================
--- seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/read_mesh_databases.f90	2013-04-08 16:39:25 UTC (rev 21763)
+++ seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/read_mesh_databases.f90	2013-04-08 16:54:30 UTC (rev 21764)
@@ -375,6 +375,15 @@
           abs_boundary_ijk(3,NGLLSQUARE,num_abs_boundary_faces), &
           abs_boundary_jacobian2Dw(NGLLSQUARE,num_abs_boundary_faces), &
           abs_boundary_normal(NDIM,NGLLSQUARE,num_abs_boundary_faces),stat=ier)
+  if (USE_VADIM) then
+     ! VM for DSM/SEM Hybrid method
+     allocate(Veloc_dsm_boundary(3,Ntime_step_dsm,NGLLSQUARE,num_abs_boundary_faces))
+     allocate(Tract_dsm_boundary(3,Ntime_step_dsm,NGLLSQUARE,num_abs_boundary_faces))
+     open(unit=IIN_veloc_dsm,file=dsmname(1:len_trim(dsmname))//'vel.bin',status='old',&
+          action='read',form='unformatted',iostat=ier)
+     open(unit=IIN_tract_dsm,file=dsmname(1:len_trim(dsmname))//'tract.bin',status='old',&
+          action='read',form='unformatted',iostat=ier)
+  end if
   if( ier /= 0 ) stop 'error allocating array abs_boundary_ispec etc.'
   if( num_abs_boundary_faces > 0 ) then
     read(27) abs_boundary_ispec

Modified: seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/specfem3D_par.f90
===================================================================
--- seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/specfem3D_par.f90	2013-04-08 16:39:25 UTC (rev 21763)
+++ seismo/3D/SPECFEM3D/branches/coupling_vadim/src/specfem3D/specfem3D_par.f90	2013-04-08 16:54:30 UTC (rev 21764)
@@ -86,6 +86,9 @@
   integer, dimension(:), allocatable :: free_surface_ispec
   integer :: num_free_surface_faces
 
+! VM for DSM/SEM Hybrid method
+  real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: Veloc_dsm_boundary,Tract_dsm_boundary
+
 ! attenuation
   integer :: NSPEC_ATTENUATION_AB,NSPEC_ATTENUATION_AB_kappa !ZN
   character(len=256) prname_Q