求助大神帮忙看看NCL脚本

无敌泡

求助大神，这是我写的一个读取WRF结果中的SWDOWN数据并输出的脚本。第10行的意思难道是输出所有时刻，坐标为（18，48）的SWDOWN吗？
如果我想输出我模式的中心点的数据，我该怎样设置呢》？多谢

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  2  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
  3
  4 f =addfile("wrfout_d01_2014-04-30_18:00:00","r")
  5
  6
  7 ; time = f->TIME
  8   swdown= wrf_user_getvar(f,"SWDOWN",-1)
  9
10  slp= f->SWDOWN(:,18,48)
11  print(slp)

longlivehj

wrf_latlon_to_ij找到中心位置对应的ij，然后替换程序第10行中的（18，48）即可。

andrewsoong

这里不是坐标。。。。。。。沙发说得对。。。。

无敌泡

longlivehj 发表于 2014-7-11 13:09
wrf_latlon_to_ij找到中心位置对应的ij，然后替换程序第10行中的（18，48）即可。

NCL真的是太牛了，谢谢您！