关于meteoinfo的一维图

黄小仙儿

最近研究meteoinfo，发现它的功能真的很强大，画图比grads方便，漂亮，能读取的数据也比较全面。但是在摸索的时候发现了一个问题，就是在用netcdf数据画一维图的时候
1、它的高度值是从0到1的，这个还好理解，可能是按照eat层算的？（不知道是不是我想的这样）
但是经纬度的数值却不知道怎么换算了，比如经度的选择范围是从-218182.7143986一直到219594.62560134，纬度也差不多。感觉这数字完全没有规律可循。不知道大家有没有注意到这个问题？
2、在选择w垂直速度时，只有surface的数据。不知道是我的数据有问题还是meteoinfo的bug,不是很懂。

MeteoInfo

1、高度是从数据里读出来的，具体表示什么需要看数据的信息。气象数据并非只是等经纬度的，如果是投影后的数据，格点坐标通常是以米为单位。
2、w到底是几层可以通过查看数据信息来确定，在Meteo Data对话框中打开数据后点击工具栏第二个按钮（Show Data Information）可以查看数据的详细信息。

黄小仙儿

MeteoInfo 发表于 2014-3-27 10:49
1、高度是从数据里读出来的，具体表示什么需要看数据的信息。气象数据并非只是等经纬度的，如果是投影后的 ...

那个投影后的经纬度数据，有没有什么换算公式？我的是兰伯特投影的。
还有刚刚看了下数据信息（再一次感受了meteoinfo的强大功能），关于垂直速度的描述是这样的：
NC_FLOAT W(Time,bottom_top_stag,south_north,west_east);
                W: FieldType = 104;
                W: MemoryOrder = "XYZ";
                W: description = "z-wind component";
                W: units = "m s-1";
                W: stagger = "Z";
                W: coordinates = "XLONG XLAT";
和U,V的表述信息完全一样：但是其他两个就有所有层的数据，W就没有
NC_FLOAT U(Time,bottom_top,south_north,west_east_stag);
                U: FieldType = 104;
                U: MemoryOrder = "XYZ";
                U: description = "x-wind component";
                U: units = "m s-1";
                U: stagger = "X";
                U: coordinates = "XLONG_U XLAT_U";
        NC_FLOAT V(Time,bottom_top,south_north_stag,west_east);
                V: FieldType = 104;
                V: MemoryOrder = "XYZ";
                V: description = "y-wind component";
                V: units = "m s-1";
                V: stagger = "Y";
                V: coordinates = "XLONG_V XLAT_V";

MeteoInfo

黄小仙儿 发表于 2014-3-27 10:59
那个投影后的经纬度数据，有没有什么换算公式？我的是兰伯特投影的。
还有刚刚看了下数据信息（再一次感 ...

wrf out数据不止一个高度维，之前对这种情况没考虑。更新了程序处理这种问题，需要下载MeteoInfo最新文件（MeteoInfo最新文件http://bbs.06climate.com/forum.php?mod=viewthread&tid=1126&fromuid=106
(出处: 气象家园)
）

黄小仙儿

MeteoInfo 发表于 2014-3-27 14:32
wrf out数据不止一个高度维，之前对这种情况没考虑。更新了程序处理这种问题，需要下载MeteoInfo最新文件 ...

哇，王老师你太给力了。
中午估计都没好好吃饭，辛苦了。相信meteoinfo会有越来越多人用的

MeteoInfo

黄小仙儿 发表于 2014-3-27 14:38
哇，王老师你太给力了。
中午估计都没好好吃饭，辛苦了。相信meteoinfo会有越来越多人用的

还是有些问题，wrf out的x, y, z坐标都是两套，比如x坐标有west_east和west_east_stag，后者比前者多一个数。你前面贴出来的U, V变量事实上x, y坐标是不一致的。我有一些问题希望有人能解答：
1、wrf为什么会这样处理？
2、能不能这样理解：west_east坐标指的是x方向格点的中心坐标，west_east_stag指的是x方向格点边界坐标？
3、实际绘图等处理中是不是应该先将stag坐标插值为非stag坐标（貌似NCL是这样处理的）？

前面帖子提到里的U, V变量的信息：
       NC_FLOAT U(Time,bottom_top,south_north,west_east_stag);
                U: FieldType = 104;
                U: MemoryOrder = "XYZ";
                U: description = "x-wind component";
                U: units = "m s-1";
                U: stagger = "X";
                U: coordinates = "XLONG_U XLAT_U";
        NC_FLOAT V(Time,bottom_top,south_north_stag,west_east);
                V: FieldType = 104;
                V: MemoryOrder = "XYZ";
                V: description = "y-wind component";
                V: units = "m s-1";
                V: stagger = "Y";
                V: coordinates = "XLONG_V XLAT_V";

MeteoInfo

基本明白了，这个网页有很好的解释：http://www.openwfm.org/wiki/How_to_interpret_WRF_variables

黄小仙儿

MeteoInfo 发表于 2014-3-27 15:23
基本明白了，这个网页有很好的解释：http://www.openwfm.org/wiki/How_to_interpret_WRF_variables

才看完，又学到了很多东西~

OOOO1

学习中！！！！！