请教各位大侠，我的程序有两个错误，是Error: Constants and expressions are invalid in read-only I/O lists.   [FDATA]，不知道该怎么解决，请兄弟姐妹们帮帮我好吗，我把程序复制到下面了，附件里也有的
           program inte1
    use imsl
    
    integer iweigh,ndata,n
    parameter (ndata=300,n=1)
    real a,b,omega,errabs,errrel,result,errest,f !积分子函数程序的参数
    
    real fdata(ndata),xdata(ndata),float,value,xvec !插值程序的参数
    external f
    common /part1/rset
    common/part2/xdata
    common/part3/fdata

    open(unit=1,file='22.den',status='old')
    a=0.
    b=12.
    iweigh=2
    omega=0.5
    errabs=0.0
    errrel=0.001
    do 10 i=1,300
         rset=float(i)*0.04
10    continue
    call qdawo(f,a,b,iweigh,omega,errabs,errrel,result,errest)
    write(*,*)'result',result    
    end
    real function f(rset)        
    do 20 i=1,300
    read(1,*) xdata(i),fdata(i)
20    continue
    xvec=rset
    call csiez(ndata,xdata,fdata,n,xvec,value)
    f=value
    end