Datagrid Mapping to a Viewport

import numpy as np
import matplotlib.pyplot as plt
import s3dlib.surface as s3d

#.. Mapping a Datagrid to a Viewport 

# 1. Define function to examine ....................................

Z=np.load('data/jacksboro_fault_dem.npz')['elevation']
datagrid = np.flip(Z,0)

# 2. Setup and map surfaces .........................................
rez=6
vp = [0.15,0.3,0.35,0.70]

fault = s3d.SphericalSurface.grid(250,500, 'r',facecolor='honeydew')
fault.map_cmap_from_datagrid(datagrid,viewport=vp)
fault.shade(direction=[0,1,1],contrast=0.5)

# 3. Construct figure, add surfaces, and plot ......................

fig = plt.figure(figsize=plt.figaspect(1) )
fig.text(0.975,0.975,str(fault), ha='right', va='top', fontsize='smaller', multialignment='right')
ax = plt.axes(projection='3d', aspect='equal')
minmax = (-0.75,0.75)
ax.set(xlim=minmax, ylim=minmax, zlim=minmax)
ax.set_title('Jacksboro Fault')
ax.set_axis_off()
ax.view_init(0,90)

ax.add_collection3d(fault)

fig.tight_layout(pad=2.5)
plt.show()