ncl 绘制中国区域、青藏高原区域的月平均温度

ncl 绘制中国区域、青藏高原区域的月平均温度

针对月平均温度数据：
Regional.nc

提取中国区域和青藏高原区域的数据：

// ncl代码：
// HOME_PATH : DATA PATH
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shapefile_mask_data.ncl"

begin
dir="$HOME_PATH/"
dir1="$HOME_PATH/tpmap/"
dir2="$HOME_PATH/chinamap/"

file1="Temperature.nc"   ; Output: China  and  TP  regional Total Temperature
system("rm -rf "+ dir +file1)
fout =addfile(dir+file1, "c")

f1=addfile(dir+"Regional.nc","r")

lat=f1->lat;
lon=f1->lon;
data=f1->Temperture

nlat=dimsizes(lat)
nlon=dimsizes(lon)
printVarSummary(data)
; data 为2维  月平均的数据

data!0      = "lat"
data!1      = "lon"
data&lat    = lat
data&lon    = lon

; 青藏高原区域
shp_filename = dir1 + "DBATP_Polygon.shp"
data_tp = shapefile_mask_data(data,shp_filename,True)
data_tp!0      = "lat"
data_tp!1      = "lon"
data_tp&lat    = lat
data_tp&lon    = lon

; 中国区域
shp_filename1 = dir2 + "bou2_4l.shp"
data_China = shapefile_mask_data(data,shp_filename,True)
data_China!0      = "lat"
data_China!1      = "lon"
data_China&lat    = lat
data_China&lon    = lon

fout->time=f1->time
fout->TP_temp=data_tp
fout->China_temp=data_China
fout->Regional=data
end

结果：

青藏高原区域:

最后附上所需文件：
链接：https://pan.baidu.com/s/1B0QQbaH_5cHoOF3skq9lhQ
提取码：rrqq

祝好运！！！