MeteoInfoLab如何把GSMaP网站获取的GrADS数据转换为netCDF数据？

hlli

GSMaP网站获取的0.1度日降雨数据为GrADS数据，打算用MeteoInfoLab将其转换为netCDF数据，可是报错。请问如何解决？谢谢！

run script...
Loading classes...
F:\data\GSMap(dat)\2016\201601\
sdate: 2016-01-01 00:00:00
F:\data\GSMap(dat)\2016\201601\gsmap_nrt.20160101.0.1d.daily.00Z-23Z.ctl
Open GrADS data file...
Create output NC file: F:\test\GSMap\result\201601\gsmap_nrt.20160101.0.1d.daily.00Z-23Z.nc
Add dimensions...
Add global attributes...
Write x variable...
Write y variable...
Write time variable...
Write data variable...
Time: 1
Traceback (most recent call last):
  File "<iostream>", line 202, in <module>
    at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.getFilePath_Template(GrADSDataInfo.java:1241)
    at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.readGrADSData_Grid_LonLat(GrADSDataInfo.java:1429)
    at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.getGridData_LonLat(GrADSDataInfo.java:1395)
    at org.meteoinfo.data.meteodata.MeteoDataInfo.getGridData(MeteoDataInfo.java:914)
    at org.meteoinfo.data.meteodata.MeteoDataInfo.getGridData(MeteoDataInfo.java:894)
    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)

java.lang.NullPointerException: java.lang.NullPointerException

控制文件内容如下：
DSET ^gsmap_nrt.%y4%m2%d2.0.1d.daily.00Z-23Z.dat
TITLE GSMaP_NRT 0.1deg Daily (00:00Z-23:59Z)
UNDEF -999.9
OPTIONS yrev little_endian TEMPLATE
XDEF   3600 LINEAR  0.05 0.1
YDEF   1200  LINEAR -59.95 0.1
zdef 1 levels 1
tdef 3650 linear 00:00z1dec2007 1dy
VARS 1
precip          0 99 daily averaged precip(mm/hr)
ENDVARS

hlli

用的脚本是王老师2015年的MeteoInfo脚本示例：GrDAS to netCDF

MeteoInfo

参照此贴：http://bbs.06climate.com/forum.p ... &extra=page%3D1

hlli

MeteoInfo 发表于 2018-8-29 12:50
参照此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=43980&extra=page%3D1

run script...
Loading classes...
F:\data\GSMap(dat)\2016\201601\
sdate: 2016-01-01 00:00:00
F:\data\GSMap(dat)\2016\201601\gsmap_nrt.20160101.0.1d.daily.00Z-23Z.ctl
Open GrADS data file...
Variable: precip
2016-01-01 00:00
Variable: precip
Traceback (most recent call last):
  File "<iostream>", line 165, in <module>
  File "C:\MeteoInfo\pylib\mipylib\dimvariable.py", line 155, in __getitem__
    rr = self.dataset.read(self.name, origin, size, stride).reduce()
  File "C:\MeteoInfo\pylib\mipylib\dimdatafile.py", line 86, in read
    return self.dataset.read(varname, origin, size, stride)
        at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.getFilePath_Template(GrADSDataInfo.java:1241)
        at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.readXY(GrADSDataInfo.java:1333)
        at org.meteoinfo.data.meteodata.grads.GrADSDataInfo.read(GrADSDataInfo.java:1310)
        at org.meteoinfo.data.meteodata.MeteoDataInfo.read(MeteoDataInfo.java:829)
        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)

java.lang.NullPointerException: java.lang.NullPointerException

还是报错，MeteoInfo用的1.3.6版，是否需要升级到新版本才能正确运行？

MeteoInfo

hlli 发表于 2018-8-29 15:14
run script...
Loading classes...
F:\data\GSMap(dat)\2016\201601\

用最新版本试试

hlli

MeteoInfo 发表于 2018-8-29 16:41
用最新版本试试

下载网址打不开了，

MeteoInfo

hlli 发表于 2018-8-30 11:08
下载网址打不开了，

可以加入MeteoInfo QQ群，里面有。

hlli

好的，谢谢！