| 查看: 342 | 回复: 0 | ||
[求助]
各位大神给看看fortran语言编写的双温方程怎么理解
|
|
DO 1000 CLOCK=1,Ntype !************************************************************************* !!!!!!!!!!!!!!!//////////TTM Model\\\\\\\\\\\!!!!!!!!!!!!!!!!!!!!!!!!!!!!! DO 50 j=1,mm2 GVV(j)=0.0 50 CONTINUE DO 100 i=((CLOCK-1)*10+1),CLOCK*10 !FD步长是MD步长的10倍 if (CLOCK.LT.1200 .AND. CLOCK.GT.1000) then DO 110 j=1,mm Ce(j)=Ce0*Te(j) ok(j)=Ke0*B*Te(j)/(A*Te(j)**2+B*Tl(j+50)) !修改 f(j)=ok(j)*dt1/(Ce(j)*dx*dx) g(j)=g0*Tl(j+50)/Te(j) !temperature dependence of e-ph factor alf(j)=g(j)**2*dt1*dt1/(Cl*Ce(j)*(1+g(j)*dt1/Cl))-1-2*f(j)-g(j)*dt1/Ce(j) alf(j)=alf(j)/f(j) blf(j)=-Te(j)-(Wa/Ce(j))*exp(-alpha*j*dx)*exp(-4*0.69315*((i-11000)/1000.0)**2)-g(j)*dt1*Tl(j+50)/(Ce(j)*(1+g(j)*dt1/Cl)) !修改 blf(j)=blf(j)/f(j) |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
