Notebook to Check NonLinLoc Grids

Daniel Trugman, 2024

from nllgrid import NLLGrid

Verification of the velocity grids and timetables

# Model grid
grdfile = "./model/delph_poly2_lcc.P.mod"

# instantiate
grd = NLLGrid(grdfile)

# show
print(grd)

# plot - notice slow_len units
grd.plot(slice_index='max') #line 1224

basename: ./model/delph_poly2_lcc.P.mod
nx: 249 ny: 249 nz: 105
x_orig: -124.0 y_orig: -124.0 z_orig: -4.0
dx: 1.0 dy: 1.0 dz: 1.0
grid_type: SLOW_LEN
float_type: FLOAT
transform: TRANSFORM  LAMBERT RefEllipsoid WGS-84  LatOrig 48.000000  LongOrig -121.300000  FirstStdParal 46.900000  SecondStdParal 49.100000  RotCW 0.000000

# Time grid
grdfile = "./time/delph_poly2_lcc.P.ALKI.time.buf"

# instantiate
grd = NLLGrid(grdfile)

# show
print(grd)

# plot - notice slow_len units
grd.plot()

basename: ./time/delph_poly2_lcc.P.ALKI.time
nx: 249 ny: 249 nz: 105
x_orig: -124.0 y_orig: -124.0 z_orig: -4.0
dx: 1.0 dy: 1.0 dz: 1.0
grid_type: TIME
float_type: FLOAT
station: ALKI sta_x: -84.086821 sta_y: -46.649995 sta_z: -0.001
transform: TRANSFORM  LAMBERT RefEllipsoid WGS-84  LatOrig 48.000000  LongOrig -121.300000  FirstStdParal 46.900000  SecondStdParal 49.100000  RotCW 0.000000