爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5523|回复: 4

参考清风大侠的《fortran+grads站点作图详解》,多时次站点数据的问题

[复制链接]

新浪微博达人勋

发表于 2014-8-31 16:13:40 | 显示全部楼层 |阅读模式
GrADS
系统平台:
问题截图:
问题概况: 程序:
Program jiangshui
parameter (ns=160,itt=50)
!这里是程序的变量声明
character*8 stid(ns)
real lat(ns),lon(ns),rain(itt,ns),tim
integer nlev,flag,st(ns),j
!变量声明结束
!程序开始
!读取变量
open(1,file='e:\test2\3\yearrain.txt',status='old')

do j=1,itt
do i=1,ns

20 read(1,*,end=30)stid(i),lat(i),lon(i),rain(j,i)
goto 20
30 continue


enddo
end do
close(1)
!读取结束
!写入数据
open(2,file='e:\test2\3\yearrain.grd',form='binary')

do j=1,itt
tim=0.0
nlev=1
flag=1
do i=1,ns

stid(i)=char(st(i))
write(2)stid(i),lat(i),lon(i),tim,nlev,flag,rain(j,i)
enddo
nlev=0
flag=0
write(2)stid(ns),lat(ns),lon(ns),tim,nlev,flag
enddo
close(2)
!写入结束
!程序结束

end
ctl文件:
DSET e:\test2\3\yearrain.grd
Dtype station
stnmap e:\test2\3\yearrain.map
undef -9999.0
title station data sample
tdef 50 linear jun1958 1yr
vars 1
rain 0 99 rain data
endvars

并且也生成了map文件,我想这是不是没有读入站点信息,但是在上面的程序里没有找出错误,求大神指点!
我看过提问的智慧: 看过
自己思考时长(天): 1

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

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

x

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-31 16:47:22 | 显示全部楼层
问题说的不明不白的,你直接给张截图是什么意思···你是画站点分布,降水等值线,还是什么?grads有没有错误提示?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-31 16:48:36 | 显示全部楼层
不要flag=0,flag=0表示没有地面变量~~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-9-7 11:31:29 | 显示全部楼层
http://bbs.06climate.com/forum.php?mod=viewthread&tid=27685
不要脸的推荐下本人这个帖子,可能对楼主有用
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-9-9 09:33:24 | 显示全部楼层
问题当时没有表述清楚,是无法读入站点信息,后已解决,谢谢大家!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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