关于wrfout在ncl上出图的方式

sgfzl

本人最近在用wrfout数据画图的时候遇到一点小疑问，希望家园的大神和版主们不吝赐教。
事情是这样的，我用三种方式画出来的图，虽然都能出图，但结果不尽相同。
法一：
我读取了wrfout变量后，对变量重新设置了coordinates,然后直接用gsn_csm_contour_map函数出图，图片和部分代码如下

                               
登录/注册后可看大图

                               
登录/注册后可看大图

法二：
同样是利用gsn_csm_contour_map函数出图，但这次没有对变量设置coordinates，而是加了res@tfDoNDCOverlay = True这个resouce,图片如下

                               
登录/注册后可看大图

法三：
既没设置coordinates，也没添加res@tfDoNDCOverlay = True，直接利用wrf_contour、wrf_map_overlays出图，图片如下

                               
登录/注册后可看大图

对比这三张图，图二和图三的投影方式虽然不同，但是位置和轮廓是相同的
而图一与图二、图三相比，却不太一样，图右下角的两个大值中心消失了，同时填色区的位置也变化了。
希望各位能帮忙解决解决，万分谢谢！

werewolf

图一的坐标没有设置正确

放逐之刃

投影不同，横纵坐标范围都不一样，当然不同了

freekiller

个人认为1和3是对的，只是用了不同的投影方式。2的问题是你用了wrf的投影数据，强制到地图上，会有位置偏移

sgfzl

放逐之刃 发表于 2020-4-24 09:31
投影不同，横纵坐标范围都不一样，当然不同了

图二我已经设置了105～125E，怎么不起作用呢？

sgfzl

freekiller 发表于 2020-4-24 10:04
个人认为1和3是对的，只是用了不同的投影方式。2的问题是你用了wrf的投影数据，强制到地图上，会有位置偏移

懂了，2应该是错的，2是直接把整个wrf的数据强制到了小的地图上。可是1和3还是有点不一样，到底哪个是对的呢，这放上来的图怎么那么模糊

sgfzl

werewolf 发表于 2020-4-24 09:14
图一的坐标没有设置正确

要怎么设置呢?我放上来的代码怎么显示不出来了。
lat = a->XLAT(0,:,0)
lon = a->XLONG(0,0,:)
lat@units = "degrees_north"
lon@units = "degrees_east"
q!0 = "time"
q!1 = "level"
q!2 = "lat"
q!3 = "lon"
q&lat = lat
q&lon = lon
上面就是我设置坐标的代码，有什么不对的地方，希望不吝赐教

werewolf

xlat = a-> XLAT(0,:,:)
xlon = a-> XLONG(0,:,:)
q@lat2d = xlat
q@lon2d = xlon

sgfzl

werewolf 发表于 2020-4-25 01:29
xlat = a-> XLAT(0,:,:)
xlon = a-> XLONG(0,:,:)
q@lat2d = xlat

这样读入的话，经纬度坐标不就是一个二维的数组吗？经纬度不应该是一维的吗？

werewolf

sgfzl 发表于 2020-4-26 12:18
这样读入的话，经纬度坐标不就是一个二维的数组吗？经纬度不应该是一维的吗？

对，二维也可以，你可以去官网看看相关文档