我的站点插值格点的程序出现draw error:invalid line coordinate的错误

sukeyangel

我是在做毕业论文，我需要将北京地区的降水站点资料插值到格点中去，但是run后满屏幕滚动draw error:invalid line coordinate，我也看了一些插值的例子，修改了相关程序，但还是无果，再检查也发现不出什么问题，非常无奈。希望科研爱好者能够帮助我解答困惑。

首先是我的原始资料转换二进制程序（这一部分经本人验证没有问题）：PROGRAM MAIN  !LSO FILE FOR ARPS
        INTEGER,PARAMETER   ::  N=11000,NUM=744

        INTEGER             ::  TOTAL
        REAL                ::  B(5)

      REAL                ::  SLAT(N),SLON(N),p(N)   
      integer             ::  SELV(N)
      CHARACTER           ::  STNAME(N)*8,name*12

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

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

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

        J=1
10      READ(11,'(a7,2f7.2,i5,f7.1)',END=20)stname(j),(B(LL),LL=2,5)
!        print*,i,j
          SLAT(J)=B(3)
          SLON(J)=B(2)
        p(j)=b(5)
               
          IF(B(4).EQ.9999.0)p(J)=9999.0
!          print*,j,slon(j),slat(j),p(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)p(j)  
                  WRITE(13,*)STNAME(j),SLAT(J),SLON(J),p(j)

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

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

接下来是我的两个ctl：
1.dset f:\wjj-data\output\r744.dat
dtype station
stnmap f:\wjj-data\output\r744.map
undef -999.0
title rain
tdef 744 linear jul2013 1hr
vars 1
p 0 99 rainfall data
ENDVARS

2.grid的ctl
dset f:\wjj-data\output\grid.grd
undef -999.0
title Sample GRID Data
xdef   301 linear 70 0.3
ydef   151 linear 15 0.3
zdef   1  linear 500 1
tdef   1 linear jul2013 1 hr
vars   1
g 0 99 grid data
endvars

3.这一段是我的gs程序：
'open f:\wjj-data\output\grid.ctl'
'open f:\wjj-data\output\p1.ctl'
'set grads off'
'set grid off'

'enable print f:\wjj-data\output\pc.png'

i=1
while(i<=744)
'c'
'set t ' i


'set lon 105 125'
'set lat 27 50'

'set mpdset cnworld'

'define a=oacres(g(t=1),p.2,1.5)'



'set gxout shaded'

'set clevs 5 8 11 14 17 20 23 26 29 32 35 38 41 44 47 50 53 56'

'run f:\wjj-data\output\color.gs'


'd a'
'cbarn 0.9 1'

'set grid off'
'set gxout contour'
'set clab auto'
'set cint 4'
'd a'
'run c:\OpenGrADS\lib\cbarn.gs'
'run f:\wjj-data\map\china.gs'
'run f:\wjj-data\map\river.gs'



'q time'
title=sublin(result,1)
title=substr(title,8,12)
'draw title ' title


'draw rec 4.3 3 6.7 5.5'



'print'

i=i+1
endwhile

'disable print'
'reinit'
;

river

没看懂，虽然提示了错误，但是最终出来图了么？你是把怎样的资料转成二进制的？转成二进制的程序呢？感觉就是提问的智慧看了和没看是一样的······

sukeyangel

river 发表于 2014-4-15 19:51
没看懂，虽然提示了错误，但是最终出来图了么？你是把怎样的资料转成二进制的？转成二进制的程序呢？感觉就 ...

不好意思，我理解能力还真是弱爆了，即使看了也没有达到把问题说清的水准。
在这里我用的就是最原始的逐一小时降水资料，然后转换为二进制的，我的二进制转换的是没有问题的，图自然没有出来，下面我把转换二进制的程序附上。

river

sukeyangel 发表于 2014-4-15 21:09
不好意思，我理解能力还真是弱爆了，即使看了也没有达到把问题说清的水准。
在这里我用的就是最原始的逐 ...

你还是看一下论坛关于grads站点作图的帖子吧，按照步骤来一遍应该接可以了。你正常生成stnmap文件了吗？

sukeyangel

river 发表于 2014-4-15 21:49
你还是看一下论坛关于grads站点作图的帖子吧，按照步骤来一遍应该接可以了。你正常生成stnmap文件了吗？

我有正常生成。我明天来研究一下论坛的帖子，谢谢了。