[cig-commits] r22739 - in seismo/2D/SPECFEM2D/trunk: UTILS src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Thu Aug 29 14:30:29 PDT 2013
Author: dkomati1
Date: 2013-08-29 14:30:28 -0700 (Thu, 29 Aug 2013)
New Revision: 22739
Modified:
seismo/2D/SPECFEM2D/trunk/UTILS/su_Seismic_Unix_seismograms_create.py
seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/write_seismograms.F90
Log:
changed .bin to .su for binary files in order to clarify that Seismic Unix can be used to process or visualize them
Modified: seismo/2D/SPECFEM2D/trunk/UTILS/su_Seismic_Unix_seismograms_create.py
===================================================================
--- seismo/2D/SPECFEM2D/trunk/UTILS/su_Seismic_Unix_seismograms_create.py 2013-08-29 16:55:41 UTC (rev 22738)
+++ seismo/2D/SPECFEM2D/trunk/UTILS/su_Seismic_Unix_seismograms_create.py 2013-08-29 21:30:28 UTC (rev 22739)
@@ -52,11 +52,11 @@
# Create headers and Su file
if seismotype==4:
- ordres=['suaddhead < pressure_file_single.bin ns='+str(nt)+' | sushw key=dt a='+str(int(deltat*1e6))+' > pressure_file.su']
+ ordres=['suaddhead < pressure_file_single.su ns='+str(nt)+' | sushw key=dt a='+str(int(deltat*1e6))+' > pressure_file.su']
ordres.append('suxwigb < pressure_file.su perc=96 '+labels+' title=" Pressure : '+title+'"&')
else:
- ordres=['suaddhead < Ux_file_single.bin ns='+str(nt)+' | sushw key=dt a='+str(int(deltat*1e6))+' > Ux_file.su']
- ordres.append('suaddhead < Uz_file_single.bin ns='+str(nt)+' | sushw key=dt a='+str(int(deltat*1e6))+' > Uz_file.su')
+ ordres=['suaddhead < Ux_file_single.su ns='+str(nt)+' | sushw key=dt a='+str(int(deltat*1e6))+' > Ux_file.su']
+ ordres.append('suaddhead < Uz_file_single.su ns='+str(nt)+' | sushw key=dt a='+str(int(deltat*1e6))+' > Uz_file.su')
ordres.append('supswigp < Ux_file.su > Ux_file.ps '+labels+' title=" Ux : '+title+'"&')
ordres.append('supswigp < Uz_file.su > Uz_file.ps '+labels+' title=" Uz : '+title+'"&')
ordres.append('suxwigb < Ux_file.su perc=96 '+labels+' title=" Ux : '+title+'"&')
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2013-08-29 16:55:41 UTC (rev 22738)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2013-08-29 21:30:28 UTC (rev 22739)
@@ -2335,10 +2335,10 @@
enddo
else
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/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')
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/write_seismograms.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/write_seismograms.F90 2013-08-29 16:55:41 UTC (rev 22738)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/write_seismograms.F90 2013-08-29 21:30:28 UTC (rev 22739)
@@ -147,28 +147,28 @@
if (save_binary_seismograms .and. myrank == 0 .and. seismo_offset == 0) then
! delete the old files
- open(unit=12,file='OUTPUT_FILES/Ux_file_single.bin',status='unknown')
+ 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.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/Ux_file_double.su',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/pressure_file_single.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/pressure_file_single.su',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/pressure_file_double.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/pressure_file_double.su',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/Uz_file_single.bin',status='unknown')
+ 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.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/Uz_file_double.su',status='unknown')
close(12,status='delete')
- open(unit=12,file='OUTPUT_FILES/Curl_file_single.bin',status='unknown')
+ open(unit=12,file='OUTPUT_FILES/Curl_file_single.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/Curl_file_double.su',status='unknown')
close(12,status='delete')
endif
@@ -178,38 +178,38 @@
! 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/pressure_file_single.su',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.su',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.su',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/pressure_file_double.su',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.su',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.su',status='unknown',access='direct',recl=8)
endif
endif
! 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.su',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.su',status='unknown',access='direct',recl=8)
endif
! 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/Curl_file_single.su',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/Curl_file_double.su',status='unknown',access='direct',recl=8)
endif
endif
More information about the CIG-COMMITS
mailing list