24小时热门版块排行榜

返回列表

凌霄12

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1
帖子: 2
在线: 10.6小时
虫号: 3517203
注册: 2014-11-04
专业: 物理电子学

[求助] 各位大神给看看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)

回复此楼

» 猜你喜欢

1楼 2016-07-07 15:10:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主凌霄12 的主题更新

返回列表