wrf_user_intrp3d 出问题

奔波尔灞

我用的tt=wrf_user_intrp3d(t2,z,"v",(/164,1200,754,1/),0,True)
指望的是从点a（164,1200）到点b（754,1）做垂直剖面

然而我printVarSUmmary（tt），得到的剖面却是
  Orientation : Cross-Section: (164,1199) to (164,1199)

仿佛只读了数组里前两个数字，但是我看一下源代码，它的确是应该读四个数字的。
这是为什么呢？？

Lighting

上程序，你只是这样描述谁能回答你的问题呢

奔波尔灞

Lighting 发表于 2016-5-11 20:40
上程序，你只是这样描述谁能回答你的问题呢

t1=nc_file1[it]->T(0,:,:,:)
  t2=t1(:23,:,:)
  tt=wrf_user_intrp3d(t2,z(:23,:,:),"v",(/164,1200,754,1/),153.82,True)
  printVarSummary(tt)

这是我的程序

Variable: tt
Type: float
Total Size: 534400 bytes
            133600 values
Number of Dimensions: 2
Dimensions and sizes:        [Vertical | 100] x [Horizontal | 1336]
Coordinates:
Number Of Attributes: 4
  _FillValue :        9.96921e+36
  units :        K
  description :        perturbation potential temperature (theta-t0)
  Orientation :        Cross-Section: (164,1200) to (164,1200)
这是printVarSummary(tt)的结果

七色彩虹

请问楼主这个问题是怎么解决的呢？

奔波尔灞

七色彩虹 发表于 2017-1-10 10:56
请问楼主这个问题是怎么解决的呢？

并没有得到解决，还是用那个一个点加一个角度的吧

zhaolinanshiwo

请问楼主wrf_user_intrp3d是怎么用的？我想得到2m处的风速，是不是需要将ua和va插值到height层上？例如：ua=wrf_user_getvar(a,"ua",-1),va=wrf_user_getvar(a,"va",-1),z=wrf_user_getvar(a,"z",-1)
u_plane = wrf_user_intrp3d(ua,z,"h",2,0.,False)???
求楼主解答~

奔波尔灞

zhaolinanshiwo 发表于 2017-11-4 17:45
请问楼主wrf_user_intrp3d是怎么用的？我想得到2m处的风速，是不是需要将ua和va插值到height层上 ...

是的，抱歉回复晚了

zhaolinanshiwo

奔波尔灞 发表于 2017-11-25 19:36
是的，抱歉回复晚了

嗯嗯，没事！能回复已经很感谢
可是我这样编写最后得到的都是缺测值，之前有人说是地形的问题，不是很理解？？

奔波尔灞

zhaolinanshiwo 发表于 2017-11-26 23:04
嗯嗯，没事！能回复已经很感谢
可是我这样编写最后得到的都是缺测值，之前有人说是地形的 ...

我也遇到过插值都是nan的情况，而且又看不了源码，也不知道问题出在哪

zhaolinanshiwo

奔波尔灞 发表于 2017-11-28 22:18
我也遇到过插值都是nan的情况，而且又看不了源码，也不知道问题出在哪

同样感谢！！
哎  插值不出来心塞塞