NCL绘图示例（六）中国区域内的站点（客观分析）图.PDFVIP

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)) ;-设定经纬度的单位 glon!0 = lon glon@long_name = longitude glon@units = degree_east glat!0 = lat glat@long_name = latitude glat@units = degree_north ;客观分析 rscan=(/5.,2.5,1/) ;不超过10 个数，通常1 至4 个数。数值依次递减，表示每 次迭代分析的半径（单位为度） val_mon := obj_anal_ic_Wrap(lon_stat,lat_stat,val_mon,glon,glat,rscan,False);creanm 插 值;creanm 插值 ;--3，图形文件设置 plot = new(1,graphic) plot2 = new(1,graphic) wks = gsn_open_wks(eps,plot-china-station) gsn_define_colormap(wks,MPL_Spectral) ;--4 ，首先绘制中国区域内的变量 res=True res@gsnDraw = False res@gsnFrame = False res@gsnAddCyclic = False res@cnFillOn = True res@cnLineColor = white res@cnLineThicknessF = 0.5 res@cnFillDrawOrder = PreDraw res@cnLineDrawOrder = PreDraw res@gsnSpreadColorEnd = -6 res@cnLevelSelectionMode = ManualLevels res@cnMinLevelValF = 8 res@cnMaxLevelValF = 32 res@cnLevelSpacingF