SuperDARN/pydarn

Discussion: Southern hemisphere plot look direction

Closed this issue · 0 comments

carleyjmartin commented

Discussion

Topic: Should the look direction for southern hemisphere geomagnetic plots be as if you are looking down through the Earth or as if you were looking upwards from below the south pole.

Description

So currently the magnetic spatial plots all have the look direction as if you're looking down through the Earth to the southern hemisphere, but the geographic plots are set up as if you're looking upwards from below the south pole. I've had a few people confused at first about why antarctica is sometimes backwards.

Do we want to:
a. keep everything as it is
b. make a standard look direction that both projections use (from above respective pole)
c. give the user an option to do either (I know sometimes it's nice for conjugate studies)

Category

  • User-friendly
  • Software Design
  • Data related
  • Capabilities
  • Clarity
  • Workflow
import matplotlib.pyplot as plt 
import datetime as dt
import pydarn

#grid_file = "/Users/carley/Documents/data/grids/20011220.8000.grid"
grid_file_s = "/Users/carley/Documents/data/grids/20150308.south.grid2"
#grid_file = "/Users/carley/Desktop/test/20200101.n.grd"
#grid_file = "/Users/carley/Desktop/sas_single.grd"
#grid_file = "/Users/carley/Documents/data/fitacf3/kod.grd"
#grid_file_s = "/Users/carley/Documents/data/grids/20210101.0001.00.zho.grd"
#grid_file_s = "/Users/carley/Documents/data/grids/20210101.0400.00.jme.grd"
grid_file = "/Users/carley/Documents/data/grids/20150301.2001.00.inv.grid"
grid_data = pydarn.SuperDARNRead(grid_file).read_grid()
grid_data_s = pydarn.SuperDARNRead(grid_file_s).read_grid()

pydarn.Grid.plot_grid(grid_data, record=30, coastline=True, lowlat=50,
                      projs=pydarn.Projs.GEO, coords=pydarn.Coords.GEOGRAPHIC)

plt.show()

pydarn.Grid.plot_grid(grid_data, record=30, coastline=True, lowlat=50)

plt.show()


pydarn.Grid.plot_grid(grid_data_s, record=10, coastline=True, lowlat=50,
                      projs=pydarn.Projs.GEO, coords=pydarn.Coords.GEOGRAPHIC)

plt.show()

pydarn.Grid.plot_grid(grid_data_s, record=10, coastline=True, lowlat=50)

plt.show()