grads站点数据插值，多层次多变量出错

miracle

看了很久没发现错在哪儿，请大家帮忙看看。
站点数据：
57749 110 27.56 925 20 19 295 3     站号   经度   维度   层次 温度 露点  风向  风俗
57749 110 27.56 850 19 17 205 13
57749 110 27.56 700 10 6 270 17
这是其中一部分，共有17个站点，8层（925,850,700,500,400,300,250,200）
下面是我的fortran程序，将站点转为二进制的：
program main
parameter (sta=17) !定义站点数
real stid(sta)
real lon(sta),lat(sta),lev(8),t(sta,8),td(sta,8),f(sta,8),v(sta,8),tim
integer nlev,flag,j,i
open(1,file='e:\探空资料\05050508.txt',status='old')
open(2,file='e:\grid\05050508.grd',status='replace',form='binary')
do i=1,sta
do j=1,nlev
read(1,*) stid(i),lon(i),lat(i),lev(j),t(i,j),td(i,j),f(i,j),v(i,j)
!10 read(1,*,end=100) stid(i),lon(i),lat(i),lev(j),t(i,j),td(i,j),f(i,j),v(i,j)
!print*,  stid(i),lon(i),lat(i),lev(j),t(i,j),td(i,j),f(i,j),v(i,j)
end do
end do
close(1)
tim=0.0
nlev=8
flag=0  !无地面变量
do i=1,sta
write(2)   stid(i),lon(i),lat(i),tim,nlev,flag
do j=1,nlev
write(2)  lev(j),t(i,j),td(i,j),f(i,j),v(i,j)
end do
end do
!goto 10
!100 continue
nlev=0
write(2) stid(i-1),lat(i-1),lon(i-1),tim,nlev,flag
end
程序可以运行，也可以生成grd文件。
ctl文件如下：DSET e:\grid\05050508.grd
DTYPE  station
STNMAP e:\grid\station.map
undef 9999
TITLE  Station Data Sample
tdef 1 linear 00Z05may2005 1mo
VARS 5
lev 8 99 presure
t 8 99 Temp
td 8 99 dewpoint temp
f 8 99 wind direction
v 8 99 wind
ENDVARS
接着在grads里运行！stnmap -i e:/grid/station.ctl
可以生成map文件，
但是提示如下错误。读不出数据。
time=1 has stn count=0     我的站数是17，怎么会读出来0呢
不知道什么原因啊，这样出来的图是空白。看了好久了，请帮帮忙。
   

mofangbao

real stid(sta)
定义错误
就一句话，没认真看懂就下手做。
请看站点作图的帖子：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903

miracle

这个我发上去就发现了，已经改过了，不是这个问题哦

miracle

这个我改过了character*8 stid(不是这个原因）。出来的图还是没有数据

miracle

我又发现了一个问题，我的ctl也写错了，但改了之后还是不催，图没有数据。我再继续看看有没有别的问题

冰点

写FORTRAN程序的时候，你先定义了NLEV，却并未给初值就开始循环读取数据 ，直到写的时候才赋值为8……按你写的程序，得出的二进制的数据文件是不是空的啊？

miracle

问题已解决，谢谢各位

whmooncool

miracle 发表于 2012-9-27 08:34

                               
登录/注册后可看大图

问题已解决，谢谢各位

请问你的问题出在哪了？我也出现了类似的问题，谢谢