- 积分
- 68
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-12-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 团团 于 2017-8-26 12:57 编辑
用ncl画全国160站的数据,代码如图,错误如图
begin
fr="e:/cyg/cygdrive/tsummer_d.txt"
r=asciiread(fr, 160, "float")
;printVarSummary(r)
;flat="e:/cyg/cygdrive/lat.txt"
;lat=asciiread(flat, 41, "float")
;flon="e:/cyg/cygdrive/lon.txt"
;lon=asciiread(flon, 71, "float")
flatlon="e:/cyg/cygdrive/lat_lon.txt"
latlon=asciiread(flatlon, (/160,2/), "float")
zlat=latlon(:,0)
zlon=latlon(:,1)
;print(zlat)
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"
olon=new(71, "float")
olat=new(41, "float")
do i=0,71
olon(i)=70+i
end do
do l=0,41
olat(l)=15+l
end do
olat!0="zlat"
olat@long_name="zlat"
olat@units="degrees_north"
olat&lat=olat
olon!0="zlon"
olon@long_name="zlon"
olon@units="degrees_east"
olon&lon=olon
data1=new((/41,71/), "float")
;cha=new((/41,71/), float)
r@_FillValue=9999
rscan=(/10,5,3/)
data1=obj_anal_ic_deprecated(zlon, zlat, r, olon ,olat, rscan, False)
|
-
|