爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4953|回复: 12

[分享资料] GRADS站点资料画图求助~~参考过有关求助答案,但是一无所获,不知道有没有大神搭救

[复制链接]

新浪微博达人勋

发表于 2014-6-7 10:34:00 | 显示全部楼层 |阅读模式

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

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

x
先贴出FORTRAN程序:


real::m(12*176),lat(12*176),lon(12*176)
character*8 stid(12*176)
integer nlev,flag
tim=0.0
nlev=1
flag=1

open(1,file='F:\m.dat')

read(1,*)m

open(2,file='F:\lat.dat')

read(2,*)lat

open(3,file='F:\lon.dat')
read(3,*)lon

close(1)
close(2)
close(3)

do i=1,176*12
if(lat(i)==0.0)then
m(i)=-999.0
endif
enddo

open(4,file='F:\station.txt')
do i=1,176*12
read(4,*)stid(i)
enddo
close(4)


open(5,file='F:\station.grd',form='binary')
do i=1,12*176
write(5,*)stid(i),lat(i),lon(i),tim,nlev,flag,m(i)
enddo
nlev=0
write(5,*)stid(i-1),lat(i-1),lon(i-1),tim,nlev,flag,m(i-1)
close(5)
end


我的CTL:
DSET  F:/station.grd
DTYPE station
STNMAP F:/grid.map
UNDEF   -999.0
TITLE DUCT Data
TDEF 1 linear 01jan2014 1mo
VARS 1
m 0 99 DUCT MAGTINUTY
ENDVARS


请大家帮帮忙,如果可以进一步帮我,加我QQ吧,我想跟大家多学点,1284569459

QQ截图20140607101821.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-7 11:21:25 | 显示全部楼层
Fortran输出的时候把*去掉。建议程序最好写的规范一点。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-7 11:45:14 | 显示全部楼层
我以前遇到过i这个问题,后来我自己重新做就解决了,不知道你的问题具体处在哪里,不行你就重头重新再做一次
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-7 16:32:14 | 显示全部楼层
lqouc 发表于 2014-6-7 11:21
Fortran输出的时候把*去掉。建议程序最好写的规范一点。

lqouc,你好!我按照你说的,已经去掉了‘*’,但是还是不行,依然出现:
forrtl: severe (256): unformatted I/O to unit open for formatted transfers, unit
1, file F:\m.dat
Image              PC        Routine            Line        Source
Text1.exe          0040A279  Unknown               Unknown  Unknown
Text1.exe          0040A0D7  Unknown               Unknown  Unknown
Text1.exe          004092B4  Unknown               Unknown  Unknown
Text1.exe          004096E9  Unknown               Unknown  Unknown
Text1.exe          0040475A  Unknown               Unknown  Unknown
Text1.exe          004010DB  Unknown               Unknown  Unknown
Text1.exe          004310D9  Unknown               Unknown  Unknown
Text1.exe          0041A4D9  Unknown               Unknown  Unknown
kernel32.dll       7C816037  Unknown               Unknown  Unknown

Incrementally linked image--PC correlation disabled.
Press any key to continue
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-7 16:34:16 | 显示全部楼层
jeffggg 发表于 2014-6-7 11:45
我以前遇到过i这个问题,后来我自己重新做就解决了,不知道你的问题具体处在哪里,不行你就重头重新再做一 ...

jeffggg,你好!我的数据时处理过后的数据,stid读取的站点是我自己制作的,但是听师兄说,这个是不影响的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-7 16:36:14 | 显示全部楼层
瀟緗繖亽 发表于 2014-6-7 16:32
lqouc,你好!我按照你说的,已经去掉了‘*’,但是还是不行,依然出现:
forrtl: severe (256): unform ...

按照报错把你的1号文件格式声明补充完整,用循环读取。
程序规范是很重要的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-7 16:54:26 | 显示全部楼层
lqouc 发表于 2014-6-7 16:36
按照报错把你的1号文件格式声明补充完整,用循环读取。
程序规范是很重要的。

恩,我先试试,有问题再找你,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-10 19:09:55 | 显示全部楼层
你的dat资料是二进制的还是十进制的啊,二进制的话,读的时候也不需要 * 的啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-23 08:28:44 | 显示全部楼层
river 发表于 2014-6-10 19:09
你的dat资料是二进制的还是十进制的啊,二进制的话,读的时候也不需要 * 的啊

能用MATLAB直接打开,应该是十进制
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-23 12:08:06 | 显示全部楼层
瀟緗繖亽 发表于 2014-6-23 08:28
能用MATLAB直接打开,应该是十进制

没有这么说的,matlab也可以打开二进制的文件。你用文本打开你的资料看看是不是乱码
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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