Plotting Spatial Map of 2-m Temperature

This example plots the spatial distribution of 2-m temperature from UFS WM outputs.

import Ngl, Nio
import pandas as pd
import numpy as np
import xarray as xr
import netCDF4 as nc
from netCDF4 import Dataset
import matplotlib.pyplot as plt

# Read in NetCDF data
fv3=nc.MFDataset('GFSv16veta/GFSPRS.GrbF96.nc')
tmp2m=fv3["TMP_2maboveground"][0,:,:]
tmp2mF=(tmp2m-273.15)*9/5+32
lat=fv3["latitude"]
lon=fv3["longitude"]

# Read in NetCDF data
wks_type = "png"
wks = Ngl.open_wks(wks_type,"FV3_T2_20200121_1200_GFSv16beta")
cnres                 = Ngl.Resources()

# Contour resources
cnres.cnFillOn        = True
cnres.cnFillPalette   = "NCV_bright_white"
cnres.cnLinesOn       = False
cnres.cnLineLabelsOn  = False
cnres.cnLevelSelectionMode = "ManualLevels"
cnres.cnMinLevelValF = -30
cnres.cnMaxLevelValF = 100
cnres.nglMaximize = True
cnres.tiXAxisString = "Lon"
cnres.tiYAxisString = "Lat"
cnres.lbOrientation   = "horizontal"

# Scalar field resources
cnres.sfXArray        = lon[:]
cnres.sfYArray        = lat[:]

# Map resources
cnres.mpFillOn               = False
cnres.mpGridAndLimbOn = False
cnres.mpLandFillColor        = "Transparent"
cnres.mpOceanFillColor       = "Transparent"
cnres.mpInlandWaterFillColor = "Transparent"
cnres.mpLimitMode        = "LatLon"

# Map region
cnres.mpMaxLonF = -65
cnres.mpMaxLatF = 55
cnres.mpMinLonF = -125
cnres.mpMinLatF = 20
cnres.mpDataBaseVersion     = "MediumRes"
cnres.mpOutlineBoundarySets = "USStates"

# Plot resources
cnres0=cnres
cnres0.tiMainString = "FV3 GFSv16beta: 2m Temp (F) FV3 ~C~Initilized 12Z 17 Jan 2020 valid 12Z 21 Jan 2020"#"{} ({}) at {:.2f} hPa with fast_sat".format(clwmr.long_name,clwmr.units,vpfull)
cnres0.tiMainFontHeightF=0.012
contour = Ngl.contour_map(wks,tmp2mF[:,:],cnres0)

# Draw the plot
Ngl.draw(contour)
Ngl.end()