MeteoInfo脚本示例：GRIB to ARL

MeteoInfo · 发表于 2015-1-13 15:08:55

llg2002 发表于 2015-1-13 15:06
配上ctl，用aDataInfo = GrADSDataInfo()，aDataInfo.ReadDataInfo(infile)应该能读取数据

嗯，可以不用GRADSDataInfo类，如下：
mydata = MeteoDataInfo()
mydata.openGrADSData(infile)

MeteoInfo · 发表于 2015-1-15 09:55:34

llg2002 发表于 2015-1-13 15:06
配上ctl，用aDataInfo = GrADSDataInfo()，aDataInfo.ReadDataInfo(infile)应该能读取数据

收到了你的Email，解决方案如下：
1、你的数据是GRIB1格式，不需要配ctl和idx文件，MeteoInfo可以直接读GRIB1和GRIB2格式数据。
2、请下载最新的MeteoInfo Java 1.1.9R3版本。
3、可以用MeteoInfo软件打开你的数据，看到所有的变量名和相关说明。示例脚本程序：

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2014-10-24
# Purpose: Convert GRIB 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/test/data"
#---- Set output data file
outFile = os.path.join(dataDir, 'fnl.arl')
#---- Read a GRIB data file
mydata = MeteoDataInfo()
infile = os.path.join(dataDir, 'fnl_20110101_00_00_c')
print infile
mydata.openNetCDFData(infile)
print 'GRIB file has been opened...'
#---- Set output ARL data info
arlDI = ARLDataInfo()
#---- Set variable and level list
gvars = ['Pressure_surface','Temperature_height_above_ground',\
'u-component_of_wind_height_above_ground','v-component_of_wind_height_above_ground',\
'Geopotential_height_isobaric','Temperature_isobaric',\
'u-component_of_wind_isobaric','v-component_of_wind_isobaric','Vertical_velocity_pressure_isobaric',\
'Relative_humidity_isobaric']
avars = ['PRSS','T02M','U10M','V10M','HGTS','TEMP','UWND','VWND','WWND','RELH']
levels = [0,200,300,400,500,600,700,800,\
850,900,925,950,975,1000]
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(), 'FNL1', 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*1.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 GRIB 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/test/data" 
 
#---- Set output data file 
outFile = os.path.join(dataDir, 'fnl.arl') 
 
#---- Read a GRIB data file 
mydata = MeteoDataInfo() 
infile = os.path.join(dataDir, 'fnl_20110101_00_00_c') 
print infile 
mydata.openNetCDFData(infile) 
print 'GRIB file has been opened...' 
 
#---- Set output ARL data info 
arlDI = ARLDataInfo() 
 
#---- Set variable and level list 
gvars = ['Pressure_surface','Temperature_height_above_ground',\ 
    'u-component_of_wind_height_above_ground','v-component_of_wind_height_above_ground',\ 
    'Geopotential_height_isobaric','Temperature_isobaric',\ 
    'u-component_of_wind_isobaric','v-component_of_wind_isobaric','Vertical_velocity_pressure_isobaric',\ 
    'Relative_humidity_isobaric'] 
avars = ['PRSS','T02M','U10M','V10M','HGTS','TEMP','UWND','VWND','WWND','RELH'] 
levels = [0,200,300,400,500,600,700,800,\ 
    850,900,925,950,975,1000] 
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(), 'FNL1', 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*1.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!'

llg2002 · 发表于 2015-1-15 11:08:59

谢谢王老师，gvars的命名方式和ctl的不一致，要用MeteoInfo软件查看变量名

sandytiny · 发表于 2015-3-31 12:43:27

你好，王老师，我的任务是在把GRADS数据转为ARL数据，但是在运行程序过程中（MeteoInfo Java版），出现下面提示，不知是哪里出了问题。烦请王老师指点。
Traceback (most recent call last):
File "<iostream>", line 71, in <module>
at org.meteoinfo.data.meteodata.arl.ARLDataInfo.packARLGridData(ARLDataInfo.java:1505)
at org.meteoinfo.data.meteodata.arl.ARLDataInfo.writeGridData(ARLDataInfo.java:1457)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
at java.lang.reflect.Method.invoke(Unknown Source)

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2014-10-24
# Purpose: Convert GRAD 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/test.arl')
#---- Read a GRADS data file
mydata = MeteoDataInfo()
infile = os.path.join(dataDir, 'grads/gzmm15021412.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 = ['u','v','t','h','Qv','t2m','q2m','u10m','v10m','psl','rain']
avars = ['UWND','VWND','TEMP','HGTS','RELH','T02M','RH2M','U10M','V10M','MSLP','TPP1']
levels = [0,10,20,30,50,70,100,150,200,250,300,400,500,600,700,850,925,1000]
for l in levels:
arlDI.levels.add(l)
if l == 0:
arlDI.LevelVarList.add(['T02M','RH2M','U10M','V10M','MSLP','TPP1'])
else:
arlDI.LevelVarList.add(['UWND','VWND','TEMP','HGTS','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 == 'TPP1':
gData = gData.div(1000)
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!'

复制代码

sandytiny · 发表于 2015-3-31 12:46:48

你好，王老师，我的任务是在把GRADS数据转为ARL数据，但是在运行程序过程中（MeteoInfo Java版），出现下面提示，不知是哪里出了问题。烦请王老师指点。
严重: null
java.io.IOException: Negative seek offset
at java.io.RandomAccessFile.seek(Unknown Source)
at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.readGrADSData_Grid_LonLat(GrADSDataInfo.java:1179)
at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.getGridData_LonLat(GrADSDataInfo.java:1135)
at org.meteoinfo.data.meteodata.MeteoDataInfo.getGridData(MeteoDataInfo.java:691)
at org.meteoinfo.data.meteodata.MeteoDataInfo.getGridData(MeteoDataInfo.java:671)
at sun.reflect.GeneratedMethodAccessor18.invoke(Unknown Source)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
at java.lang.reflect.Method.invoke(Unknown Source)
at org.python.core.PyReflectedFunction.__call__(PyReflectedFunction.java:186)
at org.python.core.PyReflectedFunction.__call__(PyReflectedFunction.java:204)
at org.python.core.PyObject.__call__(PyObject.java:420)
at org.python.core.PyObject.__call__(PyObject.java:424)
at org.python.core.PyMethod.__call__(PyMethod.java:136)
at org.python.pycode._pyx7.f$0(<iostream>:76)
at org.python.pycode._pyx7.call_function(<iostream>)
at org.python.core.PyTableCode.call(PyTableCode.java:165)
at org.python.core.PyCode.call(PyCode.java:18)
at org.python.core.Py.runCode(Py.java:1302)
at org.python.util.PythonInterpreter.execfile(PythonInterpreter.java:235)
at org.python.util.PythonInterpreter.execfile(PythonInterpreter.java:230)
at org.meteoinfo.desktop.forms.FrmTextEditor$22.doInBackground(FrmTextEditor.java:668)
at org.meteoinfo.desktop.forms.FrmTextEditor$22.doInBackground(FrmTextEditor.java:607)
at javax.swing.SwingWorker$1.call(Unknown Source)
at java.util.concurrent.FutureTask.run(Unknown Source)
at javax.swing.SwingWorker.run(Unknown Source)
at java.util.concurrent.ThreadPoolExecutor.runWorker(Unknown Source)
at java.util.concurrent.ThreadPoolExecutor$Worker.run(Unknown Source)
at java.lang.Thread.run(Unknown Source)
Traceback (most recent call last):
File "<iostream>", line 71, in <module>
at org.meteoinfo.data.meteodata.arl.ARLDataInfo.packARLGridData(ARLDataInfo.java:1505)
at org.meteoinfo.data.meteodata.arl.ARLDataInfo.writeGridData(ARLDataInfo.java:1457)
at sun.reflect.GeneratedMethodAccessor23.invoke(Unknown Source)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
at java.lang.reflect.Method.invoke(Unknown Source)
java.lang.NullPointerException: java.lang.NullPointerException

sandytiny · 发表于 2015-3-31 12:50:32

我的ctl文件如附件。

MeteoInfo · 发表于 2015-4-7 13:01:20

sandytiny 发表于 2015-3-31 12:50
我的ctl文件如附件。

如果想让我帮你测试，需要提供ctl和数据文件以及Jython脚本。

a8453414 · 发表于 2015-4-9 11:49:11

{:5_213:}

良辰 · 发表于 2015-6-18 14:37:27

sandytiny 发表于 2015-3-31 12:46
你好，王老师，我的任务是在把GRADS数据转为ARL数据，但是在运行程序过程中（MeteoInfo Java版），出现下面 ...

大神，您的问题解决了吗？能指导一下吗？

良辰 · 发表于 2015-6-24 22:42:38

王老师，我在进行grads转arl时，遇到输出结果hysplt不认识的问题。一开始提示需要地面气压和高度变量，我把输出的avars 修改为它要的名字，hyspit不报这个错误了，但是还是不能使用，想请教一下您，知道问题出在哪吗？hsplit需要的arl格式资料对要素名有要求吗？
另外，我在进行nc转arl时，遇到nc的高度是地形跟随坐标，该作何处理那？输出的要素(hyslit需要使用的)是如何选取的哪? 给王老师添麻烦了，的确这两个问题困扰很久。

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