作者psallen (psallen)
看板Fortran
标题[问题] 程式码请教
时间Wed Feb 15 00:01:50 2012
小弟是fortran新手,最近在看一段程式码要将其改成c语言,
但目前有一段一直看不懂想请教版上的高手能否帮忙看看,
感激不尽!
撷取片段程式码如下:
PARAMETER (M1=18,N1=2**M1)
implicit real*8 (a-h,o-z)
parameter (mtot=1000000)
REAL*8 X(mtot),vvz(mtot,6),xout(mtot,6),xx(n1),xall(n1,6)
complex*16 xfl(n1)
c**
c** nn=number of input responses
c** T=period
c**
write(*,*) ' [1]dt (in inp1) [2] Nterms (1-6) [3] M (13=16s)'
read (*,*) edt,nterm,M
N=2**M
if (edt>10.0) edt=1.0/edt
do i=1,mtot
read (1,*,iostat=kk) ajj,(vvz(i,j),j=1,nterm)
if (kk.ne.0) exit
enddo
write(*,*) ' dt=',edt
看不懂的地方是倒数第四行read那行代表的意思?
因为我执行过第一次的write和read会显示出来,而第四行的read却不会显示出来
要我输入东西,之前留下的执行档事会在我输入完edt,nterm,M的值时要我输入input档,
但我撷取出来这段产生的执行档却输入完edt,nterm,M的值时就跳到最後一个write了,
麻烦大家了
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.121.203.142
※ 编辑: psallen 来自: 122.121.203.142 (02/15 00:07)