站点资料处理为格点资料的问题，有没有大神可以告诉我这种格式的资料要怎么写循环啊！

DGT

有谁能帮帮我解决这个问题，妹纸一定感恩戴德，十万火急的事情，大家帮帮忙，帮帮忙啊！有组织里的人更要帮帮了！天王盖地虎啊！

程序运行出来就是这样的

shuju.txt (1.19 MB, 下载次数: 2)

2014-4-30 15:16 上传

点击文件名下载附件
我需要处理的资料就是这种格式的，共27个站，50年的气温资料，想把它转成二进制，但琢磨了很久，实在是无能 ... zhandian.f90 (370 Bytes, 下载次数: 1)

2014-4-30 15:20 上传

点击文件名下载附件
这个是用过的第一个程序 shuju.f90 (543 Bytes, 下载次数: 1)

2014-4-30 15:21 上传

点击文件名下载附件
这是第二个，但都不可以
第一个程序
real lon,lat,vis(6)
real tim
integer nlev,nflag
character stationid*8
tim=0.0
open(88,file='test.dat',form='binary')
open(11,file='shuju.txt')
read(11,*)
nlev=1
nflag=1
do while(.not.eof(11))
read(11,*)lat,lon,vis
write(88)stationId,lat,lon,tim,nlev,nflag,(vis(i),i=1,6)
end do
nlev=0
write(88)stationId,lat,lon,tim,nlev,nflag
close(88)
close(11)
end
第二个程序
real vec(27)
open(1,file='d\shuju.txt',status='old')
read(1,*) (vec(i),i=1,27)
close(1)
call stntogrd(vec)
end
subroutine stntogrd(x)
real lat(16200),lon(27),x(27)
character*8 stid(27)
open(2,file='d\test.dat')
do 20 k=1,16200
20 read(2,*)lat(k),lon(k)
close(2)
do 2 i=1,27
2 stid(i)=char(i)
open(3,file='d:\data\wj.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do 40 i=1,27
write(3) stid(i),lat(i),lon(i)
# ,tim,nlev,nflag,x(i)
40 continue
nlev=0
write(3) stid(i-1),lat(i-i),lon(i-1),tim,nlev,nflag
close(3)
return
end
数据也截上来一段：
站号 纬度 经度 年份 月份 月平均气温 月均最高气温 月均最低气温
51886 3813 9031 1961 1 -14.1 -4.2 -23.1
51886 3813 9031 1961 2 -10.7 -1.3 -20.8
51886 3813 9031 1961 3 -4.8 5.1 -14.8
51886 3813 9031 1961 4 4.2 13.9 -7.1
51886 3813 9031 1961 5 8.2 18 -5.2
51886 3813 9031 1961 6 13.7 21.8 2.9
51886 3813 9031 1961 7 15.7 24.7 5.4
51886 3813 9031 1961 8 14 23.3 3.8
51886 3813 9031 1961 9 8.5 19.2 -2.3
51886 3813 9031 1961 10 1.1 11 -8.2
51886 3813 9031 1961 11 -9.3 1.2 -19
51886 3813 9031 1961 12 -18.2 -7.2 -27.4
51886 3813 9031 1962 1 -16.2 -5.5 -26.2

lqouc

第一，这个时候就不要卖萌了，这里不是糗百，你把自己炖了蘑菇也没人能帮你。
第二，请下次一定把程序直接贴上来，问题说清楚，不要都写在附件的题目上，这是对人对己的负责。这次我给你贴上来了。

报错是说没找到你的数据文件，还有在你第二个程序的子程序里打开了一个不知道是什么的文件，反正你没贴上来。你的问题其实没什么，主要就是程序写的不好，至于怎么写论坛有很详细的说明，去搜索一下清风的帖子，grads站点数据的保存方式。

DGT

lqouc 发表于 2014-4-30 15:50
第一，这个时候就不要卖萌了，这里不是糗百，你把自己炖了蘑菇也没人能帮你。
第二，请下次一定把程序直接 ...

很犀利啊！不过还是谢啦！因为第一次发帖子，有些问题不是很清楚！

lqouc

DGT 发表于 2014-4-30 16:21
很犀利啊！不过还是谢啦！因为第一次发帖子，有些问题不是很清楚！

不清楚没关系，这不是告诉你了么。
犀利算不上吧，因为这里不是糗百（不歧视糗百，我也喜欢看，只是看看），所以你那么问可能得不到帮助，这里有这里的方式，不停喊叫帮忙没有用，不如把问题描述清楚才会得到帮助，或者在论坛找找，初学者的问题这里基本都有帖子进行总结了。
那么~欢迎常来。

DGT

lqouc 发表于 2014-4-30 16:38
不清楚没关系，这不是告诉你了么。
犀利算不上吧，因为这里不是糗百（不歧视糗百，我也喜欢看，只是看看 ...

谢啦！刚刚有搜清风的《fortran+grads站点作图详解》，但是找不到，是原来发在站子里的吗？

lqouc

DGT 发表于 2014-4-30 16:40
谢啦！刚刚有搜清风的《fortran+grads站点作图详解》，但是找不到，是原来发在站子里的吗？

可能是论坛自带的搜索引擎的问题，我直接把附件给你好了，先看看。有问题在来问吧。 GrADS站点文件作图详细解决方案.pdf (1.35 MB, 下载次数: 90)

2014-4-30 17:02 上传

点击文件名下载附件

xsjvn_

你的资料是不是几十年的气候资料类型的？我正在处理的资料格式和楼主你的有点类似，建议你先把文件的读取写好，然后文件用循环批处理实现，耐心研究还是不难的，。。

DGT

xsjvn_ 发表于 2014-5-1 21:50
你的资料是不是几十年的气候资料类型的？我正在处理的资料格式和楼主你的有点类似，建议你先把文件的读取写 ...

好的，谢谢你！