- 积分
- 671
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-7-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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)这个变量出现负数?
|
|