Based on this NCL

library(eixport)
library(raster)
#> Loading required package: sp
library(stars)
#> Loading required package: abind
#> Loading required package: sf
#> Linking to GEOS 3.9.0, GDAL 3.2.2, PROJ 7.2.1
library(cptcity)
library(sf)

Reading Temperature

wrfo <- "/media/sergio/ext41/wrfo_sebrazil/wrfout_d02_2014-11-03_000000"
t2 <- wrf_get(wrfo, "T2", as_raster = T)
#> Warning in wrf_get(wrfo, "T2", as_raster = T): The option as_raster will be
#> deprecated, see eixport::wrf_raster()
t2 <- t2$T2_2014.11.03_00.00.00 -273.15# we select one

Find colour palette for temperature

find_cpt("temperature")
#> [1] "arendal_temperature"    "idv_temperature"        "jjg_misc_temperature"  
#> [4] "kst_03_red_temperature"

Based on NCL:

plot(t2, main = "Temperature using plot", col = cpt("arendal_temperature"))
contour(t2, add = T)

#spplot(t2, main = "Temperature 2m using spplot", scales=list(draw = TRUE),
#       col.regions = cpt("idv_temperature"),
#       sp.layout = list("sp.lines", as_Spatial(cl), col = "black"))