关于二进制文件的读入问题？

evanustc

有点激动，整了快2天的时间，对Fortran理解不少，太有意思了

evanustc

我把 DIMENSION WORK(300)改为INTEGER  WORK(300) 能够读入74 10 22 45但是后面小数就不行了，读的全是整数。如果WORK为REAL型，为啥不能把74当做74.0读入呢？类似C++里面的读入处理

evanustc

necrohan 发表于 2020-5-17 12:49
整数和浮点数在二进制文件里表示方法不一样，不能通用，否则读出来数据是错的。
你附件那个二进制文件按6 ...

是的，二进制读入没有问题。我明白您的意思，只是有点疑惑的就是小工具中关于二进制代码读取的部分：
      DIMENSION WORK(1)
      OPEN(UNIT=IOPDS, FILE=FILNAM, ERR=100, STATUS='UNKNOWN',
     &     ACCESS='SEQUENTIAL', FORM='UNFORMATTED',IOSTAT=IOS)
      READ(UNIT=IOPDS,ERR=300,IOSTAT=IOS)  LENG,(WORK(I),I=1,LENG)
      CLOSE(UNIT=IOPDS, ERR=200, STATUS='KEEP',IOSTAT=IOS)
      RETURN74
WORK声明的是数组，也没有指定类型，在UBUNTU下编译后，转出的10进制文件是正确的。我推断WORK中即读入了整型也有实型，不然打印成十进制不可能对的上。这部分我在win下实现后就出现了我刚才问的问题，屏幕输出的前4个本应该是74 10 22 45却变为了神秘小数。。

evanustc

li913 发表于 2020-5-17 16:13
研究了两天，发现这个代码是比较奇葩的，典型的国外开源代码风格，故意让人看不懂（也可能是代码太老了）。 ...

李老师，请收下我的膝盖orz....太厉害了。在结合工具代码的基础上，阅读您的代码，我明白了意思。在您的代码里，有一项我不太明白原理，就是a(1)读进来的时候还是REAL型的小数，但是在最后的循环中，使用了 b=a(1),然后整形n就成了74了。。。  equivalence(b,n)这一句，实型b和整型n共用同一块内存地址，是由于整型和实型的原因导致取同一块地址里的值时导致不一样的吗？为啥不用n=b转换呢

evanustc

我思考了下，是因为共用一块地址，地址里的值是二进制的，如010110，那么分别用不同类型的变量取值这段二进制的值时，就会翻译为相应的数据类型值。而n=b这种是数据上的强制转换

evanustc

li913 发表于 2020-5-17 20:15
1、equivalence(b,n), 二者的二进制表示形式一样(同一块内存)，但值不一样；
2、n=b, 值一样(或近似，存在 ...

谢谢李老师，彻底明白了！感谢您的指点！继续前行，加油!