关于NCL脚本问题

BDRUSH

最近在做站点降水图，然后以http://blog.sina.com.cn/s/blog_44f7c0840100on44.html为例，对其中的几句并不是太懂，麻烦各位大神讲解下
首先创建存放差值后生成数据的数组 根据亚洲截图的经纬度而定北纬17-57，东经72-138这个矩形框内插值
  olon = new(66,"float");
  olat = new(40,"float");
  data1 = new((/40,66/),"float")
  do i=0,65
     olon(i) =72+i
  end do
  do l=0,39
     olat(l) = 17+l
  end do
和最后调用插值函数
  rscan = (/10,5,3/)   ;连续的有效半径大小，最大为10，依次递减
  data1 = obj_anal_ic_deprecated(lon,lat,R,olon,olat,rscan, False)  ;Creanm插值
这两段具体是什么意思？第一段为什么这么做？一直不太理解

longlivehj

第一部分就是生成输出的网格，实际上可以简单一点：
olon = ispan(72, 138, 1)
olat = ispan(17, 57, 1)

第二部分就是调用插值函数，按声明一个个参数往里放。

BDRUSH

longlivehj 发表于 2014-3-26 16:08
第一部分就是生成输出的网格，实际上可以简单一点：
olon = ispan(72, 138, 1)
olat = ispan(17, 57, 1)
...

这两部分自己看着看着就看懂了...不过也谢谢你的回答，我还有一个疑问就是
lines = asciiread(argu(2)+argu(3),-1,"string") 中argu92）+argu(3)又是什么意思呢？

longlivehj

BDRUSH 发表于 2014-3-26 17:02
这两部分自己看着看着就看懂了...不过也谢谢你的回答，我还有一个疑问就是
lines = asciiread(argu(2)+a ...

就是把config.txt中的第三行何第四行组合成一个字符串，应该就是数据文件的完整路径。
你可以把config.txt的四行复制保存到本地，然后在ncl里面一行一行的执行代码，并打印相关记过，就容易明白程序的目的了。

BDRUSH

longlivehj 发表于 2014-3-26 17:12
就是把config.txt中的第三行何第四行组合成一个字符串，应该就是数据文件的完整路径。
你可以把config.t ...

我没有源文件，就自己找了个数据文件一步一步进行print，然后运行到这步就不能出结果了……

BDRUSH

额，两个小财神卡，给力

BDRUSH

longlivehj 发表于 2014-3-26 17:12
就是把config.txt中的第三行何第四行组合成一个字符串，应该就是数据文件的完整路径。
你可以把config.t ...

应该差不多是像你所说的那样理解，3Q拉

言墨

求问插值半径有什么含义？好像插值半径不同  画出来的图也不太一样