| ²é¿´: 204 | »Ø¸´: 0 | |||
215215Òø³æ (СÓÐÃûÆø)
|
[½»Á÷]
¡¾ÇóÖú¡¿´ó¼Ò¿´ÏÂÕâ¸ö³ÌÐòÊÇ˵Ã÷ʲôÎÊÌâµÄ£¿
|
|
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 ] |
» ²ÂÄãϲ»¶
265Çóµ÷¼Á
ÒѾÓÐ9È˻ظ´
085600²ÄÁÏÓ뻯¹¤µ÷¼Á
ÒѾÓÐ20È˻ظ´
ר˶ 351 086100 Ò²ÊÇ¿¼µÄ²Ä¿Æ»ù ±¾¿ÆÒ²ÊDzÄÁÏ
ÒѾÓÐ6È˻ظ´
085600ר˶²ÄÁÏÓ뻯¹¤348·ÖÇóµ÷¼Á
ÒѾÓÐ10È˻ظ´
085600 295·ÖÇóµ÷¼Á
ÒѾÓÐ21È˻ظ´
285Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
Ò»Ö¾Ô¸°²»Õ´óѧ0817»¯Ñ§¹¤³ÌÓë¼¼Êõ£¬Çóµ÷¼Á
ÒѾÓÐ9È˻ظ´
Ò»Ö¾Ô¸0817»¯Ñ§¹¤³ÌÓë¼¼Êõ£¬Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
271·ÖÇóµ÷¼ÁѧУ
ÒѾÓÐ3È˻ظ´
ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£©
ÒѾÓÐ7È˻ظ´















»Ø¸´´ËÂ¥