WRF输出的带ctl的dat文件的经纬度坐标提取

543267158

王老师的不错

MeteoInfo

当未来 发表于 2016-3-24 14:12
太棒了，感谢王老师。。这个转好的经纬度如何输出呢？我看只有nc的

1. #Add file
   
2. fn = 'D:/Temp/grads/wrfout_d02.ctl'
   
3. f = addfile(fn)
   
4. #Get a variable
   
5. v = f['T2']
   
6. #Get X/Y dimension value
   
7. x = v.dimvalue(2)
   
8. y = v.dimvalue(1)
   
9. #Project X/Y to Lon/Lat
   
10. toproj = projinfo()
    
11. xx, yy = meshgrid(x, y)
    
12. lon, lat = project(xx, yy, fromproj=f.proj, toproj=toproj)
    
13. #Write lon/lat
    
14. lon.savegrid(x, y, 'C:/Temp/test/lon.txt')
    
15. lat.savegrid(x, y, 'C:/Temp/test/lat.txt')
    
16. #Plot for test
    
17. scatter(lon[::4,::4], lat[::4,::4], edge=False, size=3, color='b')

#Add file<br /> fn = 'D:/Temp/grads/wrfout_d02.ctl'<br /> f = addfile(fn)<br /> #Get a variable<br /> v = f['T2']<br /> #Get X/Y dimension value<br /> x = v.dimvalue(2)<br /> y = v.dimvalue(1)<br /> #Project X/Y to Lon/Lat<br /> toproj = projinfo()<br /> xx, yy = meshgrid(x, y)<br /> lon, lat = project(xx, yy, fromproj=f.proj, toproj=toproj)<br /> #Write lon/lat<br /> lon.savegrid(x, y, 'C:/Temp/test/lon.txt')<br /> lat.savegrid(x, y, 'C:/Temp/test/lat.txt')<br /> #Plot for test<br /> scatter(lon[::4,::4], lat[::4,::4], edge=False, size=3, color='b')