Merge pull request #70 from pmpowers-usgs/gail-table-rscale-69

Gail tbale fix
usgs · Oct 21, 2016 · 5495a8e · 5495a8e
2 parents a5487ed + ca151e3
commit 5495a8e
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 0 deletions.
diff --git a/src/hazFXnga13l.f b/src/hazFXnga13l.f
@@ -8210,6 +8210,8 @@ real function amean11(amag,r,rjb,Vs30,ip,jf,ka)
       if (jdl .eq. 0) write(*,*)' ERROR. CANNOT FIND DIST'
       jdu = jdl + 1
       fracd = (rl - rlog(jdl,ka)) / (rlog(jdu,ka)-rlog(jdl,ka))
+c     fracd goes negative above rmax; clamp to 0 P.Powers 10/21/16
+      if (jdflag .eq. 1) fracd = 0
 c**   fracd gives interpolation fraction in distanc
       do 40 j = jfl, jfu
       arl = gma(j,jdl,jml,ka) + fracm * (gma(j,jdl,jmu,ka)-gma(j,jdl,jml,ka))

diff --git a/src/hazgridXnga13l.f b/src/hazgridXnga13l.f
@@ -8192,6 +8192,8 @@ real function amean11(amag,r,rjb,Vs30,ip,jf,ka)
       if (jdl .eq. 0) write(*,*)' ERROR. CANNOT FIND DIST'
       jdu = jdl + 1
       fracd = (rl - rlog(jdl,ka)) / (rlog(jdu,ka)-rlog(jdl,ka))
+c     fracd goes negative above rmax; clamp to 0 P.Powers 10/21/16
+      if (jdflag .eq. 1) fracd = 0
 c**   fracd gives interpolation fraction in distanc
       do 40 j = jfl, jfu
       arl = gma(j,jdl,jml,ka) + fracm * (gma(j,jdl,jmu,ka)-gma(j,jdl,jml,ka))