计算Z指数的fortran主程序搞不定……

viloan

我编的主程序主要是为了读取TXT文档中的数据并且将计算出来的Z指数输出。
TXT中的数据格式是年月日降水量，用tab隔开的。

程序如下：
program main
implicit none
integer i,j,m
real zz
real,dimension(1:3,1:16107)::a
real,dimension(m)::zz
open(1,file='rainfall.txt')
read(1,*)((a(i,j),j=1,16107),i=1,3)
read *,m
call zindex(a(3,j),m,zz)
print (1,*)(zz(i),i=1,m)
end program main

!降水Z指数计算子程序
!aa是降水序列，zz是得到的Z指数序列，n是时间长度
!**************************************************************************   
    subroutine zindex(aa,n,zz)
    implicit none
    integer i,n
    real,dimension(n)::aa,ii,zz
    real bzc,cs,av,ls
    ls=0
    call biaozhuncha(aa,n,bzc,av)
    do i=1,n
    ls=ls+(aa(i)-av)**3.0
    end do
    cs=ls/(n*(bzc**3.0))
    do i=1,n
    ii(i)=(aa(i)-av)/bzc
    end do
    do i=1,n
    zz(i)=(ii(i)*cs/2+1)**(1.0/3)
    zz(i)=zz(i)*6/cs-6/cs+cs/6
    end do

    end
!*************************************************************************   
    subroutine biaozhuncha(aa,n,bzc,av)
    implicit none
    integer i,n
    real,dimension(n)::aa
    real bzc,av,av1
    av=0
    do i=1,n
    av=av+aa(i)
    end do
    av=av/n
    av1=0
    do i=1,n
    av1=av1+(aa(i)-av)**2.0
    end do
    av1=av1/n
    bzc=sqrt(av1)
    end

云在青山

这个程序有问题吗？看不懂啊，我也要Z指数的程序

kongfeng0824

请问降水量的Z指数是描述什么特征的呢？我只用过降水集中度和集中期。

viloan

kongfeng0824 发表于 2013-5-27 22:12

                               
登录/注册后可看大图

请问降水量的Z指数是描述什么特征的呢？我只用过降水集中度和集中期。

是描述旱涝等级的一种指标

kongfeng0824

viloan 发表于 2013-5-27 22:44

                               
登录/注册后可看大图

是描述旱涝等级的一种指标

哦呵呵谢谢

ly19850330

最近挺多人做干旱的啊

蔚蓝的天

顶一个哈哦

sweetbox21200

求程序
啊啊啊啊啊啊啊

wexuan

不太懂，纯支持

Eegle

楼主的这个程序能用么？还是有问题？