| 查看: 955 | 回复: 28 | ||
| 本帖产生 1 个 程序强帖 ，点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
格式问题
|
||
|
write(19,100) ikpt,(kpt(i),i=1,3),npw 100 FORMAT('k-point #',I3,': (',3f7.4,') npw=',I6) 这个句子哪里有问题？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
holmescn
金虫 (正式写手)
- 程序强帖: 37
- 应助: 1 (幼儿园)
- 金币: 1918.8
- 散金: 275
- 红花: 1
- 帖子: 699
- 在线: 102.6小时
- 虫号: 913482
- 注册: 2009-11-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-04-26 08:38:43
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2011-04-26 08:33:16
4楼2011-04-26 12:29:01
|
这是程序： program WAVECAR integer nkpt,nband,emax,ikpt,npw,iband real*8 A kpt dimension A(3,3),kpt(3) open(unit=12,file="WAVECAR",status="old",form="unformatted"  open(unit=19,file="wave",status="new" read(12,*) nkpt,nband,emax,((A(i,j),i=1,3),j=1,3) write(19,*) write(19,*) 'nkpt =',nkpt write(19,*) 'nband =',nband write(19,*) 'emax =',emax write(19,*) 'A=' write(19,'(3X,3(1X,f8.3))') (A(i,1),i=1,3) write(19,'(3X,3(1X,f8.3))') (A(i,2),i=1,3) write(19,'(3X,3(1X,f8.3))') (A(i,3),i=1,3) read(12,*) code do ikpt=1, nkpt read(12,*) npw, (kpt(i),i=1,3) write(19,*) write(19,100) ikpt,(kpt(i),i=1,3),npw 100 FORMAT('k-point #',I3,': ',(3f7.4),' npw=',I6) write(19,*) 'kpt =',kpt write(19,*) 'npw =',npw write(19,*) " band energy weight" c ALLOCATE(coef(npw)) ALLOCATE(ac(npw,200)) ic = 1 DO iband = 1, nband read(12,*) eval, fweight, (coef(i),i=1,npw) c write(19,'(5X,I3,5X,f8.4,5x,f8.4)') iband, dreal(eval), fweight c IF(iband >= nbandmin.and.iband.le.nbandmax) then DO i = 1,npw ac(i,ic) = coef(i) ENDDO ic = ic+1; end if ENDDO DEALLOCATE(coef) ENDDO close(unit=12) close(unit=19) c stop end c----------------------- end --------------------------- |
5楼2011-04-26 12:31:45
