[cig-commits] r8535 - in seismo/2D/SPECFEM2D/trunk: . DATA

[walter at geodynamics.org]

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 *