[cig-commits] r8535 - in seismo/2D/SPECFEM2D/trunk: . DATA
walter at geodynamics.org
walter at geodynamics.org
Fri Dec 7 15:54:24 PST 2007
Author: walter
Date: 2007-12-07 15:54:23 -0800 (Fri, 07 Dec 2007)
New Revision: 8535
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/TODO_list
seismo/2D/SPECFEM2D/trunk/meshfem2D.F90
Log:
added option to generate or not STATIONS file (useful for complex geometries) in Par_file.
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2007-06-26 00:56:08 UTC (rev 8534)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2007-12-07 23:54:23 UTC (rev 8535)
@@ -52,6 +52,7 @@
# receiver line parameters for seismograms
seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiverlines = 1 # number of receiver lines
anglerec = 0.d0 # angle to rotate components at receivers
Modified: seismo/2D/SPECFEM2D/trunk/TODO_list
===================================================================
--- seismo/2D/SPECFEM2D/trunk/TODO_list 2007-06-26 00:56:08 UTC (rev 8534)
+++ seismo/2D/SPECFEM2D/trunk/TODO_list 2007-12-07 23:54:23 UTC (rev 8535)
@@ -4,7 +4,6 @@
- partitioning using weights for load-balancing.
- choosing a way to use assign_external_model?
- adding comments.
-- adding option to create or not STATIONS file (useful for complex geometries like asteroids).
- checking for points with different normals for absorbing conditions, when the absorbing edges are not in the same elements (similar to what is done for the corners).
- scripts for translating GID/CUBIT meshes into files for xmeshfem2D.
- modifying scripts for UPPA cluster (when FS sync issues are solved and remote commands are available).
Modified: seismo/2D/SPECFEM2D/trunk/meshfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/meshfem2D.F90 2007-06-26 00:56:08 UTC (rev 8534)
+++ seismo/2D/SPECFEM2D/trunk/meshfem2D.F90 2007-12-07 23:54:23 UTC (rev 8535)
@@ -78,6 +78,7 @@
integer ixdebregion,ixfinregion,izdebregion,izfinregion
integer iregion,imaterial,nbregion,nb_materials
integer NTSTEP_BETWEEN_OUTPUT_INFO,pointsdisp,subsamp,seismotype,imagetype
+ logical generate_STATIONS
integer ngnod,nt,nx,nz,nxread,nzread,icodematread,ireceiverlines,nreceiverlines
integer, dimension(:), allocatable :: nrec
@@ -408,6 +409,7 @@
! read receiver line parameters
call read_value_integer(IIN,IGNORE_JUNK,seismotype)
+ call read_value_logical(IIN,IGNORE_JUNK,generate_STATIONS)
call read_value_integer(IIN,IGNORE_JUNK,nreceiverlines)
call read_value_double_precision(IIN,IGNORE_JUNK,anglerec)
@@ -1243,6 +1245,7 @@
!--- compute position of the receivers and write the STATIONS file
+ if ( generate_STATIONS ) then
print *
print *,'writing the DATA/STATIONS file'
print *
@@ -1293,6 +1296,7 @@
enddo
close(15)
+ end if
print *
More information about the cig-commits
mailing list