.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "generated/gallery/arcade-rtv.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

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

        :ref:`Go to the end <sphx_glr_download_generated_gallery_arcade-rtv.py>`
        to download the full example code.

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

.. _sphx_glr_generated_gallery_arcade-rtv.py:


Modeling an Arcade of Loops with the RTV Scaling Laws
=====================================================

This example shows how to model AIA emission from an arcade
of semi-circular loops who's thermal structure is modeled
using the RTV scaling laws.

.. GENERATED FROM PYTHON SOURCE LINES 9-21

.. code-block:: Python

    import astropy.time
    import astropy.units as u

    from astropy.coordinates import SkyCoord
    from sunpy.coordinates import get_earth

    import synthesizAR

    from synthesizAR.instruments.sdo import InstrumentSDOAIA
    from synthesizAR.interfaces import RTVInterface
    from synthesizAR.models import semi_circular_arcade








.. GENERATED FROM PYTHON SOURCE LINES 22-23

First, set up the coordinates for loops in the arcade.

.. GENERATED FROM PYTHON SOURCE LINES 23-27

.. code-block:: Python

    obstime = astropy.time.Time('2022-11-14T22:00:00')
    pos = SkyCoord(lon=15*u.deg, lat=25*u.deg, radius=1*u.AU, obstime=obstime, frame='heliographic_stonyhurst')
    arcade_coords = semi_circular_arcade(80*u.Mm, 5*u.deg, 50, pos, inclination=10*u.deg)








.. GENERATED FROM PYTHON SOURCE LINES 28-29

Next, assemble the arcade.

.. GENERATED FROM PYTHON SOURCE LINES 29-32

.. code-block:: Python

    strands = [synthesizAR.Strand(f'strand{i}', c) for i, c in enumerate(arcade_coords)]
    arcade = synthesizAR.Skeleton(strands)








.. GENERATED FROM PYTHON SOURCE LINES 33-35

We can make a quick plot of what these coordinates would look like as
viewed from Earth.

.. GENERATED FROM PYTHON SOURCE LINES 35-39

.. code-block:: Python

    earth_observer = get_earth(obstime)
    arcade.peek(observer=earth_observer,
                axes_limits=[(150, 300)*u.arcsec, (275, 425)*u.arcsec])




.. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_001.png
   :alt: arcade rtv
   :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_001.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 40-41

Next, model the thermal structure of each loop using the RTV scaling laws.

.. GENERATED FROM PYTHON SOURCE LINES 41-44

.. code-block:: Python

    rtv = RTVInterface(heating_rate=1e-6*u.Unit('erg cm-3 s-1'))
    arcade.load_loop_simulations(rtv)








.. GENERATED FROM PYTHON SOURCE LINES 45-46

Finally, compute the LOS integrated AIA emission.

.. GENERATED FROM PYTHON SOURCE LINES 46-49

.. code-block:: Python

    aia = InstrumentSDOAIA([0, 1]*u.s, earth_observer, pad_fov=(40, 40)*u.pixel)
    maps = aia.observe(arcade)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Files Downloaded:   0%|          | 0/1 [00:00<?, ?file/s]

    aiapy.aia_V8_all_fullinst.genx:   0%|          | 0.00/5.43M [00:00<?, ?B/s]

    aiapy.aia_V8_all_fullinst.genx:  21%|██        | 1.14M/5.43M [00:00<00:00, 11.4MB/s]

    aiapy.aia_V8_all_fullinst.genx:  92%|█████████▏| 4.98M/5.43M [00:00<00:00, 26.2MB/s]

                                                                                        
    Files Downloaded: 100%|██████████| 1/1 [00:00<00:00,  3.89file/s]
    Files Downloaded: 100%|██████████| 1/1 [00:00<00:00,  3.88file/s]
    INFO: Apparent body location accounts for 493.67 seconds of light travel time [sunpy.coordinates.ephemeris]




.. GENERATED FROM PYTHON SOURCE LINES 50-51

We can make a quick plot of what each EUV channel of AIA would look like.

.. GENERATED FROM PYTHON SOURCE LINES 51-54

.. code-block:: Python

    for k in maps:
        maps[k][0].peek()




.. rst-class:: sphx-glr-horizontal


    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_002.png
         :alt: AIA $94 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_002.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_003.png
         :alt: AIA $131 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_003.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_004.png
         :alt: AIA $171 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_004.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_005.png
         :alt: AIA $193 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_005.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_006.png
         :alt: AIA $211 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_006.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_007.png
         :alt: AIA $335 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_007.png
         :class: sphx-glr-multi-img





.. GENERATED FROM PYTHON SOURCE LINES 55-57

We can easily adjust the viewing angle to integrate the emission along
a different LOS.

.. GENERATED FROM PYTHON SOURCE LINES 57-63

.. code-block:: Python

    off_limb_observer = SkyCoord(
        lon=-70*u.deg, lat=earth_observer.lat, radius=earth_observer.radius, frame=earth_observer.frame)
    aia = InstrumentSDOAIA([0, 1]*u.s, off_limb_observer, pad_fov=(20, 20)*u.pixel,)
    maps = aia.observe(arcade)
    for k in maps:
        maps[k][0].peek(draw_limb=True)



.. rst-class:: sphx-glr-horizontal


    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_008.png
         :alt: AIA $94 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_008.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_009.png
         :alt: AIA $131 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_009.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_010.png
         :alt: AIA $171 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_010.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_011.png
         :alt: AIA $193 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_011.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_012.png
         :alt: AIA $211 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_012.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /generated/gallery/images/sphx_glr_arcade-rtv_013.png
         :alt: AIA $335 \; \mathrm{\mathring{A}}$ 2022-11-14 22:00:00
         :srcset: /generated/gallery/images/sphx_glr_arcade-rtv_013.png
         :class: sphx-glr-multi-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    INFO: Apparent body location accounts for 493.67 seconds of light travel time [sunpy.coordinates.ephemeris]





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

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


.. _sphx_glr_download_generated_gallery_arcade-rtv.py:

.. only:: html

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

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

      :download:`Download Jupyter notebook: arcade-rtv.ipynb <arcade-rtv.ipynb>`

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

      :download:`Download Python source code: arcade-rtv.py <arcade-rtv.py>`

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

      :download:`Download zipped: arcade-rtv.zip <arcade-rtv.zip>`


.. only:: html

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

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