From a620e440b47ae27bb3be69ff3caf85b608c2dc9e Mon Sep 17 00:00:00 2001 From: Frantisek Gallovic Date: Wed, 29 May 2019 13:11:14 +0200 Subject: [PATCH] Fix to ASTF calculation in Napa example --- examples/Napa/apparentSTFs/RIKstf-apparentSTFs.sh | 8 ++++++++ examples/Napa/apparentSTFs/RIKstf-apperentSTFs.sh | 12 ------------ examples/Napa/apparentSTFs/seissimul.f90 | 4 ++-- 3 files changed, 10 insertions(+), 14 deletions(-) create mode 100644 examples/Napa/apparentSTFs/RIKstf-apparentSTFs.sh delete mode 100644 examples/Napa/apparentSTFs/RIKstf-apperentSTFs.sh diff --git a/examples/Napa/apparentSTFs/RIKstf-apparentSTFs.sh b/examples/Napa/apparentSTFs/RIKstf-apparentSTFs.sh new file mode 100644 index 0000000..ef5382f --- /dev/null +++ b/examples/Napa/apparentSTFs/RIKstf-apparentSTFs.sh @@ -0,0 +1,8 @@ +#./RIKsrf2 + +ifort -openmp -oseissimul seissimul.f90 +./seissimul +python stfspec.py +python seismogramsspec.py +gnuplot specplot.gp +gnuplot timespecplot.gp diff --git a/examples/Napa/apparentSTFs/RIKstf-apperentSTFs.sh b/examples/Napa/apparentSTFs/RIKstf-apperentSTFs.sh deleted file mode 100644 index 967acbb..0000000 --- a/examples/Napa/apparentSTFs/RIKstf-apperentSTFs.sh +++ /dev/null @@ -1,12 +0,0 @@ -ifort -openmp -O3 -oRIKsrfRandomruptvel RIKsrfRandomruptvel.f90 -gcc -C -Wall -DNO_IEEE_INFINITY -O3 -c Time_2d.c -ifort -C -oRIKrandomruptvel RIKrandomruptvel.f90 Time_2d.o - -ifort -openmp -oseissimul seissimul.f90 - -./RIKstf -./seissimul -python stfspec.py -python seismogramsspec.py -gnuplot specplot.gp -gnuplot timespecplot.gp diff --git a/examples/Napa/apparentSTFs/seissimul.f90 b/examples/Napa/apparentSTFs/seissimul.f90 index 986119e..2f81f14 100644 --- a/examples/Napa/apparentSTFs/seissimul.f90 +++ b/examples/Napa/apparentSTFs/seissimul.f90 @@ -34,7 +34,7 @@ PROGRAM seissimul REAL dum INTEGER i,j,k,tshift - open(101,FILE='RIKstf.in') + open(101,FILE='RIKsrf.in') read(101,*) read(101,*)L, W read(101,*) @@ -109,7 +109,7 @@ PROGRAM seissimul open(201,FILE='sr.dat') do i=1,NSR do j=1,NT - read(201,*)dum,sr(j,i) + read(201,*)sr(j,i) enddo read(201,*) read(201,*)