求助，关于Fortran开头program+（tape)的用法

Harry · 发表于 2015-12-26 15:12:04

[Fortran] 纯文本查看 复制代码

01!Master plane
02    program Bok1(input,output,tape10=input,tape12=output,tape2,tape4,tape6,tape8)
03    dimension title(12)
04 
05    common/Contro/npoin,nelem,nnode,ndofn,ndime,nstre,ntype,ngaus,nprop,nmats,nvfix,nevab,&
06                   icase,ncase,itemp,iprob,nprob
07    common/Lgdata/coord(80,2),props(10,5),presc(40,2),asdis(160),eload(25,16),strin(4,225),&
08                   nofix(40),ifpre(40,2),lnods(25,8),matno(25)
09    common/Work/elcod(2,8),shape(8),deriv(2,8),dmatx(3,3),cartd(2,8),dbmat(3,16),bmatx(3,16),&
10                 smatx(3,16,9),posgp(3),weigp(3),gpcod(2,9),neror(24)
11 
12 
13!**********************************************************
14!  program for the solution of plane
15!   frameworks
16!**********************************************************
17!
18!  this segment controls the calling,in order, of all subroutines
19! 
20    !call Openfile
21    read(10,900) nprob
22 900 format(i5)
23    write(12,905) nprob
24 905 format(1ho,5x,23htotal no. of problems =,i5)
25    do 20 iprob=1,nprob
26    rewind 2
27    rewind 4
28    rewind 6
29    rewind 8
30    read(10,910) title
31 910 format(12a6)
32    write(12,915) iprob,title
33 915 format(/////,6x,12hproblem no. ,i3,10x,12a6)
34!
35!  call the subroutine which reads most of 
36!  the problem data
37! 
38    call Input
39!
40!  next create the element stiffness file
41!
42    call Stifps
43    do 10 icase=1,ncase
44!
45!  compute loads, after reading the relevant
46!  extra data
47!
48    call Loadps
49!
50!  merge and solve the resulting equations
51!  by the prontal solver
52! 
53    call Front
54!
55!  compute the stresses in all the elements
56!
57    call Streps
58 10 continue
59 20 continue
60    stop
61  end