ncl如何将wrfout数据插值到站点

风之牧语 · 发表于 2014-9-17 10:05:10

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

在家园上看到一篇帖子讲高斯格点插值到站点，但是在使用linint2_points_Wrap函数的时候总数报错！ MYJfiles=systemfunc("ls -1 "+ MYJ_name +"*")  MYJf=addfiles(MYJfiles,"r")
  ListSetType(MYJf,"cat")
  t2=wrf_user_getvar(MYJf[:],"T2",-1)
  lon2d = wrf_user_getvar(MYJf[:],"XLONG", 0)
  lat2d = wrf_user_getvar(MYJf[:],"XLAT", 0)
  MYJ_t2 = dim_avg_n_Wrap(t2, 0)
  hgt = wrf_user_getvar(MYJf[:],"HGT", 0)
  data  = asciiread("../aws/jieguo6.txt",-1,"float")
  argu = onedtond(data, (/2439,13/))

  aws_t2=argu(:,8)
  lat=argu(:,2)
  lon=argu(:,1)
  MYJ_hgt=linint2_points_Wrap(lon2d,lat2d,hgt, False,lon,lat, 0)
  print(MYJ_hgt)

程序是这样的，程序中lat和lon是自动站的经纬度，但是出现这样的错误！
fatal:Number of dimensions in parameter (0) of (linint2_points_Wrap) is (2), (1) dimensions were expected
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 39 in file trans.sur678.ncl

longlivehj · 发表于 2014-9-17 10:15:34

rcm2points

风之牧语 · 发表于 2014-9-17 10:24:49

longlivehj 发表于 2014-9-17 10:15
rcm2points

哦哦！是这个函数啊！我看了一些资料都说的那个函数！我再请教你一个问题，是不是在wrfout数据的格点之间的转换都是用rcm这一类的函数啊？

hxyj · 发表于 2015-10-31 19:12:25

感谢分享。

yuan/// · 发表于 2015-11-6 10:25:56

再问一个问题，先差值垂直再差值水平，和先水平再垂直，两种效果不一样，这有什么物理意义吗？

zjr · 发表于 2015-11-8 14:05:25

用插值函数就行了

Sunnystone24 · 发表于 2016-1-27 13:36:49

试试看。。。

孤独卡齐莫多 · 发表于 2016-4-13 14:47:06

zjr 发表于 2015-11-8 14:05
用插值函数就行了

请问下用什么插值函数?

礼乐王 · 发表于 2017-3-1 21:12:29

你进行过垂直方向差值没

zhajinlin · 发表于 2017-9-18 13:02:58

楼主你好，我是模式初学者，我想问一下如何将WRF模式输出的一年的数据（一天四个时次）插值到中国的700多个气象站点上。谢谢指教！！！

		自动登录	找回密码
密码			立即注册

[作图] ncl如何将wrfout数据插值到站点

登录后查看更多精彩内容~