FORTRAN作简单的统计遇到问题

咕噜

怪事，隔几个数就不统计，为0 了
具体是文件1是风速数据，文件2是0,0.1,0.2,0.3... 15的151个数据，统计观测到同一风速的次数，文件3是要写入的文件，代码如下，
PARAMETER(p1=8020,p2=151)
real wind1(p1),p(p2),num(p2),psum(p2),fx(p2)
open(1,file='d:\wind\select\m3531nj.txt')
open(2,file='d:\wind\try1.txt')
open(3,file='d:\wind\frequency\m3531nj.txt')

read(1,*)(wind1(i),i=1,p1)
read(2,*)(num(i),i=1,p2)
do i=1,p1
wind1(i)=wind1(i)*0.1
enddo


p(1)=0
do j=1,p2
do i=1,p1
if(wind1(i)==num(j))then
p(j)=p(j)+1
endif
enddo
enddo

do i=1,p2
write(3,'(1x,f10.5)'),p(i)
enddo
end

得到的结果如下，比如0.9,1.3,1.8等风速的次数就为0，而实际上是有值的，这怎么解决呢？还有很多台站要统计，也不能拿excel查找。。请帮忙看看
  483.00000
   43.00000
   65.00000
   93.00000
  106.00000
  108.00000
  152.00000
  165.00000
  202.00000
    0.00000
  223.00000
  222.00000
  187.00000
    0.00000
  195.00000
  199.00000
  201.00000
  216.00000
    0.00000
  188.00000
  192.00000
    0.00000
  186.00000
  182.00000
  162.00000
  158.00000
    0.00000
  149.00000
  140.00000
  139.00000
  140.00000
    0.00000
  126.00000
  122.00000
  106.00000
  104.00000
    0.00000
  106.00000
   82.00000
  111.00000
   78.00000
   80.00000
    0.00000
   86.00000
   62.00000
   71.00000
   70.00000
    0.00000
   40.00000
   49.00000
   54.00000
   46.00000
    0.00000
   45.00000
   41.00000
   44.00000
   41.00000
    0.00000
   41.00000
   33.00000
   28.00000
   34.00000
    0.00000
   35.00000
   26.00000
   21.00000
   18.00000
    0.00000
   19.00000
   17.00000
    8.00000
   22.00000
    0.00000
    8.00000
   13.00000
    9.00000
    9.00000
    0.00000
   11.00000
    7.00000
    7.00000
    5.00000
    9.00000
    2.00000
    0.00000
    6.00000
    3.00000
    2.00000
    4.00000
    0.00000
    1.00000
    5.00000
    5.00000
    2.00000
    0.00000
    3.00000
    0.00000
    2.00000
    2.00000
    0.00000
    1.00000
    1.00000
    0.00000
    4.00000
    0.00000
    1.00000
    1.00000
    2.00000
    1.00000
    0.00000
    1.00000
    2.00000
    0.00000
    0.00000
    0.00000
    1.00000
    0.00000
    0.00000
    1.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    1.00000
    1.00000
    0.00000
    0.00000
    0.00000
    1.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    1.00000
    1.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000

咕噜

自己弄好了，貌似是因为原来的文件风速是乘了10，然后我再程序里乘0.1，但是读取文件时没有指定读取方式，可能再这里出毛病了，所以在风速文件里实现乘以0.1，而不再程序里乘0.1，就没有出错了

zhukefu

这个太简单了吧，判断语句算法不对

咕噜

zhukefu 发表于 2014-12-8 16:01
这个太简单了吧，判断语句算法不对

哪里不对了？num（j）里存放的是0，0.1,0.2,0.3,0.4...15的风速值，只要资料wind1(i)中有等于0时，就开始计数，直到统计好资料里有多少个0风速。0.1也是，0.2也是...