[cig-commits] r18571 - in seismo/3D/CPML: tags/v1.1.3 trunk
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Thu Jun 9 15:55:55 PDT 2011
Author: dkomati1
Date: 2011-06-09 15:55:55 -0700 (Thu, 09 Jun 2011)
New Revision: 18571
Modified:
seismo/3D/CPML/tags/v1.1.3/seismic_ADEPML_2D_RK4_eighth_order.f90
seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_fourth_order.f90
seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_second_order.f90
seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_poroelastic_fourth_order.f90
seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_isotropic_MPI_OpenMP.f90
seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_viscoelastic_MPI.f90
seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_anisotropic_fourth.f90
seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_isotropic.f90
seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_3D_isotropic_OpenMP.f90
seismo/3D/CPML/trunk/seismic_ADEPML_2D_RK4_eighth_order.f90
seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_fourth_order.f90
seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_second_order.f90
seismo/3D/CPML/trunk/seismic_CPML_2D_poroelastic_fourth_order.f90
seismo/3D/CPML/trunk/seismic_CPML_3D_isotropic_MPI_OpenMP.f90
seismo/3D/CPML/trunk/seismic_CPML_3D_viscoelastic_MPI.f90
seismo/3D/CPML/trunk/seismic_PML_Collino_2D_anisotropic_fourth.f90
seismo/3D/CPML/trunk/seismic_PML_Collino_2D_isotropic.f90
seismo/3D/CPML/trunk/seismic_PML_Collino_3D_isotropic_OpenMP.f90
Log:
added sngl() to output of Gnuplot energy curves, and changed w l to w l lc to use new Gnuplot syntax
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_ADEPML_2D_RK4_eighth_order.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_ADEPML_2D_RK4_eighth_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_ADEPML_2D_RK4_eighth_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1055,8 +1055,8 @@
NPOINTS_PML,USE_PML_XMIN,USE_PML_XMAX,USE_PML_YMIN,USE_PML_YMAX,2)
open(unit=20,file='energy.dat',status='unknown')
do it2 = 1,NSTEP
- write(20,*) sngl(dble(it2-1)*DELTAT),total_energy_kinetic(it2), &
- total_energy_potential(it2),total_energy_kinetic(it2) + total_energy_potential(it2)
+ write(20,*) sngl(dble(it2-1)*DELTAT),sngl(total_energy_kinetic(it2)), &
+ sngl(total_energy_potential(it2)),sngl(total_energy_kinetic(it2) + total_energy_potential(it2))
enddo
close(20)
call write_seismograms(sisvx,sisvy,NSTEP,NREC,DELTAT)
@@ -1071,8 +1071,8 @@
! save total energy
open(unit=20,file='energy.dat',status='unknown')
do it = 1,NSTEP
- write(20,*) sngl(dble(it-1)*DELTAT),total_energy_kinetic(it), &
- total_energy_potential(it),total_energy_kinetic(it) + total_energy_potential(it)
+ write(20,*) sngl(dble(it-1)*DELTAT),sngl(total_energy_kinetic(it)), &
+ sngl(total_energy_potential(it)),sngl(total_energy_kinetic(it) + total_energy_potential(it))
enddo
close(20)
@@ -1086,8 +1086,8 @@
write(20,*)
write(20,*) 'set output "cpml_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1101,8 +1101,8 @@
write(20,*)
write(20,*) 'set output "compare_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l 1, &
- & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l 2'
+ write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l lc 1, &
+ & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l lc 2'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1117,22 +1117,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_fourth_order.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_fourth_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_fourth_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -750,8 +750,8 @@
! save total energy
open(unit=20,file='energy.dat',status='unknown')
do it = 1,NSTEP
- write(20,*) sngl(dble(it-1)*DELTAT),total_energy_kinetic(it), &
- total_energy_potential(it),total_energy_kinetic(it) + total_energy_potential(it)
+ write(20,*) sngl(dble(it-1)*DELTAT),sngl(total_energy_kinetic(it)), &
+ sngl(total_energy_potential(it)),sngl(total_energy_kinetic(it) + total_energy_potential(it))
enddo
close(20)
@@ -765,8 +765,8 @@
write(20,*)
write(20,*) 'set output "cpml_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -780,8 +780,8 @@
write(20,*)
write(20,*) 'set output "compare_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l 1, &
- & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l 2'
+ write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l lc 1, &
+ & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l lc 2'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -796,22 +796,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_second_order.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_second_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_isotropic_second_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -749,8 +749,8 @@
! save total energy
open(unit=20,file='energy.dat',status='unknown')
do it = 1,NSTEP
- write(20,*) sngl(dble(it-1)*DELTAT),total_energy_kinetic(it), &
- total_energy_potential(it),total_energy_kinetic(it) + total_energy_potential(it)
+ write(20,*) sngl(dble(it-1)*DELTAT),sngl(total_energy_kinetic(it)), &
+ sngl(total_energy_potential(it)),sngl(total_energy_kinetic(it) + total_energy_potential(it))
enddo
close(20)
@@ -764,8 +764,8 @@
write(20,*)
write(20,*) 'set output "cpml_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -779,8 +779,8 @@
write(20,*)
write(20,*) 'set output "compare_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l 1, &
- & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l 2'
+ write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l lc 1, &
+ & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l lc 2'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -795,22 +795,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_poroelastic_fourth_order.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_poroelastic_fourth_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_CPML_2D_poroelastic_fourth_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1010,8 +1010,8 @@
write(20,*) 'set ylabel "Total energy"'
write(20,*)
write(20,*) '# set output "cpml_total_energy.eps"'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1028,25 +1028,25 @@
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
write(20,*) '#set yrange [-4:4.5]'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
write(20,*) '#set yrange [-15:19]'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
write(20,*) '#set yrange [-12:16]'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
write(20,*) '#set yrange [-7:10]'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_isotropic_MPI_OpenMP.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_isotropic_MPI_OpenMP.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_isotropic_MPI_OpenMP.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1260,7 +1260,7 @@
write(20,*)
write(20,*) 'set output "CPML3D_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" t ''Total energy'' w l 1'
+ write(20,*) 'plot "energy.dat" t ''Total energy'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1275,32 +1275,32 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_001.eps"'
- write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_002.eps"'
- write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_viscoelastic_MPI.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_viscoelastic_MPI.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_CPML_3D_viscoelastic_MPI.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1430,8 +1430,8 @@
! save energy
open(unit=21,file='energy.dat',status='unknown')
do it2=1,NSTEP
- write(21,*) sngl(dble(it2-1)*DELTAT),total_energy_kinetic(it2),&
- total_energy_potential(it2),total_energy(it2)
+ write(21,*) sngl(dble(it2-1)*DELTAT),sngl(total_energy_kinetic(it2)),&
+ sngl(total_energy_potential(it2)),sngl(total_energy(it2))
enddo
close(21)
@@ -1455,11 +1455,6 @@
! save seismograms
call write_seismograms(sisvx,sisvy,NSTEP,NREC,DELTAT,t0)
-! open(unit=20,file='energy.dat',status='unknown')
-! do it = 1,NSTEP
-! write(20,*) sngl(dble(it-1)*DELTAT),total_energy(it)
-! enddo
-! close(20)
! create script for Gnuplot for total energy
open(unit=20,file='plot_energy',status='unknown')
@@ -1471,7 +1466,7 @@
write(20,*)
write(20,*) 'set output "CPML3D_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" t ''Total energy'' w l 1'
+ write(20,*) 'plot "energy.dat" t ''Total energy'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1486,32 +1481,32 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_001.eps"'
- write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_002.eps"'
- write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_anisotropic_fourth.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_anisotropic_fourth.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_anisotropic_fourth.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -646,8 +646,8 @@
write(20,*)
write(20,*) 'set output "collino_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -662,22 +662,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_isotropic.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_isotropic.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_2D_isotropic.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -578,8 +578,8 @@
write(20,*)
write(20,*) 'set output "collino_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -594,22 +594,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_3D_isotropic_OpenMP.f90
===================================================================
--- seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_3D_isotropic_OpenMP.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/tags/v1.1.3/seismic_PML_Collino_3D_isotropic_OpenMP.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -896,7 +896,7 @@
write(20,*)
write(20,*) 'set output "collino3D_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" t ''Total energy'' w l 1'
+ write(20,*) 'plot "energy.dat" t ''Total energy'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -911,32 +911,32 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_001.eps"'
- write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_002.eps"'
- write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_ADEPML_2D_RK4_eighth_order.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_ADEPML_2D_RK4_eighth_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_ADEPML_2D_RK4_eighth_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1055,8 +1055,8 @@
NPOINTS_PML,USE_PML_XMIN,USE_PML_XMAX,USE_PML_YMIN,USE_PML_YMAX,2)
open(unit=20,file='energy.dat',status='unknown')
do it2 = 1,NSTEP
- write(20,*) sngl(dble(it2-1)*DELTAT),total_energy_kinetic(it2), &
- total_energy_potential(it2),total_energy_kinetic(it2) + total_energy_potential(it2)
+ write(20,*) sngl(dble(it2-1)*DELTAT),sngl(total_energy_kinetic(it2)), &
+ sngl(total_energy_potential(it2)),sngl(total_energy_kinetic(it2) + total_energy_potential(it2))
enddo
close(20)
call write_seismograms(sisvx,sisvy,NSTEP,NREC,DELTAT)
@@ -1071,8 +1071,8 @@
! save total energy
open(unit=20,file='energy.dat',status='unknown')
do it = 1,NSTEP
- write(20,*) sngl(dble(it-1)*DELTAT),total_energy_kinetic(it), &
- total_energy_potential(it),total_energy_kinetic(it) + total_energy_potential(it)
+ write(20,*) sngl(dble(it-1)*DELTAT),sngl(total_energy_kinetic(it)), &
+ sngl(total_energy_potential(it)),sngl(total_energy_kinetic(it) + total_energy_potential(it))
enddo
close(20)
@@ -1086,8 +1086,8 @@
write(20,*)
write(20,*) 'set output "cpml_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1101,8 +1101,8 @@
write(20,*)
write(20,*) 'set output "compare_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l 1, &
- & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l 2'
+ write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l lc 1, &
+ & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l lc 2'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1117,22 +1117,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_fourth_order.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_fourth_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_fourth_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -750,8 +750,8 @@
! save total energy
open(unit=20,file='energy.dat',status='unknown')
do it = 1,NSTEP
- write(20,*) sngl(dble(it-1)*DELTAT),total_energy_kinetic(it), &
- total_energy_potential(it),total_energy_kinetic(it) + total_energy_potential(it)
+ write(20,*) sngl(dble(it-1)*DELTAT),sngl(total_energy_kinetic(it)), &
+ sngl(total_energy_potential(it)),sngl(total_energy_kinetic(it) + total_energy_potential(it))
enddo
close(20)
@@ -765,8 +765,8 @@
write(20,*)
write(20,*) 'set output "cpml_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -780,8 +780,8 @@
write(20,*)
write(20,*) 'set output "compare_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l 1, &
- & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l 2'
+ write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l lc 1, &
+ & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l lc 2'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -796,22 +796,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_second_order.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_second_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_CPML_2D_isotropic_second_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -749,8 +749,8 @@
! save total energy
open(unit=20,file='energy.dat',status='unknown')
do it = 1,NSTEP
- write(20,*) sngl(dble(it-1)*DELTAT),total_energy_kinetic(it), &
- total_energy_potential(it),total_energy_kinetic(it) + total_energy_potential(it)
+ write(20,*) sngl(dble(it-1)*DELTAT),sngl(total_energy_kinetic(it)), &
+ sngl(total_energy_potential(it)),sngl(total_energy_kinetic(it) + total_energy_potential(it))
enddo
close(20)
@@ -764,8 +764,8 @@
write(20,*)
write(20,*) 'set output "cpml_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -779,8 +779,8 @@
write(20,*)
write(20,*) 'set output "compare_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l 1, &
- & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l 2'
+ write(20,*) 'plot "energy.dat" us 1:4 t ''Total energy CPML'' w l lc 1, &
+ & "../collino/energy.dat" us 1:4 t ''Total energy Collino'' w l lc 2'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -795,22 +795,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_CPML_2D_poroelastic_fourth_order.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_CPML_2D_poroelastic_fourth_order.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_CPML_2D_poroelastic_fourth_order.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1010,8 +1010,8 @@
write(20,*) 'set ylabel "Total energy"'
write(20,*)
write(20,*) '# set output "cpml_total_energy.eps"'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1028,25 +1028,25 @@
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
write(20,*) '#set yrange [-4:4.5]'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
write(20,*) '#set yrange [-15:19]'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
write(20,*) '#set yrange [-12:16]'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
write(20,*) '#set yrange [-7:10]'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_CPML_3D_isotropic_MPI_OpenMP.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_CPML_3D_isotropic_MPI_OpenMP.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_CPML_3D_isotropic_MPI_OpenMP.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1260,7 +1260,7 @@
write(20,*)
write(20,*) 'set output "CPML3D_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" t ''Total energy'' w l 1'
+ write(20,*) 'plot "energy.dat" t ''Total energy'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1275,32 +1275,32 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_001.eps"'
- write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_002.eps"'
- write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_CPML_3D_viscoelastic_MPI.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_CPML_3D_viscoelastic_MPI.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_CPML_3D_viscoelastic_MPI.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -1430,8 +1430,8 @@
! save energy
open(unit=21,file='energy.dat',status='unknown')
do it2=1,NSTEP
- write(21,*) sngl(dble(it2-1)*DELTAT),total_energy_kinetic(it2),&
- total_energy_potential(it2),total_energy(it2)
+ write(21,*) sngl(dble(it2-1)*DELTAT),sngl(total_energy_kinetic(it2)),&
+ sngl(total_energy_potential(it2)),sngl(total_energy(it2))
enddo
close(21)
@@ -1455,11 +1455,6 @@
! save seismograms
call write_seismograms(sisvx,sisvy,NSTEP,NREC,DELTAT,t0)
-! open(unit=20,file='energy.dat',status='unknown')
-! do it = 1,NSTEP
-! write(20,*) sngl(dble(it-1)*DELTAT),total_energy(it)
-! enddo
-! close(20)
! create script for Gnuplot for total energy
open(unit=20,file='plot_energy',status='unknown')
@@ -1471,7 +1466,7 @@
write(20,*)
write(20,*) 'set output "CPML3D_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" t ''Total energy'' w l 1'
+ write(20,*) 'plot "energy.dat" t ''Total energy'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -1486,32 +1481,32 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_001.eps"'
- write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_002.eps"'
- write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_PML_Collino_2D_anisotropic_fourth.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_PML_Collino_2D_anisotropic_fourth.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_PML_Collino_2D_anisotropic_fourth.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -646,8 +646,8 @@
write(20,*)
write(20,*) 'set output "collino_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -662,22 +662,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_PML_Collino_2D_isotropic.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_PML_Collino_2D_isotropic.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_PML_Collino_2D_isotropic.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -578,8 +578,8 @@
write(20,*)
write(20,*) 'set output "collino_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l 1, "energy.dat" us 1:3 &
- & t ''Ep'' w l 3, "energy.dat" us 1:4 t ''Total energy'' w l 4'
+ write(20,*) 'plot "energy.dat" us 1:2 t ''Ec'' w l lc 1, "energy.dat" us 1:3 &
+ & t ''Ep'' w l lc 3, "energy.dat" us 1:4 t ''Total energy'' w l lc 4'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -594,22 +594,22 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
Modified: seismo/3D/CPML/trunk/seismic_PML_Collino_3D_isotropic_OpenMP.f90
===================================================================
--- seismo/3D/CPML/trunk/seismic_PML_Collino_3D_isotropic_OpenMP.f90 2011-06-09 19:40:21 UTC (rev 18570)
+++ seismo/3D/CPML/trunk/seismic_PML_Collino_3D_isotropic_OpenMP.f90 2011-06-09 22:55:55 UTC (rev 18571)
@@ -896,7 +896,7 @@
write(20,*)
write(20,*) 'set output "collino3D_total_energy_semilog.eps"'
write(20,*) 'set logscale y'
- write(20,*) 'plot "energy.dat" t ''Total energy'' w l 1'
+ write(20,*) 'plot "energy.dat" t ''Total energy'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
close(20)
@@ -911,32 +911,32 @@
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_001.eps"'
- write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_001.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_001.eps"'
- write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_001.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_001.eps"'
- write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_001.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vx_receiver_002.eps"'
- write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l 1'
+ write(20,*) 'plot "Vx_file_002.dat" t ''Vx C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vy_receiver_002.eps"'
- write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l 1'
+ write(20,*) 'plot "Vy_file_002.dat" t ''Vy C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
write(20,*) 'set output "v_sigma_Vz_receiver_002.eps"'
- write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l 1'
+ write(20,*) 'plot "Vz_file_002.dat" t ''Vz C-PML'' w l lc 1'
write(20,*) 'pause -1 "Hit any key..."'
write(20,*)
More information about the CIG-COMMITS
mailing list