计算Z指数的程序运行出错，请高人帮忙

clh914

      program main
parameter (ii=48,jj=63)
real m(ii,jj)
real s(jj)
real p(jj)
real q(jj)
real z(ii,jj)
real w=0.0
real y=0.0
real t=0.0
real u=0.0
      open(1,file='1.txt')
open(2,file='1-z.txt')
read(1,*) ((m(i,j),j=1,jj),i=1,ii)
      
do  j=1,jj
do i=1,ii
  w=w+m(i,j)
end do
s(j)=w/48.0
w=0.0
end do
    do j=1,jj
  do i=1,ii
  y=y+((m(i,j)-s(j))**2)
  t=t+((m(i,j)-s(j))**3)
       end do
p(j)=sqrt(y/48.0)
q(j)=t/(48.0*(p(j)**3))

y=0.0
t=0.0
end do
do  i=1,ii
  do j=1,jj
u=((q(j)*((m(i,j)-s(j))/p(j)))/2.0)+1
  if (u>0.0)then
z(i,j)=((6*(u**0.33333333))/q(j))-(6.0/q(j))+(q(j)/6.0)
       else if (u==0.0)then
z(i,j)=-(6.0/q(j))+(q(j)/6.0)
  else
z(i,j)=(-1)*((6*(((-1)*u)**0.33333333))/q(j))-(6.0/q(j))+(q(j)/6.0)
  end if
u=0.0
end do
      write (2,100) (z(i,j),j=1,jj)
end do
100   format(63f8.2)

end
附件是计算的数据

言深深

您好，您的程序我下载运行了一下，仅针对你的语法看了一下，算法(因为没有做过这个，所以)没有看！
我用的fortran版本是康版的6.5，运行的时候提示12个错误：
--------------------Configuration: z - Win32 Debug--------------------
Compiling Fortran...
D:\ypc2012\j\z index\z.f90
D:\ypc2012\j\z index\z.f90(8) : Error: Syntax error, found IDENTIFIER 'W' when expecting one of: ( : % . = =>
real w=0.0
-----^
D:\ypc2012\j\z index\z.f90(9) : Error: Syntax error, found IDENTIFIER 'Y' when expecting one of: ( : % . = =>
real y=0.0
-----^
D:\ypc2012\j\z index\z.f90(10) : Error: Syntax error, found IDENTIFIER 'T' when expecting one of: ( : % . = =>
real t=0.0
-----^
D:\ypc2012\j\z index\z.f90(11) : Error: Syntax error, found IDENTIFIER 'U' when expecting one of: ( : % . = =>
real u=0.0
-----^
D:\ypc2012\j\z index\z.f90(8) : Error: This variable or component must be of a derived or structure type   [REAL]
real w=0.0
^
D:\ypc2012\j\z index\z.f90(8) : Error: This is not a field name that is defined in the encompassing structure.   [W]
real w=0.0
-----^
D:\ypc2012\j\z index\z.f90(9) : Error: This variable or component must be of a derived or structure type   [REAL]
real y=0.0
^
D:\ypc2012\j\z index\z.f90(9) : Error: This is not a field name that is defined in the encompassing structure.   [Y]
real y=0.0
-----^
D:\ypc2012\j\z index\z.f90(10) : Error: This variable or component must be of a derived or structure type   [REAL]
real t=0.0
^
D:\ypc2012\j\z index\z.f90(10) : Error: This is not a field name that is defined in the encompassing structure.   [T]
real t=0.0
-----^
D:\ypc2012\j\z index\z.f90(11) : Error: This variable or component must be of a derived or structure type   [REAL]
real u=0.0
^
D:\ypc2012\j\z index\z.f90(11) : Error: This is not a field name that is defined in the encompassing structure.   [U]
real u=0.0
-----^
Error executing df.exe.

z.obj - 12 error(s), 0 warning(s)
检查发现因为8-11行的四处赋值语句出现边定义边赋值的问题，去掉定义real改为：
w=0.0
y=0.0
t=0.0
u=0.0

【此处修改时基于fortran程序的I-N规则进行，可以自行查看I-N规则】
再次运行之后得到结果 1-z.txt (23.72 KB, 下载次数: 16)

2012-5-29 10:19 上传

点击文件名下载附件

mofangbao

变量在定义的时候如果想同时初始化值的话，可以使用下面的格式：
real::var=0.0

小肥虾619

楼主，你的Z指数问题是否已解决，有木有Z指数的fortran程序，能不能给我用下，谢谢了

shuiyunjian

谢谢，分享程序，具体怎么使用还有待进一步学习

彭哈儿

你好，请问你的求Z指数的Fortran程序还保留着吗

小公子

shuiyunjian 发表于 2014-12-31 15:49
谢谢，分享程序，具体怎么使用还有待进一步学习

您好，请问Z指数的程序还保留着吗

小公子

彭哈儿 发表于 2016-4-24 09:42
你好，请问你的求Z指数的Fortran程序还保留着吗

您好，请问Z指数的程序还保留呢吗