- 积分
- 43
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-4-25
- 最后登录
- 1970-1-1
|
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
|
|