- 积分
- 12
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-9-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
用Fortran处理十进制文件转换成二进制文件
程序如下:
program dat
parameter(ms=62,mt=39)
integer is,it
integer nlev,nflag
parameter(nlev=1,nflag=1)
integer t(mt)
character*8 endstn
character*8 stid(ms)
real lat(ms),lon(ms),rain(mt,ms)
data endstn/'end_time'/
open(10,file='C:\Users\ACER-1\Desktop\graduation\data\rainfull\southernrainfull3.txt')
200 format (i4,a5,f5.2,f6.2,f6.2)
do it=1,mt
do is=1,ms
read(10,200) t(it),stid(is),lat(is),lon(is),rain(it,is)
end do
end do
close(10)
open(20,file='C:\Users\ACER-1\Desktop\graduation\data\rainfull\rainfull3.dat',form='binary',status='replace')
do it=1,mt
do is=1,ms
write(20) stid(is),lat(is),lon(is),0.0,nelv,nflag
write(20) rain(it,is)
end do
write(20) endstn,0.0,0.0,0.0,0,0
end do
close(20)
write(*,*) t(1),stid(1),lat(1),lon(1),rain(1,1)
end
相应ctl文件:
dset C:\Users\ACER-1\Desktop\graduation\data\rainfull\rainfull3.dat
dtype station
stnmap C:\Users\ACER-1\Desktop\graduation\data\rainfull\rainfull3.map
title station data
undef -999.0
tdef 39 linear mar1979 1yr
vars 1
rainfull 0 99 march rainfull data
endvars
用grads生成map映射文件报错:
invalid station hdr found in station binary file
想问问各位是哪里出错了呀?
谢谢!
|
|