多个时次的站点资料怎么写成GrADS数据啊？

wentsing

我用的是753站资料，该资料是每个站点写在一个文件（su*.txt)里，时间是从50-05年。我想该资料做一个全国降水分布的时间序列分析，采用的方法如下：
step1.把753个文件合并为一个文件：（我的系统为Linux）
cat su*>file1.txt
step2.用fortrian程序去读file1.txt,可得到如下结果：
stn(md)       站点号（站点数md=753)
lat(md) ,lon(md)        对应站点的经纬度
rain(mt,md)     站点stn(md)在mt时刻的降水
step3.写成二进制文件，fortran程序如下（部分）：
open(11,file='out.dat',form='unformatted',status='replace')
   do it=1,mt          ！mt=365为总时次
         tim=0.0
         nlev=1
         flag=1

      do id=1,md
              write(11) stn(id),lat(id),lon(id),tim,nlev,flag,rain(it,id)
      enddo

      nlev=0
      write(11) stn(id),lat(id0，lon(id),tim,nlev,flag
enddo
step4.给出out.dat文件的ctl文件：
dset  out.dat
dtype  station
stnmap  out.map
undef  9999.0
title station precipition
tdef 365 linear 01may1957 1dy
vars 1
rain 0 99 precipition
endvars
step5.生成out.map 文件（命令）
stnmap -i out.ctl
问题就出在这里，没有成功！（前面的step1和step2都没有问题）
高手都进来帮忙看看，不甚感激！
注：fortran程序的编译是（pgf90 rain.f90)

mofangbao

楼主请先搜索
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903

qxtlyf

参考网站上的资料

岳焕

请问，753站资料中观测间隔是多少？提供的只是日平均的吗？谢谢了

vera

我的资料是这样的：
站号  纬度  经度 第一个时次降水资料  第二个时次降水资料。。。。。。。。。第n个时次降水资料
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
我已转成grads的格式，也编写了ctl文件，生成map文件，问题是我想知道从第一个时次至第n个时次某一个区域（假如从北纬23至24，东经110-115范围内）平均降水量的逐时变化，不知哪位大侠可以帮我解决这个问题？