def test_declarative_global_gfs(): """Test making a global contour plot using GFS.""" data = xr.open_dataset(get_test_data('GFS_global.nc', as_file_obj=False)) cntr = ContourPlot() cntr.data = data cntr.time = datetime(2021, 1, 30, 12) cntr.field = 'Geopotential_height_isobaric' cntr.level = 300 * units.hPa cntr.contours = np.arange(0, 100000, 120).tolist() cntr.linecolor = 'darkblue' cntr.linewidth = 1 panel = MapPanel() panel.area = [-180, 180, 10, 90] panel.projection = 'ps' panel.layers = ['coastline'] panel.plots = [cntr] pc = PanelContainer() pc.size = (8, 8) pc.panels = [panel] pc.draw() return pc.figure
def test_latlon(): """Test our handling of lat/lon information.""" data = xr.open_dataset(get_test_data('irma_gfs_example.nc', as_file_obj=False)) img = ImagePlot() img.data = data img.field = 'Temperature_isobaric' img.level = 500 * units.hPa img.time = datetime(2017, 9, 5, 15, 0, 0) img.colorbar = None contour = ContourPlot() contour.data = data contour.field = 'Geopotential_height_isobaric' contour.level = img.level contour.time = img.time panel = MapPanel() panel.projection = 'lcc' panel.area = 'us' panel.plots = [img, contour] pc = PanelContainer() pc.panel = panel pc.draw() return pc.figure
def test_latlon(): """Test our handling of lat/lon information.""" data = xr.open_dataset(get_test_data('irma_gfs_example.nc', as_file_obj=False)) img = ImagePlot() img.data = data img.field = 'Temperature_isobaric' img.level = 500 * units.hPa img.time = datetime(2017, 9, 5, 15, 0, 0) contour = ContourPlot() contour.data = data contour.field = 'Geopotential_height_isobaric' contour.level = img.level contour.time = img.time panel = MapPanel() panel.projection = 'lcc' panel.area = 'us' panel.plots = [img, contour] pc = PanelContainer() pc.panel = panel pc.draw() return pc.figure
def test_declarative_contour_cam(): """Test making a contour plot with CAM data.""" data = xr.open_dataset(get_test_data('CAM_test.nc', as_file_obj=False)) contour = ContourPlot() contour.data = data contour.field = 'PN' contour.time = datetime.strptime('2020-11-29 00:00', '%Y-%m-%d %H:%M') contour.level = 1000 * units.hPa contour.linecolor = 'black' contour.contours = list(range(0, 1200, 4)) panel = MapPanel() panel.plots = [contour] panel.layout = (1, 1, 1) panel.layers = ['coastline', 'borders', 'states', 'land'] panel.plots = [contour] pc = PanelContainer() pc.panels = [panel] pc.draw() return pc.figure