Fortran Coder

标题: 计算不出结果 [打印本页]
作者: 中森天下    时间: 2015-9-1 16:46
标题: 计算不出结果
系统版本  win8.1 64位，编译环境microsoft visual studio2010      Intel.Visual.Fortran.Composer.XE.2011
是不是因为内部函数log需要另外定义一下？？求指点
[Fortran] 纯文本查看 复制代码


program ex_Ci

implicit none

    real Ls,Ci,Vs,Rl,Rv,Ht,hli,st,H0,LV

    real,intrinsic::log

    read(*,*) Ls,Vs,Rl,Rv,Ht,hli

    Ci = exp(-4.531+1.6562*(H0)+5.5496*(H0)**2-6.4695*(H0)**3+&

        (-0.474675+0.079*(H0)-1.39*(H0)**2+1.3212*(H0)**3)*log(LV)&

        +(-0.07291+0.088307*(H0)-0.49123*(H0)**2+0.43196*(H0)**3)*&

        (log(LV))**2)

    print *,"标准系数Ci:",Ci

end program

function H0(Ht,hli) 

implicit none

real Ht,hli

real H0

H0=Ht-hli

end function



function LV(Ls,Vs,Rl,Rv) 

implicit none

real Ls,Vs,Rl,Rv

real LV

LV = (Ls/Vs)*sqrt(Rl/Rv)

end function

作者: fcode    时间: 2015-9-1 17:13
H0 和 LV 作为函数，你需要给予传递参数。
作者: 中森天下    时间: 2015-9-1 17:35
怎么改呢？？
作者: fcode    时间: 2015-9-1 17:40
我猜测应该是这样：
[Fortran] 纯文本查看 复制代码
program ex_Ci

 implicit none

     real Ls,Ci,Vs,Rl,Rv,Ht,hli,st,H0,LV

     real :: h0_f , lv_f

     real,intrinsic::log

     read(*,*) Ls,Vs,Rl,Rv,Ht,hli

     h0 = h0_f( ht , hli )

     Lv = Lv_f( Ls,Vs,Rl,Rv )

     Ci = exp(-4.531+1.6562*(H0)+5.5496*(H0)**2-6.4695*(H0)**3+&

         (-0.474675+0.079*(H0)-1.39*(H0)**2+1.3212*(H0)**3)*log(LV)&

         +(-0.07291+0.088307*(H0)-0.49123*(H0)**2+0.43196*(H0)**3)*&

         (log(LV))**2)

     print *,"标准系数Ci:",Ci

 end program

 function H0_f(Ht,hli) 

 implicit none

 real Ht,hli

 real H0_f

 H0_f=Ht-hli

 end function H0_f



 function LV_f(Ls,Vs,Rl,Rv) 

 implicit none

 real Ls,Vs,Rl,Rv

 real LV_f

 LV_f = (Ls/Vs)*sqrt(Rl/Rv)

 end function LV_f
作者: 中森天下    时间: 2015-9-1 17:52
可以了  谢谢楼主！！！



欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2