Example #1
0
def latitude():
    """ latitude Profile Example """

    if show is None:
        raise ImportError("need matplotlib")

    p = ArgumentParser(description="IRI2016 latitude profile")
    p.add_argument("glon", help="geodetic  longitude (degrees)", type=float)
    p.add_argument("-glat",
                   help="geodetic latitude START STOP STEP (degrees)",
                   type=float,
                   nargs=3,
                   default=(-60, 60, 2.0))
    p.add_argument("-alt_km", help="altitude (km)", type=float, default=300.0)
    p.add_argument("-o", "--outfn", help="write data to file")
    P = p.parse_args()

    iono = iri.geoprofile(latrange=P.glat,
                          glon=P.glon,
                          altkm=P.alt_km,
                          time="2004-01-01T17")

    if P.outfn:
        outfn = Path(P.outfn).expanduser()
        print("writing", outfn)
        iono.to_netcdf(outfn)

    piri.latprofile(iono)
    show()
Example #2
0
def main():
    p = ArgumentParser()
    p.add_argument("glon", help="geodetic  longitude (degrees)", type=float)
    p.add_argument("-glat", help="geodetic latitude START STOP STEP (degrees)", type=float, nargs=3, default=(-60, 60, 2.0))
    p.add_argument("-alt_km", help="altitude (km)", type=float, default=300.0)
    p.add_argument("-o", "--outfn", help="write data to file")
    P = p.parse_args()

    iono = iri.geoprofile(latrange=P.glat, glon=P.glon, altkm=P.alt_km, time="2004-01-01T17")

    if P.outfn:
        outfn = Path(P.outfn).expanduser()
        print("writing", outfn)
        iono.to_netcdf(outfn)

    piri.latprofile(iono)
    show()
Example #3
0
def main():
    p = ArgumentParser()
    p.add_argument('glon', help='geodetic  longitude (degrees)',
                   type=float)
    p.add_argument('-glat', help='geodetic latitude START STOP STEP (degrees)',
                   type=float, nargs=3, default=(-60, 60, 2))
    p.add_argument('-alt_km', help='altitude (km)',
                   type=float, default=300)
    p.add_argument('-o', '--outfn', help='write data to file')
    p.add_argument('-q', '--quiet', help='disable plotting', action='store_true')
    P = p.parse_args()

    iono = iri.geoprofile(altkm=P.alt_km,
                          latlim=(P.glat[0], P.glat[1]),
                          dlat=P.glat[2],
                          glon=P.glon,
                          time='2004-01-01T17')

    if not P.quiet:
        piri.latprofile(iono)
        show()
Example #4
0
def main():
    p = ArgumentParser()
    p.add_argument('glon', help='geodetic  longitude (degrees)', type=float)
    p.add_argument('-glat', help='geodetic latitude START STOP STEP (degrees)',
                   type=float, nargs=3, default=(-60, 60, 2.))
    p.add_argument('-alt_km', help='altitude (km)', type=float, default=300.)
    p.add_argument('-o', '--outfn', help='write data to file')
    P = p.parse_args()

    iono = iri.geoprofile(latrange=P.glat,
                          glon=P.glon,
                          altkm=P.alt_km,
                          time='2004-01-01T17')

    if P.outfn:
        outfn = Path(P.outfn).expanduser()
        print('writing', outfn)
        iono.to_netcdf(outfn)

    piri.latprofile(iono)
    show()