ncl mask中国区域的数据
这里用到中国地图的文件:bou2_4p.shp,shapefile_mask_data.ncl etc.
链接:https://pan.baidu.com/s/1O9xr1eIgj80bduxRftPtnw
提取码:u9yn
下面是ncl 脚本
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/gsn_csm.ncl”
load “$ NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl”
load “./China/shapefile_mask_data.ncl”
begin
;—Generate some dummy data over map
dir ="/China/"
dir1=“China/chinamap/”
f1=addfile(dir+“Tempature.nc”,“r”)
lat=f1->lat
lon=f1->lon
file1=“mask.nc”
system("rm -rf "+ dir +file1)
fout =addfile(dir+file1, “c”)
nlat=dimsizes(lat)
nlon=dimsizes(lon)
data=generate_2d_array(10, 10, 0., 100., 0, (/nlat,nlon/))
printVarSummary(data)
data!0 = “lat”
data!1 = “lon”
data&lat = lat
data&lon = lon
shp_filename = dir1 + “bou2_4p.shp”
data_mask = shapefile_mask_data(data,shp_filename,True)
fout->data_mask=data_mask
end
挺实用的code,祝 好运
文件失效可留言获取