grads中站点资料转换为格点资料时，制作站点映射文件时出现错误

sukeyangel

我现在做的是一个将站点资料转换为格点资料的工作，资料是一个逐一小时的闪电资料，从2013070100时开始（共37个时次），目前是在做“ga>！ stnmap”这一步上出现了问题，grads报错如下：

grads报错

如下是我的ctl：
dset f:\wjj-data\output\lighting\light2.dat
dtype station
stnmap f:\wjj-data\output\lighting\light2.map
undef -999.0
title lighting
tdef 37 linear jul2013 1hr
vars 1
l 0 99 lighting power（l是回流闪电强度）
ENDVARS


=====================================================================================


如下是我的转换二进制的fortran程序：
     PROGRAM MAIN  !LSO FILE FOR ARPS
         INTEGER,PARAMETER   ::  N=10000,NUM=37

         INTEGER             ::  TOTAL
         REAL                ::  B(5)

      REAL                ::  SLAT(N),SLON(N),l(N)   
      integer             ::  SELV(N)
      CHARACTER           ::  STNAME(N)*6,name*14

      open(12,file='f:\wjj-data\output\lighting\light2.dat',form='binary')
          open(13,file='f:\wjj-data\output\lighting\light2.txt')                  
      OPEN(15,FILE='f:\wjj-data\namel.txt')

        DO I=1,num
      print *,i
          READ(15,*)NAME

          OPEN(11,FILE='f:\wjj-data\201307lighting\'//NAME)
        do kk=1,4
          READ(11,*)
          enddo

        J=1
10      READ(11,'(a5,3f10.4)',END=20)stname(j),(B(LL),LL=2,4)
!        print*,i,j
          SLAT(J)=B(3)
          SLON(J)=B(2)
        l(j)=b(4)
               
          IF(B(4).EQ.9999.0)l(J)=9999.0
!          print*,j,slon(j),slat(j),l(j)
          J=J+1
          GOTO  10
!       print *,j
20      CONTINUE
        TOTAL=J-1
          CLOSE(11)

          DO  J=1,TOTAL
            WRITE(12)STNAME(j),SLAT(J),SLON(J),0.0,1,1
          WRITE(12)l(j)  
                  WRITE(13,*)STNAME(j),SLAT(J),SLON(J),l(j)

            ENDDO
            WRITE(12)STNAME(1),SLAT(1),SLON(1),0.0,0,1

            
      ENDDO
         CLOSE(12)
    CLOSE(13)
        close(15)
        END PROGRAM

lqouc

sukeyangel 发表于 2014-4-20 15:52
谢谢提醒，我已经贴在主题里了

程序的输出部分有错误，站号要求8位字符，输出的数据要和前面的各种信息写在一行，不能分开写，楼主可以看一看论坛里清风写的教程。

lqouc

可以把f文件也贴出来看看。

sukeyangel

lqouc 发表于 2014-4-20 15:29
可以把f文件也贴出来看看。

谢谢提醒，我已经贴在主题里了

yuppieflu

是不是应该把路径中的‘\’ 改成‘/’

sukeyangel

lqouc 发表于 2014-4-20 16:04
程序的输出部分有错误，站号要求8位字符，输出的数据要和前面的各种信息写在一行，不能分开写，楼主可以 ...

谢谢啊，已经解决了问题了，感谢建议的教程。

立志天下

我今天也遇到这个问题了看到 清风的帖子后 就解决了