您所在位置：主题：求助 高手 找错
00illusion00
[专家分：0] 发布于 2010-11-05 13:36:00
找错，整个程序有3000行 这是开始的一部分，我郁闷坏了
以前的程序，现在用6.6打开出错
dimension grain(200,200)
real AVG_MISORIENTION0(200,200),AVG_MISORIENTION(200,200),MAX_MISORIENTION(200,200),l(200),num(300)
REAL Node_STRAIN(200,200),STORED_ENERGY(200,200),delE(6),TEMP_ENERGY(6),AVG_D0(200,200),AVG_D(200,200),MAX_D(200,200),AVG_BOUNDARY_ENERGY(200,200)
integer f(5),Nuc(200,200),NODE_DIRECTION2(5000,200),MEDIUM_NODE_DIRECTION(18)
character f0*3,f11*1,f1*1,f2*1,f3*1,f4*1,f5*4,ff*12
integer i,j,k,Q,n,nn,DIRECTION_NUMBER,num1,cnuc
integer T,xxd,yyd,num2,xxd2,yyd2,directionm,un_rec_node,MCS
real ENG0,r,Boltz,jb,E,D,dd,Du,m,max_strain,MOBILITY,MOBILITY0,MENG,MEAN_STRAIN,MEDIUM,RATIO_ANGDF
real MEDIUM_D0,MEDIUM_MISORIENTION0,MEDIUM_Node_STRAIN,NUCLEATION_NUM,SUBGRAIN_NUM,angdf,angdf0
real tt,DELttmin,k3,k4
common /c1/NEAR_NODE_DIRECTION(18),NODE_DIRECTION(200,200)
common /c2/NEAR_NODE_DIRECTION1(6),NODE_DIRECTION1(1000,1000)
common /c3/grain,DIS_MEAN_DENSITY
common /c4/avgsize,MCS,dd
common /c5/xd(18),yd(18)
jb=3.956E-9
Q=180
T=653
BOLTZ=1.38E-23
angdf0=15.0 !微观取向参数
RENG=0.324 !J/M^2
BURGERS_VECTOR=1.65E-10 !柏氏向量AL:1.65E-10;IF:2.03E-10
print*,'MEAN_STRAIN=0.8'
read*,MEAN_STRAIN !平均应变
!平均储存能
MEAN_STORED_ENERGY=ENG_STRAIN(MEAN_STRAIN)
!根据平均储存能确定的网格常数
D_CRYm=4*RENG/MEAN_STORED_ENERGY
dd=D_CRYm
print*,MEAN_STORED_ENERGY
print*,"网格常数dd=",dd
Du0=1.71E-4 !AL扩散系数(m2/s)
!Du0= If钢的扩散系数
DIFFUSE_ACTIVATION_EBERGY=1.089E-19 !Al扩散激活能(J)
!DIFFUSE_ACTIVATION_EBERGY=6.56E4 !扩散激活能(J/mol)
du=du0*exp(-DIFFUSE_ACTIVATION_EBERGY/(BOLTZ*T))
MOBILITY0=(BURGERS_VECTOR*BURGERS_VECTOR*Du)/(6*BOLTZ*T)
print*,du,MOBILITY0
C1=1/18.0 ! 是晶界处储存能的量度
KT=7.21
goto 10000
numnuc=0
!给每个节点赋予取向
do 210 i=1,500
do 220 j=1,500
Call random_number(ran1)
DIRECTION_NUMBER=int(ran1*Q)+1
NODE_DIRECTION1(2*i-1,2*j-1)=DIRECTION_NUMBER
Call random_number(ran2)
DIRECTION_NUMBER=int(ran2*Q)+1
NODE_DIRECTION1(2*i,2*j)=DIRECTION_NUMBER
NODE_DIRECTION1(2*i-1,2*j)=0
NODE_DIRECTION1(2*i,2*j-1)=0
220 continue
以前的程序，现在用6.6打开出错
dimension grain(200,200)
real AVG_MISORIENTION0(200,200),AVG_MISORIENTION(200,200),MAX_MISORIENTION(200,200),l(200),num(300)
REAL Node_STRAIN(200,200),STORED_ENERGY(200,200),delE(6),TEMP_ENERGY(6),AVG_D0(200,200),AVG_D(200,200),MAX_D(200,200),AVG_BOUNDARY_ENERGY(200,200)
integer f(5),Nuc(200,200),NODE_DIRECTION2(5000,200),MEDIUM_NODE_DIRECTION(18)
character f0*3,f11*1,f1*1,f2*1,f3*1,f4*1,f5*4,ff*12
integer i,j,k,Q,n,nn,DIRECTION_NUMBER,num1,cnuc
integer T,xxd,yyd,num2,xxd2,yyd2,directionm,un_rec_node,MCS
real ENG0,r,Boltz,jb,E,D,dd,Du,m,max_strain,MOBILITY,MOBILITY0,MENG,MEAN_STRAIN,MEDIUM,RATIO_ANGDF
real MEDIUM_D0,MEDIUM_MISORIENTION0,MEDIUM_Node_STRAIN,NUCLEATION_NUM,SUBGRAIN_NUM,angdf,angdf0
real tt,DELttmin,k3,k4
common /c1/NEAR_NODE_DIRECTION(18),NODE_DIRECTION(200,200)
common /c2/NEAR_NODE_DIRECTION1(6),NODE_DIRECTION1(1000,1000)
common /c3/grain,DIS_MEAN_DENSITY
common /c4/avgsize,MCS,dd
common /c5/xd(18),yd(18)
jb=3.956E-9
Q=180
T=653
BOLTZ=1.38E-23
angdf0=15.0 !微观取向参数
RENG=0.324 !J/M^2
BURGERS_VECTOR=1.65E-10 !柏氏向量AL:1.65E-10;IF:2.03E-10
print*,'MEAN_STRAIN=0.8'
read*,MEAN_STRAIN !平均应变
!平均储存能
MEAN_STORED_ENERGY=ENG_STRAIN(MEAN_STRAIN)
!根据平均储存能确定的网格常数
D_CRYm=4*RENG/MEAN_STORED_ENERGY
dd=D_CRYm
print*,MEAN_STORED_ENERGY
print*,"网格常数dd=",dd
Du0=1.71E-4 !AL扩散系数(m2/s)
!Du0= If钢的扩散系数
DIFFUSE_ACTIVATION_EBERGY=1.089E-19 !Al扩散激活能(J)
!DIFFUSE_ACTIVATION_EBERGY=6.56E4 !扩散激活能(J/mol)
du=du0*exp(-DIFFUSE_ACTIVATION_EBERGY/(BOLTZ*T))
MOBILITY0=(BURGERS_VECTOR*BURGERS_VECTOR*Du)/(6*BOLTZ*T)
print*,du,MOBILITY0
C1=1/18.0 ! 是晶界处储存能的量度
KT=7.21
goto 10000
numnuc=0
!给每个节点赋予取向
do 210 i=1,500
do 220 j=1,500
Call random_number(ran1)
DIRECTION_NUMBER=int(ran1*Q)+1
NODE_DIRECTION1(2*i-1,2*j-1)=DIRECTION_NUMBER
Call random_number(ran2)
DIRECTION_NUMBER=int(ran2*Q)+1
NODE_DIRECTION1(2*i,2*j)=DIRECTION_NUMBER
NODE_DIRECTION1(2*i-1,2*j)=0
NODE_DIRECTION1(2*i,2*j-1)=0
220 continue