站点插值成格点出现错误：warning from oacres:less than two stations

边缘人糊涂乐

站点转格点程序：
!转化成格点
program main
parameter(m=22)
real jd(m),wd(m),a(m),b(m),c(m)
character*8 sta(m)
!-----read data
open(1,file='e:\f\egvt.txt')
do i=1,m
read(1,*) sta(i),jd(i),wd(i),a(i),b(i),c(i)
end do
close(1)
!第一模态
open (1,file='e:\f\1.grd',form='binary')
tim=0.0
nlev=1
nflag=1
do j=1,m
write(1) sta(j),wd(j),jd(j),tim,nlev,nflag,a(j)
print*,sta(j),wd(j),jd(j),tim,nlev,nflag,a(j)
end do
tim=0.0
nlev=0
nflag=1
write(1) sta(j-1),wd(j-1),jd(j-1),tim,nlev,nflag
close (1)
end


背景格点程序：
program sx
parameter (nx=200,ny=200)
real lat(ny),lon(nx),s(nx,ny)
open(1,file='e:\f\grid.grd',form='binary')
lat(1)=35
lon(1)=118
do j=1,ny-1
lat(j+1)=lat(j)+0.1
end do
do i=1,nx-1
lon(i+1)=lon(i)+0.1
end do
print *, 'ok'
do i=1,nx
   do j=1,ny
      s(i,j)=1
   end do
end do
write(1)s
end

两个ctl如下：
dset e:\f\1.grd
dtype station
stnmap e:\f\1.map
undef -999.0
title rain
tdef 1  linear may1981 1yr
vars 1
p  0 99  rain data
endvars

dset e:\f\grid.grd
title  grid
undef -999.0
xdef 200 linear 118  0.1
ydef 200 linear 35  0.1
zdef 1  levels 1000   1
tdef 1  linear may1981 1yr
vars 1
g  0 99 grid data
endvars


gs如下：'reinit'
'open e:\f\grid.ctl'
'open e:\f\1.ctl'
'set grid off'
'set grads off'
'set mpdset  cnworld'
'set lat 35 54'
'set lon 118 137'
'define a=oacres(g,p.2)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set gxout contour'
'set clab forced'
'set cint 0.05'
'd aa'
'enable print e:\f\1.gmf'
'print'
'disable print'
;

出现错误：




求大神解答！万分感谢。



来自群组: 成于大气 信达天下

勇敢的悲伤

应该是缺测值的问题吧~你查一查你的缺测值是多少，和你的CTL文件一样吗

aa22556688

直接用grads打开1.grd的ctl，看是否可以显示出站点的具体数值，如果显示不出来应该是程序的问题，目前没有看出来问题

四叶草

同意@aa22556688 说法，如果站点信息都不能显示，那就是读写时出错了

小预报员

我也遇到过这个问题啊

涛声依旧否

我的是站点能显示出现，但是画等值线就不行了，就出现：OACRES:less than two stations

cch513

请问这个问题是怎么解决的？我的站点资料能绘填图，但是插值格点时出现楼上同样的错误提示，能共享一下解决方案吗？

hello呀喵

求问楼主这个问题解决没有呢？？

紫竹

我今天画图时也出现同样的问题，单独在站点资料显示填图时，能显示，但用oacres插值时就出错，请问楼主台何解决的？

紫竹

我今天画站点资料时，将站点数据写成grd格式，再写ctl，然后用此ctl显示站点资料时，能正常显示。但用oacres插值时就出错，错误提示warning from oacres:less than two stations ，我试着调整网格格点间隔，仍不能解决。 请问楼主你怎么解决的？谢谢！