- 积分
- 7026
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-7-29
- 最后登录
- 1970-1-1
|
NCL
系统平台: |
|
问题截图: |
- |
问题概况: |
批量读入4个fnl文件,都能正常读入,我要查找500hpa上位势高度最小值所在的经纬度,结果却只输出第一个文件的计算结果(52,117),怎么没有输出四组数据? |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
7 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
diri = "/home/liluck/ncl/data/"
fils = systemfunc ("ls "+diri+"*.grib")
f= addfiles (fils, "r")
;print(f)
ListSetType (f, "cat")
hgt=f[:]->HGT_3_ISBL({500},{35:60},{115:145})
dims = dimsizes(hgt)
x1d = ndtooned(hgt) ; convert 2D array to 1D for use in minind
inds = ind_resolve(minind (x1d), dims) ; convert 1D array back to 2D
ilat = inds(0,0)
ilon = inds(0,1)
lat_min = hgt&lat_3(ilat)
lon_min = hgt&lon_3(ilon)
print("Minimum value located at "+lat_min+", "+lon_min)
end
|
|