不知哪一步错了，请大神们帮我看看

小王子殿下

for i=1:1:180
     lat=90.5-i;
     for j=1:1:360
         lon=180.5-j;
         for z=1:1:13350
             if CO(z,606)==lat && CO(z,607)==lon
                 for x=1:1:605
                     C(i,j,x)=CO(z,x);
                 end
             else
                 for x=1:1:605
                     C(i,j,x)=NaN;
                 end
             end
         end
     end
end


CO是一个13350*607的数组，606和607行是它的经纬度信息，前605行是这个经纬度点的时间序列，然后我建立一个180*360的循环，如果是这个经纬度就把值赋给它，如果不在数据里就设为NaN，希望建立一个180*360*605的。
但是我这段程序运行完之后，最后的结果全是NaN，不知道哪里出了问题，麻烦大神们帮我看看哪里出了问题。
。

x6xin

很不错的东西，谢谢！

二爷名声在外

为什么没有结果的截图？

qingqing0917

不是很懂Fortran编程，不过你这个lat和lon不是浮点型吗，浮点型能==这么比较吗？C里面是不能这么比的

天狼007

很不错，谢谢

Adam-Lee

不是很明白为什么要转换成3维的嘞

hulihutu

CO(z,606)==lat && CO(z,607)==lon
有可能是数据中没有数据符合这个条件。小数位数的原因，你试试CO(1,606)==-54.4，看返回的是1还是0

xiazaimatlab

转置一下试试

小王子殿下

二爷名声在外 发表于 2018-2-7 21:24
为什么没有结果的截图？

结果就是180*360*605全部都是nan，我就比如说查看C(:,:,5)，就全是NaN，按照道理来说应该有的点上是有值的呀

小王子殿下

Adam-Lee 发表于 2018-2-24 17:27
不是很明白为什么要转换成3维的嘞

就是转换成180*360然后再乘时间，这样的话就对应每个经纬度，方便后续的处理