生成.map是出问题

lightxx

看了之前的一些相关帖子，可是我还是不知道我这个的问题出在哪里了，求哪位大神帮忙看下哦
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

求帮忙解答一下，谢谢

orient

可能是读入的站号数据有问题，最好对照输出一下站号值

福禄寿囍财

问题解决了吗，我也遇到类似的问题

小小书童

应该是站号写入错误了，do i=1,160
                                   stid(i)=char(i)
                               end do
之后才能写入。

手心里的阳光

不知道楼主问题解决了吗？我也是遇到类似问题，看发帖子时间久远，希望还是把问题解决后的帖子~~

Hatter

同问

goodmood

求解决的办法

goodmood

Hatter 发表于 2017-9-23 12:51
同问

请问你的问题有解决吗？？？？