[cig-commits] r13918 - in seismo/2D/SPECFEM2D/branches/BIOT: . DATA

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-01-22 06:28:48 -0800 (Thu, 22 Jan 2009)
New Revision: 13918

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION
   seismo/2D/SPECFEM2D/branches/BIOT/DATA/Par_file
   seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90
   seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
Log:
BIOT2D: modif input files


Modified: seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION	2009-01-22 05:22:55 UTC (rev 13917)
+++ seismo/2D/SPECFEM2D/branches/BIOT/DATA/CMTSOLUTION	2009-01-22 14:28:48 UTC (rev 13918)
@@ -11,16 +11,3 @@
 Mzz                             = 1.             # Mzz component (for a moment tensor source only)
 Mxz                             = 0.             # Mxz component (for a moment tensor source only)
 factor                          = 1.d10          # amplification factor
-#source 2
-source_surf                     = .false.        # source inside the medium or at the surface
-xs                              = 1800.          # source location x in meters
-zs                              = 2900.          # source location z in meters
-source_type                     = 1              # elastic force or acoustic pressure = 1 or moment tensor = 2
-time_function_type              = 1              # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
-f0                              = 15.0           # dominant source frequency (Hz) if not Dirac or Heaviside
-t0                              = 0.0            # time shift when multi sources (if one source, must be zero)
-angleforce                      = 0.             # angle of the source (for a force only)
-Mxx                             = 1.             # Mxx component (for a moment tensor source only)
-Mzz                             = 1.             # Mzz component (for a moment tensor source only)
-Mxz                             = 0.             # Mxz component (for a moment tensor source only)
-factor                          = 1.d10          # amplification factor

Modified: seismo/2D/SPECFEM2D/branches/BIOT/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/DATA/Par_file	2009-01-22 05:22:55 UTC (rev 13917)
+++ seismo/2D/SPECFEM2D/branches/BIOT/DATA/Par_file	2009-01-22 14:28:48 UTC (rev 13918)
@@ -10,7 +10,7 @@
 materials_file                  = ./DATA/yangluo_mesh_overthrust/materials_file   # file containing the material number for each element
 free_surface_file               = ./DATA/yangluo_mesh_overthrust/free_surface_file   # file containing the free surface
 absorbing_surface_file          = ./DATA/yangluo_mesh_overthrust/absorbing_surface_file   # file containing the absorbing surface
-receivers_file                  = ./DATA/yangluo_mesh_overthrust/receivers_file   # file containing the absorbing surface
+receivers_file                  = ./DATA/yangluo_mesh_overthrust/receivers_file   # file containing the receivers coordinates 
 
 # parameters concerning partitionning
 nproc                           = 1              # number of processes

Modified: seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90	2009-01-22 05:22:55 UTC (rev 13917)
+++ seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90	2009-01-22 14:28:48 UTC (rev 13918)
@@ -787,8 +787,11 @@
 
         ! check if we are in the last layer, which contains topography,
         ! and modify the position of the source accordingly if it is located exactly at the surface
-        if(source_surf(1) .and. ilayer == number_of_layers) & !yang use first source
-             zs = value_spline(xs(1),xinterface_top,zinterface_top,coefs_interface_top,npoints_interface_top)
+        do i_source =1,NSOURCE
+        if(source_surf(i_source) .and. ilayer == number_of_layers) then 
+             zs(i_source) = value_spline(xs(i_source),xinterface_top,zinterface_top,coefs_interface_top,npoints_interface_top)
+        endif
+        enddo
 
         ! compute the offset of this layer in terms of number of spectral elements below along Z
         if(ilayer > 1) then

Modified: seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2009-01-22 05:22:55 UTC (rev 13917)
+++ seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2009-01-22 14:28:48 UTC (rev 13918)
@@ -5540,7 +5540,7 @@
     write(IOUT,*) 'drawing image of vertical component of displacement vector...'
     endif
 
-    call compute_vector_whole_medium(b_potential_acoustic,b_displ_elastic,b_displs_poroelastic,&
+    call compute_vector_whole_medium(potential_acoustic,displ_elastic,displs_poroelastic,&
           elastic,poroelastic,vector_field_display, &
           xix,xiz,gammax,gammaz,ibool,hprime_xx,hprime_zz,nspec,npoin)