如何控制wrfout变量输出维度是经纬度而不是south_north和west_east

15104562864

大家好：
本人菜鸟，我只能用最直接的方式表达我的问题，非专业之处请大家见谅，也拜托大家不要打击我。

1）我的wrfout文件中所有变量维度都是(Time, south_north, west_east)，我最近在学NCL，系统总是不能识别这种坐标（错误信息见图1）。后面我发现所有的官网脚本中变量都是（Time, lat, lon）degrees_east、degrees_north的形式。咨询如何将(Time, south_north, west_east)变成（Time, lat, lon）？可以控制wrfout变量维度的输出么？
2）如何制作自己的目标nc文件，举例子，比如我想提取压力P（一维）、XLAT中的纬度（一维）、XLONG中的经度（一维），合成新的二维nc文件P（纬度，经度），保留时间（一维）、纬度（一维）、经度（一维）三个单独变量，形成一个具有四个变量的nc文件（见图二）。
3）我的XLAT变量也很诡异，坐标轴的单位是km，官网案例中文件的单位都是度（见图三）。

因在qq群中提问讽刺挖苦的比较多，完事还不给说怎么解决，还有就是让我多看教程，大家也可以给我推荐一些相关的教程。帮帮忙吧！

croton

(1)楼主需要对sst进行处理：
sst@lat2d = lat
sst@lon2d = lon
其中的lat,lon就是wrfout输出的XLAT,XLONG。
（2）写入nc文件，可详细看一下addfile函数。
（3）没有看清楼主描述问题的图片，所以不能轻易给出答案。

15104562864

croton 发表于 2019-1-3 17:55
(1)楼主需要对sst进行处理：
sst@lat2d = lat
sst@lon2d = lon

嗯嗯，我大概明白意思，先试下，再问问题，感谢！

chickai1109

ncl有专门针对wrfout的绘图函数啊，建议楼主用处理wrfout的函数绘图

xhb3256789

wrfout输出格点都是有经纬度坐标的，是通过km换算的，获取经纬度坐标可以用wrf_user_getvar(a,"lon",0)得到。你可能理解错了south_north和west_east,这只是度数的单位，实际上就是经纬度。

15104562864

xhb3256789 发表于 2019-1-4 17:20
wrfout输出格点都是有经纬度坐标的，是通过km换算的，获取经纬度坐标可以用wrf_user_getvar(a,"lon",0)得到 ...

刚试了一下，真的是我想要的结果。但是有了一维经纬度，怎么和变量组合起来呢，变成新的nc文件？

xhb3256789

15104562864 发表于 2019-1-4 20:04
刚试了一下，真的是我想要的结果。但是有了一维经纬度，怎么和变量组合起来呢，变成新的nc文件？

https://blog.csdn.net/weixin_42762673/article/details/81198436   参考这篇文章，因为你直接从wrfout提取，所以前面的定义经纬度属性就不需要了，直接参考后半部分就好

15104562864

xhb3256789 发表于 2019-1-5 13:52
https://blog.csdn.net/weixin_42762673/article/details/81198436   参考这篇文章，因为你直接从wrfout ...

我的天 好多