早教吧作业答案频道 -->其他-->
急求fortran大神帮忙,编译出现Theshapesofthearrayexpressionsdonotconform.[U]ubroutineInit(U,dx,dy)implicitdoubleprecision(a-h,o-z)common/Gdef/Gama,Pi,Jx,Jy,dLx,dLy,TT,Sfdimension::U(Jx+1,Jy+1,4)u=3.0v=0.0rou=1.0p=0.71429dx=dLx/
题目详情
急求fortran大神帮忙,编译出现The shapes of the array expressions do not conform. [U]
ubroutine Init(U,dx,dy)
implicit double precision (a-h,o-z)
common/G_def/Gama,Pi,Jx,Jy,dLx,dLy,TT,Sf
dimension ::U(Jx+1,Jy+1,4)
u=3.0
v=0.0
rou=1.0
p=0.71429
dx=dLx/Jx
dy=dLy/Jy
do i=2,Jx/4
do j=1,Jy
U(i,j,1)=rou
U(i,j,2)=rou*u
U(i,j,3)=rou*v
U(i,j,4)=p/(Gama-1.0)+rou*(u*u+v*v)/2.0
end do
end do
do i=Jx/4+1,Jx
do j=Jy/5+1,Jy
U(i,j,1)=rou
U(i,j,2)=rou*u
U(i,j,3)=rou*v
U(i,j,4)=p/(Gama-1.0)+rou*(u*u+v*v)/2.0
end do
end do
ubroutine Init(U,dx,dy)
implicit double precision (a-h,o-z)
common/G_def/Gama,Pi,Jx,Jy,dLx,dLy,TT,Sf
dimension ::U(Jx+1,Jy+1,4)
u=3.0
v=0.0
rou=1.0
p=0.71429
dx=dLx/Jx
dy=dLy/Jy
do i=2,Jx/4
do j=1,Jy
U(i,j,1)=rou
U(i,j,2)=rou*u
U(i,j,3)=rou*v
U(i,j,4)=p/(Gama-1.0)+rou*(u*u+v*v)/2.0
end do
end do
do i=Jx/4+1,Jx
do j=Jy/5+1,Jy
U(i,j,1)=rou
U(i,j,2)=rou*u
U(i,j,3)=rou*v
U(i,j,4)=p/(Gama-1.0)+rou*(u*u+v*v)/2.0
end do
end do
▼优质解答
答案和解析
U和u一样,fortran不分大小写的,最后几行都是数等于数组,所以错了
看了急求fortran大神帮忙,编...的网友还看了以下:
还是matlab画图问题,要输入一个变量的clc; clear all; h=input('h=' 2020-05-16 …
怎样用MATLAB算函数:y=h*cos(pi*x/L)/2+h/2的曲率半径;其中L=20;h= 2020-05-16 …
已知a+b+c=H a+b+e=J a+d+e=K b+c+d=M c+d+e=N 求a=?b=? 2020-05-16 …
v=h*pi*r^2=>h=v/(pi*r^2)中的*手写是什么数字 2020-05-16 …
求一数学几何证明题、今天二测的.三角形ABC,角C为钝角,点E、H是边AB上的点,点K和M分别是边 2020-06-03 …
如图,△ABC中,AD平分∠BAC,EF⊥AD交AB于点E,交AC于点F,交BC的延长线于点H.求 2020-06-27 …
空间四边形ABCD中,E,F是AB,AD的中点,G,H在BC,DC上,且BG:GC=DH:HC=1 2020-06-27 …
求证E,F,G,M,N,H六点共面已知正方体ABCD-A1B1C1D1中,点E,F,G,H,M,N 2020-07-09 …
A+B+C=84,D+E+F=111,H+I+J=138,A+D+H=124,B+E+I=148,C 2020-12-14 …
点O是以四边形abcd的外接园和内切园的园心,内切园与四边形各边分别相切于点e,f,g,h.求证:四 2020-12-25 …