- 积分
- 16
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-3-10
- 最后登录
- 1970-1-1
|
发表于 2014-4-2 14:04:28
|
显示全部楼层
大哥,在不 有问题想请教啊
看过你的几个帖子 还是搞不懂哪出错了
这是我的fortran
Program sta2grd
parameter (ns=2109)
!这里是程序的变量声明
character*8 stid(ns)
real lat(ns),lon(ns),var(ns),rain(ns),tim
integer nlev,flag
!变量声明结束
!程序开始
!读取变量
open(1,file='F:\date\MicapsRain\R6\13081308.txt',form='formatted',status='old')
do i=1,ns
read(1,*)stid(i),lat(i),lon(i),var(i),rain(i)
end do
close(1)
!读取结束
!写入数据
open(2,file='F:\date\MicapsRain\R6\13081308.dat',status='replace',form='binary')
tim=0.0
nlev=1
flag=1
do i=1,ns
write(2)stid(i),lat(i),lon(i),var(i),tim,nlev,flag,rain(i)
end do
nlev=0
flag=0
write(2)stid(ns),lat(ns),lon(ns),tim,nlev,flag
close(2)
!写入结束
!程序结束
end
ctl:
dset F:\date\MicapsRain\R6\13081308.dat
dtype station
STNMAP F:\date\MicapsRain\R6\13081308.map
undef -999.0
title 13081308 Rain Data Sample
tdef 1 linear 08Z13Aug2013 6hr
vars 1
rain 0 99 Rain Data
endvars
运行了之后是这样的:
ga-> ! stnmap -i F:/date/MicapsRain/R6/13081308.ctl
Name of binary data set: F:\date\MicapsRain\R6\13081308.da
Number of times in the data set: 1
Number of surface variables: 1
Number of level dependent variables: 0
Starting scan of station data binary file.
Binary data file open: F:\date\MicapsRain\R6\13081308.dat
Processing time step 1
Invalid station hdr found in station binary file
Possible causes: Invalid level count in hdr
Descriptor file mismatch
File not station data
Invalid relative time
levs = 0 flag = 1 time = 101
求帮忙解答一下,谢谢 |
|