[cig-commits] r17984 - in seismo/3D/SPECFEM3D/trunk: src/specfem3D utils/seis_process
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Fri Feb 25 18:19:14 PST 2011
Author: dkomati1
Date: 2011-02-25 18:19:14 -0800 (Fri, 25 Feb 2011)
New Revision: 17984
Added:
seismo/3D/SPECFEM3D/trunk/utils/seis_process/ascii_2_sep.f90
Removed:
seismo/3D/SPECFEM3D/trunk/src/specfem3D/ascii_2_sep.f90
Log:
moved ascii_2_sep.f90 to "utils"
Deleted: seismo/3D/SPECFEM3D/trunk/src/specfem3D/ascii_2_sep.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/ascii_2_sep.f90 2011-02-26 02:16:28 UTC (rev 17983)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/ascii_2_sep.f90 2011-02-26 02:19:14 UTC (rev 17984)
@@ -1,41 +0,0 @@
-
- program ascii_2_sep
-
-! to convert ASCII seismograms to SEP binary format
-
- implicit none
-
-! Parameters to be defined to set the dimension of the arrays:
-! - nr = maximum number of receivers.
-! - ntime = maximum number of points of each seismogram.
- integer, parameter :: nr = 101
- integer, parameter :: ntime = 512
-
- real(kind=4), dimension(ntime,nr) :: sy
-
- integer :: it,ir,i
-
-! read seismogram component in ASCII text format
- print *,'Reading seismograms in text format'
-
- open(unit=11,file='U_file.txt',status='unknown')
- do ir=1,nr
- do it=1,ntime
- read(11,*) sy(it,ir)
-
-! invert sign of vertical component if needed depending on the reference frame convention
-! sy(it,ir) = - sy(it,ir)
-
- enddo
- enddo
- close(11)
-
-! write seismogram component in binary SEP format
- print *,'Saving seismograms in SEP format'
-
- open(unit=11,file='U_file.sep',status='unknown',access='direct',recl=ntime*nr*4)
- write(11,rec=1) (sy(i,1),i=1,ntime*nr)
- close(11)
-
- end program ascii_2_sep
-
Copied: seismo/3D/SPECFEM3D/trunk/utils/seis_process/ascii_2_sep.f90 (from rev 17969, seismo/3D/SPECFEM3D/trunk/src/specfem3D/ascii_2_sep.f90)
===================================================================
--- seismo/3D/SPECFEM3D/trunk/utils/seis_process/ascii_2_sep.f90 (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/utils/seis_process/ascii_2_sep.f90 2011-02-26 02:19:14 UTC (rev 17984)
@@ -0,0 +1,41 @@
+
+ program ascii_2_sep
+
+! to convert ASCII seismograms to SEP binary format
+
+ implicit none
+
+! Parameters to be defined to set the dimension of the arrays:
+! - nr = maximum number of receivers.
+! - ntime = maximum number of points of each seismogram.
+ integer, parameter :: nr = 101
+ integer, parameter :: ntime = 512
+
+ real(kind=4), dimension(ntime,nr) :: sy
+
+ integer :: it,ir,i
+
+! read seismogram component in ASCII text format
+ print *,'Reading seismograms in text format'
+
+ open(unit=11,file='U_file.txt',status='unknown')
+ do ir=1,nr
+ do it=1,ntime
+ read(11,*) sy(it,ir)
+
+! invert sign of vertical component if needed depending on the reference frame convention
+! sy(it,ir) = - sy(it,ir)
+
+ enddo
+ enddo
+ close(11)
+
+! write seismogram component in binary SEP format
+ print *,'Saving seismograms in SEP format'
+
+ open(unit=11,file='U_file.sep',status='unknown',access='direct',recl=ntime*nr*4)
+ write(11,rec=1) (sy(i,1),i=1,ntime*nr)
+ close(11)
+
+ end program ascii_2_sep
+
More information about the CIG-COMMITS
mailing list