爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 11789|回复: 11

用GRADS作图时出现的问题,请大家帮帮忙,很急~~谢谢

[复制链接]

新浪微博达人勋

发表于 2014-6-7 10:24:06 | 显示全部楼层 |阅读模式
GrADS
系统平台: XP
问题截图:
问题概况: 我将自己出来出来的m值和经纬度值用MATLAB读取出来,变成了dat文件,然后准备用这些数据绘制填色图,但是在!stnmap出现了问题:
Starting X server under /Xming
Starting GRADSD~1 under f:\grads\GrADS19\win32 ...

Grid Analysis and Display System (GrADS) Version 1.9.0-rc1
Copyright (c) 1988-2007 by Brian Doty and IGES
Center for Ocean-Land-Atmosphere Studies (COLA)
Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information

Config: v1.9.0-rc1 32-bit little-endian readline sdf/xdf netcdf dap dap-stn lats
athena printim

Issue 'q config' command for more information.

Landscape mode? (no for portrait):
GX Package Initialization: Size = 11 8.5
ga-> !stnmap -i F:/station.ctl
Name of binary data set: F:/station.grd
Number of times in the data set: 1
Number of surface variables: 1
Number of level dependent variables: 0

Starting scan of station data binary file.
Binary data file open: F:/station.grd
Processing time = 1
Invalid station hdr found in station binary file
Possible causes: Invalid level count in hdr
Descriptor file mismatch
File not station data
Invalid relative time
levs = 808463173 flag = 774905888 time = 6.40969e-10
我看过提问的智慧: 看过
自己思考时长(天): 15

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

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

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.grdDTYPE station STNMAP F:/grid.mapUNDEF   -999.0TITLE DUCT Data TDEF 1 linear 01jan2014 1moVARS 1m 0 99 DUCT MAGTINUTYENDVARS请大家帮帮忙,谢谢~~有进一步讨论的好心人,可以加我QQ:1284569459

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

新浪微博达人勋

发表于 2014-6-12 13:24:15 | 显示全部楼层
建议楼主把你的数据文件也上传一份,这样他人也好对你的程序进行测试,从而更快的找到解决问题的方法。
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-6-7 16:14:22 | 显示全部楼层
把open(5,file='F:\station.grd',form='binary')这句改成
open(5,file='F:\station.grd',form='unformatted', recordtype='stream')
再试试看
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-7 16:26:58 | 显示全部楼层
linikiku 发表于 2014-6-7 16:14
把open(5,file='F:\station.grd',form='binary')这句改成
open(5,file='F:\station.grd',form='unformatte ...

不行,C:\Documents and Settings\Administrator\桌面\QQ截图20140607162356.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-7 16:28:04 | 显示全部楼层
linikiku 发表于 2014-6-7 16:14
把open(5,file='F:\station.grd',form='binary')这句改成
open(5,file='F:\station.grd',form='unformatte ...

forrtl: severe (257): formatted I/O to unit open for unformatted transfers, unit
5, file F:\station.grd
Image              PC        Routine            Line        Source
Text1.exe          0040A9F9  Unknown               Unknown  Unknown
Text1.exe          0040A857  Unknown               Unknown  Unknown
Text1.exe          00409A34  Unknown               Unknown  Unknown
Text1.exe          00409E69  Unknown               Unknown  Unknown
Text1.exe          00406121  Unknown               Unknown  Unknown
Text1.exe          00448B2D  Unknown               Unknown  Unknown
Text1.exe          0043F809  Unknown               Unknown  Unknown
Text1.exe          00428BD9  Unknown               Unknown  Unknown
kernel32.dll       7C816037  Unknown               Unknown  Unknown

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

新浪微博达人勋

发表于 2014-6-8 10:42:59 | 显示全部楼层
瀟緗繖亽 发表于 2014-6-7 16:28
forrtl: severe (257): formatted I/O to unit open for unformatted transfers, unit
5, file F:\stat ...

这样的话,我就不明白了。请高人来解答吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-11 13:44:44 | 显示全部楼层
解决没,楼主。我也遇到了,一模一样
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-11 14:47:50 | 显示全部楼层
character*8 stid(12*176)
我的是这句有问题,改完就好了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-23 08:26:49 | 显示全部楼层
nuister 发表于 2014-6-11 14:47
character*8 stid(12*176)
我的是这句有问题,改完就好了。

我把程序重新启动了一下就好了。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-6-23 08:27:11 | 显示全部楼层
四叶草 发表于 2014-6-12 13:24
建议楼主把你的数据文件也上传一份,这样他人也好对你的程序进行测试,从而更快的找到解决问题的方法。

谢谢,已经解决
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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