MeteoInfoLab脚本汇总贴

MeteoInfo

ruanjinshuai 发表于 2016-12-6 23:59
王老师，我的数据是2009年GLADS(grib1格式)全年的三小时一次的数据，固有2920个grb文档，我想分析其时间序 ...

你给我发几个数据文件做例子我先试试

ruanjinshuai

我先把数据传上去，这是一天的数据，麻烦王老师了！

ruanjinshuai

王老师，看到你的回复说通过dos命令来合并grib数据http://bbs.06climate.com/forum.p ... id=19271&page=2，我尝试了一下，合并出来的数据大小和原始数据总的大小一致，然而我画图发现同一时刻的T分布是不一致的。以下是我的图片，是我处理错了，还是本来dos这种处理方式有问题？图1为未合成是8.7.00时分布，图2为合成得到数据后的8.7：00时分布，

MeteoInfo

ruanjinshuai 发表于 2016-12-7 13:20
我先把数据传上去，这是一天的数据，麻烦王老师了！

目前MeteoInfoLab还没有类似NCL的addfiles函数同时读入多个数据文件（以后会考虑加上此项功能），可以用循环来处理。

1. datadir = 'E:/Temp'
   
2. st = datetime.datetime(2009,1,1,0)
   
3. et = datetime.datetime(2009,1,1,21)
   
4. tn = 8
   
5. data = None
   
6. i = 0
   
7. while st <= et:
   
8.     fn = os.path.join(datadir, 'GLDAS_NOAH10_3H.A' + st.strftime('%Y%m%d.%H') + \
   
9.         '00.020.grb')
   
10.     print fn
    
11.     f = addfile(fn)
    
12.     ws = f['Wind_Speed_surface'][0,:,:]
    
13.     if i == 0:
    
14.         yn = ws.dimlen(0)
    
15.         xn = ws.dimlen(1)
    
16.         data = zeros((tn,yn,xn))
    
17.     data[i,:,:] = ws.array
    
18.     i += 1
    
19.     st = st + datetime.timedelta(hours=3)
    
20. data = mean(data, axis=0)
    
21. lat = ws.dimvalue(0)
    
22. lon = ws.dimvalue(1)
    
23. #Plot
    
24. axesm()
    
25. lworld = shaperead('D:/Temp/map/country1.shp')
    
26. geoshow(lworld)
    
27. layer = imshowm(lon, lat, data)
    
28. colorbar(layer)
    
29. title('Mean wind speed (' + et.strftime('%Y-%m-%d') + ')')

datadir = 'E:/Temp'<br /> st = datetime.datetime(2009,1,1,0)<br /> et = datetime.datetime(2009,1,1,21)<br /> tn = 8<br /> data = None<br /> i = 0<br /> while st <= et:<br />     fn = os.path.join(datadir, 'GLDAS_NOAH10_3H.A' + st.strftime('%Y%m%d.%H') + \<br />         '00.020.grb')<br />     print fn<br />     f = addfile(fn)<br />     ws = f['Wind_Speed_surface'][0,:,:]<br />     if i == 0:<br />         yn = ws.dimlen(0)<br />         xn = ws.dimlen(1) <br />         data = zeros((tn,yn,xn))<br />     data[i,:,:] = ws.array<br />     i += 1<br />     st = st + datetime.timedelta(hours=3)<br /> data = mean(data, axis=0)<br /> lat = ws.dimvalue(0)<br /> lon = ws.dimvalue(1)<br /> #Plot<br /> axesm()<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld)<br /> layer = imshowm(lon, lat, data)<br /> colorbar(layer)<br /> title('Mean wind speed (' + et.strftime('%Y-%m-%d') + ')')

ruanjinshuai

谢谢王老师！

MeteoInfo

ruanjinshuai 发表于 2016-12-7 14:02
王老师，看到你的回复说通过dos命令来合并grib数据http://bbs.06climate.com/forum.php?mod=viewthread&tid ...

我试了用dos合并grib文件，没问题呀

半颗麦芽糖

王老师您好，我有一个数据，第一列是高度，后面是一天中一些时刻所对应的数据值，我想请问一下能不能画出我下面附上的图？附部分数据和图，请老师指导。
all_.txt (7.79 KB, 下载次数: 5)

2016-12-13 22:02 上传

点击文件名下载附件

MeteoInfo

半颗麦芽糖 发表于 2016-12-13 22:06
王老师您好，我有一个数据，第一列是高度，后面是一天中一些时刻所对应的数据值，我想请问一下能不能画出我 ...

我抽空看看。

有新的问题最好开一个新帖子讨论。

MeteoInfo

半颗麦芽糖 发表于 2016-12-13 22:06
王老师您好，我有一个数据，第一列是高度，后面是一天中一些时刻所对应的数据值，我想请问一下能不能画出我 ...

数据格式有说明吗？看不明白

半颗麦芽糖

MeteoInfo 发表于 2016-12-14 10:32
数据格式有说明吗？看不明白

王老师，数据是自己处理合并的，30、60、90……表示高度，后面那些表示一天内一些时刻对应的数据值，数据有些混乱，我重新整理一下再问王老师吧