search

SuperDARN / pydarn

Python library for visualizing SuperDARN Data
GNU Lesser General Public License v3.0
31 stars 11 forks source link

Discussion: Southern hemisphere plot look direction #362

Closed carleyjmartin closed 3 months ago

carleyjmartin commented 8 months ago

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

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()