Since this method is deprecated, how can you plot the same using xarray, cartopy plotly etc:

import warnings
warnings.filterwarnings("ignore")
import matplotlib.pyplot as plt
import numpy as np
import cartopy.crs as ccrs
import cartopy.feature as feat
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
import xarray as xr
import subprocess

# URL of the data file
url = "https://psl.noaa.gov/thredds/fileServer/Datasets/interp_OLR/olr.mon.mean.nc"

# Path to save the downloaded file
output_path = "olr.mon.mean.nc"

# Download the file using curl
subprocess.run(["curl", "-o", output_path, url])

# Open the downloaded file using xarray
ds = xr.open_dataset(output_path)

def align_coords(ds):
    ds.coords['lon'] = (ds.coords['lon'] + 180) % 360 - 180
    ds = ds.sortby(["lat", "lon"])
    return ds

ds = align_coords(ds)
olr = ds['olr'][0]
lat = ds.lat.values
lon = ds.lon.values

Input 1, line 1👈 This is a link! Replace this text with a note about this line...

This is excellent! Thanks for the good description.