- 积分
- 2682
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-3-11
- 最后登录
- 1970-1-1
|
GrADS
系统平台: |
|
问题截图: |
|
问题概况: |
程序:
Program jiangshui
parameter (ns=160,itt=50)
!这里是程序的变量声明
character*8 stid(ns)
real lat(ns),lon(ns),rain(itt,ns),tim
integer nlev,flag,st(ns),j
!变量声明结束
!程序开始
!读取变量
open(1,file='e:\test2\3\yearrain.txt',status='old')
do j=1,itt
do i=1,ns
20 read(1,*,end=30)stid(i),lat(i),lon(i),rain(j,i)
goto 20
30 continue
enddo
end do
close(1)
!读取结束
!写入数据
open(2,file='e:\test2\3\yearrain.grd',form='binary')
do j=1,itt
tim=0.0
nlev=1
flag=1
do i=1,ns
stid(i)=char(st(i))
write(2)stid(i),lat(i),lon(i),tim,nlev,flag,rain(j,i)
enddo
nlev=0
flag=0
write(2)stid(ns),lat(ns),lon(ns),tim,nlev,flag
enddo
close(2)
!写入结束
!程序结束
end
ctl文件:
DSET e:\test2\3\yearrain.grd
Dtype station
stnmap e:\test2\3\yearrain.map
undef -9999.0
title station data sample
tdef 50 linear jun1958 1yr
vars 1
rain 0 99 rain data
endvars
;
并且也生成了map文件,我想这是不是没有读入站点信息,但是在上面的程序里没有找出错误,求大神指点! |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
1 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
|
|