def _test(vs30file=None):
print('Testing creation of Sites object with Vs30 file of %s...' % vs30file)
cx = -118.2
cy = 34.1
dx = 0.0083
dy = 0.0083
xspan = 3.0
yspan = 3.0
mysite = Sites.createFromCenter(cx,cy,xspan,yspan,dx,dy,vs30File=vs30file,
padding=True,resample=False)
sc = mysite.getSitesContext()
cx = -118.2
cy = 83
dx = 0.0083
dy = 0.0083
xspan = 3.0
yspan = 3.0
mysite = Sites.createFromCenter(cx,cy,xspan,yspan,dx,dy,vs30File=vs30file,padding=True,resample=False)
xmin = 116.234
xmax = 120.876
ymin = 20.12345
ymax = 24.75435
dx = 0.0083
dy = 0.0083
mysite = Sites.createFromBounds(xmin,xmax,ymin,ymax,dx,dy,vs30File=vs30file,padding=False,resample=False)
print('Passed creation of Sites object with Vs30 file of %s' % vs30file)
def _test(vs30file=None):
print('Testing creation of Sites object with Vs30 file of %s...' %
vs30file)
cx = -118.2
cy = 34.1
dx = 0.0083
dy = 0.0083
xspan = 3.0
yspan = 3.0
mysite = Sites.createFromCenter(cx,
cy,
xspan,
yspan,
dx,
dy,
vs30File=vs30file,
padding=True,
resample=False)
sc = mysite.getSitesContext()
cx = -118.2
cy = 83
dx = 0.0083
dy = 0.0083
xspan = 3.0
yspan = 3.0
mysite = Sites.createFromCenter(cx,
cy,
xspan,
yspan,
dx,
dy,
vs30File=vs30file,
padding=True,
resample=False)
xmin = 116.234
xmax = 120.876
ymin = 20.12345
ymax = 24.75435
dx = 0.0083
dy = 0.0083
mysite = Sites.createFromBounds(xmin,
xmax,
ymin,
ymax,
dx,
dy,
vs30File=vs30file,
padding=False,
resample=False)
print('Passed creation of Sites object with Vs30 file of %s' % vs30file)
def test(vs30file=None):
vs30file = os.path.join(shakedir, 'data/Vs30_test.grd')
cx = -118.2
cy = 34.1
dx = 0.0083
dy = 0.0083
xspan = 0.0083 * 5
yspan = 0.0083 * 5
mysite = Sites.createFromCenter(cx, cy, xspan, yspan, dx, dy,
vs30File=vs30file, padding=True,
resample=False)
grd = mysite.getVs30Grid().getData()
grd_target = np.array(
[[ 409.14199829, 428.24468994, 395.90026855, 426.31658936,
427.53845215, 433.40481567, 436.59918213, 440.43664551,
413.60662842],
[ 423.60855103, 426.64892578, 398.89712524, 428.88549805,
428.78335571, 428.58578491, 430.54354858, 433.59750366,
415.18811035],
[ 391.67398071, 426.20635986, 425.57946777, 428.21954346,
426.06726074, 421.86233521, 423.53192139, 426.25296021,
365.48190308],
[ 385.00396729, 428.14602661, 430.05944824, 429.3427124 ,
426.13626099, 409.76391602, 383.07299805, 372.39117432,
426.48486328],
[ 356.68756104, 432.64077759, 434.55209351, 432.21600342,
395.53771973, 419.31866455, 421.67749023, 426.23449707,
427.32818604],
[ 337.0546875 , 345.14605713, 403.78097534, 385.49118042,
413.04779053, 428.22869873, 427.00268555, 426.8951416 ,
425.62023926],
[ 403.78884888, 336.48217773, 347.82220459, 425.96798706,
432.0640564 , 429.40097046, 427.74179077, 427.00006104,
419.5062561 ],
[ 420.07458496, 330.57504272, 392.33255005, 430.33862305,
432.01391602, 429.43969727, 427.30435181, 425.96151733,
426.15856934]])
np.testing.assert_allclose(grd, grd_target)
sc = mysite.getSitesContext()
scsamp = mysite.sampleFromSites(np.array([34.1, 34.111]),
np.array([-118.2, -118.222]))
xmin = -118.2
xmax = -118.12
ymin = 34.05
ymax = 34.1
dx = 0.0083
dy = 0.0083
mysite = Sites.createFromBounds(xmin, xmax, ymin, ymax, dx, dy,
vs30File=vs30file, padding=False,
resample=False)
grd = mysite.getVs30Grid().getData()
grd_target = np.array(
[[ 426.13626099, 409.76391602, 383.07299805, 372.39117432,
426.48486328, 422.45065308, 417.21685791, 410.43737793,
405.33862305, 400.5239563 , 396.03030396],
[ 395.53771973, 419.31866455, 421.67749023, 426.23449707,
427.32818604, 424.84860229, 415.17590332, 407.5760498 ,
404.54788208, 402.28500366, 398.95913696],
[ 413.04779053, 428.22869873, 427.00268555, 426.8951416 ,
425.62023926, 419.60952759, 411.7961731 , 407.53509521,
406.22122192, 405.31622314, 402.92575073],
[ 432.0640564 , 429.40097046, 427.74179077, 427.00006104,
419.5062561 , 411.10830688, 407.48901367, 406.53305054,
406.59658813, 406.24887085, 406.17044067],
[ 432.01391602, 429.43969727, 427.30435181, 425.96151733,
426.15856934, 427.56121826, 397.67102051, 399.21054077,
404.54968262, 407.18515015, 410.00823975],
[ 426.76913452, 428.86270142, 425.99606323, 423.5692749 ,
423.59835815, 425.92758179, 408.44885254, 406.55810547,
409.06945801, 413.75210571, 417.71130371],
[ 423.1892395 , 427.91104126, 424.53796387, 419.47485352,
418.177948 , 424.14065552, 428.57913208, 432.95300293,
427.77731323, 431.46524048, 442.75540161],
[ 410.56729126, 423.15557861, 424.48355103, 419.27658081,
418.60211182, 423.86721802, 428.06176758, 432.42089844,
438.54446411, 448.37237549, 462.25509644],
[ 405.86877441, 414.27923584, 419.53878784, 396.18017578,
414.54452515, 419.80047607, 425.19824219, 424.48919678,
424.72061157, 431.52270508, 432.84347534]])
np.testing.assert_allclose(grd, grd_target)
def test(vs30file=None):
vs30file = os.path.join(shakedir, 'tests/data/Vs30_test.grd')
cx = -118.2
cy = 34.1
dx = 0.0083
dy = 0.0083
xspan = 0.0083 * 5
yspan = 0.0083 * 5
mysite = Sites.createFromCenter(cx, cy, xspan, yspan, dx, dy,
vs30File=vs30file, padding=True,
resample=False)
grd = mysite.getVs30Grid().getData()
grd_target = np.array(
[[ 409.14199829, 428.24468994, 395.90026855, 426.31658936,
427.53845215, 433.40481567, 436.59918213, 440.43664551,
413.60662842],
[ 423.60855103, 426.64892578, 398.89712524, 428.88549805,
428.78335571, 428.58578491, 430.54354858, 433.59750366,
415.18811035],
[ 391.67398071, 426.20635986, 425.57946777, 428.21954346,
426.06726074, 421.86233521, 423.53192139, 426.25296021,
365.48190308],
[ 385.00396729, 428.14602661, 430.05944824, 429.3427124 ,
426.13626099, 409.76391602, 383.07299805, 372.39117432,
426.48486328],
[ 356.68756104, 432.64077759, 434.55209351, 432.21600342,
395.53771973, 419.31866455, 421.67749023, 426.23449707,
427.32818604],
[ 337.0546875 , 345.14605713, 403.78097534, 385.49118042,
413.04779053, 428.22869873, 427.00268555, 426.8951416 ,
425.62023926],
[ 403.78884888, 336.48217773, 347.82220459, 425.96798706,
432.0640564 , 429.40097046, 427.74179077, 427.00006104,
419.5062561 ],
[ 420.07458496, 330.57504272, 392.33255005, 430.33862305,
432.01391602, 429.43969727, 427.30435181, 425.96151733,
426.15856934]])
np.testing.assert_allclose(grd, grd_target)
sc = mysite.getSitesContext()
scsamp = mysite.sampleFromSites(np.array([34.1, 34.111]),
np.array([-118.2, -118.222]))
xmin = -118.2
xmax = -118.12
ymin = 34.05
ymax = 34.1
dx = 0.0083
dy = 0.0083
mysite = Sites.createFromBounds(xmin, xmax, ymin, ymax, dx, dy,
vs30File=vs30file, padding=False,
resample=False)
grd = mysite.getVs30Grid().getData()
grd_target = np.array(
[[ 426.13626099, 409.76391602, 383.07299805, 372.39117432,
426.48486328, 422.45065308, 417.21685791, 410.43737793,
405.33862305, 400.5239563 , 396.03030396],
[ 395.53771973, 419.31866455, 421.67749023, 426.23449707,
427.32818604, 424.84860229, 415.17590332, 407.5760498 ,
404.54788208, 402.28500366, 398.95913696],
[ 413.04779053, 428.22869873, 427.00268555, 426.8951416 ,
425.62023926, 419.60952759, 411.7961731 , 407.53509521,
406.22122192, 405.31622314, 402.92575073],
[ 432.0640564 , 429.40097046, 427.74179077, 427.00006104,
419.5062561 , 411.10830688, 407.48901367, 406.53305054,
406.59658813, 406.24887085, 406.17044067],
[ 432.01391602, 429.43969727, 427.30435181, 425.96151733,
426.15856934, 427.56121826, 397.67102051, 399.21054077,
404.54968262, 407.18515015, 410.00823975],
[ 426.76913452, 428.86270142, 425.99606323, 423.5692749 ,
423.59835815, 425.92758179, 408.44885254, 406.55810547,
409.06945801, 413.75210571, 417.71130371],
[ 423.1892395 , 427.91104126, 424.53796387, 419.47485352,
418.177948 , 424.14065552, 428.57913208, 432.95300293,
427.77731323, 431.46524048, 442.75540161],
[ 410.56729126, 423.15557861, 424.48355103, 419.27658081,
418.60211182, 423.86721802, 428.06176758, 432.42089844,
438.54446411, 448.37237549, 462.25509644],
[ 405.86877441, 414.27923584, 419.53878784, 396.18017578,
414.54452515, 419.80047607, 425.19824219, 424.48919678,
424.72061157, 431.52270508, 432.84347534]])
np.testing.assert_allclose(grd, grd_target)