- 积分
- 641
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-6-24
- 最后登录
- 1970-1-1
|
发表于 2012-10-22 13:06:08
|
显示全部楼层
hzwjy 发表于 2012-10-22 11:51
我不怎么用QQ,有事可以在论坛聊,交流的经验别人也能看到,这样不是更好?
那挺好的。呵呵
do while(date_start .lt. "20101231 23" )
year=systemfunc("date +%Y -d '"+date_start+"'")
month=systemfunc("date +%m -d '"+date_start+"'")
day=systemfunc("date +%d -d '"+date_start+"'")
hour=systemfunc("date +%H -d '"+date_start+"'")
precip=cbinread("FY2E_FDI_ALL_NOM_"+year+""+month+""+day+"_"+hour+"00.rad",(/601,1001/),"double")
res@tiMainString = "GuangXi Radiation 1km Data"
res@gsnLeftString=date_start
res@gsnRightString="W/m~S~2~N~"
plot=gsn_csm_contour_map(wks,precip,res)
draw(plot)
system("rm -f 1km_gx*.png")
frame(wks)
system("mv 1km_gx*.png "+year+month+day+hour+".png")
date_start=systemfunc("date +%Y%m%d' '%H -d ' "+date_start+" 1 hour'")
print(date_start+"")
end do
这个已经可以画图,但是因为有的数据没有,中间就会停止画图,还得人工改时间才能进行画图。有没有一个更好的方法来进行批量画图,希望大神能知道一下,谢谢
这个问题 |
|