请教关于wrfout出来的文件用ncl转换成直线投影的问题

hardshield

在处理wrfout数据时想把兰伯特投影转换成直线投影，查阅资料发现该用rcm2rgrid_Wrap函数，但是对xgrd  = rcm2rgrid_Wrap(lat2d,lon2d,x,lat,lon,0)中的lat和lon不是很理解。官网上似乎是给的全球坐标，而翻阅其他帖子很多前辈给的是自己设定的坐标，因为wrf转出来的数据是某一区域的，这个区域的经纬度范围不知道如何dump出来。想请问一下各位前辈，这个lon和lat需要怎样根据自身的wrfout设置？

淡定冷静

编个程序吧，lat=f->xlat  print(lon(0,1,:,:))类似的，可能稍微麻烦些~

hardshield

淡定冷静 发表于 2020-8-17 10:21
编个程序吧，lat=f->xlat  print(lon(0,1,:,:))类似的，可能稍微麻烦些~

谢谢您的回复，请问如果我想插值成1*1的数据，那么就应该按照最接近于print出来的经纬度的整数来设计插值时的lat，lon么？

淡定冷静

hardshield 发表于 2020-8-17 15:57
谢谢您的回复，请问如果我想插值成1*1的数据，那么就应该按照最接近于print出来的经纬度的整数来设计插值 ...

是的。但一定要取范围区间内的最接近的

Mocca

我前不久试着处理wrfout，画出的是等经纬度投影。读取数据时用了下列的函数：
lat = wrf_user_getvar(f,"lat",it)
lon = wrf_user_getvar(f,"lon",it)
p = wrf_user_getvar(f,"pressure",it)
z = wrf_user_getvar(f, "z",it)     

z_925=  wrf_user_intrp3d(z,p,"h",925,0.,False)   ;select z on 925hPa
z_925@lat2d = lat
z_925@lon2d = lon

后面就是按一维的lat和lon的形式来画啦，不需要用插值。供参考~