brendanjmeade / celeri

Next generation earthquake cycle kinematics
BSD 3-Clause "New" or "Revised" License
24 stars 7 forks source link

Plotting is ignoring x and y limits for Japan model #79

Closed brendanjmeade closed 2 years ago

brendanjmeade commented 2 years ago

Both celeri.plot_input_summary and celeri.plot_estimation_summary are not respecting the lon_range and lat_range tuples that are passed to them. I've verified that they are being passed correctly and have the right values when xlim and ylim are called. Somehow this is ignored for the Japan model but still followed for the Western North America model.

brendanjmeade commented 2 years ago

However, this works fine from within celeri_direct_dense_japan.ipynb with:

lon_range=(130.0, 145.0)
lat_range=(30.0, 50.0)
for i in range(len(segment)):
    if segment.dip[i] == 90.0:
        plt.plot(
            [segment.lon1[i], segment.lon2[i]],
            [segment.lat1[i], segment.lat2[i]],
            "-k",
            linewidth=0.5,
        )
    else:
        plt.plot(
            [segment.lon1[i], segment.lon2[i]],
            [segment.lat1[i], segment.lat2[i]],
            "-r",
            linewidth=0.5,
        )
plt.xlim([lon_range[0], lon_range[1]])
plt.ylim([lat_range[0], lat_range[1]])
plt.gca().set_aspect("equal", adjustable="box")