[cig-commits] [commit] devel: changes output format of seismogram names to **network**.**station**.., also for adjoint source filenames (78316b2)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Tue Dec 23 05:06:47 PST 2014
Repository : https://github.com/geodynamics/specfem3d
On branch : devel
Link : https://github.com/geodynamics/specfem3d/compare/c92aeae6613ca38c29e96dd169c292160dc53e21...c3b9c7ef07feb7d057fa71da18323d53b9c27233
>---------------------------------------------------------------
commit 78316b2310efc7e8d15ac2bebd48f1375b4fdefd
Author: daniel peter <peterda at ethz.ch>
Date: Tue Dec 23 10:39:30 2014 +0100
changes output format of seismogram names to **network**.**station**.., also for adjoint source filenames
>---------------------------------------------------------------
78316b2310efc7e8d15ac2bebd48f1375b4fdefd
src/specfem3D/compute_add_sources_acoustic.f90 | 8 +--
src/specfem3D/compute_add_sources_poroelastic.f90 | 4 +-
src/specfem3D/compute_add_sources_viscoelastic.f90 | 8 +--
src/specfem3D/locate_receivers.f90 | 66 +++++++++++-----------
src/specfem3D/setup_sources_receivers.f90 | 4 +-
src/specfem3D/write_seismograms.f90 | 5 +-
6 files changed, 48 insertions(+), 47 deletions(-)
diff --git a/src/specfem3D/compute_add_sources_acoustic.f90 b/src/specfem3D/compute_add_sources_acoustic.f90
index d26c617..88baf48 100644
--- a/src/specfem3D/compute_add_sources_acoustic.f90
+++ b/src/specfem3D/compute_add_sources_acoustic.f90
@@ -257,8 +257,8 @@
if (myrank == islice_selected_rec(irec)) then
irec_local = irec_local + 1
- ! reads in **sta**.**net**.**LH**.adj files
- adj_source_file = trim(station_name(irec))//'.'//trim(network_name(irec))
+ ! reads in **net**.**sta**.**BH**.adj files
+ adj_source_file = trim(network_name(irec))//'.'//trim(station_name(irec))
call compute_arrays_adjoint_source(myrank,adj_source_file, &
xi_receiver(irec),eta_receiver(irec),gamma_receiver(irec), &
adj_sourcearray, xigll,yigll,zigll, &
@@ -711,8 +711,8 @@
if (myrank == islice_selected_rec(irec)) then
irec_local = irec_local + 1
- ! reads in **sta**.**net**.**LH**.adj files
- adj_source_file = trim(station_name(irec))//'.'//trim(network_name(irec))
+ ! reads in **net**.**sta**.**BH**.adj files
+ adj_source_file = trim(network_name(irec))//'.'//trim(station_name(irec))
call compute_arrays_adjoint_source(myrank,adj_source_file, &
xi_receiver(irec),eta_receiver(irec),gamma_receiver(irec), &
adj_sourcearray, xigll,yigll,zigll, &
diff --git a/src/specfem3D/compute_add_sources_poroelastic.f90 b/src/specfem3D/compute_add_sources_poroelastic.f90
index 8f4f7ef..4ef49bd 100644
--- a/src/specfem3D/compute_add_sources_poroelastic.f90
+++ b/src/specfem3D/compute_add_sources_poroelastic.f90
@@ -261,8 +261,8 @@
! compute source arrays
if (myrank == islice_selected_rec(irec)) then
irec_local = irec_local + 1
- ! reads in **sta**.**net**.**LH**.adj files
- adj_source_file = trim(station_name(irec))//'.'//trim(network_name(irec))
+ ! reads in **net**.**sta**.**BH**.adj files
+ adj_source_file = trim(network_name(irec))//'.'//trim(station_name(irec))
call compute_arrays_adjoint_source(myrank,adj_source_file, &
xi_receiver(irec),eta_receiver(irec),gamma_receiver(irec), &
adj_sourcearray, xigll,yigll,zigll, &
diff --git a/src/specfem3D/compute_add_sources_viscoelastic.f90 b/src/specfem3D/compute_add_sources_viscoelastic.f90
index 47ce893..4f9c470 100644
--- a/src/specfem3D/compute_add_sources_viscoelastic.f90
+++ b/src/specfem3D/compute_add_sources_viscoelastic.f90
@@ -246,8 +246,8 @@
! compute source arrays
if (myrank == islice_selected_rec(irec)) then
irec_local = irec_local + 1
- ! reads in **sta**.**net**.**LH**.adj files
- adj_source_file = trim(station_name(irec))//'.'//trim(network_name(irec))
+ ! reads in **net**.**sta**.**BH**.adj files
+ adj_source_file = trim(network_name(irec))//'.'//trim(station_name(irec))
call compute_arrays_adjoint_source(myrank,adj_source_file, &
xi_receiver(irec),eta_receiver(irec),gamma_receiver(irec), &
adj_sourcearray, xigll,yigll,zigll, &
@@ -750,8 +750,8 @@
! compute source arrays
if (myrank == islice_selected_rec(irec)) then
irec_local = irec_local + 1
- ! reads in **sta**.**net**.**LH**.adj files
- adj_source_file = trim(station_name(irec))//'.'//trim(network_name(irec))
+ ! reads in **net**.**sta**.**BH**.adj files
+ adj_source_file = trim(network_name(irec))//'.'//trim(station_name(irec))
call compute_arrays_adjoint_source(myrank,adj_source_file, &
xi_receiver(irec),eta_receiver(irec),gamma_receiver(irec), &
adj_sourcearray, xigll,yigll,zigll, &
diff --git a/src/specfem3D/locate_receivers.f90 b/src/specfem3D/locate_receivers.f90
index 3e2c32f..a83acfe 100644
--- a/src/specfem3D/locate_receivers.f90
+++ b/src/specfem3D/locate_receivers.f90
@@ -180,7 +180,7 @@
! define topology of the control element
call usual_hex_nodes(NGNOD,iaddx,iaddy,iaddz)
- ! opens STATIONS file
+ ! opens STATIONS or STATIONS_ADJOINT file
open(unit=IIN,file=trim(rec_filename),status='old',action='read',iostat=ier)
if (ier /= 0) call exit_mpi(myrank,'error opening file '//trim(rec_filename))
@@ -250,39 +250,39 @@
! allocate memory for arrays using number of stations
allocate(stlat(nrec), &
- stlon(nrec), &
- stele(nrec), &
- stbur(nrec), &
- stutm_x(nrec), &
- stutm_y(nrec), &
- horiz_dist(nrec), &
- elevation(nrec), &
- ix_initial_guess(nrec), &
- iy_initial_guess(nrec), &
- iz_initial_guess(nrec), &
- x_target(nrec), &
- y_target(nrec), &
- z_target(nrec), &
- x_found(nrec), &
- y_found(nrec), &
- z_found(nrec), &
- final_distance(nrec), &
- ispec_selected_rec_all(nrec), &
- xi_receiver_all(nrec), &
- eta_receiver_all(nrec), &
- gamma_receiver_all(nrec), &
- x_found_all(nrec), &
- y_found_all(nrec), &
- z_found_all(nrec), &
- final_distance_all(nrec), &
- nu_all(3,3,nrec),stat=ier)
+ stlon(nrec), &
+ stele(nrec), &
+ stbur(nrec), &
+ stutm_x(nrec), &
+ stutm_y(nrec), &
+ horiz_dist(nrec), &
+ elevation(nrec), &
+ ix_initial_guess(nrec), &
+ iy_initial_guess(nrec), &
+ iz_initial_guess(nrec), &
+ x_target(nrec), &
+ y_target(nrec), &
+ z_target(nrec), &
+ x_found(nrec), &
+ y_found(nrec), &
+ z_found(nrec), &
+ final_distance(nrec), &
+ ispec_selected_rec_all(nrec), &
+ xi_receiver_all(nrec), &
+ eta_receiver_all(nrec), &
+ gamma_receiver_all(nrec), &
+ x_found_all(nrec), &
+ y_found_all(nrec), &
+ z_found_all(nrec), &
+ final_distance_all(nrec), &
+ nu_all(3,3,nrec),stat=ier)
if (ier /= 0) stop 'error allocating arrays for locating receivers'
! loop on all the stations
do irec=1,nrec
read(IIN,*,iostat=ier) station_name(irec),network_name(irec), &
- stlat(irec),stlon(irec),stele(irec),stbur(irec)
+ stlat(irec),stlon(irec),stele(irec),stbur(irec)
if (ier /= 0) call exit_mpi(myrank, 'Error reading station file '//trim(rec_filename))
@@ -298,9 +298,9 @@
if (myrank == 0) then
! limits user output if too many receivers
if (nrec < 1000 .and. (.not. SU_FORMAT)) then
- write(IMAIN,*) 'Station #',irec,': ',station_name(irec)(1:len_trim(station_name(irec))), &
- '.',network_name(irec)(1:len_trim(network_name(irec))), &
- ' horizontal distance: ',sngl(horiz_dist(irec)),' km'
+ write(IMAIN,*) 'Station #',irec,': ', &
+ network_name(irec)(1:len_trim(network_name(irec)))//'.'//station_name(irec)(1:len_trim(station_name(irec))), &
+ ' horizontal distance: ',sngl(horiz_dist(irec)),' km'
endif
endif
@@ -820,7 +820,7 @@
! checks stations location
if (final_distance(irec) == HUGEVAL) then
- write(IMAIN,*) 'error locating station # ',irec,' ',station_name(irec),network_name(irec)
+ write(IMAIN,*) 'error locating station # ',irec,' ',trim(network_name(irec)),' ',trim(station_name(irec))
call exit_MPI(myrank,'error locating receiver')
endif
@@ -828,7 +828,7 @@
if (nrec < 1000 .and. (.not. SU_FORMAT )) then
write(IMAIN,*)
- write(IMAIN,*) 'station # ',irec,' ',station_name(irec),network_name(irec)
+ write(IMAIN,*) 'station # ',irec,' ',trim(network_name(irec)),' ',trim(station_name(irec))
write(IMAIN,*) ' original latitude: ',sngl(stlat(irec))
write(IMAIN,*) ' original longitude: ',sngl(stlon(irec))
diff --git a/src/specfem3D/setup_sources_receivers.f90 b/src/specfem3D/setup_sources_receivers.f90
index c9cc703..ae80ec9 100644
--- a/src/specfem3D/setup_sources_receivers.f90
+++ b/src/specfem3D/setup_sources_receivers.f90
@@ -708,8 +708,8 @@
! updates counter
nadj_rec_local = nadj_rec_local + 1
- ! checks **sta**.**net**.**BH**.adj files for correct number of time steps
- adj_source_file = trim(station_name(irec))//'.'//trim(network_name(irec))
+ ! checks **net**.**sta**.**BH**.adj files for correct number of time steps
+ adj_source_file = trim(network_name(irec))//'.'//trim(station_name(irec))
do icomp = 1,NDIM
filename = OUTPUT_FILES_PATH(1:len_trim(OUTPUT_FILES_PATH))// &
'/../SEM/'//trim(adj_source_file) // '.'// comp(icomp) // '.adj'
diff --git a/src/specfem3D/write_seismograms.f90 b/src/specfem3D/write_seismograms.f90
index 116de34..90ff69d 100644
--- a/src/specfem3D/write_seismograms.f90
+++ b/src/specfem3D/write_seismograms.f90
@@ -522,8 +522,9 @@
if (length_network_name < 1 .or. length_network_name > MAX_LENGTH_NETWORK_NAME) &
call exit_MPI(myrank,'wrong length of network name')
- write(sisname,"(a,'.',a,'.',a3,'.sem',a1)") station_name(irec)(1:length_station_name),&
- network_name(irec)(1:length_network_name),channel,component
+ ! writes out **net**.**sta**.**BH**.sem* files
+ write(sisname,"(a,'.',a,'.',a3,'.sem',a1)") network_name(irec)(1:length_network_name), &
+ station_name(irec)(1:length_station_name),channel,component
! directory to store seismograms
final_LOCAL_PATH = OUTPUT_FILES_PATH(1:len_trim(OUTPUT_FILES_PATH)) // '/'
More information about the CIG-COMMITS
mailing list