24小时热门版块排行榜

返回列表

215215

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 388.2
散金: 134
红花: 1
帖子: 125
在线: 53.5小时
虫号: 860062
注册: 2009-09-29
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】大家看下这个程序是说明什么问题的？

implicit real*8 (a-h,o-z)
   dimension dump(20),num(2),x(100),y(100),z(100),chg(100,100,800)
   dimension chgt(1000),a(3,3)

   write(*,*) 'Choose input file (LOCPOT=1, CHGCAR=2):'
   read(*,*) input
   if (input.eq.1) then
   open(15,file='LOCPOT')
   write(6,*) ' Enter fermi energy: '
   read(*,*) ef
   elseif (input.eq.2) then
   open(15,file='CHGCAR')
   ef=0.0
   elseif ((input.ne.1).and.(input.ne.2)) then
   write(*,*) ' INPUT ERROR, input must equal to 1 or 2 '
   stop
   endif

   write(*,*) 'Spin polarized calculation? (no=1, yes=2):'
   read (*,*) ispin
   if (ispin.eq.1) then
   if (input.eq.1) open(16,file='workfn.dat')
   if (input.eq.2) open(16,file='chgave.dat')

[ Last edited by 215215 on 2010-8-27 at 13:03 ]

回复此楼

» 猜你喜欢

nevergiveup!

1楼 2010-07-20 20:25:27

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

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

返回列表