- 积分
- 58
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-4-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我按照GRADS教材中91页上编写的将站点数据转换成二进制的程序,提示有错,但菜鸟看不懂,求高手指导。
一下是程序:
program main
real vec(4997)
open(1,file='F:\decimaldata\SK\2014050220.txt',status='old')
read(1,*)(vec(i),i=1,4997)
close(1)
call stntogrd(vec)
end
subroutine stntogrd(x)
real lat(4997),lon(4997),x(4997)
character*8 stid(4997)
open(2,file='F:\decimaldata\SK\ZD.txt')
do 2 k=1,4997
20 read(2,*) lat(k),lon(k)
close(2)
do 2 i=1,4997
2 stid(i)=char(i)
open(3,file='2014050220.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do 40 i=1,4997
write(3) stid(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
|
|