求助关于FORTRAN生成的二进制文件显示报错，Contouring: nan to nan interval inf

nuistcc

利用FRTRAN编程，读取二进制文件中每日00时刻温度数据，再进行日平均运算，运算结果保存为二进制文件。程序运行正常，读取的日数据也没为题，日平均运算后生成的二进制文件在grads中能打开，但无法显示，出现 Contouring: nan to nan interval inf

不知道出了什么问题，是ctl文件缺省值设置问题吗？缺省值现在是 undef 1.e30

求助各位高手，急！！谢谢

mofangbao

只有说仔细检查了，ctl的缺省值需要和你运算完之后的缺省值一致，既然你描述的过程认为都是正确的，那看不出来有什么问题，把你认为正确的再仔细检查

nuistcc

mofangbao 发表于 2012-4-25 11:54

                               
登录/注册后可看大图

只有说仔细检查了，ctl的缺省值需要和你运算完之后的缺省值一致，既然你描述的过程认为都是正确的，那看不出 ...

谢谢清风。程序只是对二进制数据中每天00时刻的温度进行了日平均计算，将计算结果保存为一个新的二进制文件。Contouring: nan to nan interval inf 网上搜了一下说可能是缺省值有问题。结果的ctl文件是将日平均运算处理前的dat文件的ctl改成的，没改缺省值啊~

兰溪之水

检查一下你的二进制文件是否输出的是实型数据

godenflame135

你的计算代码中，是不是有变量或数组没有赋初值，（如果没有，fortran默认赋值为0），于是被当做分母，就会出现nan。

godenflame135

你的计算代码中，是不是有变量或数组没有赋初值，（如果没有，fortran默认赋值为0），于是被当做分母，就会出现nan。

残月

错误应该是楼上所讲，完全同意

nuistcc

godenflame135 发表于 2012-4-25 13:26

                               
登录/注册后可看大图

你的计算代码中，是不是有变量或数组没有赋初值，（如果没有，fortran默认赋值为0），于是被当做分母，就会 ...

谢谢斑竹。。。代码中的确定义了个浮点型数组real,allocatable:: T2_ave(:,:)用来存放平均值计算结果，没赋初值。allocate(T2_ave(93,84))。
代码中涉及结果文件值计算的只有下面几行。计算时没涉及数组元素做分母的情况呀。
open(200,file="E:\WRF_study\data_test\T2_beijing_d03_2010-10-01-05_ave00.dat",status='REPLACE',form='unformatted',access='direct',recl=93*84)
do j=1,84
    do i=1,93
     T2_ave(i,j)=(T2temp1(i,j)+T2temp2(i,j)+T2temp3(i,j)+T2temp4(i,j)+T2temp5(i,j))/5.0
   end do
  end do
write(200,rec=1) T2_ave(:,:)
close(200)

而且我将T2temp1(i,j)、T2temp2(i,j)、T2temp3(i,j)、T2temp4(i,j)、T2temp5(i,j)这几个从原始文件中读取的数组分别写成二进制文件，用grads画图都没问题，就只有计算结果T2_ave(i,j)写成的二进制出现了以上问题。

冰枫

你算平均的时候应该考虑缺测点的情况吧，不然T2_ave(i,j)=(T2temp1(i,j)+T2temp2(i,j)+T2temp3(i,j)+T2temp4(i,j)+T2temp5(i,j))/5.0
若有个点为缺测的，你算出来的值就不对了呀！即不是缺测也不是正确的值。个人觉得这段应该加if语句，考虑某点缺测时的情况

godenflame135

冰枫 发表于 2012-4-26 08:24

                               
登录/注册后可看大图

你算平均的时候应该考虑缺测点的情况吧，不然T2_ave(i,j)=(T2temp1(i,j)+T2temp2(i,j)+T2temp3(i,j)+T2temp ...

确实应该考虑，一般缺测值会是个很大的值，5个缺测值除以5.0，还是很大的。你的那5个T*temp的缺测值从二进制文件中读取出来的缺测值是多大啊？编译时要不要考虑高低字节互换问题?