以下给出利用ncl实现：1980-2010年1月500 hPa 环流平均图，1980年1月高度距平图，1980年1月的环流纬偏图

由于缺少数据文件，我无法直接为您提供完整的代码。但是，我可以给出一个实现的大体思路和一些相关的 NCL 函数，您可以根据自己的数据和需求进行修改和调整。

1. 绘制1980-2010年1月500 hPa 环流平均图

可以使用 NCL 中的 climo 函数计算出时间序列的月平均数据，然后再使用平均数据绘制环流图。具体步骤如下：

; 读取数据
f = addfile("data.nc", "r")
u = f->U(0,:,:,:)   ; 第一个维度为时间，选择第一层和全球范围

; 计算时间序列的月平均数据
u_climo = climo(u, "time", "jan")   ; jan 表示计算一月份的平均值

; 绘制环流图
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

res = True
res@cnFillOn = True
res@cnFillMode = "RasterFill"
res@cnLineLabelsOn = False
res@cnRasterSmoothingOn = True
res@cnInfoLabelOn = False
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = -15
res@cnMaxLevelValF = 15
res@cnLevelSpacingF = 1
res@lbTitleString = "500 hPa 环流平均图 (1980-2010年1月)"
res@gsnDraw = False
res@gsnFrame = False

plot = gsn_csm_contour_map(wks, u_climo, res)

2. 绘制1980年1月高度距平图

可以使用 NCL 中的 climo 函数计算出时间序列的月平均数据和对应的月平均高度场数据，然后再计算出 1980 年 1 月的高度距平，并绘制高度距平图。具体步骤如下：

; 读取数据
f = addfile("data.nc", "r")
u = f->U(0,:,:,:)   ; 第一个维度为时间，选择第一层和全球范围
z = f->Z(0,:,:,:)

; 计算时间序列的月平均数据和月平均高度场数据
u_climo = climo(u, "time", "jan")
z_climo = climo(z, "time", "jan")

; 计算1980年1月的高度距平
z_jan1980 = z_climo(1,:,:,:)
z_jan1980_anom = z_jan1980 - z_climo(0,:,:,:)

; 绘制高度距平图
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

res = True
res@cnFillOn = True
res@cnFillMode = "RasterFill"
res@cnLineLabelsOn = False
res@cnRasterSmoothingOn = True
res@cnInfoLabelOn = False
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = -100
res@cnMaxLevelValF = 100
res@cnLevelSpacingF = 10
res@lbTitleString = "1980年1月高度距平图"
res@gsnDraw = False
res@gsnFrame = False

plot = gsn_csm_contour_map(wks, z_jan1980_anom, res)

3. 绘制1980年1月的环流纬偏图

可以使用 NCL 中的 addfile 函数读取数据文件，然后选择经度为全球范围，纬度为 20-90 度之间的数据。接着，可以使用 NCL 中的 gsn_csm_contour_map 函数绘制环流纬偏图。具体步骤如下：

; 读取数据
f = addfile("data.nc", "r")
v = f->V(0,:,:,20:88)   ; 第一个维度为时间，选择第一层和纬度为 20-88 度之间的数据

; 绘制环流纬偏图
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

res = True
res@cnFillOn = True
res@cnFillMode = "RasterFill"
res@cnLineLabelsOn = False
res@cnRasterSmoothingOn = True
res@cnInfoLabelOn = False
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = -15
res@cnMaxLevelValF = 15
res@cnLevelSpacingF = 1
res@lbTitleString = "1980年1月的环流纬偏图"
res@gsnDraw = False
res@gsnFrame = False

plot = gsn_csm_contour_map(wks, v, res)