[cig-commits] r8452 - in seismo/2D/SPECFEM2D/trunk: . DATA
walter at geodynamics.org
walter at geodynamics.org
Fri Dec 7 15:47:11 PST 2007
Author: walter
Date: 2007-12-07 15:47:10 -0800 (Fri, 07 Dec 2007)
New Revision: 8452
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA
seismo/2D/SPECFEM2D/trunk/meshfem2D.f90
seismo/2D/SPECFEM2D/trunk/specfem2D.f90
seismo/2D/SPECFEM2D/trunk/write_seismograms.f90
Log:
cleaned SEP binary seismograms in 2D code and now writes all the output files
to subdirectory OUTPUT_FILES. Also copied final clean Par_files to DATA/
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2005-01-25 20:58:50 UTC (rev 8451)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2007-12-07 23:47:10 UTC (rev 8452)
@@ -7,7 +7,7 @@
#
# title of job, and file that contains interface data
#
-title = Test Sinusoide
+title = Test sinusoide elastique
interfacesfile = interface_sinus.dat
#
# geometry of the model (origin lower-left corner = 0,0) and mesh description
@@ -42,8 +42,11 @@
source_type = 1 ! force = 1 or explosion = 2
time_function_type = 1 ! Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4
f0 = 10.0 ! dominant source frequency (Hz) if not Dirac
-angle = 0. ! angle of the source (for a force only)
-factor = 1.d4 ! amplification factor
+angleforce = 0. ! angle of the source (for a force only)
+Mxx = 1. ! Mxx component (for a moment tensor source only)
+Mzz = 1. ! Mzz component (for a moment tensor source only)
+Mxz = 0. ! Mxz component (for a moment tensor source only)
+factor = 1.d12 ! amplification factor
#
# receiver line parameters
#
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic 2005-01-25 20:58:50 UTC (rev 8451)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic 2007-12-07 23:47:10 UTC (rev 8452)
@@ -7,7 +7,7 @@
#
# title of job, and file that contains interface data
#
-title = Test Sinusoide
+title = Test sinusoide acoustique
interfacesfile = interface_sinus.dat
#
# geometry of the model (origin lower-left corner = 0,0) and mesh description
@@ -42,7 +42,10 @@
source_type = 1 ! force = 1 or explosion = 2
time_function_type = 2 ! Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4
f0 = 10.0 ! dominant source frequency (Hz) if not Dirac
-angle = 0. ! angle of the source (for a force only)
+angleforce = 0. ! angle of the source (for a force only)
+Mxx = 1. ! Mxx component (for a moment tensor source only)
+Mzz = 1. ! Mzz component (for a moment tensor source only)
+Mxz = 0. ! Mxz component (for a moment tensor source only)
factor = 1.d4 ! amplification factor
#
# receiver line parameters
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic 2005-01-25 20:58:50 UTC (rev 8451)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic 2007-12-07 23:47:10 UTC (rev 8452)
@@ -7,7 +7,7 @@
#
# title of job, and file that contains interface data
#
-title = Test Sinusoide
+title = Test sinusoide elastique
interfacesfile = interface_sinus.dat
#
# geometry of the model (origin lower-left corner = 0,0) and mesh description
@@ -42,8 +42,11 @@
source_type = 1 ! force = 1 or explosion = 2
time_function_type = 1 ! Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4
f0 = 10.0 ! dominant source frequency (Hz) if not Dirac
-angle = 0. ! angle of the source (for a force only)
-factor = 1.d4 ! amplification factor
+angleforce = 0. ! angle of the source (for a force only)
+Mxx = 1. ! Mxx component (for a moment tensor source only)
+Mzz = 1. ! Mzz component (for a moment tensor source only)
+Mxz = 0. ! Mxz component (for a moment tensor source only)
+factor = 1.d12 ! amplification factor
#
# receiver line parameters
#
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA 2005-01-25 20:58:50 UTC (rev 8451)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA 2007-12-07 23:47:10 UTC (rev 8452)
@@ -42,7 +42,10 @@
source_type = 1 ! force = 1 or explosion = 2
time_function_type = 1 ! Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4
f0 = 10.0 ! dominant source frequency (Hz) if not Dirac
-angle = 0. ! angle of the source (for a force only)
+angleforce = 0. ! angle of the source (for a force only)
+Mxx = 1. ! Mxx component (for a moment tensor source only)
+Mzz = 1. ! Mzz component (for a moment tensor source only)
+Mxz = 0. ! Mxz component (for a moment tensor source only)
factor = 1.d10 ! amplification factor
#
# receiver line parameters
Modified: seismo/2D/SPECFEM2D/trunk/meshfem2D.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/meshfem2D.f90 2005-01-25 20:58:50 UTC (rev 8451)
+++ seismo/2D/SPECFEM2D/trunk/meshfem2D.f90 2007-12-07 23:47:10 UTC (rev 8452)
@@ -533,7 +533,7 @@
print *
print *,'Ecriture de la grille format Gnuplot...'
- open(unit=20,file='gridfile.gnu',status='unknown')
+ open(unit=20,file='OUTPUT_FILES/gridfile.gnu',status='unknown')
! dessin des lignes horizontales de la grille
print *,'Ecriture lignes horizontales'
@@ -570,7 +570,7 @@
close(20)
! cree le script de dessin pour gnuplot
- open(unit=20,file='plotgnu',status='unknown')
+ open(unit=20,file='OUTPUT_FILES/plotgnu',status='unknown')
write(20,*) '#set term postscript landscape monochrome solid "Helvetica" 22'
write(20,*) '#set output "grille.ps"'
write(20,*) 'plot "gridfile.gnu" title "Macrobloc mesh" w l'
@@ -582,7 +582,7 @@
! *** generation de la base de donnees
- open(unit=15,file='Database',status='unknown')
+ open(unit=15,file='OUTPUT_FILES/Database',status='unknown')
write(15,*) '#'
write(15,*) '# Database for SPECFEM2D'
Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.f90 2005-01-25 20:58:50 UTC (rev 8451)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.f90 2007-12-07 23:47:10 UTC (rev 8452)
@@ -191,7 +191,7 @@
!
!***********************************************************************
- open (IIN,file='Database')
+ open (IIN,file='OUTPUT_FILES/Database')
! uncomment this to write to file instead of standard output
! open (IOUT,file='results_simulation.txt')
Modified: seismo/2D/SPECFEM2D/trunk/write_seismograms.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/write_seismograms.f90 2005-01-25 20:58:50 UTC (rev 8451)
+++ seismo/2D/SPECFEM2D/trunk/write_seismograms.f90 2007-12-07 23:47:10 UTC (rev 8452)
@@ -22,22 +22,27 @@
integer nrec,NSTEP,it,sismostype
double precision t0,deltat
- double precision,dimension(NSTEP,nrec) :: sisux,sisuz
+ double precision, dimension(NSTEP,nrec) :: sisux,sisuz
double precision st_xval(nrec)
character(len=MAX_LENGTH_STATION_NAME), dimension(nrec) :: station_name
character(len=MAX_LENGTH_NETWORK_NAME), dimension(nrec) :: network_name
- integer irec,length_station_name,length_network_name,iorientation,isample
+ integer irec,irecord,length_station_name,length_network_name,iorientation,isample
character(len=4) chn
character(len=1) component
character(len=150) sisname
+! to write seismograms in single precision SEP binary format
+ real(kind=4), dimension(NSTEP*nrec) :: buffer_SEP_binary
+
! scaling factor for Seismic Unix xsu dislay
double precision, parameter :: FACTORXSU = 1.d0
+!----
+
! write seismograms in ASCII format
! save displacement, velocity or acceleration
@@ -101,20 +106,32 @@
!----
-! write seismograms in single precision binary format
+! write seismograms in single precision SEP binary format
! X component
- open(unit=11,file='OUTPUT_FILES/Ux_file.bin',status='unknown',access='direct',recl=NSTEP*4)
+ irecord = 0
do irec=1,nrec
- write(11,rec=irec) (sngl(sisux(isample,irec)),isample=1,NSTEP)
+ do isample=1,NSTEP
+ irecord = irecord + 1
+ buffer_SEP_binary(irecord) = sngl(sisux(isample,irec))
+ enddo
enddo
+
+ open(unit=11,file='OUTPUT_FILES/Ux_file.bin',status='unknown',access='direct',recl=4*NSTEP*nrec)
+ write(11,rec=1) buffer_SEP_binary
close(11)
! Z component
- open(unit=11,file='OUTPUT_FILES/Uz_file.bin',status='unknown',access='direct',recl=NSTEP*4)
+ irecord = 0
do irec=1,nrec
- write(11,rec=irec) (sngl(sisuz(isample,irec)),isample=1,NSTEP)
+ do isample=1,NSTEP
+ irecord = irecord + 1
+ buffer_SEP_binary(irecord) = sngl(sisuz(isample,irec))
+ enddo
enddo
+
+ open(unit=11,file='OUTPUT_FILES/Uz_file.bin',status='unknown',access='direct',recl=4*NSTEP*nrec)
+ write(11,rec=1) buffer_SEP_binary
close(11)
!----
More information about the cig-commits
mailing list