求助！读取MPL数据做气溶胶消光系数，编译没错，运行时出错

liutaoz21

201105300600.csv (6.98 MB, 下载次数: 5)

2013-5-8 13:28 上传

点击文件名下载附件

如题，提示错误：forrtl：severe<24>:end-of-file during read，unit22，file d:\fs\data\201105300200.csv

不明白到底怎么出错的，还请各位大大予以指导
附件是所用到的数据之一
程序如下：
        PROGRAM MAIN !用MPL求气溶胶消光系数xa
        IMPLICIT NONE
        INTEGER i,j,k,cn
        INTEGER,PARAMETER::Z=200 !设参照高度hc=6km，即计算至200行
        REAL,ALLOCATABLE::p(:,:) !动态数组，回波信号p
        REAL::pa(Z)=0,xa(Z)=0,xm(Z)=0
        !pa为每小时的平均回波信号;xa为气溶胶消光系数;xm为大气消光系数
        !REAL,DIMENSION(0:h)::h,D,jf
        REAL,DIMENSION(0:Z)::h
        REAL,DIMENSION(0:Z)::D
        REAL,DIMENSION(0:Z)::jf
        !h为高度,MPL每0.03km一个信号;D为回波信号的距离平方修正函数D=p*h*h;jf为(p*h*h)从h到h的积分
        CHARACTER*10 fname !字符串长度10要跟filelist里每一行文件名字符长度一致
        CHARACTER*1 sign,C0(1:3601)
        !3601第1行中包含分隔符在内的字符个数，估计值，用于保证完整地读取一行
        sign=','!设置分隔符“,”
        OPEN(15,file='d:\fs\filelist.txt')!打开列有csv文件名的文件
        DO WHILE(.not.EOF(15))!只要文件没有结尾就进行循环
        READ(15,*)fname !读取文件名
      PRINT*,fname
        cn=1
        !!!!!!!!!!!!!!!!!!!计算csv文件的列数(开始)!!!!!!!!!!!!!!!!!!!!!
        OPEN(12,file='d:\fs\data\'//fname//'00.csv',ACTION='READ')
        !打开刚读取名字的文件
        READ(12,100) C0
100        FORMAT(<3601>A1)
        !PRINT*,C0 !显示第一行数据
        DO i=1,3600,1
          IF( (C0(i)==sign) .AND. (C0(i+1)/=sign) ) cn=cn+1
        ENDDO
        CLOSE(12)
        PRINT*,cn !列数
        !!!!!!!!!!!!!!!!!!!计算csv文件的列数(结束)!!!!!!!!!!!!!!!!!!!!!
        OPEN(22,file='d:\fs\data\'//fname//'00.csv')
        ALLOCATE(p(Z,cn)) !动态数组分配内存
        DO i=2,Z
                DO j=1,cn
                READ(22,*),p(i,j) !将文件数据读入数组中，第一列数据为高度，舍去
                ENDDO
        ENDDO
        CLOSE(22)
        DO i=1,Z
                DO j=2,cn
                pa(i)=pa(i)+p(i,j)
                ENDDO
        pa(i)=pa(i)/(cn-1) !计算每行每小时的平均回波信号
        ENDDO
        DO i=1,Z
        h(i)=h(i)+0.03 !设定高度数据
        D(i)=pa(i)*h(i)*h(i)
        xm(i)=8*3.1415926/3*(1.54e-6)*exp(-h(i)/7000) !标准大气模型，确定h高度处的大气小光系数
        ENDDO
        OPEN(20,file='d:\fs\xishudata\'//fname//'.txt',form='formatted')!打开有格式的txt文件
        WRITE(20,200) 'h','Pave','D','xa','xm'
200        FORMAT(1X,A3,4A6)
        h(0)=h(1)
        D(0)=D(1)
        DO i=1,Z
                DO k=i,Z
                jf(k)=jf(k-1)+(D(k)+D(k-1))*0.015 !0.015为高度步长0.03km的一般，梯形面积法
                ENDDO
        xa(i)=D(i)/(D(Z)/xm(Z)+2*jf(Z))-xm(i) !klett法计算气溶胶消光系数xa
        WRITE(20,300) h(i),pa(i),D(i),xa(i),xm(i)
300        FORMAT(1X,F3.2,4F6.4)
        ENDDO
        CLOSE(20)                        
        ENDDO !打开flielist的do while的结尾
        CLOSE(15)
        ENDPROGRAM

liutaoz21

求助啊{:eb303:}{:eb303:}

悠然

也想做这个 请问楼主问题解决了吗