Fortran 读取SEGY文件提示（95）浮点数转换失败

kuangshenzk · 发表于 2021-6-16 15:19:15

本帖最后由 kuangshenzk 于 2021-6-16 15:53 编辑

读取IBM格式的SEGY数据时，浮点数提示转换失败，使用的是Open语句中 convert ibm，1、这个有什么解决方案，修改这个错误 2、经验证，转换失败处的浮点数应该为0，属于无效值，但是因为存在这个问题所以读取整个数据时间花费太多，有没有什么方案可以提速

kuangshenzk · 发表于 2021-6-16 15:30:45

program main
!use IFWin
!
!aa=call CreateFile("cc.dat",0,0,0,0,0,0)
integer:: aa
integer(kind=1):: a(3600)
integer(kind=1):: b(240)
!character(len=5):: aa
real:: m(16)
!aa="11"//"bb"
open(100,file="BCO3_3D.sgy",form="unformatted",access="stream",convert="ibm")
read(100)a
!do i=1,1000,1
!read(100)b
!read(100,IOSTAT = aa)m
!end do
read(100)b
read(100)m

end program main

kuangshenzk · 发表于 2021-6-16 15:32:21

提示的错误：
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed
forrtl: info (95): floating point conversion failed