WRF运算结果画图不知道哪里出错？

tianmeng

风之牧语 发表于 2014-10-29 22:04
我觉得你还是先去看看wrf_user_intrp3d的用法吧！你这样用是不对的！

我看了，这个函数既可以插值高度，也可以插值气压

风之牧语

tianmeng 发表于 2014-10-29 22:14
我看了，这个函数既可以插值高度，也可以插值气压

嗯嗯！不好意思，是我弄错了，既然楼主的数据是wrfout数据，是不能用gsn画图的，除非你进行了插值！是要用  contour = wrf_contour(f[0],wks,result_uv,res)
  plots = wrf_map_overlays(f[0],wks,(/contour/),plres,mpres)进行画图的，我用你帖子里面的程序就改了一下这个地方，图就出来了，没有什么问题！不过你1500m的图哈！细节问题你自己调吧！

tianmeng

风之牧语 发表于 2014-10-30 09:12
嗯嗯！不好意思，是我弄错了，既然楼主的数据是wrfout数据，是不能用gsn画图的，除非你进行了插值！是要 ...

多谢，不过我注意到你的图也是1500m以上的图就画不出来了，如果是从地形高度开始1500米要怎么画呢？

风之牧语

tianmeng 发表于 2014-10-30 18:26
多谢，不过我注意到你的图也是1500m以上的图就画不出来了，如果是从地形高度开始1500米要怎么画呢？

wrfout数据里面是有地形高度变量的，你直接加上1500m作为你的z就行啦！

tianmeng

风之牧语 发表于 2014-10-30 21:27
wrfout数据里面是有地形高度变量的，你直接加上1500m作为你的z就行啦！

我加入了地形
a =addfile("./geo_em.d01.nc","r")
  ter = a->HGT_M(0,:,:)
result_u  = wrf_user_intrp3d( avg_u,z,"h",ter+70.,0.,False)
  result_v  = wrf_user_intrp3d( avg_v,z,"h",ter+70.,0.,False)
但是报错,应该是维度对不上

tianmeng

不知道什么原因，加入地形也不可以

longlivehj

tianmeng 发表于 2014-10-30 23:01
我加入了地形
a =addfile("./geo_em.d01.nc","r")
  ter = a->HGT_M(0,:,:)

维数不同，就用conform处理一下。楼主试试下面的代码，看看行不行。
a =addfile("./geo_em.d01.nc","r")
ter = conform(z, a->HGT_M(0,:,:), (/1, 2/))
result_u  = wrf_user_intrp3d( avg_u, z - ter,"h",70.,0.,False)
result_v  = wrf_user_intrp3d( avg_v, z - ter,"h",70.,0.,False)

tianmeng

longlivehj 发表于 2014-11-3 17:39
维数不同，就用conform处理一下。楼主试试下面的代码，看看行不行。
a =addfile("./geo_em.d01.nc","r") ...

大牛跪拜呀！我按你说的，就调整好了！但是好像还是有点点问题，感觉结果有点不可信

tianmeng

tianmeng

tianmeng 发表于 2014-11-4 11:07
大牛跪拜呀！我按你说的，就调整好了！但是好像还是有点点问题，感觉结果有点不可信

这是我的脚本