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()
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()
def cli():
p = ArgumentParser(description="IRI2016 latitude profile")
p.add_argument("time", help="time of simulation")
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 = main(P.time, P.alt_km, P.glat, P.glon, P.outfn)
try:
from matplotlib.pyplot import show
import iri2016.plots as piri
piri.latprofile(iono)
show()
except ImportError:
pass
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()
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()