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

星雨

“想要写出对应的经纬度数组”----------这个问题你可以在grads中实现，如果你不是非得要nc格式输出的。你用fwrite命令就可以提数据啊。

兰北

找到了，是sdfwrite

洋咩咩

感觉重新装WRF是比较保险的，说不定比你想其他办法还要快一些

当未来

543267158 发表于 2016-3-24 08:11
你最终的目的是什么？要转成nc格式的数据？还是用数据画图？画图的话没必要非转成nc格式再画。
ctl中有经 ...

最终目的是要用这个数据计算，不是画图；

当未来

兰北 发表于 2016-3-24 08:37
我记得Grads中有个直接输出为nc格式的命令（似乎是writenc，记不清楚了），你可以试一下。

grads的这个ctl已经写了pdef插值，我想要原始的经纬度数据

当未来

星雨 发表于 2016-3-24 08:47
“想要写出对应的经纬度数组”----------这个问题你可以在grads中实现，如果你不是非得要nc格式输出的。你 ...

你好，我的ctl里面有pdef插值了，能否去掉pdef直接输出原始的经纬度呢？因为原始的数据是351*252的，而pdef插值成了809×532的

星雨

当未来 发表于 2016-3-24 09:55
你好，我的ctl里面有pdef插值了，能否去掉pdef直接输出原始的经纬度呢？因为原始的数据是351*252的，而pd ...

你好！你是说直接删掉pdef吗？这个貌似不可用，直接删掉的话位置信息会变。你就把资料正常打开，然后设置时间循环、层次循环，set x 1——*，set y 1——*，写出来就行了。

当未来

星雨 发表于 2016-3-24 10:08
你好！你是说直接删掉pdef吗？这个貌似不可用，直接删掉的话位置信息会变。你就把资料正常打开，然后设置 ...

你好，我明白你的意思。但是grads打开的时候已经经过pdef插值了， x 1 y 1 即为xdef和ydef定义的第一个。
但是我想要原始的那个351×252的经纬度（兰伯特投影的二维经纬度矩阵）

MeteoInfo

当未来 发表于 2016-3-24 11:23
你好，我明白你的意思。但是grads打开的时候已经经过pdef插值了， x 1 y 1 即为xdef和ydef定义的第一个。 ...

可以写一个MeteoInfoLab脚本获得相关经纬度数据。

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. #Plot for test
    
14. 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 /> #Plot for test<br /> scatter(lon[::4,::4], lat[::4,::4], edge=False, size=3, color='b')

当未来

MeteoInfo 发表于 2016-3-24 11:33
可以写一个MeteoInfoLab脚本获得相关经纬度数据。

太棒了，感谢王老师。。这个转好的经纬度如何输出呢？我看只有nc的