I think it's worth to extend global slice example to better show obsts mask (examples/slcf)/slcf_example_global.py).
For less experienced users it can be helpful:
def main():
sim = fds.Simulation("./fds_multimesh")
# Get the first slice
slc = sim.slices[1]
# Set fill value larger than vmax below
data = slc.to_global(masked=True, fill=50)
# Set colormap
cmap = cm.get_cmap('rainbow')
# Set obsts color
cmap.set_over('white')
# Plot the slice
plt.imshow(data[-1].T,
vmin=0,
vmax=2,
cmap=cmap,
origin="lower")
plt.colorbar()
plt.show()
if __name__ == "__main__":
main()
Hi,
I think it's worth to extend global slice example to better show obsts mask (examples/slcf)/slcf_example_global.py). For less experienced users it can be helpful: