- 积分
- 6328
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-11-12
- 最后登录
- 1970-1-1
![[SilenceLiu-] 粉丝数:68 微博数:85 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 liuxiaoyue 于 2017-12-29 20:25 编辑
最近尝试了一下用NCL中system函数调用wget,批量下载中央气象台的天气分析图,上面是下载亚欧500hPa的,把代码分享给大家!~
附件中一个是下天气图的,另一个是下GRAPES产品的,在使用时只需更改时间即可(year, mon, day)
begin
; date
year = "2017"
mon = "12"
day = "29"
; 500 hPa weather map, EA
folder = "http://image.nmc.cn/product" + "/" + year + "/" + mon + "/" + day + "/"+ "WESA/medium/"
fprefix = "SEVP_NMC_WESA_SFER_EGH_ACHN_L50_P9_"
time = year + mon + day + sprinti("%0.2i", ispan(0, 12, 12)) + "0000000.bmp"
URL = folder + fprefix + time
print(URL)
; write URL to file
diri = "./500/EA" + "/" + mon + day + "/"
system("mkdir " + diri)
asciiwrite (diri + "fig.txt", URL)
; download with wget
system("wget -N -P " + diri + " -i " + diri + "fig.txt")
delete(time)
delete(URL)
end
|
|