爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3313|回复: 2

[秀图] 求助,关于读取grads地图数据的问题

[复制链接]

新浪微博达人勋

发表于 2012-11-22 17:20:14 | 显示全部楼层 |阅读模式

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

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

x
fortran读取grads地形数据的时候出现head为负值。想知道这是什么原因?
program main
integer*1:: ihead(3),ilon(4),ilat(4)
integer lon,lat
real alon(255),alat(255)
equivalence(lon,ilon),(lat,ilat)
open (12,file='yunnan.txt')
write(12,'("diamond 9")')
write(12,*)'0 0 0 0 0 0 0 0 0 0 '
open (11,file='yunnan',form='binary')
do while (.not.eof(11))
read (11) ihead
print *,(ihead(i),i=1,3)
!if(ihead(3)<=0) cycle
do j=1,ihead(3)
read (11) (ilon(i),i=3,1,-1)
read (11) (ilat(i),i=3,1,-1)
alon(j)=float(lon)/1e4
alat(j)=float(lat)/1e4-90.0
enddo
write(12,*) ihead(3),'   map 276 1 65535'
write(12,'(255(f7.3,1x,f7.3,1x))') (alon(j),alat(j),j=1,ihead(3))
enddo
close(11)
close(12)
end


ihead(3)这个变量出现负数?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-11-23 09:32:57 | 显示全部楼层
大侠们,给我点希望吧?不胜感激
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-1 12:59:33 | 显示全部楼层
同样的问题,是不是数据存储方式的差异
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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