| 查看: 832 | 回复: 28 | |||
| 本帖产生 1 个 程序强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
格式问题
|
|||
|
write(19,100) ikpt,(kpt(i),i=1,3),npw 100 FORMAT('k-point #',I3,': (',3f7.4,') npw=',I6) 这个句子哪里有问题？ |
回复此楼» 猜你喜欢
请教限项目规定
已经有3人回复
-
拟解决的关键科学问题还要不要写
已经有7人回复
-
存款400万可以在学校里躺平吗
已经有15人回复
-
Materials Today Chemistry审稿周期
已经有6人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
推荐一本书
已经有13人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
|
这是程序： 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
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2011-04-26 08:33:16
holmescn
金虫 (正式写手)
- 程序强帖: 37
- 应助: 1 (幼儿园)
- 金币: 1918.8
- 散金: 275
- 红花: 1
- 帖子: 699
- 在线: 102.6小时
- 虫号: 913482
- 注册: 2009-11-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-04-26 08:38:43
4楼2011-04-26 12:29:01