[cig-commits] [commit] devel: changed seismogram filename extensions, as discussed (ac61542)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Dec 11 11:44:53 PST 2014
Repository : https://github.com/geodynamics/specfem2d
On branch : devel
Link : https://github.com/geodynamics/specfem2d/compare/9b5755d28b601d380cb365841af5f95863815e3d...ac615429b347b44208cf08988e81863e4b7887e2
>---------------------------------------------------------------
commit ac615429b347b44208cf08988e81863e4b7887e2
Author: rmodrak <rmodrak at princeton.edu>
Date: Thu Dec 11 14:30:44 2014 -0500
changed seismogram filename extensions, as discussed
>---------------------------------------------------------------
ac615429b347b44208cf08988e81863e4b7887e2
src/specfem2D/setup_sources_receivers.F90 | 6 +--
src/specfem2D/write_seismograms.F90 | 84 ++++++++++++++++++++++++-------
2 files changed, 69 insertions(+), 21 deletions(-)
diff --git a/src/specfem2D/setup_sources_receivers.F90 b/src/specfem2D/setup_sources_receivers.F90
index d68dafd..437b751 100644
--- a/src/specfem2D/setup_sources_receivers.F90
+++ b/src/specfem2D/setup_sources_receivers.F90
@@ -184,13 +184,13 @@ subroutine add_adjoint_sources_SU
integer :: ios
irec_local = 0
- write(filename, "('./SEM/Ux_file_single.bin.adj')")
+ write(filename, "('./SEM/Ux_file_single.su.adj')")
open(111,file=trim(filename),access='direct',recl=240+4*NSTEP,iostat = ios)
if (ios /= 0) call exit_MPI(' file '//trim(filename)//'does not exist')
- write(filename, "('./SEM/Uy_file_single.bin.adj')")
+ write(filename, "('./SEM/Uy_file_single.su.adj')")
open(112,file=trim(filename),access='direct',recl=240+4*NSTEP,iostat = ios)
if (ios /= 0) call exit_MPI(' file '//trim(filename)//'does not exist')
- write(filename, "('./SEM/Uz_file_single.bin.adj')")
+ write(filename, "('./SEM/Uz_file_single.su.adj')")
open(113,file=trim(filename),access='direct',recl=240+4*NSTEP,iostat = ios)
if (ios /= 0) call exit_MPI(' file '//trim(filename)//'does not exist')
diff --git a/src/specfem2D/write_seismograms.F90 b/src/specfem2D/write_seismograms.F90
index 9d3e8a5..eb4aa85 100644
--- a/src/specfem2D/write_seismograms.F90
+++ b/src/specfem2D/write_seismograms.F90
@@ -258,6 +258,7 @@
character(len=4) chn
character(len=1) component
+ character(len=4) sfx
character(len=150) sisname
! to write seismograms in single precision SEP and double precision binary format
@@ -310,21 +311,29 @@
number_of_components = NDIM
endif
+
+! filename extension
+ if (.not. SU_FORMAT) then
+ sfx = '.bin'
+ else
+ sfx = '.su'
+ endif
+
allocate(buffer_binary(number_of_components,NSTEP_BETWEEN_OUTPUT_SEISMOS/subsamp_seismos,nrec))
if (save_binary_seismograms .and. myrank == 0 .and. seismo_offset == 0) then
-! delete the old files
+! delete old files
open(unit=12,file='OUTPUT_FILES/Ux_file_single.bin',status='unknown')
close(12,status='delete')
open(unit=12,file='OUTPUT_FILES/Ux_file_double.bin',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/pressure_file_single.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/Uy_file_single.bin',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/pressure_file_double.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/Uy_file_double.bin',status='unknown')
close(12,status='delete')
open(unit=12,file='OUTPUT_FILES/Uz_file_single.bin',status='unknown')
@@ -333,51 +342,90 @@
open(unit=12,file='OUTPUT_FILES/Uz_file_double.bin',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/Curl_file_single.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/Up_file_single.bin',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Up_file_double.bin',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Uc_file_single.bin',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Uc_file_double.bin',status='unknown')
+ close(12,status='delete')
+
+
+! delete old Seismic Unix files
+ open(unit=12,file='OUTPUT_FILES/Ux_file_single.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Ux_file_double.su',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/Curl_file_double.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/Uy_file_single.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Uy_file_double.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Uz_file_single.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Uz_file_double.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Up_file_single.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Up_file_double.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Uc_file_single.su',status='unknown')
+ close(12,status='delete')
+
+ open(unit=12,file='OUTPUT_FILES/Uc_file_double.su',status='unknown')
close(12,status='delete')
endif
+
if (save_binary_seismograms .and. myrank == 0) then
-! write the new files
+ ! write the new files
if(save_binary_seismograms_single) then
if(seismotype == 4 .or. seismotype == 6) then
- open(unit=12,file='OUTPUT_FILES/pressure_file_single.bin',status='unknown',access='direct',recl=4)
+ open(unit=12,file='OUTPUT_FILES/Up_file_single'//sfx,status='unknown',access='direct',recl=4)
else if(.not.p_sv) then
- open(unit=12,file='OUTPUT_FILES/Uy_file_single.bin',status='unknown',access='direct',recl=4)
+ open(unit=12,file='OUTPUT_FILES/Uy_file_single'//sfx,status='unknown',access='direct',recl=4)
else
- open(unit=12,file='OUTPUT_FILES/Ux_file_single.bin',status='unknown',access='direct',recl=4)
+ open(unit=12,file='OUTPUT_FILES/Ux_file_single'//sfx,status='unknown',access='direct',recl=4)
endif
endif
if(save_binary_seismograms_double) then
if(seismotype == 4 .or. seismotype == 6) then
- open(unit=13,file='OUTPUT_FILES/pressure_file_double.bin',status='unknown',access='direct',recl=8)
+ open(unit=13,file='OUTPUT_FILES/Up_file_double'//sfx,status='unknown',access='direct',recl=8)
else if(.not.p_sv) then
- open(unit=13,file='OUTPUT_FILES/Uz_file_double.bin',status='unknown',access='direct',recl=8)
+ open(unit=13,file='OUTPUT_FILES/Uz_file_double'//sfx,status='unknown',access='direct',recl=8)
else
- open(unit=13,file='OUTPUT_FILES/Ux_file_double.bin',status='unknown',access='direct',recl=8)
+ open(unit=13,file='OUTPUT_FILES/Ux_file_double'//sfx,status='unknown',access='direct',recl=8)
endif
endif
-! no Z component seismogram if pressure
+ ! no Z component seismogram if pressure
if(seismotype /= 4 .and. seismotype /= 6 .and. p_sv) then
if(save_binary_seismograms_single) &
- open(unit=14,file='OUTPUT_FILES/Uz_file_single.bin',status='unknown',access='direct',recl=4)
+ open(unit=14,file='OUTPUT_FILES/Uz_file_single'//sfx,status='unknown',access='direct',recl=4)
if(save_binary_seismograms_double) &
- open(unit=15,file='OUTPUT_FILES/Uz_file_double.bin',status='unknown',access='direct',recl=8)
+ open(unit=15,file='OUTPUT_FILES/Uz_file_double'//sfx,status='unknown',access='direct',recl=8)
endif
-! curl output
+ ! curl output
if(seismotype == 5) then
if(save_binary_seismograms_single) &
- open(unit=16,file='OUTPUT_FILES/Curl_file_single.bin',status='unknown',access='direct',recl=4)
+ open(unit=16,file='OUTPUT_FILES/Uc_file_single'//sfx,status='unknown',access='direct',recl=4)
if(save_binary_seismograms_double) &
- open(unit=17,file='OUTPUT_FILES/Curl_file_double.bin',status='unknown',access='direct',recl=8)
+ open(unit=17,file='OUTPUT_FILES/Uc_file_double'//sfx,status='unknown',access='direct',recl=8)
endif
endif
More information about the CIG-COMMITS
mailing list