求助帖关于fortran

吴阿好大人

将站点资料处理为二进制文件时出错，就是GrADS教材91页的程序。编译时没有出现错误，链接时却报错，有没有大神知道原因啊，谢谢！
program main
implicit none
integer i
real vec(37)
open(1,file='1.txt',status='old')
read(1,*)(vec(i),i=1,37)
close(1)
!将文本记录转换为GrADS所支持的二进制记录
call stntogrd(vec)
end
!子程序
subroutine stntogrd(x)
real lat(37),lon(37),x(37)
character*8 stid(37)
!从文件china.dat中读入经纬度
open(2,file='china.txt')
do 20 k=1,37
20 read(2,*)lat(k),lon(k)
close(2)
!给出各站的站号
do 2 i=1,37
2 stid(i)=char(i)
!按排序向二进制文件1.grd写入记录
open(3,file='1.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do 40 i=1,37
write(3) tid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
40 continue
!在文件最后给出一个特别记录表示这个时间组的记录结束
nlev=0
write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
close(3)
return
end

river

我的电脑上也是这个问题，具体不知道什么原因了。建议你换一个程序试试吧

言深深

write(3) tid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
写错了，tid应该是stid吧？？？

小虫小虫

我好像也错过。哈哈。{:eb513:}

吴阿好大人

言深深 发表于 2015-4-5 20:20
write(3) tid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
写错了，tid应该是stid吧？？？

谢谢！！真的啊！

尼罗河星星

我好想也错过

尼罗河星星

怎么这么像是我们寝室的人

sdzxlk

{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}