- 积分
- 352
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-9-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
以下是我的NCL脚本,我用它读取全国1212个站的资料,进行插值,但是发现,插值后的最大值比起站点数据小太多
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"
begin
f=asciiread("E:/program_plot/NCL/NCL_script/data.000", -1,"string")
;print(f)
;if(False)
;;;;;;;;;;;;;;;;;;;;;;;经度纬度起始行;;;;;;;;;;;;;;;;;;;;;
latlon_line_start_num=0
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;经度纬度结束行;;;;;;;;;;;;;;;;;;;;;
latlon_line_end_num=1212
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;经度所在列;;;;;;;;;;;;;;;;;;;;;;;;;
lon_colome_num=1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;维度所在列;;;;;;;;;;;;;;;;;;;;;;;;
lat_colome_num=2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;变量所在列;;;;;;;;;;;;;;;;;;;;;;;;
var_colome_num=3
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
lon= tofloat(str_get_field(f,lon_colome_num," "));获取全国1212个站点经度
lon@units="degrees_east"
lat= tofloat(str_get_field(f,lat_colome_num," "));获取全国1212个站点纬度
lat@units="degrees_north"
var=tofloat(str_get_field(f, var_colome_num," "));获取全国1212个站点气象要素(24h降水量)
olon=fspan(72,137,132) ;我要插值的经度
olon@units="degrees_east" ;
olat=fspan(17,56,80) ;我要插值的维度
olat@units="degrees_north"
rscan=(/10,5,3/) ;插值半径
var_interp=obj_anal_ic_deprecated(lon, lat, var, olon, olat, rscan, False) ;cressman插值
;var_interp=natgrid_Wrap(lon,lat,var,olon,olat)
var_interp@_FillValue = 9.96921e+36
print(max(var));输出站点资料最大值
print(max(var_interp));输出插值后格点资料最大值
;end if
end
输出:
(0) 45.7
(0) 28.86712
为什么插值后的最大值( 45.7)比起站点数据(28.86712) 最大值小太多 ?
|
|