求助，提取753站站点资料问题~~~

勇敢的悲伤

这次的资料是08年725个站点逐日资料，数据排列如下图所示：

预览前100条数据记录
v01000	V04001	V04002	V04003	v10004	v12001	v12052	v12053	v13003	v13201	v11002	v14032
台站	年	月	日	平均本站气压	平均气温	日最高气温	日最低气温	平均相对湿度	20-20时降水量	平均风速	日照时数

我用fortran读出来了，可是为什么建立站点CTL的时候不能生成MAP文件？高手们帮帮我啊~~

fortran程序：
integer,parameter::m=725,n=366      
real lat(m),lon(m),yr(m,n),mo(m,n),day(m,n),p(m,n),t(m,n),tmax(m,n),tmin(m,n),rh(m,n),r(m,n),u(m,n),rz(m,n),tim,nlev,flag
character*8 stid(m),stid1(m,n)
open(1,file='f:\huatu\743stagai.txt',status='old')
do i=1,725
read(1,*) stid(i),lon(i),lat(i)
enddo
open(2,file='f:\huatu\08.txt',status='old')
open(3,file='f:\huatu\08.grd',form='binary')
read(2,1000)
1000 format(12i7)
do j=1,366
tim=0.0
nlev=1
flag=1
do i=1,725
read(2,*) stid1(i,j),yr(i,j),mo(i,j),day(i,j),p(i,j),t(i,j),tmax(i,j),tmin(i,j),rh(i,j),r(i,j),u(i,j),rz(i,j)
write(3) stid(i),lat(i),lon(i),p(i,j),t(i,j),tmax(i,j),tmin(i,j),rh(i,j),r(i,j),u(i,j),rz(i,j)
enddo
nlev=0
Write(3)stid,lat,lon,tim,nlev,flag
Enddo
end



站点CTL如下：
DSET   f:\huatu\08.grd
DTYPE  station  
STNMAP f:\huatu\08.map
undef -9.99e+33
TITLE  Station Data Sample
TDEF  366 linear 1jan2008 1dy
VARS 8
p        0  99  Surface Pressure
t        0  99   Temperature
tmax    0  99  max Temperature
tmin     0  99  min Temperature
rh       0  99  xiangduishidu
r        0  99  yuliang
u        0  99  u-wind
rz       0  99  rizhao
ENDVARS

mofangbao

你还有很多要先看会再动手做：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903

river

楼上正解，那么多时次那么多变量没那么简单的

做个霸气的木头

grads做站点数据需要很多工作的，一步到位，需要好好理解清风的东西。再者这类数据的降水有分类说明，像降雪，露什么的，都有算法，需要提前处理下，不然后面都是问题。

陶蜜云

太赞二楼大大分享的教程，对于小白太重要！

Snowyea

多谢二楼大神分享~

Foradslab

能否告知753站降水资料在哪下啊？谢谢~

勇敢的悲伤

Foradslab 发表于 2018-1-2 20:22
能否告知753站降水资料在哪下啊？谢谢~

这个是当时导师给的资料