用grads读取两个文件的问题

stone

现有一个txt文件，里面有时间，根据这个时间提取nc数据，现在就是不知道读取这个txt的时间，怎么才能和nc文件的匹配呢？

mofangbao

没明白楼主的意思

stone

就是txt里面有一列是时间，该怎么设置才能让grads读为时间

冬日情愫+

读了txt中的时间以后，用这个时间设定grads里面的时间维度就可以使两者对应了。

stone

非常感谢你，按照你的方法确实可以实现，但是对于单个文件结果是对的，对于多个文件确是错的，gs如下

'reinit'
'open W:/work1/fnl_201007.ctl'
*'set gxout fwrite'
*'set fwrite W:/work1/ss.grd'

**---------循环-----------

i=1
while(i<2)
aa=read('W:/work1/rain/jg2007.txt')
aa1=sublin(aa,2)
station=subwrd(aa1,1)
say station
lonn=subwrd(aa1,2)
say lon
latn=subwrd(aa1,3)
say lat
rain=subwrd(aa1,5)

res=subwrd(aa1,4)
say res
year=substr(res,9,12)
say year

month=substr(res,6,3)
say month

day=substr(res,4,2)
say day   

hour=substr(res,1,2)
say hour


tt=1
'set grads off'
'set grid off'
'set xlopts 1 8 0.25'
'set ylopts 1 8 0.25'
'set clopts 1 8 0.25'
'set cthick 16'
'set mpdset cnworld'
while(tt<41)
'set t 'res''

   'set lat 'latn''
   'set lon 'lonn''
'grads2ascii VVELprs W:/work1/rain/w.txt a %g 1'



tt=tt+1

endwhile





i=i+1
endwhile

ff=close('W:/work1/rain/jg2007.txt')

return
;