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

当未来

WRF3.2输出的华南区域的3km分辨率产品（351×252），元数据3km分辨率，lambert投影，想转成nc格式，但经纬度信息无法对应。
1.因为ctl中有pdef属性，cdo直接转nc     ——失败！
2.中心点（22.878°N，113.323°E）为x方向第176个点，y方向上第126.5个点。x和y方向都为3km间隔，想要写出对应的经纬度数组。这个没考虑出来   ——失败

目前的话，不能重新转wrf输出nc，但是还是想将ctl的dat文件转为nc，不知道各位大神有没有办法？

andrewsoong

wrf输出的不就可以是netcdf格式么？

543267158

wrf输出的就是nc格式。你可以去网上找找fortran读和写netcdf，将wrfout文件中需要的变量信息和经纬度信息提取出来，写成nc格式。

当未来

andrewsoong 发表于 2016-3-23 21:34
wrf输出的不就可以是netcdf格式么？

可以选择grads的，我现在只有grads的这个版本。。想转成nc的

当未来

543267158 发表于 2016-3-23 21:41
wrf输出的就是nc格式。你可以去网上找找fortran读和写netcdf，将wrfout文件中需要的变量信息和经纬度信息提 ...

默认是nc的啊，我现在只有grads的，没机会再跑一遍了

andrewsoong

当未来 发表于 2016-3-23 21:43
可以选择grads的，我现在只有grads的这个版本。。想转成nc的

wrf输出的时候选择nc个输出就好了啊。你是说你只会用grads？

当未来

andrewsoong 发表于 2016-3-23 21:45
wrf输出的时候选择nc个输出就好了啊。你是说你只会用grads？

这是之前别人跑的，没办法再跑一次了。。现在只有grads版本的

andrewsoong

当未来 发表于 2016-3-23 22:03
这是之前别人跑的，没办法再跑一次了。。现在只有grads版本的

那我帮不了你了

543267158

你最终的目的是什么？要转成nc格式的数据？还是用数据画图？画图的话没必要非转成nc格式再画。
ctl中有经纬度信息，xdef一行表示有809个点，起始经度是107.84959，间隔经度是0.01351351；同样ydef是纬度的信息。不过你要画图的话可以直接用grads画啊。画出来边界有缺省值。

兰北

我记得Grads中有个直接输出为nc格式的命令（似乎是writenc，记不清楚了），你可以试一下。