[cig-commits] r15115 - seismo/3D/SPECFEM3D_SESAME/trunk
pieyre at geodynamics.org
pieyre at geodynamics.org
Thu Jun 4 02:19:36 PDT 2009
Author: pieyre
Date: 2009-06-04 02:19:36 -0700 (Thu, 04 Jun 2009)
New Revision: 15115
Modified:
seismo/3D/SPECFEM3D_SESAME/trunk/locate_receivers.f90
seismo/3D/SPECFEM3D_SESAME/trunk/specfem3D.f90
Log:
for clarity : modified required structure of the DATA/STATIONS_ADJOINT file for adjoint simulations, the first line doesn't have to contain the number of stations any more
Modified: seismo/3D/SPECFEM3D_SESAME/trunk/locate_receivers.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/locate_receivers.f90 2009-06-03 21:08:01 UTC (rev 15114)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/locate_receivers.f90 2009-06-04 09:19:36 UTC (rev 15115)
@@ -824,12 +824,20 @@
double precision stlat,stlon,stele,stbur
character(len=MAX_LENGTH_STATION_NAME) station_name
character(len=MAX_LENGTH_NETWORK_NAME) network_name
+ character(len=150) dummystring
+ nrec = 0
nrec_filtered = 0
open(unit=IIN, file=trim(filename), status = 'old', iostat = ios)
if (ios /= 0) call exit_mpi(myrank, 'No file '//trim(filename)//', exit')
- read(IIN, *) nrec
+ do while(ios == 0)
+ read(IIN,"(a)",iostat = ios) dummystring
+ if(ios == 0) nrec = nrec + 1
+ enddo
+ close(IIN)
+
+ open(unit=IIN, file=trim(filename), status = 'old', iostat = ios)
do irec = 1, nrec
read(IIN, *) station_name, network_name, stlat, stlon, stele, stbur
if(stlat >= LATITUDE_MIN .and. stlat <= LATITUDE_MAX .and. stlon >= LONGITUDE_MIN .and. stlon <= LONGITUDE_MAX) &
Modified: seismo/3D/SPECFEM3D_SESAME/trunk/specfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/specfem3D.f90 2009-06-03 21:08:01 UTC (rev 15114)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/specfem3D.f90 2009-06-04 09:19:36 UTC (rev 15115)
@@ -1180,7 +1180,7 @@
call get_value_string(filtered_rec_filename, 'solver.STATIONS_FILTERED', 'DATA/STATIONS_ADJOINT_FILTERED')
call station_filter(myrank,rec_filename,filtered_rec_filename,nrec, &
LATITUDE_MIN, LATITUDE_MAX, LONGITUDE_MIN, LONGITUDE_MAX)
- if (nrec < 1) call exit_MPI(myrank, 'adjoint simulation needs at least one source')
+ if (nrec < 1) call exit_MPI(myrank, 'adjoint simulation needs at least one receiver')
call sync_all()
endif
More information about the CIG-COMMITS
mailing list