MeteoInfo脚本示例：MICAPS 1类数据转为GrADS站点数据

MeteoInfo

这个脚本演示了如何将MICAPS第1类数据（地面全要素站点观测）转为GrADS站点数据格式（需要MeteoInfo最新文件，见置顶帖子）。

转换后的GrADS站点格式数据绘图：




脚本代码如下：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2012-11-24                                               
   
4. # Purpose: Write GrADS station data  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------
   
7. import clr
   
8. from System import *
   
9. from System.Collections.Generic import *
   
10. clr.AddReference("MeteoInfoC.dll")
    
11. from MeteoInfoC import *
    
12. from MeteoInfoC.Data import *
    
13. from MeteoInfoC.Data.MeteoData import *
    
14. 
    
15. #---- Set directories
    
16. dataDir = "D:\\Temp\"
    
17. 
    
18. #---- Set output data file
    
19. outCtlFile = dataDir + 'grads\\m1test.ctl'
    
20. outDataFile = dataDir + 'grads\\m1test.dat'
    
21. 
    
22. #---- Set time
    
23. sTime = DateTime.Parse("2010-10-14 14:00")
    
24. eTime = DateTime.Parse("2010-10-14 20:00")
    
25. tNum = 2
    
26. 
    
27. #---- Read a micaps data file
    
28. mydata = MeteoDataInfo()
    
29. infile = sTime.ToString("yyMMddHH") + ".000"
    
30. mydata.OpenMICAPSData(dataDir + "micaps\" + infile)
    
31. 
    
32. #---- Define and set output GrADS data info
    
33. gDataInfo = GrADSDataInfo()
    
34. gDataInfo.FileName = outCtlFile
    
35. gDataInfo.DSET = outDataFile
    
36. gDataInfo.TITLE = 'MICAPS data'
    
37. gDataInfo.DTYPE = 'STATION'
    
38. gDataInfo.UNDEF = 9999
    
39. gDataInfo.TDEF.Type = 'LINEAR'
    
40. gDataInfo.TDEF.TNum = tNum
    
41. gDataInfo.TDEF.STime = sTime
    
42. gDataInfo.TDEF.TDelt = '6hr'
    
43. for v in mydata.GetVariableNames():
    
44.         aVar = VAR()
    
45.         aVar.VName = v
    
46.         aVar.LevelNum = 0
    
47.         aVar.Units = '99'
    
48.         aVar.Description = v
    
49.         gDataInfo.VARDEF.AddVar(aVar)
    
50. 
    
51. #---- Write GrADS data file
    
52. gDataInfo.CreateDataFile(outDataFile)
    
53. aTime = sTime
    
54. for t in range(0, tNum):
    
55.         infile = aTime.ToString("yyMMddHH") + ".000"
    
56.         mydata.OpenMICAPSData(dataDir + "micaps\" + infile)
    
57.         stInfoData = mydata.GetStationInfoData()
    
58.         gDataInfo.WriteStationData(stInfoData)
    
59.         aTime = aTime.AddHours(6)
    
60. 
    
61. gDataInfo.CloseDataFile()
    
62. 
    
63. #---- Write GrADS control file
    
64. gDataInfo.WriteGrADSCTLFile()
    
65. 
    
66. print 'Finished!'

复制代码

rhythm9806

用这个脚本出错，求救！！
Traceback (most recent call last):
  File "<string>", line 38, in <module>
AttributeError: 'GrADSDataInfo' object has no attribute 'UNDEF'

MeteoInfo

改为MissingValue

federa123

LZ厉害，ＬＺ加油

风荷1

gDataInfo.VARDEF.AddVar(aVar）这句脚本出错
typeerror:expected variable,got var

MeteoInfo

风荷1 发表于 2015-7-30 15:45
gDataInfo.VARDEF.AddVar(aVar）这句脚本出错
typeerror:expected variable,got var

这是很久以前的脚本了，有些东西改变了。

你尝试把
aVar = VAR()
改为
aVar = Variable()
试试

风荷1

А

风荷1

感谢在百忙之中回复，谢谢

正能量

正能量

求JAVA 版本的代码谢谢了