爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 4845|回复: 0

NCL批量读入4个文件后,怎么只输出一个文件的计算结果

[复制链接]
发表于 2014-12-2 18:31:28 | 显示全部楼层 |阅读模式
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

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表