ncl如何添加线shp文件_NCL绘图示例（六）中国区域内的站点（客观分析）图.PDF...

NCL绘图示例(六)中国区域内的站点(客观分析)图.PDF

NCL 绘图示例(六)：

中国区域内的站点(客观分析)图

施宁

(南京信息工程大学 大气科学学院)

begin

;--1, 读入2016 年7 月中国839 站逐日地表气温

temp=asciiread("./temp-839-201607.txt",-1,"integer") ;全部读入为1 维数组

var = reshape(temp,(/839,31,13/)) ;转换为 839 站X31 天X13 列

var@_FillValue = 32766

val_mon = dim_avg_n(var(:,:,7),1)/10. ;计算839 站的月平均气温

lat_stat = var(:,0,1)/100. ;每个站点的纬度

lon_stat = var(:,0,2)/100. ;每个站点的经度

;--2 进行客观分析，则将台站资料插值至等经纬度网格网格点上

delta = 0.5 ; 网格距

lon_s = 70

lon_e = 140

lat_s = 15

lat_e = 55

nlat = ceil((lat_e-lat_s)/delta) + 1 ;29

nlon = ceil((lon_e-lon_s)/delta) + 1 ;21

glat = fspan(lat_s,lat_e,toint(nlat))

glon = fspan(lon_s,lon_e,toint(nlon))

;-设定经纬度的单位

<