| 查看: 213 | 回复: 0 | |||
[交流]
【求助】大家看下这个程序是说明什么问题的？
|
|
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 ] |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
