Fortran Coder

查看: 118|回复: 2

[数值问题] 我这个小程序哪儿有问题啊,谢谢

[复制链接]

17

帖子

5

主题

0

精华

入门

F 币
84 元
贡献
51 点
发表于 2020-7-10 10:18:41 | 显示全部楼层 |阅读模式
[Fortran] 纯文本查看 复制代码
program tyt20
implicit none
integer i,j,tmp
real*8::Z,Zi,u1,u2,Zf1,Zf2,R
real,parameter::pi=3.14159
real*8 A,B,C,D
real*8 A2,A4,A01,A02,b1,b2,n1,n2
real*8 E
real*8 F
real*8 H1,H2
real*8 t         !!反射系数的平方用t表示
character(len=80)::filename="data.dat"
integer,parameter::fileid=10 
open(10,file='data.dat')

R=2
H1=80
do i=1,int(40/0.01)
    H1=H1+0.01
    !n2=n2+0.001
    !R=R+1   
H2=80
do j=1,int(40/0.01)
     H2=H2+0.01
u1=1+1200.0/H1
u2=1+1200.0/H2
  n1=0.45
  n2=0.3
  A01=0.5
  A02=0.3
  Zi=-3                 
  Z=50.0                 
  A2=0.375    
  A4=0.35       
  Zf1=sqrt(u1/14.0)  
  Zf2=sqrt(u2/14.0)  
b1=(Zf1*tan(sqrt(u1*14.0)*n1*2.0*pi*A01)+Z*tan((1-n1)*2.0*pi*A01))/(Z-Zf1*tan(sqrt(u1*14.0)*n1*2.0*pi*A01)*tan((1-n1)*2.0*pi*A01))
b2=(Zf2*tan(sqrt(u2*14.0)*n2*2.0*pi*A02)+Z*tan((1-n2)*2.0*pi*A02))/(Z-Zf2*tan(sqrt(u2*14.0)*n2*2.0*pi*A02)*tan((1-n2)*2.0*pi*A02))
A=(cos(2*pi*A2)+sin(2*pi*A2)/b2)*cos(2*pi*A4)-sin(2*pi*A2)*sin(2*pi*A4)
B=((cos(2*pi*A2)+sin(2*pi*A2)/b2)*sin(2*pi*A4)+sin(2*pi*A2)*cos(2*pi*A4))*Z
C=(((-1/b1)*cos(2*pi*A2)+sin(2*pi*A2)-((1/b1)*sin(2*pi*A2)+cos(2*pi*A2))/b2)*cos(2*pi*A4)+((1/b1)*sin(2*pi*A2)+cos(2*pi*A2))*sin(2*pi*A4))/Z
D=((1/b1)*cos(2*pi*A2)-sin(2*pi*A2)+((1/b1)*sin(2*pi*A2)+cos(2*pi*A2))/b2)*sin(2*pi*A4)+((1/b1)*sin(2*pi*A2)+cos(2*pi*A2))*cos(2*pi*A4)
E=(A*R*(D-C*Zi)+C*R*(B+A*Zi))/(((D-C*Zi)**2+(R*C)**2)*Z)
F=((B+A*Zi)*(D-C*Zi)-(A*C)*(R**2))/(((D-C*Zi)**2+(R*C)**2)*Z)
t=100*((E**2-1+F**2)**2+4*(F**2))/(((E+1)**2+F**2)**2)
  tmp=t
if (t(i,j)<tmp)then
    tmp=t(i,j)
end if
end do
write(*,"('H1:'F10.4,'H2:'F10.4,'t:'F11.4)")H1,H2,t
write(10,"(F10.4,' 'F10.4,''F11.4)")H1,H2,t
end do
close(10)
stop
end
回复

使用道具 举报

17

帖子

5

主题

0

精华

入门

F 币
84 元
贡献
51 点
 楼主| 发表于 2020-7-10 10:19:17 | 显示全部楼层
显示错误 错误        1         error #6410: This name has not been declared as an array or a function.   [T]        D:\\2019\tyt20\tyt20\tyt20.f90       

58

帖子

0

主题

0

精华

专家

F 币
418 元
贡献
184 点

规矩勋章元老勋章新人勋章

发表于 2020-7-10 13:12:12 | 显示全部楼层
T 没有被定义为数组
天之道,损有余而补不足
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2020-8-6 23:58

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表