- 积分
- 6169
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-1-5
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2015-1-27 16:28:26
|
显示全部楼层
$debug
real lon(3500),lat(3500),x(3500),num(3500),h(3500),y(2)
character*8 stid(3500)
n=1314
open(1,file='D:\cvf65\MyProjects\rain\s.000',status='old')
do 10 i=1,n
READ(1,*)num(i),lon(i),lat(i),h(i),x(i)
10 continue
close(1)
do 30 i=1,n
30 stid(i)=char(i)
open(2,file='D:\grads\win32e\by\rain.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do 40 i=1,n
write(2)stid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
40 continue
nlev=0
write(2)stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
close(2)
end |
|