.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/03-projectiles_and_scaling/3.1-plot_transition_diameter.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_auto_examples_03-projectiles_and_scaling_3.1-plot_transition_diameter.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_auto_examples_03-projectiles_and_scaling_3.1-plot_transition_diameter.py:


The simple-to-complex transition diameter
=========================================

.. rubric:: By David Minton

This example demonstrates how to use the :ref:`Scaling <ug-scaling>` class to calculate the simple-to-complex transition diameter as a function of surface gravity bodies using the Monte Carlo scaling model. Because the scaling model is non-deterministic, we sample each body 100 times to get a distribution of transition diameters. We plot the mean and standard deviation of the transition diameter is plotted for each body, showing results that resemble Fig. 7 of Schenk et al. (2021) [#]_.

References
----------

.. [#] Schenk, P., Castillo-Rogez, J., Otto, K.A., Marchi, S., O’Brien, D., Bland, M., Hughson, K., Schmidt, B., Scully, J., Buczkowski, D., Krohn, K., Hoogenboom, T., Kramer, G., Bray, V., Neesemann, A., Hiesinger, H., Platz, T., De Sanctis, M.C., Schroeder, S., Le Corre, L., McFadden, L., Sykes, M., Raymond, C., Russell, C.T., (2021) Compositional control on impact crater formation on mid-sized planetary bodies: Dawn at Ceres and Vesta, Cassini at Saturn. Icarus 359, 114343. `doi: 10.1016/j.icarus.2021.114343 <https://doi.org/10.1016/j.icarus.2021.114343>`_

.. GENERATED FROM PYTHON SOURCE LINES 15-92



.. image-sg:: /auto_examples/03-projectiles_and_scaling/images/sphx_glr_3.1-plot_transition_diameter_001.png
   :alt: 3.1 plot transition diameter
   :srcset: /auto_examples/03-projectiles_and_scaling/images/sphx_glr_3.1-plot_transition_diameter_001.png
   :class: sphx-glr-single-img





.. code-block:: Python


    import matplotlib.pyplot as plt
    import numpy as np

    from cratermaker import Scaling

    silicate_bodies = ["Vesta", "Moon", "Mercury", "Mars", "Earth", "Venus"]
    ice_bodies = [
        "Mimas",
        "Miranda",
        "Tethys",
        "Dione",
        "Rhea",
        "Ceres",
        "Ganymede",
        "Europa",
    ]

    fig, ax = plt.subplots(figsize=(5, 5))
    ax.set_xscale("log")
    ax.set_yscale("log")
    ax.set_xlabel("Surface Gravity [cm/s²]")
    ax.set_ylabel("Transition Diameter [km]")
    ax.set_xlim((4, 2000))
    ax.set_ylim((0.6, 200))
    color_index = 0
    label_color_map = {}
    nsamples = 100
    bodies = silicate_bodies + ice_bodies
    for body in bodies:
        scaling = Scaling.maker(target=body)
        Dt_list = []
        g_list = []
        for _ in range(nsamples):
            Dt_list.append(scaling.transition_diameter * 1e-3)
            g_list.append(scaling.target.gravity * 100)
            scaling.recompute()

        Dt_arr = np.array(Dt_list)
        g_arr = np.array(g_list)
        g_mean = np.mean(g_arr)
        Dt_mean = np.mean(Dt_arr)
        Dt_std = np.std(Dt_arr)
        if body == "Europa":
            label_offset = 0.6 * Dt_mean
        elif body == "Ceres":
            label_offset = 0.57 * Dt_mean
        elif body == "Rhea":
            label_offset = 0.53 * Dt_mean
        elif body == "Mercury" or body == "Venus":
            label_offset = 0.75 * Dt_mean
        else:
            label_offset = 0.5 * Dt_mean

        if body in silicate_bodies:
            marker = "o"
            va = "bottom"
            y_text = Dt_mean + label_offset
        else:
            marker = "^"
            va = "top"
            y_text = Dt_mean - label_offset

        ax.plot(g_arr, Dt_arr, marker, markersize=2, alpha=0.5, label=f"{body} samples")
        ax.errorbar(
            g_mean,
            Dt_mean,
            yerr=Dt_std,
            fmt=marker,
            color="k",
            capsize=3,
            label=f"{body} mean",
        )
        ax.text(g_mean, y_text, body, fontsize=9, ha="center", va=va)

    plt.tight_layout()
    plt.show()


.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 3.622 seconds)


.. _sphx_glr_download_auto_examples_03-projectiles_and_scaling_3.1-plot_transition_diameter.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: 3.1-plot_transition_diameter.ipynb <3.1-plot_transition_diameter.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: 3.1-plot_transition_diameter.py <3.1-plot_transition_diameter.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: 3.1-plot_transition_diameter.zip <3.1-plot_transition_diameter.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_