brendanjmeade

import matplotlib.pyplot as plt
import cartopy.io.shapereader as shpreader
from shapely.geometry import Polygon, MultiPolygon, box
import matplotlib.patches as mpatches
import matplotlib.collections as mc
from shapely.geometry import LineString, MultiLineString, GeometryCollection


def plot_land(lon_min, lat_min, lon_max, lat_max):
    """Plot filled gray land within the given extent using plt."""

brendanjmeade

brendanjmeade

brendanjmeade

"""
Inertial vs. non-inertial DEMs (gravitational settling)
This file implements a Discrete Element Method (DEM) simulation that
compares two fundamentally different approaches to particle dynamics:

1. An inertial model (F = ma)
Uses Newton's second law with Verlet integration, tracking
accelerations and velocities over time, allowing particles to
accumulate momentum and exhibit dynamic settling behavior with
some kinetic energy dissipation