fortran读取数据问题

qxtlyf

好像用MATLAB可以读出，NaN是缺省值

kongfeng0824

lqouc 发表于 2013-3-15 08:24

                               
登录/注册后可看大图

确实数字会变成字符型，不过要转换回来也很简单，但是整形或者实型不能够识别NAN吧，不像是matlab比较智能 ...

那楼主一般如何判断NaN，并转换成-999呢？敬请指点。

lqouc

kongfeng0824 发表于 2013-3-15 21:49

                               
登录/注册后可看大图

那楼主一般如何判断NaN，并转换成-999呢？敬请指点。

就把数据以字符型读出来啊，读一个判断一个，然后运算之前全部转换成数字就好了，在fortran里面非常好实现的。read，if，write。就这么简单。

kongfeng0824

lqouc 发表于 2013-3-15 22:10

                               
登录/注册后可看大图

就把数据以字符型读出来啊，读一个判断一个，然后运算之前全部转换成数字就好了，在fortran里面非常好实现 ...

请问转换字符到数字应该使用real直接限制就可以吧？我是这么理解的。我没用过这种方法。

lqouc

kongfeng0824 发表于 2013-3-15 22:24

                               
登录/注册后可看大图

请问转换字符到数字应该使用real直接限制就可以吧？我是这么理解的。我没用过这种方法。

你说的我不太清楚，应该是不可以的，用read语句可以转换的。每个人都有自己的方法，不用和别人一样，只是我数据经常这样作比较方便。

godenflame135

楼主：
你的问题描述不是很全面。我有问题想问下。
1 你的数据文件是二进制？是否有效正确读入？
2 如果正确读入，那么在程序中应该加判断，在重新输出到新的数组或文件中。

F平方

lqouc 发表于 2013-3-12 19:44

                               
登录/注册后可看大图

你读取的时候吧数组设置成字符型的，然后读出来每一个数据后进行判断，如果是NAN就赋-999，如果不是的话就转 ...

fortran90/95 有isnan函数

huaner1112

嗯恩，学习了

cgccq

可以借鉴。非常好，学习了！

winne

lqouc 发表于 2013-3-12 21:04
其实这个我不该帮你的，应该你自己来写程序，网上是很容易查到这些方法的，而且我的方法也不一定好，以后尽 ...

您好，我也遇到类似的问题，CSV文件含nan，顺序读取，借鉴了你的方法，代码如下：
open(100,file=trim(fname1))
         do j=1,nlon
            do i=1,nlat
            read(100,*,iostat=error)  dat
            if (dat == 'NaN') then
               dmsc_2d(i,j)=-9999
            else
               read(dat(1:14),*) dmsc_2d(i,j)
            end if
            end do
         end do

但是显示错误：“PGFIO-F-225/list-directed read/internal file/lexical error-- unknown token type”
求指导，非常感谢！