在线急等：专业人员帮我看看这段gs

texas

我写了一段gs，在gs中要调用一个自定义的函数，如代码：
'reinit'
'clear'
'sdfopen d:\hgt.2009.nc'
'set mproj lambert'
'set grads off'
'set lon 60 150'
'set lat 10  80'
'set t 1 365'
'set lev 500'
'enable print d:\fu.gmf'
i=1
while(i<=3)
'set t 'i''
'set grid off'
'set gxout contour'
'set vpage 0 8.5 0 11'
'set parea 0 8   0 10'
'set mpdset mres guangdong'
'set poli on'
'set map color 1 5'
'set cint 2'
'd hgt/10'
date=getDate(2009,10)            -----------这里调用函数
'draw title hgt500 'date''            -----------date变量始终没有值
'set grads off'
'print'
'clear'
'set grads off'
i=i+1
endwhile
'disable print'

function getDate(cyear,day)
return "op"

程序运行后，画的图中date部分为空啊。请问我这错在哪里？

texas

这个。。没有办法吗？

mofangbao

没看懂你这function干什么的，只返回一个“op”？

watertea

试着屏幕输出date变量看看有没数值。。。
另外，你的自定义函数是干嘛用的？如果只是get时间，用query time不就可以了么？

冬日情愫+

是不是返回值定义有问题，你只返回了”op",op是什么？

WWW

'set t 'i''
'query time'
resA=subwrd(result,3)
timeA=substr(resA,1,12)
'draw title hgt500 'timeA''            

只是获取nc文件时间的话，类似这样改改第三四行也许就可以满足你需求了，不需要另外自定义函数。