Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

RTP terminator #409

Open
wants to merge 5 commits into
base: develop
Choose a base branch
from
Open

RTP terminator #409

wants to merge 5 commits into from

Conversation

RemingtonRohel
Copy link
Collaborator

Scope

Addition of terminator plotting (shading the nightside cells) on RTP plots.

Screenshot from 2024-10-11 19-06-10

issue: closes #321

Approval

Number of approvals: 1

Test

matplotlib version: 3.9.15

Try out different options for RTP.plot_range_time(), RTP.plot_coord_time() functions with terminator=True.

E.g.

import dmap
from pydarn import RTP, RangeEstimation, Coords
import matplotlib.pyplot as plt


if __name__ == '__main__':
    data = dmap.read_fitacf('fitacf/20211102.0000.00.rkn.a.fitacf')

    fig, axes = plt.subplots(1, 3, figsize=(10,6))
    RTP.plot_range_time(data, ax=axes[0], beam_num=3, range_estimation=RangeEstimation.RANGE_GATE, terminator=True, groundscatter=True, colorbar=False)
    RTP.plot_range_time(data, ax=axes[1], beam_num=3, range_estimation=RangeEstimation.SLANT_RANGE, terminator=True, groundscatter=True, colorbar=False)
    RTP.plot_range_time(data, ax=axes[2], beam_num=3, range_estimation=RangeEstimation.GSMR, terminator=True, groundscatter=True)
    axes[0].set_ylabel("Range Gate")
    axes[1].set_ylabel("Slant Range")
    axes[2].set_ylabel("GSMR")
    plt.show()
    plt.close()

    fig, axes = plt.subplots(2, 2, figsize=(10, 10))
    RTP.plot_coord_time(data, ax=axes[0, 0], beam_num=3, range_estimation=RangeEstimation.GSMR, coords=Coords.GEOGRAPHIC, terminator=True, latlon='lat')
    RTP.plot_coord_time(data, ax=axes[0, 1], beam_num=3, range_estimation=RangeEstimation.SLANT_RANGE, coords=Coords.GEOGRAPHIC, terminator=True, latlon='lat')
    RTP.plot_coord_time(data, ax=axes[1, 0], beam_num=3, range_estimation=RangeEstimation.GSMR, coords=Coords.AACGM, terminator=True, latlon='lon')
    RTP.plot_coord_time(data, ax=axes[1, 1], beam_num=3, range_estimation=RangeEstimation.SLANT_RANGE, coords=Coords.AACGM, terminator=True, latlon='lon')
    plt.show()
    plt.close()

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

New Feature: Terminator in RTPs
1 participant