24小时热门版块排行榜

返回列表

cquccc

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 17
散金: 385
红花: 1
帖子: 597
在线: 294.8小时
虫号: 1085708
注册: 2010-08-31
性别: GG
专业: 理论和计算化学

[求助] Fortran 产生新文件的文件名长度改变方法

用FORTRAN 编译成功一程序，40行代码以内的那种，发现新生成的文件名的长度不够，但是总是修改不成功，请大神看一眼。附代码和拟转化数据。

   PROGRAM GROUP
   IMPLICIT DOUBLE PRECISION (A-H,O-Z)
   INTEGER B,F,G,H,PENT
   CHARACTER*3 PGROUP,filin*12,filen*20,base*5
   INTEGER ORDER
C -----------------------------------------------------------------
C    FORTRAN 77 WITH NAG LIBRARY ROUTINES F02ABF AND F04AAF -----------------------------------------------------------------
   PARAMETER (NMAX = 120)
   PARAMETER (MMAX = NMAX/2+2)
   DIMENSION A(NMAX,NMAX),B(NMAX,3),C(MMAX,MMAX)
   DIMENSION F(MMAX),G(3*NMAX),H(MMAX,6)
   DIMENSION ATOM(3,NMAX),BOND(3,3*NMAX),FACE(3,NMAX)
   DIMENSION AXIS(3,5*NMAX)
   DIMENSION PENT (31)
C
   read(5,2084) filin
2084 format(a12)
   open(unit=8,status='old',file=filin)

   inumb=0

1  READ (8,*) N
   M = N/2+2

   read(8,*) pent
   DO 24 I=1,M
   iii=pent(i)
24    f(i)=iii

C
   inumb=inumb+1
C
   ICONST=ICHAR('0')
   IC1=inumb/100+ICONST ! WHY /100
   IR=MOD(inumb,100) ! WHY TO GET MOD
   IC2=IR/10+ICONST
   IC3=MOD(IR,10)+ICONST
   FILEN=filin(1:5)//CHAR(IC1)//CHAR(IC2)//CHAR(IC3)//'.gjf'
   open(unit=7,file=filen,status='new',form='formatted')

   GO TO 1
   END

附件中question为操作说明和问题。
附件中fort为拟转化文件。
附件中transfer为编译后的可执行程序。

回复此楼