作者kc5566 (56专业文)
看板Fortran
标题[请益]主程式的变数与副程式共用?
时间Wed Jan 9 10:55:40 2013
program
implicit none
integer::num,i,j
integer,parameter::ROW=100,COL=100 !矩阵大小
real::U(ROW)
real::Matrix_Bij(ROW,COL)
real::T(ROW) !数值解
real::Y(ROW) !理论解
real::b,a,G,E,pr,n,L,f,&
xi,xj,bij,bij1,bij2,bij3,bij4,&
uyi,PI,a1,a2,a3,a4,Tyy
open (9,file='input.txt')
read(9,*) f,b,num,PI,E,Pr,uyi
open(10,file='output.txt')
G=E/(2*(1+pr))
L=b*f
write(10,*)"f=",f
write(10,*)"num=",num
write(10,*)"E=",E
write(10,*)"Pr=",pr
write(10,*)"G=",G
write(10,*)"b=",b
write(10,*)"L=",L*f
write(10,*)"uyi:",uyi
call slove_Bij
(num,i,j,xi,xj,a,bij,bij1,bij2,bij3,bij4,pr,G,Matrix_Bij,U,uyi,L,b,a1,a2,a3,a4)
call gauss(U,Matrix_Bij,T,num,uyi)
stop
end program
!--------------------slove_Bij------------------------
subroutine slove_Bij
(num,i,j,xi,xj,a,bij,bij1,bij2,bij3,bij4,pr,G,Matrix_Bij,U,uyi,L,b,a1,a2,a3,a4)
implicit none
integer,parameter::ROW=100,COL=100 !矩阵大小
integer::num,i,j
real::U(ROW)
real::Matrix_Bij(ROW,COL)
real::T(ROW)
real::a,b,G,xi,xj,pr,n,L,bij,bij1,bij2,bij3,bij4,uyi,PI,E,a1,a2,a3,a4
PI=3.14159
E=5000
pr=0.25
G=E/(2*(1+pr))
uyi=-0.04
上面这些参数已经在主程式的时候开启了，
如何让这些参数在副程式的时候直接读取不用在定义一次
如果要用COMMON这个指令是要在主程式就先宣告吗?
请版上各位高手指教一下小弟了~
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.120.28.103
1F：→ kkarthur:最方便的方法就是利用module 01/09 11:12