Fix to ASTF calculation in Napa example

elifo · May 29, 2019 · a620e44 · a620e44
1 parent f18aa2e
commit a620e44
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 14 deletions.
diff --git a/examples/Napa/apparentSTFs/RIKstf-apparentSTFs.sh 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
diff --git 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,*)