MeteoInfo脚本示例：GRIB to ARL

MeteoInfo · 发表于 2015-6-24 23:26:41

良辰发表于 2015-6-24 22:42
王老师，我在进行grads转arl时，遇到输出结果hysplt不认识的问题。一开始提示需要地面气压和高度变量 ...

可以参考1楼示例脚本里的变量。

良辰 · 发表于 2015-6-25 16:12:26

MeteoInfo 发表于 2015-6-24 23:26
可以参考1楼示例脚本里的变量。

我再试试看，非常感谢王老师的帮助，辛苦了

良辰 · 发表于 2015-6-25 16:58:32

本帖最后由良辰于 2015-6-25 16:59 编辑

GRADS转arl后，MI绘图为空，为什么？请王老师指点一下，非常感谢！

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2014-10-24
# Purpose: Convert Grads data to ARL data
# Note: Sample
#-----------------------------------------------------------
from org.meteoinfo.data.meteodata import MeteoDataInfo
from org.meteoinfo.data.meteodata.arl import ARLDataInfo
from org.meteoinfo.data.meteodata.arl import DataLabel
import os
#---- Set directories
dataDir = "D:/Temp"
#---- Set output data file
outFile = os.path.join(dataDir, 'arl/test2.arl')
#---- Read a GRADS data file
mydata = MeteoDataInfo()
infile = os.path.join(dataDir, 'grads/wrfout_d02.ctl')
print infile
mydata.openGrADSData(infile)
print 'GRADS file has been opened...'
#---- Set output ARL data info
arlDI = ARLDataInfo()
#---- Set variable and level list
gvars = ['PSFC','T2','U10M','V10M','Z','TK','UMET','VMET','W','RH']
avars = ['PRSS','T02M','U10M','V10M','HGTS','TEMP','UWND','VWND','WWND','RELH']
levels = [0,1000,925,850,700,500,300,200]
for l in levels:
arlDI.levels.add(l)
if l == 0:
arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M'])
else:
arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND','WWND','RELH'])
#---- Write ARL data file
dataInfo = mydata.getDataInfo()
arlDI.createDataFile(outFile)
arlDI.X = dataInfo.getXDimension().getValues()
arlDI.Y = dataInfo.getYDimension().getValues()
variables = dataInfo.getVariables()
tNum = dataInfo.getTimeNum()
for t in range(0, tNum):
mydata.setTimeIndex(t)
atime = dataInfo.getTimes().get(t)
aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'GRPS', 2)
arlDI.writeIndexRecord(atime, aDH)
lidx = 0
for l in arlDI.levels:
print l
for v in arlDI.LevelVarList[lidx]:
vName = gvars[avars.index(v)]
print vName
if lidx == 0:
mydata.setLevelIndex(lidx)
else:
variable = dataInfo.getVariable(vName)
nlidx = variable.getZDimension().getDimValue().indexOf(l*100.0)
mydata.setLevelIndex(nlidx)
gData = mydata.getGridData(vName)
if v == 'PRSS' or v == 'WWND':
gData = gData.div(100)
aDL = DataLabel(atime)
aDL.setLevel(lidx)
aDL.setVarName(v)
aDL.setGrid(99)
aDL.setForecast(0)
arlDI.writeGridData(aDL, gData)
lidx += 1
arlDI.closeDataFile()
print 'Finished!'

#--------------------------------------------------------         
# Author: Yaqiang Wang                                           
# Date: 2014-10-24                                                
# Purpose: Convert Grads data to ARL data   
# Note: Sample                                                    
#----------------------------------------------------------- 
from org.meteoinfo.data.meteodata import MeteoDataInfo 
from org.meteoinfo.data.meteodata.arl import ARLDataInfo 
from org.meteoinfo.data.meteodata.arl import DataLabel 
import os 
#---- Set directories 
dataDir = "D:/Temp" 
#---- Set output data file 
outFile = os.path.join(dataDir, 'arl/test2.arl') 
#---- Read a GRADS data file 
mydata = MeteoDataInfo() 
infile = os.path.join(dataDir, 'grads/wrfout_d02.ctl') 
print infile 
mydata.openGrADSData(infile) 
print 'GRADS file has been opened...' 
#---- Set output ARL data info 
arlDI = ARLDataInfo() 
#---- Set variable and level list 
gvars = ['PSFC','T2','U10M','V10M','Z','TK','UMET','VMET','W','RH'] 
avars = ['PRSS','T02M','U10M','V10M','HGTS','TEMP','UWND','VWND','WWND','RELH'] 
levels = [0,1000,925,850,700,500,300,200] 
for l in levels: 
  arlDI.levels.add(l) 
  if l == 0: 
    arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M']) 
  else: 
    arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND','WWND','RELH']) 
#---- Write ARL data file 
dataInfo = mydata.getDataInfo() 
arlDI.createDataFile(outFile) 
arlDI.X = dataInfo.getXDimension().getValues() 
arlDI.Y = dataInfo.getYDimension().getValues() 
variables = dataInfo.getVariables() 
tNum = dataInfo.getTimeNum() 
for t in range(0, tNum): 
  mydata.setTimeIndex(t) 
  atime = dataInfo.getTimes().get(t) 
  aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'GRPS', 2) 
  arlDI.writeIndexRecord(atime, aDH) 
  lidx = 0 
  for l in arlDI.levels: 
    print l 
    for v in arlDI.LevelVarList[lidx]: 
      vName = gvars[avars.index(v)] 
      print vName 
      if lidx == 0: 
        mydata.setLevelIndex(lidx) 
      else: 
        variable = dataInfo.getVariable(vName) 
        nlidx = variable.getZDimension().getDimValue().indexOf(l*100.0) 
        mydata.setLevelIndex(nlidx) 
      gData = mydata.getGridData(vName) 
      if v == 'PRSS' or v == 'WWND': 
        gData = gData.div(100)       
        aDL = DataLabel(atime) 
        aDL.setLevel(lidx) 
        aDL.setVarName(v) 
        aDL.setGrid(99) 
        aDL.setForecast(0) 
        arlDI.writeGridData(aDL, gData) 
    lidx += 1 
arlDI.closeDataFile() 
print 'Finished!'

MeteoInfo · 发表于 2015-6-25 17:16:18

良辰发表于 2015-6-25 16:58
GRADS转arl后，MI绘图为空，为什么？请王老师指点一下，非常感谢！

把数据贴出来吧，要不然没办法找问题。

良辰 · 发表于 2015-6-26 12:46:59

MeteoInfo 发表于 2015-6-25 17:16
把数据贴出来吧，要不然没办法找问题。

dset ^wrfout_d02.dat
undef 1.e35
pdef 312 285 lcc 39.511 121.599 1.000 1.000 60. 30. 126.500 3000. 3000.
xdef 1098 linear 119.9 0.01351351
ydef 748 linear 38.4 0.01351351
zdef 7 levels
1000.00000
925.00000
850.00000
700.00000
500.00000
300.00000
200.00000
tdef 85 linear 12z17jun2015 1hr
vars 27
UMET 7 0 U Compoment of wind - rotated (diagnostic)
VMET 7 0 V Component of wind - rotated (diagnostic)
W 7 0 W Component of wind
THETA 7 0 Theta
TK 7 0 Temperature in K
Z 7 0 Height (m)
QVAPOR 7 0 Vapor (kg/kg)
TD 7 0 Dewpoint Temperature in C (diagnostic)
RH 7 0 Relative Humidity (diagnostic)
RAINC 0 0 ACCUMULATED TOTAL CUMULUS PRECIPITATION
RAINNC 0 0 ACCUMULATED TOTAL GRID SCALE PRECIPITATION
slvl 0 0 sea level pressure
T2 0 0 TEMP at 2 M
U10M 0 0 U at 10 M - rotated
V10M 0 0 V at 10 M - rotated
XLAT 0 0 LATITUDE, SOUTH IS NEGATIVE
XLONG 0 0 LONGITUDE, WEST IS NEGATIVE
PSFC 0 0 SURFACE PRESSURE(PA)
RH2 0 0 Relative humidity(%) at 2m (diagnostic)
TD2 0 0 Dewpoint temperature(K) at 2m (diagnostic)
T2MAX 0 0 Maximum temperature(K) at 2m (diagnostic)
T2MIN 0 0 Minimum temperature(K) at 2m (diagnostic)
HIGHCLD 0 0 High cloud fraction[%] (diagnostic)
MIDCLD 0 0 Middle cloud fraction[%] (diagnostic)
LOWCLD 0 0 Low cloud fraction[%] (diagnostic)
TOTCLD 0 0 Total cloud fraction[%] (diagnostic)
VIS 0 0 Visibility in km (diagnostic)
endvars

复制代码

数据太大，上传不了呀。

MeteoInfo · 发表于 2015-6-26 13:53:49

良辰发表于 2015-6-26 12:46
数据太大，上传不了呀。

可以用百度云盘之类的工具共享数据。

良辰 · 发表于 2015-6-26 19:41:31

MeteoInfo 发表于 2015-6-26 13:53
可以用百度云盘之类的工具共享数据。

http://pan.baidu.com/share/link? ... 5&uk=2704002612

MeteoInfo · 发表于 2015-6-27 09:01:56

良辰发表于 2015-6-26 19:41
http://pan.baidu.com/share/link?shareid=2675874425&uk=2704002612

你的模式结果数据不适合做轨迹分析，底层的风场等变量有大片的缺测值，会导致无法计算轨迹。事实上ARL格式的气象数据里没有缺测值的处理，也就是说不允许有缺测值。

良辰 · 发表于 2015-6-27 09:43:17

MeteoInfo 发表于 2015-6-27 09:01
你的模式结果数据不适合做轨迹分析，底层的风场等变量有大片的缺测值，会导致无法计算轨迹。事实上ARL格 ...

非常非常感谢王老师！如果必须用这个做，有什么办法吗？王老师。这个数据是nc转的grads,那么这个nc数据也不适合做轨迹了吗？再次感谢王老师，辛苦辛苦，那么大数据下载都要一阵子的。

良辰 · 发表于 2015-6-27 09:50:24

本帖最后由良辰于 2015-6-27 09:57 编辑

良辰发表于 2015-6-27 09:43
非常非常感谢王老师！如果必须用这个做，有什么办法吗？王老师。这个数据是nc转的grads,那么这个nc数据也 ...

王老师，刚才我查看了一下原始的nc数据，里面数据没有缺失，但是垂直高度上，U\V一致，W比UV多。如果我要用nc转arl，数据为地形追随坐标，levels如何设置？另外aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'GRPS', 2)中的2修改为1就可以了吗？

		自动登录	找回密码
密码			立即注册