def METHOD_NAME():
    # Arrange
    settings = Settings(n_sd=10, dt=1 * si.s, n_substep=5)
    simulation = Simulation(settings)
    # Act
    output = simulation.run()
    # Assert
    np.testing.assert_allclose(output["z"][-1], (1.2 + 0.1) * si.km, rtol=0.005)
    np.testing.assert_allclose(
        output["liquid water mixing ratio"][-1], 2.17, rtol=0.02
    )