wrfout数据如何得到以p为垂直坐标的三维数据

601337653

各位好，请问我想将wrfout数据的垂直坐标转换成气压坐标，该怎么实现呢？我知道有个函数是wrf_user_intrp3d，但是出来的结果貌似是水平或垂直剖面的二维数据，我希望得到的结果是以p为垂直坐标的三维数据，请问应该怎么做？谢谢

ljchen1989

官网有一个p_interp的程序~

601337653

ljchen1989 发表于 2017-1-18 11:21
官网有一个p_interp的程序~

谢谢你！请问是在wrf官网吗？还是ncl官网？

601337653

ljchen1989 发表于 2017-1-18 11:21
官网有一个p_interp的程序~

我在ncl官网上找到了这个wrf_user_vert_interp，先试一下，谢谢哈~

ljchen1989

601337653 发表于 2017-1-18 11:53
我在ncl官网上找到了这个wrf_user_vert_interp，先试一下，谢谢哈~

http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Tools/p_interp.htm

打翻醋瓶的猫

楼主很赞！

TrialBlazer

601337653 发表于 2017-1-18 11:53
我在ncl官网上找到了这个wrf_user_vert_interp，先试一下，谢谢哈~

你用这个函数好使吗？

601337653

TrialBlazer 发表于 2018-4-11 15:52
你用这个函数好使吗？

挺好使的，插值也比较方便，层数自己设置

TrialBlazer

601337653 发表于 2018-4-12 12:01
挺好使的，插值也比较方便，层数自己设置

我也用这个函数插值但是报错了，能帮我看看是什么原因吗？
我准备先用这个函数插到气压层上，然后用另一个插到经纬度网格上，最后画剖面，但是运行这个函数的时候报错，83行就是wrf_user_vert_interp所在的行
fatal:wrf_vintrp: The pres and field arrays must have the same dimensionality
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 2348 in file /usr/local/ncl6.4.0/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 83 in file wrf_cross_q_pl2.ncl

1. lon = fspan(110,125,15)
   
2. lat = fspan(25,55,30)
   
3. level = ispan(200,1000,100)
   
4. lon@units = "degrees_east"
   
5. lat@units = "degrees_noeth"
   
6. level@units = "hPa"
   
7. a = addfile("wrfout_d01","r")
   
8. times = wrf_user_getvar(a,"times",-1)
   
9. lat2d = wrf_user_getvar(a,"XLAT",0)
   
10. lon2d = wrf_user_getvar(a,"XLONG",0)
    
11. opts = True
    
12. do it = 0,ntimes-1,12
    
13. qrain = wrf_user_getvar(a,"QRAIN",it)
    
14. qrain2 = wrf_user_vert_interp(a,qrain,"pressure",level,opts)
    
15. qrain3 = rcm2rgrid_Wrap(lat2d,lon2d,qrain2,lat,lon,1)
    
16. qrain3!0 = "pressure"
    
17. qrain3!1 = "lat"
    
18. qrain3!2 = "lon"
    
19. qrain3&pressure = level
    
20. qrain3&lat = lat
    
21. qrain3&lon = lon
    
22. plot(ii) = gsn_csm_pres_hgt(wks,qrain3(:,:,{115}),res)
    
23. end do

复制代码

TrialBlazer

TrialBlazer 发表于 2018-4-13 16:09
我也用这个函数插值但是报错了，能帮我看看是什么原因吗？
我准备先用这个函数插到气压层上，然后用另一 ...

弄明白了，我运算的qrain是三维的，但是这个函数读取的气压是四维的，需要在opts@time中设置一下