请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3001|回复: 7

[分享资料] grads站点数据插值,多层次多变量出错

[复制链接]

新浪微博达人勋

发表于 2012-9-24 21:59:33 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
看了很久没发现错在哪儿,请大家帮忙看看。
站点数据:
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

新浪微博达人勋

0
早起挑战累计收入
发表于 2012-9-24 22:11:24 | 显示全部楼层
real stid(sta)
定义错误
就一句话,没认真看懂就下手做。
请看站点作图的帖子:
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-9-24 22:28:02 | 显示全部楼层
这个我发上去就发现了,已经改过了,不是这个问题哦
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-9-24 22:34:30 | 显示全部楼层
这个我改过了character*8 stid(不是这个原因)。出来的图还是没有数据
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-9-24 22:50:25 | 显示全部楼层
我又发现了一个问题,我的ctl也写错了,但改了之后还是不催,图没有数据。我再继续看看有没有别的问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-9-25 09:05:37 | 显示全部楼层
写FORTRAN程序的时候,你先定义了NLEV,却并未给初值就开始循环读取数据 ,直到写的时候才赋值为8……按你写的程序,得出的二进制的数据文件是不是空的啊?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-9-27 08:34:58 | 显示全部楼层
问题已解决,谢谢各位
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-18 16:30:10 | 显示全部楼层
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表